sigh of log2FC values in DESeq2
1
0
Entering edit mode
yampolsk • 0
@194b0276
Last seen 4 weeks ago
United States

Can someone explain to me how is the sign of log2FoldChange is set in the results of DEseq2? I was pretty sure that it is calculated as log2(Counts_treatment/Count_reference), where reference is determined alphabetically (unless specified by user).

It does appear to be the case - in most, but not all cases. Sometime the sign is the opposite to what you would expect by just looking at the read counts data. Moreover, sometimes the sign changes depending on whether shrinkage is use or not and logFC vs. LRT tests.

So it seems that under certain circumstances either the estimated fold change can switch signs (seems unlikely?) or the order of what is reference and what is treatment is reversed.

Could someone comment on this?

Thanks!

sign FoldChange DESeq2 reference • 185 views
ADD COMMENT
0
Entering edit mode

Well, no, the software isn't going to randomly flip flop. Though maybe shrinkage would change a fold change that was marginally positive to marginally negative. It would probably help if you provided normalized counts and colData for a few genes where you think the software has erred.

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

Agree with swbarnes comment. No the sign doesn’t randomly flip. Just look at the sign of significant genes using plotCounts.

ADD COMMENT

Login before adding your answer.

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