## create the object
dgeList <- DGEList(counts,group)
## calculate the sie factors
dgeList <- calcNormFactors(dgeList)
## plot them
apply(combn(rownames(dgeList$samples),2),
2,
function(co,obj){plotNormalizationFactors(obj,co[1],co[2])},dgeList)
Run the code above in your browser using DataLab