Entering edit mode

6.2 years ago

mannoulag1
▴
120

Hi biostars,

I did a pearson correlation to my data (expression matrix), and I keep only the correlation >0.8 . How can I obtain the sub expression matrix of only these highly correlated genes. Thank you

```
data<-t(matrix)
cor = cor(data, use="pairwise.complete.obs", method="pearson")
cor<-cor[abs(cor)>0.8]
```

Thank you Jean-Karim, I do this :

Then I have to remove the duplicated genes from 'correlated.genes' ?

This was just to give you quick pointer. What I think you want is to get unique indices. Something like:

Thank you Jean-karim