This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Finding most co-expressed genes of input gene based on Pearson correlation

I have a matrix of microarray (around 30000 genes X 800 samples) as a back-end data set. In the workflow, if someone input a gene of their interest, I need to find it's highly co-expressed genes based on back-end data set. For that, I need to calculate pairwise Pearson correlation of input gene against data set (in this case, 30000 pairs). I have tried psych package before and for this pairwise calculation, there might be some other better method. Also,I assume this process take a little long time. Can any experienced person can suggest on this? Thanks

pearson-correlation microarray
• 3,146 views • link
updated by Ram Ram • written by kozhaki.seq kozhaki.seq

So you want a faster method for calculating the correlation matrix? Although you don't need WGCNA, I do believe the people in WGCNA has a faster implementation for the calculation of correlation. you can have a look

• 0 views • link
written by Sam Sam

Yes, since my matrix dimension is high, I am concern about the time..Yes, I will look into it @Sam

• 0 views • link
written by kozhaki.seq kozhaki.seq

1 answer

amap::Dist, a R function, can calculate the distance or correlation in parallel. See manual http://www.inside-r.org/packages/cran/amap/docs/Dist

• 0 views • link
updated by Ram Ram • written by Zhilong Jia Zhilong Jia

Log in to answer this question.

More posts like this