Question: Enrichr analysis using ranked gene list
0
gravatar for halo22
22 months ago by
halo22120
Indianapolis, IN
halo22120 wrote:

Hello,

I am trying to perform enrichment analysis and it seems that enrichr is a better option to do this. I want to understand if enrichr performs enrichment using hypergeometric method or the GSEA (Kolmogorov-Smirnov test) method, in case of hypergeometric method is it possible to provide a list of background genes and interesting/significant genes? Can enrichr handle a ranked list, rank could be anything like the t-statistic, fold-change etc?

Thanks

ADD COMMENTlink modified 15 months ago by adnbps0 • written 22 months ago by halo22120
3
gravatar for Santosh Anand
22 months ago by
Santosh Anand4.7k
Santosh Anand4.7k wrote:

Did you see their publication? https://academic.oup.com/nar/article-lookup/doi/10.1093/nar/gkw377

Computing enrichment

Enrichr implements three approaches to compute enrichment. The first one is a standard method implemented within most enrichment analysis tools: the Fisher exact test. This is a proportion test that assumes a binomial distribution and independence for probability of any gene belonging to any set. The second test is a correction to the Fisher exact test that we developed based on intuition. We first compute enrichment using the Fisher exact test for many random input gene lists in order to compute a mean rank and standard deviation from the expected rank for each term in each gene-set library. Then, using a lookup table of expected ranks with their variances, we compute a z-score for deviation from this expected rank, this can be a new corrected score for ranking terms. Alternatively, we combined the p-value computed using the Fisher exact test with the z-score of the deviation from the expected rank by multiplying these two numbers as follows:

And I'd say NO to all of your queries above.

ADD COMMENTlink written 22 months ago by Santosh Anand4.7k

Thank you Santosh. I do have one more questions for you. What does the old.P.value in the enrichR output mean? How is it different from P.value. I am doing a KEGG analysis using enrichR.

ADD REPLYlink modified 18 months ago • written 18 months ago by halo22120

I think they changed the P-value calculation in the algorithm and to be compatible with old results, they still provide those values. IMO, you can safely ignore that.

ADD REPLYlink written 18 months ago by Santosh Anand4.7k

Great. Thank you very much.

ADD REPLYlink written 18 months ago by halo22120
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.3.0
Traffic: 1778 users visited in the last hour