Finding most co-expressed genes of input gene based on Pearson correlation
1
1
Entering edit mode
6.1 years ago
kozhaki.seq ▴ 60

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 pairwise • 1.9k views
ADD COMMENT
0
Entering edit mode

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

ADD REPLY
0
Entering edit mode

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

ADD REPLY
2
Entering edit mode
6.1 years ago
Zhilong Jia ★ 2.0k

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

ADD COMMENT

Login before adding your answer.

Traffic: 1824 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6