deseq: all padj values equal to 1
2
0
Entering edit mode
ccheung • 0
@ccheung-7248
Last seen 9.1 years ago
European Union

Hi Everybody,

I performed a deseq analysis on my samples (11 controls and 33 experimental samples). It is true that there are huge variations within each group, but although the pval are very low, the padj are always 1.  Plus when I put in the command:

resSig = res[ res$padj <.05,]

the list that it gives me includes all the transcripts although the padj are 1.

I would really appreciate any advice or help. Thanx in advance.

 

carol

deseq padj • 4.7k views
ADD COMMENT
0
Entering edit mode

Please post the complete code that you used including some output to show snippets of the result tables, and maybe also the output of DESeq2' diagnostic plots.

And, in case you were using DESeq rather than DESeq2, please consider starting over with DESeq2. DESeq is outdated.

ADD REPLY
1
Entering edit mode
@mikelove
Last seen 11 hours ago
United States

This is just a natural consequence of the Benjamini Hochberg correction. It means none of the p-values were small enough to survive multiple test correction. Here's an example:

> p.adjust(c(.8,.9,1), method="BH")

[1] 1 1 1

 

 

ADD COMMENT
0
Entering edit mode
@15150688610-7694
Last seen 9.6 years ago
China

Hello,I am encounting the same question with you.

Could you please tell that have you ever found out the reasons?

If so,could you please give me some advice.

ADD COMMENT

Login before adding your answer.

Traffic: 455 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