Entering edit mode
2.6 years ago
biology_inform
▴
50
Hi
data <- readDGE(files)
data <- head(data, -5) group <- c(rep("wt", 3),rep("res", 3))
dge = DGEList(counts=data, group=group)
dge <-calcNormFactors(dge)
dge <- estimateDisp(dge)
et <- exactTest(dge, pair=c("wt", "res"))
write.csv(et$table,"dacomitinib-edgeR.csv")
I am analyisng my sample rnaseq data with EdgeR and with this code above; I wonder that am I comparing wt vs res or res vs wt? Where can I understand this thanks:)