Question: design for DEseq2: LRT vs Wald & interactions
gravatar for cintapq
3 months ago by
cintapq10 wrote:

I would like advise about how to make a design for my experiment using DEseq2. In principle my question is straightforward, I'd simply like to find differentially expressed genes between two conditions. However, I've noticed that I have a large batch effect (I have two batches) and also large differences due to gender. Thus, I think I should include the two factors in my design:

> dds <- DESeqDataSetFromTximport(txi, sampleTable, ~BATCH + GENDER + Condition)
> resultsNames(dds)
[1] "Intercept"                                            "Batch_2_vs_1"
[3] "GENDER_MALE_vs_FEMALE"                                "Condition_T_vs_N"

Q1: if my question is which are the DE genes between conditions taking into account the batch and gender effects, which is the better option:

dds <-DESeq(dds)
res <- results(dds)
summary(res, name="Condition_T_vs_N")


dds_LRT <- DESeq(dds, test="LRT", reduced=~BATCH + GENDER)
res_LRT <- results(dds_LRT)

Q2: if I add an interaction term to the design, to which question am I answering?

dds <- DESeqDataSetFromTximport(txi, sampleTable, ~BATCH + GENDER + Condition:GENDER)

Many thanks for your comments, Cinta

deseq2 • 108 views
ADD COMMENTlink modified 3 months ago by Michael Love23k • written 3 months ago by cintapq10
Answer: design for DEseq2: LRT vs Wald & interactions
gravatar for Michael Love
3 months ago by
Michael Love23k
United States
Michael Love23k wrote:

The Wald test and the LRT should give similar (not exactly the same) results. Either is correct as you have it above.

The second design provides an baseline for sex and then two different condition effects, one for M and one for F.

ADD COMMENTlink written 3 months ago by Michael Love23k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 141 users visited in the last hour