blob: fcbf6e1ee3aceca94c35c5e1fd756e8ab2699e5d (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
|
library(edgeR)
# It's easier to read environment variables than parse command line arguments
counts <- read.delim(Sys.getenv("counts"), row.names=1)
out <- Sys.getenv("out")
group <- rep(c("A", "B"), each = 2)
design <- model.matrix(~group)
dge <- DGEList(counts=counts)
keep <- filterByExpr(dge, design)
dge <- dge[keep,,keep.lib.sizes=FALSE]
dge <- calcNormFactors(dge)
logCPM <- cpm(dge, log=TRUE, prior.count=3)
fit <- lmFit(logCPM, design)
fit <- eBayes(fit, trend=TRUE)
write.table(topTable(fit, coef=ncol(design)), file = out)
|
