| Previous changeset 81:f86cb6b724ee (2014-02-26) Next changeset 83:3bdf17623d87 (2014-02-26) |
|
Commit message:
Uploaded |
|
modified:
rgedgeRpaired_nocamera.xml |
| b |
| diff -r f86cb6b724ee -r b21cd7bec41f rgedgeRpaired_nocamera.xml --- a/rgedgeRpaired_nocamera.xml Wed Feb 26 01:50:06 2014 -0500 +++ b/rgedgeRpaired_nocamera.xml Wed Feb 26 21:49:45 2014 -0500 |
| b |
| @@ -758,9 +758,10 @@ DGEList = estimateGLMTrendedDisp(DGEList,mydesign) DGEList = estimateGLMTagwiseDisp(DGEList,mydesign) } - norm.factor = calcNormFactors(DGEList) + calcNormFactors(DGEList) # -> object$samples$norm.factors + ls = colSums(DGEList$counts) * DGEList$samples$norm.factors pdf("VOOM_mean_variance_plot.pdf") - dat.voomed = voom(DGEList, mydesign, plot = TRUE, lib.size = colSums(workCM) * norm.factor) + dat.voomed = voom(DGEList, mydesign, plot = TRUE, lib.size = ls) dev.off() # Use limma to fit data fit = lmFit(dat.voomed, mydesign) |