Displaying 1 result from an estimated 1 matches for "y_gp123".
2010 May 05
1
rcorr p-values for pearson's correlation coefficients
...ne1=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))
{
if(y_cor_p123$r[i,j]>0.8)
{
gene1=rbind(gene1,rownames(y_cor_p123$r[i,j,drop=F]))
gene2=rbind(gene2,colnames(y_cor_p123$r[i,j,drop=F]))
pcc=rbind(pcc,y_cor_p123$r[i,j])
}
}
}
y_gp123=cbind(gene1,gene2,pcc)
colnames(y_gp123)=c("Gene1","Gene2","PCC")
y_gp123=y_gp123[-1,]
#Selecting p-values for gene pairs with pcc above 0.8
z=matrix(0,nrow(y_gp123),1)
for(i in 1:nrow(y_gp123))
{rowno=grep(y_gp123[i,1],rownames(y_cor_p123$P))
colno=grep(y_gp123[i,2],...