use strict;
use warnings;
use Acme::Tools;

my %data;
while (<DATA>) {
    my ($test, @vals) = (split)[1..4];
    my $col = 3;
    for my $val (@vals) {
        push @{ $data{$test}{$col} }, $val;
        $col++;
    }
}
#use Data::Dumper; print Dumper(\%data);

for my $test (sort keys %data) {
    for my $col (sort keys %{ $data{$test} }) {
        my $med = median(@{ $data{$test}{$col} });
        print "test=$test, col=$col, med=$med\n";
    }
}

__DATA__
contig1 test1 1e-28 28 55
contig1 test2 1e-10 22 54
contig2 test1 1e-10 24 78
contig3 test2 10 78 57
contig4 test3 1e-5 200 55
contig4 test2 10 100 43

##</code><code>##

test=test1, col=3, med=5e-11
test=test1, col=4, med=26
test=test1, col=5, med=66.5
test=test2, col=3, med=10
test=test2, col=4, med=78
test=test2, col=5, med=54
test=test3, col=3, med=1e-5
test=test3, col=4, med=200
test=test3, col=5, med=55