矩阵转数据框
#!perl
use warnings;
use strict;
open FA, $ARGV[0] or die $!;
my $head = <FA>;
chomp($head);
my @out = (split /\s+/, $head)[1..12];
print "family\tsample\tvalue\n";
my %hash;
while(<FA>)
{
chomp;
my @tmp = split;
my $n = 0;
foreach my $i(@tmp[1..12])
{
$n ++ if($i == 0);
}
next if($n > $ARGV[1]);
my $s = join "\t", @tmp[0..12];
$hash{$s} = 0;
}
foreach(keys %hash)
{
my ($f, @v) = split;
for(my $i = 0; $i < @v; $i ++)
{
print "$f\t$out[$i]\t$v[$i]\n";
}
}
#!perl
use warnings;
use strict;
open FA, $ARGV[0] or die $!;
my $head = <FA>;
chomp($head);
my @ids = (split /\t/, $head)[1..2];
my @id = map {$_ =~ s/-RPKM//; $_} @ids;
my ($inner_product, $vectorA, $vectorB);
while(<FA>)
{
chomp;
my @tmp = split;
if($tmp[1] == 0.001 and $tmp[2] == 0.001)
{
next;
}else{
$inner_product += $tmp[1] * $tmp[2];
$vectorA += $tmp[1] * $tmp[1];
$vectorB += $tmp[2] * $tmp[2];
}
}
$vectorA = sqrt($vectorA);
$vectorB = sqrt($vectorB);
my $cos = $inner_product / ($vectorA * $vectorB);
print $cos . "\n";