Question: Arguments for camera (gene set analysis)
1
2.9 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?

Claire Levy

University of Washington/Fred Hutchinson Cancer Research Center

Seattle, WA USA

gsea camera • 1.4k views
modified 2.9 years ago by Steve Lianoglou12k • written 2.9 years ago by claire.n.levy10
Answer: Arguments for camera (gene set analysis)
1
2.9 years ago by
Gordon Smyth36k
Walter and Eliza Hall Institute of Medical Research, Melbourne, Australia
Gordon Smyth36k 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.

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

Answer: Arguments for camera (gene set analysis)
0
2.9 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.