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

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

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.