When subsetting an
EList object as output by
voomWithQualityWeights() there is a list element
When subsetting using the single bracket method, this element returns the initial vector, instead of one that matches the subset.
This doesn't seem to be the correct behaviour.
Is this able to be fixed?
y <- matrix(rnbinom(10000,mu=5,size=2),ncol=4) d <- DGEList(counts=y, group=rep(1:2,each=2)) v <- voomWithQualityWeights(d) v1 <- v[,1:2] length(v$sample.weights) == ncol(v) length(v1$sample.weights) == ncol(v1)