Displaying 1 result from an estimated 1 matches for "y_cor_p123".
2010 May 05
1
rcorr p-values for pearson's correlation coefficients
...tion, I wanted to select those pairs whose
PCC's are above 0.8 cut-off and then for those gene pairs their
corresponding p-values. I did following:
library(Hmisc)
#loading expression matrix
load("ratio_exp123.RData")
#calculation correlation and p-values of genes using rcorr function
y_cor_p123=rcorr(t(ratio_exp123),type="pearson")
save(y_cor_p123,file="y_cor_p123.RData")
diag(y_cor_p123$r)=0
#Selecting gene pairs with pcc value greater than 0.8
gene1=matrix(0,1,1)
gene2=matrix(0,1,1)
pcc=matrix(0,1,1)
for(i in 1:nrow(y_cor_p123$r))
{
for(j in 1:nrow(y_cor_p123$r))
{...