edgeR comparison of 2 samples
1
0
Entering edit mode
2.0 years ago

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 • 553 views
ADD COMMENT
0
Entering edit mode
2.0 years ago

A stated in the manual,

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

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

ADD COMMENT

Login before adding your answer.

Traffic: 2399 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6