# HG changeset patch # User fubar # Date 1393474952 18000 # Node ID 3bdf17623d87aacb5bd89be9eabd84543025b927 # Parent b21cd7bec41fc0852244c5a09ab91ed768edd36a Uploaded diff -r b21cd7bec41f -r 3bdf17623d87 rgedgeRpaired_nocamera.xml --- a/rgedgeRpaired_nocamera.xml Wed Feb 26 21:49:45 2014 -0500 +++ b/rgedgeRpaired_nocamera.xml Wed Feb 26 23:22:32 2014 -0500 @@ -758,8 +758,8 @@ DGEList = estimateGLMTrendedDisp(DGEList,mydesign) DGEList = estimateGLMTagwiseDisp(DGEList,mydesign) } - calcNormFactors(DGEList) # -> object$samples$norm.factors - ls = colSums(DGEList$counts) * DGEList$samples$norm.factors + calcNormFactors(DGEList) + ls = colSums(DGEList\$counts) * DGEList\$samples\$norm.factors pdf("VOOM_mean_variance_plot.pdf") dat.voomed = voom(DGEList, mydesign, plot = TRUE, lib.size = ls) dev.off()