Random effects with limma's camera()
1
0
Entering edit mode
@stephen-hoang-6156
Last seen 3.8 years ago

Hi All,

Is there a way to incorporate the consensus correlation calculated by duplicateCorrelation() into the gene set enrichment test performed by camera()? In short, I would like to treat one of my factors as a random effect via duplicateCorrelation() (this is what I did for differential expression testing). As far as I can tell, camera() is only aware of fixed effects. I don't have the option of treating my random effect as a fixed effect because I end up with a design matrix that is not of full rank. It seems like my only option is to run camera() using only fixed effects, and treat the results as overly-conservative.

Thanks in advance,

Steve 

limma camera duplicatecorrelation • 1.1k views
ADD COMMENT
1
Entering edit mode
@gordon-smyth
Last seen 1 hour ago
WEHI, Melbourne, Australia

Sorry, you can't use duplicateCorrelation() with camera(). It has been on my to-do list for a while, but I haven't managed to find time to develop and implement it.

You can use duplicateCorrelation() with roast(), mroast() or romer().

ADD COMMENT
2
Entering edit mode

Can this now be addressed by using the recent pre-ranked cameraPR() function? i.e. run lmFit() with block and correlation, then eBayes() (or treat() if appropriate?) and use the t-statistics.

ADD REPLY
0
Entering edit mode

Yes, that's right (as long as you're happy with a preset value for the inter-gene correlation).

ADD REPLY
0
Entering edit mode

Thanks Gordon.

ADD REPLY

Login before adding your answer.

Traffic: 190 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