Displaying 2 results from an estimated 2 matches for "aaex_b3242_12".
2010 Apr 30
0
extracting pairs from correlation matrix and p-value matrix
...values using rcorr function in library Hmisc. Now I wish to
filter out those pairs who's PCC value is above 0.8 cut off and
p-value is less than 0.05.
>library(Hmisc)
>mat_cor=rcorr(t(mat),type="pearson")
>head(mat_cor)
aaeA_b3241_14 aaeB_b3240_15 aaeR_b3243_15 aaeX_b3242_12
aas_b2836_14 aat_b0885_14 abgA_b1338_14 abgB_b1337_15 abgR_b1339_15
abgT_b1336_15
aaeA_b3241_14 1.00 0.12 0.64 0.21
0.10 -0.68 -0.61 -0.62 -0.66
-0.67
aaeB_b3240_15 0.12 1.00 -0.57 0.26...
2009 Oct 27
2
column names of a correlation matrix
...ish
to obtain pairs with highest correlation values. So, I did this
> b=matrix(data=NA,nrow=4217,ncol=1)
> rownames(b)=rownames(cor_expN)
> for(i in 1:4217){b[i,]=max(cor_expN[i,])}
> head(b)
[,1]
aaeA_b3241_14 0.7181912
aaeB_b3240_15 0.7513084
aaeR_b3243_15 0.7681684
aaeX_b3242_12 0.5230587
aas_b2836_14 0.6615927
aat_b0885_14 0.6344144
Now I want the corresponding columns for the above values. For that I tried
this
> c=matrix(data=NA,nrow=4217,ncol=1)
> for(i in 1:4217){b[i,]=colnames(max(cor_expN[i,]))}
And got the following error:
Error in b[i, ] = colnames(ma...