Fold Change and statistical significance
1
0
Entering edit mode
g.atla ▴ 10
@gatla-9491
Last seen 4.6 years ago

I am using DESeq2 to find the DE genes in a data set with 21 replicates ( 21 + 21, paired-design). I get the DE genes but with very small log2 fold changes (less than +/- 0.5) , as you can see in the plot below:

> summary(resdds, 0.05)
out of 22025 with nonzero total read count
adjusted p-value < 0.05
LFC > 0 (up)     : 2479, 11%
LFC < 0 (down)   : 3127, 14%
outliers [1]     : 0, 0%
low counts [2]   : 854, 3.9%
(mean count < 4)

Unshrunken LFC:

I could see that the mean expression is good and the sample size is large. I would like to know if the results makes sense and are trustable.

P.S: The GO analysis looks good as well.

deseq2 logfoldchange • 1.5k views
ADD COMMENT
2
Entering edit mode
@mikelove
Last seen 23 hours ago
United States

"I would like to know if the results makes sense and are trustable."

Yes they are. This is just a consequence of having true differences and sufficient read depth and sample size to detect them.

We discuss this in the section "Hypothesis tests with thresholds on effect size" of the DESeq2 paper.

ADD COMMENT

Login before adding your answer.

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