DESeq2 factors and ANOVA
1
2
Entering edit mode
tony.fox2016 ▴ 40
@tonyfox2016-9796
Last seen 8.5 years ago

Hi,

I have a couple of questions about DESeq2:

1) in the DESeq2 workflow, to test the dex treatment effect, the following command was used:

dds <- DESeqDataSet(se, design = ~ cell + dex)

And I guess the only factor here is "dex" and it might be more straightforward to use:

dds <- DESeqDataSet(se, design = ~ dex)

2)If I am performing an ANOVA-like comparison of 3 groups A, B and C:

design(dds) = ~ group
dds = DESeq(dds, test = "LRT", reduced = ~ 1)
results(dds)

To extract the DE genes between B-vs-A and C-vs-B. We should use the following commands and no post-hoc tests are needed. Is it correct?

res.B.vs.A <-results(dds, contrast=c("group","A","B"))

res.C.vs.B <-results(dds, contrast=c("group","B","C"))

Many thanks for Michael Love and colleagues for contributing such an excellent software and their answers in the forum!

deseq2 • 6.8k views
ADD COMMENT
0
Entering edit mode
@mikelove
Last seen 2 days ago
United States

So hopefully this thread

C: DESeq2 ANODEV with 3 sample groups

answered this question

ADD COMMENT

Login before adding your answer.

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