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
|
#!/usr/bin/env perl
use strict;
use warnings;
my $usage = "usage: $0 tpm.matrix\n\n";
my $matrix_file = $ARGV[0] or die $usage;
open (my $fh, $matrix_file) or die $!;
my $header = <$fh>;
my @tpms;
while (<$fh>) {
chomp;
my @x = split(/\t/);
shift @x; # gene accession
my $max_tpm = shift @x;
while (@x) {
my $tpm = shift @x;
if ($tpm > $max_tpm) {
$max_tpm = $tpm;
}
}
push (@tpms, $max_tpm);
}
@tpms = reverse sort {$a<=>$b} @tpms;
my $min_tpm_thresh = int($tpms[0] + 0.5);
my $num_features = 1;
print "neg_min_tpm\tnum_features\n";
shift @tpms;
while (@tpms) {
my $tpm = shift @tpms;
$tpm = int($tpm+0.5);
if ($tpm < $min_tpm_thresh) {
print "" . (-1*$min_tpm_thresh) . "\t$num_features\n";
$min_tpm_thresh = $tpm;
}
$num_features++;
}
print "$min_tpm_thresh\t$num_features\n";
exit(0);
|
