Adjusting DESeq2 p-values through IHW using pathway information as covariate
Entering edit mode
Last seen 2.9 years ago
University of Copenhagen

Greetings everyone,

I have been reading about IHW as a method to adjust p-values for multiple testing by assigning weights to each hypothesis based on a prior (informative) knowledge. My rationale here is that given a list of differentially expressed genes (as returned by DESeq2, but also other methods), I would like to adjust their p-values while accounting for the prior knowledge of which pathway the genes belong to. The idea behind this is that I would expect genes belonging to the same pathways to be more or less co-regulated.

That said, I have immediately faced the harsh truth that most genes are annotated for more then one and sometimes dozens of pathways and so the only way I could come up to pass this information as a vector-like object is to generate a binary string as long as the total number of observed pathways for my list of genes, where each position is either a zero (not belonging) or a one (the gene belongs to that pathway).

I am currently trying this out. Does anyone have a better idea?

Thanks in advance!

IHW deseq2 Pathway • 685 views
Entering edit mode
Entering edit mode
Last seen 14 hours ago
United States

This is novel territory, but you may want to see what the IHW developers have to say. I was thinking earlier that you had a particular pathway of interest.

One way to deal with overlapping categories is if you can order them. This happens in genome annotation, for example, where a region can fall under multiple categories, e.g. intron of one gene vs promoter of another. If you order the categories, you can then come up with a partition.


Login before adding your answer.

Traffic: 867 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6