矩阵转数据框和余弦定理

矩阵转数据框

#!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";