Question: DESeq2 with 2 factors: problem to fit the negative binomial and weird p-value distributions
gravatar for Aurelie MLB
5.4 years ago by
Aurelie MLB340
United Kingdom
Aurelie MLB340 wrote:


I have RNA-seq data from 27 samples to analyse. I have 2 differents treatments (Compound1 and compound 2) + control (DMSO) and I have 3 times. I have 3 replicates for each combination.

I am trying to use DESeq2 to find the differentially expressed genes. First, may I double check with you that the formula looks OK: ~treatment*time ?? (I have fount ~cell + dex on and I am confused. I am wondering why you would use this formula. Is there something special I should now about DESeq2 formulas?)

So, using a pooled analysis with all the samples at once and the design ~treatment*time:

- I tried to plot the dispersion versus the mean normalized counts (using plotDispEsts(dds)). I am not sure what is the red line in this plot because I thought the parameters were calculated per gene using the replicates... But for the low mean normalized counts (where the dispersion is bigger), this red line is well above the black and shrinked blue points. Is this the sign of a problem?

-I also looked at the p-value plots for different comparisons and some of them are not flat at all. The pvalues seem to slowly increase for the low p-values to the high pvalues with even some spike at 0.4 or 0.8. Is this classical for DESeq2 or is this the sign of a problem?

I thought that there weird observations were coming from a high variability between times so I tried to do a similar analysis but per time (formula: ~treatment ). My problem is that doing this I obtain very different lists of differentially expressed genes...

I am not sure which one is the best analysis really... any suggestion please?

Many thanks



rna-seq deseq2 • 2.5k views
ADD COMMENTlink modified 5.4 years ago by Michael Love2.1k • written 5.4 years ago by Aurelie MLB340
gravatar for Michael Love
5.4 years ago by
Michael Love2.1k
United States
Michael Love2.1k wrote:
We discourage cross posting the same question on multiple sites because it duplicates everyone's effort in answering your questions. At the least, please link to the Bioc support site posts.
ADD COMMENTlink written 5.4 years ago by Michael Love2.1k

Apologies! I did not know. I posted here first and then saw that the Bioconductor support website was recommended in your documentation so I thought it would be more appropriate to post there finally. All I can say now is that it will not happen again...

The bioconductor post is there with your answer is there:

ADD REPLYlink written 5.4 years ago by Aurelie MLB340
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1785 users visited in the last hour