betweenLaneNormalization and calcNormFactors in RUVs
0
0
Entering edit mode
6.8 years ago

I'm analyzing the RNAseq data using the RUVs tutorial. In this tutorial they use two functions,

uq <- betweenLaneNormalization(filtered, which = "upper") 
s <- RUVs(uq, 1:nrow(uq), k = 2, groups)
design <- model.matrix(~x + s$W)
y <- DGEList(counts = filtered, group = x)
y <- calcNormFactors(y, method = "upperquartile")
y <- estimateGLMCommonDisp(y, design)
y <- estimateGLMTagwiseDisp(y, design)
fit <- glmFit(y, design)

First question: Can I use the "median" instead of "upper" in betweenLaneNormalisation because it is giving me a RLE plot better than "upper"?

Second question: If I'm using "median" in above question then is it appropriate to calculate the normalization factor using the "upperquartile"method? If not then what is better option?

Thanks!

Differential-Gene-Expression RUVs RNA-Seq R • 2.1k views
ADD COMMENT

Login before adding your answer.

Traffic: 1705 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6