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.
1 answer
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"))
•
written
by
asta.tranholm •
Log in to answer this question.
-
P-value Threshold Consideration in Multi-Sample RNA-Seq Experiment
written by Netanel •
-
DESeq2 extract desired result
written by beans •
-
Multi factor analysis in DESeq2 and design formula
written by nik.kraemer •
-
Interaction between the treatment and the drug
written by wmsalsah •
-
What will the the classic approach is to make pairwise comparisons between the groups in edgeR?
written by Info.shi •
-
Specifying unequal replicates in deseq2 with htseq-count input
written by kousi31 •
-
DESEQ2 application on control vs treatment
written by bioinfonerd •
-
Deseq2 Contrast question
written by disulfidebond •
-
DESeq2 model design: dose and time effect
written by CuriusScientist
-
RNA-Seq analysis with DESeq2. Confusing results.
written by Armin •
