-
-type importstatsplot struct{}
-
-func (cmd *importstatsplot) RunCommand(prog string, args []string, stdin io.Reader, stdout, stderr io.Writer) int {
- err := cmd.Plot(stdin, stdout)
- if err != nil {
- log.Errorf("%s", err)
- return 1
- }
- return 0
-}
-
-func (cmd *importstatsplot) Plot(stdin io.Reader, stdout io.Writer) error {
- var stats []importStats
- err := json.NewDecoder(stdin).Decode(&stats)
- if err != nil {
- return err
- }
-
- p, err := plot.New()
- if err != nil {
- return err
- }
- p.Title.Text = "coverage preserved by import (excl X<0.65)"
- p.X.Label.Text = "input base calls ÷ sequence length"
- p.Y.Label.Text = "output base calls ÷ input base calls"
- p.Add(plotter.NewGrid())
-
- data := map[string]plotter.XYs{}
- for _, stat := range stats {
- data[stat.InputLabel] = append(data[stat.InputLabel], plotter.XY{
- X: float64(stat.InputCoverage) / float64(stat.InputLength),
- Y: float64(stat.TileCoverage) / float64(stat.InputCoverage),
- })
- }
-
- labels := []string{}
- for label := range data {
- labels = append(labels, label)
- }
- sort.Strings(labels)
- palette, err := colorful.SoftPalette(len(labels))
- if err != nil {
- return err
- }
- nextInPalette := 0
- for idx, label := range labels {
- s, err := plotter.NewScatter(data[label])
- if err != nil {
- return err
- }
- s.GlyphStyle.Color = palette[idx]
- s.GlyphStyle.Radius = vg.Millimeter / 2
- s.GlyphStyle.Shape = draw.CrossGlyph{}
- nextInPalette += 7
- p.Add(s)
- if false {
- p.Legend.Add(label, s)
- }
- }
- p.X.Min = 0.65
- p.X.Max = 1
-
- w, err := p.WriterTo(8*vg.Inch, 6*vg.Inch, "svg")
- if err != nil {
- return err
- }
- _, err = w.WriteTo(stdout)
- return err
-}