Random effects with limma's camera()

Question: Random effects with limma's camera()

3.8 years ago by

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 • 911 views

ADD COMMENT • link •

modified 3.8 years ago by Gordon Smyth ♦ 39k • written 3.8 years ago by Stephen Hoang • 30

Answer: Random effects with limma's camera()

3.8 years ago by

Gordon Smyth ♦ 39k

Walter and Eliza Hall Institute of Medical Research, Melbourne, Australia

Gordon Smyth ♦ 39k wrote:

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 • link modified 3.8 years ago • written 3.8 years ago by Gordon Smyth ♦ 39k

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 • link written 2.1 years ago by jamie.gearing • 40

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

ADD REPLY • link written 2.1 years ago by Gordon Smyth ♦ 39k

Thanks Gordon.

ADD REPLY • link written 2.1 years ago by jamie.gearing • 40

Similar posts • Search »