Displaying 1 result from an estimated 1 matches for "colno".
Did you mean:
coln
2010 May 05
1
rcorr p-values for pearson's correlation coefficients
...23$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],colnames(y_cor_p123$P))
z[i,]=y_cor_p123$P[rowno,colno]}
y_gp123=cbind(y_gp123,z)
> head(y_gp123)
Gene1 Gene2 PCC P-Value
[1,] "10003_f_at" "10166_at" "0.816270470619202" "0"
[2,] "10003_f_at&quo...