How to do a hypergeometric test to identify the overall enrichment of an element containing genes among differentially expressed genes
0
0
Entering edit mode
6.9 years ago
pubudu • 0

I have a differntially expressed geneset. Some genes have an specific element in their promoters. I want to perform a hypergeometric test to see whether genes that have this element of their promoters are more likely to be differentially expressed than genes that do not have this element. I found that below one should be the relevent R function.

phyper(success-in-sample-1, success-in-bkgd, failure-in-bkgd, sample-size, lower.tail= FALSE);

What are the relevent parameters of my dataset. I can extract all the information from the dataset but dont understand what are they. For an example let's say I have 5000 differentially expressed genes. I found 4000 genes that have that element. So what are other parameters?

overrepresentation analysis hypergeometric test • 3.8k views
ADD COMMENT
0
Entering edit mode

Typing ?phyper in R will bring up the help information of that function.

ADD REPLY
0
Entering edit mode

Thank you. I know that. But I dont know how to relate those paramters in my data set. For an example let's say I have 5000 differentially expressed genes. I found 4000 genes that have that element. So what are other parameters?

ADD REPLY
1
Entering edit mode

The other parameters to a hypergeometric test usually involve background. In other words, you have a set of all genes in your genome/universe ("background"). You measure how many of all genes have some characteristic of interest ("successes" in background). Then you have a subset of genes: genes of interest ("sample"). You measure how many of those have the same characteristic of interest (successes in sample). Those are your four parameters.

ADD REPLY
0
Entering edit mode

To add to that: the background are all the genes you could detect with your technology. For expression study that would be all genes which are expressed (optionally above a certain cutoff).

ADD REPLY

Login before adding your answer.

Traffic: 2806 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6