Revision 4816f2a300d7bbae1084b1940959ea35973104e1 authored by Mingzhi Lin on 23 August 2017, 14:12:59 UTC, committed by Mingzhi Lin on 23 August 2017, 14:12:59 UTC
1 parent 4453664
collect.go
package mcorr
import (
"encoding/json"
"fmt"
"os"
)
// PipeOutCorrResults pipe the the channel of CorrResults out to a file.
func PipeOutCorrResults(corrResChan chan CorrResults, outFile string) chan CorrResults {
c := make(chan CorrResults)
go func() {
defer close(c)
f, err := os.Create(outFile)
if err != nil {
panic(err)
}
defer f.Close()
encoder := json.NewEncoder(f)
for res := range corrResChan {
if err := encoder.Encode(res); err != nil {
panic(err)
}
c <- res
}
}()
return c
}
// Collect feed correlation results into boostrappers and return them.
func Collect(corrResChan chan CorrResults, numBoot int) []*Bootstrap {
// prepare bootstrappers.
bootstraps := []*Bootstrap{}
notBootstrap := NewBootstrap("all", 1.0)
notBootstrap.SetRandom(false)
bootstraps = append(bootstraps, notBootstrap)
for i := 0; i < numBoot; i++ {
id := fmt.Sprintf("boot_%d", i)
sampleRatio := 1.0
bootstraps = append(bootstraps, NewBootstrap(id, sampleRatio))
}
for corrResults := range corrResChan {
for _, bs := range bootstraps {
bs.Add(corrResults)
}
}
return bootstraps
}
// CollectWrite collects and writes the correlation results.
func CollectWrite(corrResChan chan CorrResults, outFile string, numBoot int) {
bootstraps := Collect(corrResChan, numBoot)
w, err := os.Create(outFile)
if err != nil {
panic(err)
}
defer w.Close()
w.WriteString("l,m,v,n,t,b\n")
for _, bs := range bootstraps {
results := bs.Results()
for _, res := range results {
w.WriteString(fmt.Sprintf("%d,%g,%g,%d,%s,%s\n", res.Lag, res.Mean, res.Variance, res.N, res.Type, bs.ID))
}
}
}

Computing file changes ...