Log fold change from rlog

Question

[Arindam Ghosh]

How to calculate the logged fold-change from rlog values (obtatined from DESeq2)? I have a list of genes with rlog values in two conditions (no other information). Is lfc = rlog_A / rlog_B or fc = rlog_A - rlog_B?

Answer 1

[swbarnes2]

But note that DESeq doesn't use the rlog values when computing fold changes. The rlog values are just nice for some kinds of visualizations.

[Arindam Ghosh]

Yes, I am aware of that. In my situation I have only the rlog values only, unfortunately. How do I calculate the fold-change from this information?

[ATpoint]

I would at least put them into limma then to get somewhat reliable stats rather than doing that by hand.