Question: limma barcodeplot(): calculation of enrichment score
gravatar for khandarius
4 weeks ago by
khandarius20 wrote:



I've performed gene set enrichment analysis with the limma-package's camera -function. I'm using limma's barcodeplot -function to visualize my results (one index set, no gene weights). I would like to know how the enrichment score on the plot (y-axis) is calculated - I can't find any mention of it in the help files or limma's user guide. An example of a barcode plot can be found here (figure 9):


Best regards,


ADD COMMENTlink modified 4 weeks ago by Aaron Lun21k • written 4 weeks ago by khandarius20
gravatar for Aaron Lun
4 weeks ago by
Aaron Lun21k
Cambridge, United Kingdom
Aaron Lun21k wrote:

It's done using the function tricubeMovingAverage, with the span defined with span.worm. Imagine a sliding window that is applied to the set of all genes (ordered by log-fold change for your contrast of interest). In each window, you calculate the "density" of genes in the specified gene set, where the contribution of each gene in the window is weighted according to its distance from the center. The computed value for each position of the window is the enrichment score. I would say that it's mostly for visual purposes and should not be used for any quantitative decision-making.

ADD COMMENTlink written 4 weeks ago by Aaron Lun21k

See also How barcode-plot enrichment is calculated?

ADD REPLYlink written 4 weeks ago by Gordon Smyth35k

Thank you Gordon. I don't know how I missed this post...

ADD REPLYlink written 4 weeks ago by khandarius20
Please log in to add an answer.


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