Question: Is it fair to compare the significance of overlap of 2 different sized gene sets with a third gene set using Fisher's exact test?
0
gravatar for ebrudermanver
4 months ago by
ebrudermanver40 wrote:

Say I have 3 sets A, B and C of genes chosen from an entire set of 20K genes, and I want to see whether set A has a more significant gene overlap with set C than set B. That's, I want to compute a p-value p associated with the overlap of A and C, and another p-value q associated with the overlap of B and C, and I want to check whether p < q. I know that the p-values p and q in this problem can be computed using Fisher's exact test. My question is as follows: If the set A is bigger than set B (e.g. 150 genes in A vs. 30 genes in B), is it still a fair comparison when I compare the overlaps of each set with C based on the p-values computed using Fisher's exact test?

statistics • 331 views
ADD COMMENTlink modified 4 months ago by Nicolas Rosewick5.9k • written 4 months ago by ebrudermanver40

My initial thought is that you might have to compare proportional overlaps and not the raw numbers themselves. And even them I'm not sure this makes a lot of sense in edge cases, e.g., comparing two set memberships which have lets say 10 and 10000 members respectively.

ADD REPLYlink written 4 months ago by mforde841.0k
1
gravatar for Jean-Karim Heriche
4 months ago by
EMBL Heidelberg, Germany
Jean-Karim Heriche13k wrote:

The overlap between two sets follows a hypergeometric distribution. When doing a hypergeometric test, the various sizes are taken into account so because Fisher's exact test is based on the hypergeometric distribution, it also takes sample size into consideration. Note that for this (and the Chi-squared test), you need to use counts not proportions.
So if you want to know if A∩C is more likely than B∩C, you can compare the p-values.

ADD COMMENTlink modified 4 months ago • written 4 months ago by Jean-Karim Heriche13k
0
gravatar for Nicolas Rosewick
4 months ago by
Belgium, Brussels, Université Libre de Bruxelles / Université de Liège
Nicolas Rosewick5.9k wrote:

You can use superExactTest R package for this : https://cran.r-project.org/web/packages/SuperExactTest/vignettes/set_html.html

ADD COMMENTlink written 4 months ago by Nicolas Rosewick5.9k
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: 748 users visited in the last hour