Question: RUVseq uses the the edgeR differential analysis
gravatar for tonja.r
3.7 years ago by
tonja.r460 wrote:

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) 

R • 1.4k views
ADD COMMENTlink modified 3.7 years ago by Devon Ryan90k • written 3.7 years ago by tonja.r460
gravatar for Devon Ryan
3.7 years ago by
Devon Ryan90k
Freiburg, Germany
Devon Ryan90k wrote:

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.

ADD COMMENTlink written 3.7 years ago by Devon Ryan90k
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.3.0
Traffic: 1645 users visited in the last hour