This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Hypergeometric Test On Gene Set

Hi!

This is not the first time this question is being asked, but I am confused from the previous post.

I have say two lists. List1 has 598 genes and List2 has 5500 genes and the total genes available in the pool from which these two are drawn is of size 23000 (say).

Now, if I have to make compute whether the overlap between the two list which is of 89 genes is significant or not.

I have two formulas:

method 1

phyper=(overlap-1,list1,PopSize-list1,list2,lower.tail = FALSE, log.p = FALSE)
phyper=(88,598,23000-598,5500,lower.tail = FALSE, log.p = FALSE)

method 2

phyper=(overlap,list1,PopSize,list2,lower.tail = FALSE, log.p = FALSE)
phyper=(89,598,23000,5500,lower.tail = FALSE, log.p = FALSE)

Now which method shall I use and why?

I am really confused.

Thank you

r statistics
• 27,415 views • link
updated by Ram Ram • written by ChIP ChIP

2 answers

Your method1 looks like the correct one. AFAIK, in phyper=(q,m,n,k)

n should be PopSize-list1

you can check this stackoverflow thread as well

• 236 views • link
written by Sudeep Sudeep

Can anyone explain about the q-1, why or why not?

• 0 views • link
written by Madelaine Gogol Madelaine Gogol

answer is here http://stats.stackexchange.com/questions/16247/calculating-the-probability-of-gene-list-overlap-between-an-rna-seq-and-a-chip-c

phyper(x, m, n, k) gives the probability of getting x or fewer, so phyper(x, m, n, k) is the same as sum(dhyper(0:x, m, n, k)).

The lower.tail=FALSE is a bit confusing. phyper(x, m, n, k, lower.tail=FALSE) is the same as 1-phyper(x, m, n, k), and so is the probability of x+1 or more

• 358 views • link
written by arno.guille arno.guille

I think method 1 is the correct one, because.

Method 1 gives the same result as this site: https://www.geneprof.org/GeneProf/tools/hypergeometric.jsp

Method 2 gives a different result.

• 1 views • link
written by Alejandro Jimenez Sanchez Alejandro Jimenez Sanchez

Log in to answer this question.

More posts like this