RNA-Seq differential gene expression workflow using edgeR
edgeR: Differential Expression AnalysisedgeR is a Bioconductor package for RNA-seq differential expression analysis. It models count data using the negative binomial distribution and uses empirical Bayes methods to estimate dispersion.
library(edgeR)
y <- DGEList(counts)
keep <- filterByExpr(y)
y <- y[keep,,keep.lib.sizes=FALSE]
y <- normLibSizes(y)
design <- model.matrix(~group)
y <- estimateDisp(y,design)
fit <- glmQLFit(y,design)
qlf <- glmQLFTest(fit)
Robinson MD, McCarthy DJ, Smyth GK (2010).
edgeR: a Bioconductor package for differential expression analysis of digital gene expression data.
Bioinformatics 26(1):139-140.
doi:10.1093/bioinformatics/btp616