1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
|
#!/usr/bin/perl
use strict;
use Statistics::Descriptive;
use SVG::Graph;
use SVG::Graph::Data;
use SVG::Graph::Data::Datum;
my %stat;
my $graph = SVG::Graph->new(width => 1000, height => 800, margin => 40);
my $parentgroup = $graph->add_group;
my $statgroup = $parentgroup->add_group;
my $bargroup = $parentgroup->add_group;
my @data;
my $ymax = undef;
my $stat = Statistics::Descriptive::Full->new;
while(<>){
chomp;
if(/^#stat:/){
process_stat($_);
next;
}
my($x,$y) = split /\s+/;
$ymax = $y > $ymax ? $y : $ymax;
$stat->add_data($x);
push @data, SVG::Graph::Data::Datum->new(x => $x , y => $y);
}
$parentgroup->add_glyph('axis');
warn $stat{mean};
my @statdata = (
SVG::Graph::Data::Datum->new(x => $stat->min, y => 0),
# SVG::Graph::Data::Datum->new(x => $stat{mean}, y => $ymax),
# SVG::Graph::Data::Datum->new(x => $stat{mean} - $stat{standard_deviation}, y => $ymax),
# SVG::Graph::Data::Datum->new(x => $stat{mean} + $stat{standard_deviation}, y => 0),
SVG::Graph::Data::Datum->new(x => $stat{median}, y => $ymax),
SVG::Graph::Data::Datum->new(x => $stat{quartile1}, y => $ymax),
SVG::Graph::Data::Datum->new(x => $stat{quartile3}, y => 0),
SVG::Graph::Data::Datum->new(x => $stat->max, y => 0),
);
$statgroup->add_data(SVG::Graph::Data->new(data => \@statdata));
$statgroup->add_glyph('barflex', 'stroke' => 'black', 'fill' => 'black', 'stroke-opacity' => 1.0, 'fill-opacity' => 0.4 );
$bargroup->add_data(SVG::Graph::Data->new(data => \@data));
$bargroup->add_glyph('bar', 'stroke' => 'red', 'fill' => 'red', 'stroke-opacity' => 1.0, 'fill-opacity' => 0.8 );
print $graph->draw;
sub process_stat {
my $line = shift;
my($name,$value) = $line =~ /^#stat:(\S+)\t(\S+)$/;
$stat{$name} = $value;
}
|
