It looks like calcNormFactors wants a matrix, not a DGEList. Did I do something wrong? d=DGEList(counts=TotalReads, group=treatments,genes=rownames(TotalReads)) Calculating library sizes from column totals. dNorm=calcNormFactors(d) Error in calcNormFactors(d) : 'dataMatrix' needs to be a matrix Naomi