This is a test version of Biostars. For the public version, visit https://www.biostars.org.
edgeR comparison of 2 samples

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:)

edger rna-seq

1 answer

A stated in the manual,

exactTest(dge, pair=c("wt", "res")) 

will find genes differentially expressed (DE) in "res" vs "wt".

Log in to answer this question.