my $c=Statistics::RankCorrelation->new($probehash{$ProbenameArray[$i]},$probehash{$ProbenameArray[$j]},sorted=>1);
my $corr_value=$c->spearman;
my $ZValue = 0.5*log((1.0+$corr_value)/(1.0-$corr_value));
$ZValue = $ZValue*sqrt(26-3);
my $corrPValue = 2.0*(1.0 - pnorm(abs($ZValue)));
No comments:
Post a Comment