Question: EdgeR GLM following voom conversion of counts
gravatar for adam.p.taranto
2.7 years ago by
adam.p.taranto0 wrote:

Documentation for limma-voom analysis of RNA-seq data uses eBayes to perform pairwise differential expression tests. For more complex experimental designs, would it be reasonable to use edgeR's GLM method to call DE genes from dge objects containing voom transformed counts?

ADD COMMENTlink modified 2.7 years ago by Ryan C. Thompson7.0k • written 2.7 years ago by adam.p.taranto0
gravatar for Ryan C. Thompson
2.7 years ago by
The Scripps Research Institute, La Jolla, CA
Ryan C. Thompson7.0k wrote:

LImma can already accomadate any design matrix that edgeR can. It definitely doesn't make sense to use edgeR on voom-transformed counts, because edgeR models values on a raw count scale, and voom performs a transformation to logCPM.

ADD COMMENTlink written 2.7 years ago by Ryan C. Thompson7.0k

Indeed, voom will happily generate negative log-expression values when the counts are low relative to the library size. It would be completely nonsensical to try to treat such values as counts in an edgeR analysis.

ADD REPLYlink written 2.7 years ago by Aaron Lun21k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 233 users visited in the last hour