构建edgeR的DGEList对象,并归一化,拟合模型 dge <- DGEList(counts=PCVSCP_exprSet) dge <- calcNormFactors(dge) v <- voom(dge,design,plot=TRUE, normalize="quantile") fit <- lmFit(v, design) fit2 <- contrasts.fit(fit,cont.matrix) fit2 <- eBayes(fit2) tmp <- topTable(fit2, coef=...