Question: Bias in hypergeometric test
0
gravatar for kakukeshi
2.0 years ago by
kakukeshi50
kakukeshi50 wrote:

Dear colleagues,

I calculated the significance of the overlap between drug-targeted genes and some disease genes using a hypergeometric test. The problem is that I noticed that the drugs with a lower p-value are the ones that are targeting more genes in general. What does this mean? Is my result is biased by the current knowledge about drug-target interactions? there's any way to normalize for this?.

P.S I'm using STITCH for the drug-target information

Thanks

ADD COMMENTlink modified 2.0 years ago by ivivek_ngs4.8k • written 2.0 years ago by kakukeshi50
1
gravatar for Devon Ryan
2.0 years ago by
Devon Ryan91k
Freiburg, Germany
Devon Ryan91k wrote:

You've chosen a test that works by checking the number of genes changed in a group. That a drug that changes more things is associated with more group enrichments is completely expected. Whether this is biologically correct is a different question and not one that can be answered with simple statistics (rather, you'd need to have a good understanding about how changes in protein level will affect activity elsewhere in the pathway that isn't evidenced by protein-level changes).

ADD COMMENTlink written 2.0 years ago by Devon Ryan91k
1
gravatar for ivivek_ngs
2.0 years ago by
ivivek_ngs4.8k
Seattle,WA, USA
ivivek_ngs4.8k wrote:

The drug-target gene interactions cannot be done based on the hypergeometric and this cannot be biologically correlated since such interactions depend on many factors which will be variable attributes that will be evidence for the interaction as also stated by @Devon Ryan. There are tools which do this specifically like Open-targets or DGIdb. Take a look at them for your results and the have combinatorial p-value and association score statistics. They also employ evidence from various NGS datasets and pretty well also use ExAC database and GWAS catalogue for associativity to get druggable hit in a genome. This should be more insightful for your data.

ADD COMMENTlink written 2.0 years ago by ivivek_ngs4.8k
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: 1387 users visited in the last hour