DESEQ2 Results log2fold and pvalue changes after removing genes
1
0
Entering edit mode
Bine ▴ 40
@bine-23912
Last seen 3 months ago
UK

Dear all,

we have been discussing internally but couldnt find an answer, maybe you could help us. The following happened: I run my DESEQ2 analysis with the full dataset and got genes with a certain log2fold change and pvalue via results(dds0). I then run the analysis again after removing some genes. We expected that the log2fold change and pvalue would not change and that only the padj value would change. But all values changed.. Can anyone explain why this is the case?

Thank you, Bine

DESeq2 • 855 views
ADD COMMENT
1
Entering edit mode

I suppose it depends how many genes you removed, but removing genes early in the process will change the library normalization and dispersion estimates.

ADD REPLY
0
Entering edit mode

Thank you for clarifying.

ADD REPLY
2
Entering edit mode
@mikelove
Last seen 13 hours ago
United States

Swbarnes is correct. The pvalue depends on the dispersion which borrows information across genes. Also the size factors could differ.

ADD COMMENT
0
Entering edit mode

Thank you!

ADD REPLY

Login before adding your answer.

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