This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Log fold change from rlog

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?

logarithms r deseq

1 answer

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.

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?

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

Log in to answer this question.