Question: What doest this p.value distribution mean?
0
gravatar for chipolino
11 months ago by
chipolino0
chipolino0 wrote:

Hi everyone,

I am doing an RNA-seq analysis with DESeq2. I have 3 patients with disease (affected) and 3 healthy patients (unaffected). After following standard DESeq2 pipeline (as described here), I got my table with results and p.value column. I decided to check the distribution of p.values and here it is:

I do know that we expect to see something similar to uniform distribution, apparently we don't see it here... Does it mean that there is a problem with the data or an experiment design (I used simple formula, which was taking into account genders of the patients: ~ Gender + Condition)?

Thank you

deseq2 pvalue • 359 views
ADD COMMENTlink modified 11 months ago by Michael Love23k • written 11 months ago by chipolino0
1

Can you show a PCA plot? This sometimes comes from batch effects, or other related correlations or differences in the variance than from the model (e.g. if one group had much larger dispersion, or an outlier).

ADD REPLYlink written 11 months ago by Michael Love23k

PCA plot is not great either... Zeros are unaffected patients, Ones are affected. PC2 divides them by gender.

ADD REPLYlink written 11 months ago by chipolino0
Answer: What doest this p.value distribution mean?
3
gravatar for Michael Love
11 months ago by
Michael Love23k
United States
Michael Love23k wrote:

I guess the non-uniformity could have to do with the partial confounding of sex and condition. This seems like a good case where you could try to alleviate the distribution using fdrtool using Bernd Klaus' tutorial:

http://www-huber.embl.de/users/klaus/Teaching/DESeq2Predoc2014.html#inspection-and-correction-of-pvalues

ADD COMMENTlink written 11 months ago by Michael Love23k

Thank you, I will check it

ADD REPLYlink written 11 months ago by chipolino0
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 16.09
Traffic: 377 users visited in the last hour