Search
Question: Random effects with limma's camera()
0
gravatar for Stephen Hoang
2.3 years ago by
Stephen Hoang30 wrote:

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 

ADD COMMENTlink modified 2.3 years ago by Gordon Smyth33k • written 2.3 years ago by Stephen Hoang30
1
gravatar for Gordon Smyth
2.3 years ago by
Gordon Smyth33k
Walter and Eliza Hall Institute of Medical Research, Melbourne, Australia
Gordon Smyth33k 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 COMMENTlink modified 2.3 years ago • written 2.3 years ago by Gordon Smyth33k
2

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 REPLYlink written 7 months ago by jamie.gearing20

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

ADD REPLYlink written 7 months ago by Gordon Smyth33k

Thanks Gordon.

ADD REPLYlink written 7 months ago by jamie.gearing20
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 424 users visited in the last hour