[BioC] DESeq / DESeq2 for translation efficiency in polysome profiling or ribsome profiling experiments
gilgi.friedlander at weizmann.ac.il
Tue Dec 24 15:57:44 CET 2013
Thanks a lot!
From: Simon Anders [mailto:anders at embl.de]
Sent: Tuesday, December 24, 2013 4:11 PM
To: Gilgi Friedlander; 'Michael Love'
Cc: Gilgi [guest]; bioconductor at r-project.org; DESeq Maintainer
Subject: Re: [BioC] DESeq / DESeq2 for translation efficiency in polysome profiling or ribsome profiling experiments
On 24/12/13 15:08, Simon Anders wrote:
> For the interaction, you get not only a p value but also a log2 fold
> change in the results data frame. The latter is an estimate of the
> log2 of the double ratio you mention. Using this estimate is
> preferrable to simply using the ratio itself because it has undergone
> empirical-Bayes shrinkage, which gets rid of the exaggerated values
> that the naive ratio tends to exhibit when counts are low.
Just noticed that you were using DESeq1, not DESeq2. Shrinkage is offered only in DESeq2 (and it might be useful for you to switch for this reason). In DESeq1, the reported log fold change is the unshrunken GLM estimate, which, in case of a balanced design, tends to be pretty close to the actual log ratio.
More information about the Bioconductor