Displaying 1 result from an estimated 1 matches for "ratio_exp123".
2010 May 05
1
rcorr p-values for pearson's correlation coefficients
...the correlation matrix (9275 X 9275) and pvalue matrix (9275 X 9275)
obtained using rcorr function, 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...