DESeq2: design and/or results for 2 groups and 2 treatments
1
0
Entering edit mode
20 months ago

I have RNAseq data from 4 samples with 3 replicates each:

Treated A, Vehicle A, Treated B, Vehicle B,

I have defined sampleCondition as T or V and sampleGroup as A or B and created a sampleTable accordingly. I am interested in obtaining results for the differentially expressed genes due to treatment in group A and B (TA vs VA and TB vs VB), as well as the differentially expressed genes in the groups, treatment aside (VA vs VB and TA vs TB).

I'm unsure how to define design so I can obtain these comparisons? I suspect it may be design=~group + condition + condition:group?

Additionally, I'm unsure how I should write results() in order to obtain the specific data comparisons.

DESeq2 design results interactions contrast • 473 views
ADD COMMENT
1
Entering edit mode
19 months ago

For anyone looking for this answer: I decided to define them differently so they had unique identifiers (TA, VA, TS, VS). Therefore I could define the comparisons I wanted to obtain that way by listing design as design=~group and results as follows:

res_TVA<- results(dds, contrast=c("group","TA","VA")) 
res_TVS<- results(dds, contrast=c("group","TS","VS")) 
res_VAS<- results(dds, contrast=c("group","VA","VS")) 
res_TAS<- results(dds, contrast=c("group","TA","TS"))
ADD COMMENT

Login before adding your answer.

Traffic: 3112 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