DESeq2 Log2 fold change does not make sense given change in TPM
1
0
Entering edit mode
dadca596 • 0
@8a6484d7
Last seen 21 months ago
Australia

Suppose that for three replicates of treatment 1, the average TPM for a particular gene is 80 (and not much difference between the three replicates).

And for three replicates of treatment 2, the average TPM for the same gene is 990 (again, not much difference between replicates).

80 / 990 = 0.08 so doesn't that mean that the Log2 fold change should be approximately -3.64 in treatment 1? I am getting a Log2 fold change of only -0.75, which seems way way off.

What's going on here?

DESeq2 TPM Foldchange • 1.9k views
ADD COMMENT
0
Entering edit mode

Please add code and the actual counts (not the TPMs, the normalized counts from DESeq2) and a section of your results object for people to comment. Textual descriptions are hard to debug.

ADD REPLY
0
Entering edit mode
@mikelove
Last seen 1 hour ago
United States

I recommend plotCounts. It will help you to see the measurements used for LFC calculation.

ADD COMMENT

Login before adding your answer.

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