RUVseq uses the the edgeR differential analysis
I am investigating the RUVseq right now and have noticed that they are using glm model of the edgeR to estimate the common and tag wise dispersion. From edgeR manual I read that one would use GLM model for general experiments with multiple factors. However, in RUVseq tutorial in the example they have only one factor control and treated. Why do they still use GLM model and not normal estimateCommonDisp() and estimateTagwiseDisp()?
design <- model.matrix(~x, data=pData(set))
y <- DGEList(counts=counts(set), group=x)
y <- calcNormFactors(y, method="upperquartile")
y <- estimateGLMCommonDisp(y, design)
y <- estimateGLMTagwiseDisp(y, design)
1 answer
They do this for simplicity so they don't need to implement both methods. The results should be similar enough that I can't really blame them.
• 0 views
•
link
written
by
Devon Ryan
Log in to answer this question.
More posts like this
-
RUVg correction
written by Shaimaa Gamal •
-
How to improve batch correction
written by Shaimaa Gamal •
-
Extract Deviance Residuals from DESeq2
written by chasem •
-
betweenLaneNormalization and calcNormFactors in RUVs
written by manish20072013 •
-
edgeR: no residual df: Cannot estimate dispersion
written by dec986
-
Which model to choose in edgeR analysis?
written by biolab
-
edgeR estamateDisp give other results when broken into two sub function
written by tonja.r
-
dispersion estimation from DESeq and edgerR
written by tonja.r
-
edgeR / deSeq parallel usage
written by russhh
-
edgeR: Correct pipeline for DE analysis with multiple conditions and batches
written by MT •
