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
|
#!/usr/bin/env perl
use strict;
use warnings;
my $ALLOWED_DELTA_PCT = 5;
main: {
&verify("edgeR_genes/diffExpr.P0.001_C2.matrix", 26);
&verify("edgeR_genes/Trinity_genes.counts.matrix.heatshock_vs_log_growth.edgeR.DE_results.P0.001_C2.DE.subset.GOseq.enriched", 112);
&verify("edgeR_genes/Trinity_genes.counts.matrix.heatshock_vs_log_growth.edgeR.DE_results.P0.001_C2.log_growth-UP.subset.GOseq.enriched", 46);
exit(0);
}
####
sub verify {
my ($filename, $target_count) = @_;
my $count = `cat $filename | wc -l`;
$count =~ s/\s+//g;
my $delta = abs($count - $target_count);
my $pct_delta = $delta / $target_count * 100;
if ($pct_delta <= $ALLOWED_DELTA_PCT) {
print STDERR "$filename, count: $count DE features OK.\n";
}
else {
die "Error, $filename has count: $count, $pct_delta % different from expected: $target_count";
}
return;
}
|