DEseq2: comparison the same when reversing factor order?
1
0
Entering edit mode
3.5 years ago
gt ▴ 30

Say I want to compare a treatment group to a control group. Will treatment vs. control give me the same results as control vs. treatment? Just a slightly different interpretation? I know that log(A/B) = -1*log(B/A), so I would think so. Is this true for any method I use that computes a log fold change? What if I was using FCROS or edgeR instead?

DESeq2 FCROS RNA-seq edgeR • 1.3k views
ADD COMMENT
0
Entering edit mode
3.5 years ago

If you specify with contrast what to compare to what, the factor order doesn't matter, just the order in the contrast.

ADD COMMENT
0
Entering edit mode

I think I wasn't clear. If I do test vs. control in the contrast and get log fold change values for the genes. Then I do control vs. test in the contrast and get log fold change values again. Are the log fold change values for the two separate contrasts going to be the same? As in... log2(treatment/control) = -1*log2(control/treatment)?

ADD REPLY
0
Entering edit mode

Well, why don't you try it and look for yourself? With DESeq, I think the LFCShrink operation will slightly change the results depending on the chosen control.

ADD REPLY
0
Entering edit mode

Yep, will do. Was just wondering if this would be true for any package that calculates a log fold change... such as edgeR and FCROS as well.

ADD REPLY

Login before adding your answer.

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