File: average.pl

package info (click to toggle)
trinityrnaseq 2.2.0%2Bdfsg-2
  • links: PTS, VCS
  • area: main
  • in suites: stretch
  • size: 212,452 kB
  • ctags: 5,067
  • sloc: perl: 45,552; cpp: 19,678; java: 11,865; sh: 1,485; makefile: 613; ansic: 427; python: 313; xml: 83
file content (39 lines) | stat: -rwxr-xr-x 804 bytes parent folder | download | duplicates (4)
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
#!/usr/bin/env perl

use strict;
use warnings;
use FindBin;
use lib ("$FindBin::RealBin/../../PerlLib");
use BHStats;

my $count = 0;
my $sum = 0;
my @values;

while (<STDIN>) {
    chomp;
    $sum += $_;
    $count++;
    push (@values, $_);
}

@values = sort {$a<=>$b} @values;

if ($count) {
    my $average = ($sum/$count);
    my $median_pos = int ($count/2);
    print "\n"; 
    print "MIN: " . $values[0] . "\n";
    print "MAX: " . $values[$#values] . "\n";
    print "Sum: $sum\n";
    printf ("Average: %.2f\n", $average);
    print "Median: " . BHStats::median(@values) . "\n";
    my $stdev = BHStats::stDev(@values);
    printf ("stDev from Average: %.2f\n", $stdev);

    my $geoMean = &BHStats::geometric_mean(@values);
    if ($geoMean) {
	printf ("geoMean: %.2f\n", $geoMean);
    }
}