Question: Arguments for camera (gene set analysis)
1
gravatar for claire.n.levy
3.3 years ago by
United States
claire.n.levy10 wrote:

Hi Bionconductors,

I've analyzed my microarray data for DE probes using limma and I'd like to use camera for gene set analysis, but I'm not quite sure I'm putting the right information in to the arguments. 

I am interested in whether genes in MSigDB Human C7 collection sets (downloaded from http://bioinf.wehi.edu.au/software/MSigDB/ ) are more DE compared to genes not in that set. Most of the DE probes that I found were in just one of the contrasts, so I'm just looking at that contrast.

With regard to the camera arguments, my questions are:

Is y supposed to be the exprs from the Eset that I used to do my DE analysis, i.e. the whole thing, NOT just exprs of the probes that were found to be DE in my contrast of interest?

For the index argument, if I am using ids2indices function, are the identifiers supposed to be the Entrez IDs of the probes that were DE in my contrast of interest?

Finally, in the result, does NGenes represent the number of genes from my input (index) that are also in that particular set?

Thanks in advance for your sage advice!

Claire Levy

University of Washington/Fred Hutchinson Cancer Research Center

Seattle, WA USA

 

camera gsea • 1.6k views
ADD COMMENTlink modified 3.3 years ago by Steve Lianoglou12k • written 3.3 years ago by claire.n.levy10
Answer: Arguments for camera (gene set analysis)
1
gravatar for Gordon Smyth
3.3 years ago by
Gordon Smyth37k
Walter and Eliza Hall Institute of Medical Research, Melbourne, Australia
Gordon Smyth37k wrote:

Any gene set testing analysis has to be done independently of the regular DE analysis. At no stage can you use information about which genes were individually DE, otherwise you will be double-dipping and the gene set p-values will be invalid.

'y' should be the same Eset that you passed to lmFit and eBayes for the regular DE analysis.

'index' simply identifies which genes are in each set. The 'identifers' are the gene IDs from 'y'.

'NGenes' is the number of genes in the set for which you have data, i.e., which correspond to a row of 'y'. It's the total size of the set.

ADD COMMENTlink modified 3.3 years ago • written 3.3 years ago by Gordon Smyth37k

Great, thanks for the clarification and also for smoothing the way with the ids2indices function.

ADD REPLYlink written 3.3 years ago by claire.n.levy10
Answer: Arguments for camera (gene set analysis)
0
gravatar for Steve Lianoglou
3.3 years ago by
Denali
Steve Lianoglou12k wrote:

Also, because we're talking about camera parameters here, I'll drop a link to this conversation about the inter.gene.cor parameter which was added somewhat recently:

C: Do the outputs of Limma's competitive gene set methods (camera, romer) require a

Which may be of use to the OP, but also to the next traveler who lands on this page by way of google.

ADD COMMENTlink written 3.3 years ago by Steve Lianoglou12k

Thanks, I need to add that thread to my reading list!

ADD REPLYlink written 3.3 years ago by claire.n.levy10
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 16.09
Traffic: 273 users visited in the last hour