I have several lists of genes, say A, B, C. I want to test if the overlap between A and B; A and C are significant or not. I am planning on doing a hypergeometric test using R.
Please help me to verify if this is correct or no:
phyper(q, m, n, k, lower.tail = FALSE, log.p = FALSE)
q: # of overlap between A and B/C - 1
m: # of genes in A
n: total # of genes in sample - # of genes in A
k: # of genes in B/C
Also, what if there is no overlap? Do I still do the same thing but with q=0 or q=-1?
Thanks a lot for the help!
•
written
by
wstla27 •
1 answer
So if I get this straight. You want to see if there is an enrichment of genes from A in either B or C. Overall it looks like your setup is good.
Your universe (n) is all the genes in across your samples.
The q-value depends on the question you want to ask. In your case you seem to be making a decision on asking the following:
P(Observed less than q overlaps) using use q-1
Another option however would be to ask:
P(Observed q or less overlaps) using q
As to your last question: The probability of observing less than 0 overlaps is 0 so doing that test wouldn't be very informative.
Log in to answer this question.
-
How to perform enrichment p-value for a motif
written by Apex92 •
-
Enricher - Hyper Geometric Test Details
written by Manuel Sokolov Ravasqueira •
-
SNP‐based enrichment analysis for WGCNA modules using an hypergeometric test
written by andres.firrincieli
-
Hypergeometric test in R using (phyper)
written by Rimma •
-
Is the bioconductor RRHO R package p-value computation for two.sided completely wrong ?
written by Anthony •
-
Advice for gene overlap with direction between two sets.
written by simplitia •
-
hypergeometric test on outlier SNPs
written by Ana
-
Fisher's exact test gives p-value 0
written by Adrian Pelin
-
Hypergeometric Test On Gene Set
written by ChIP
-
Probability of gene list overlap
written by Nasir •
Is there anyway to do this test for mutliple lists. Like lets say that I want to test whether gene lists for seven tests are more similar than expected by chance. I know that I can do 7+6+5...1 pairwise tests. But it would seem more elegant to do a global test for overlap.
What you're asking is unclear. What do you mean by "to test whether gene lists for seven tests are more similar than expected by chance.".
The test is for deciding how likely/unlikely the overlap between two sets is, there's no notion of similarity involved. Also consider posting as a new question providing a link to bioinformatics otherwise this is probably a purely statistical question best addressed on Cross Validated.