Two small corrections: Treated <- rep(c(1,-1),5) design <- cbind(Dye=1,Treated=Treated) fit <- lmFit(MA, design) fit <- eBayes(fit) topTable(fit, coef="Treated") Gordon