Saturday, September 20, 2008

P value of spearman

                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: