Question: Enrichr analysis using ranked gene list
0
gravatar for halo22
3.5 years ago by
halo22180
Indianapolis, IN
halo22180 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 2.9 years ago by adnbps10 • written 3.5 years ago by halo22180
4
gravatar for Santosh Anand
3.5 years ago by
Santosh Anand5.2k
Santosh Anand5.2k 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 3.5 years ago by Santosh Anand5.2k

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 3.1 years ago • written 3.1 years ago by halo22180

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 3.1 years ago by Santosh Anand5.2k

Great. Thank you very much.

ADD REPLYlink written 3.1 years ago by halo22180
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: 1186 users visited in the last hour