How to perform DESEQ2 analysis when there are three treatments (one of them being a reference)?
1
1
Entering edit mode
7 weeks ago
mrj ▴ 50

I am trying to analyze a gene expression dataset with three treatments (A, B, and control). There should not be any interaction effect. I will be thankful to you if can briefly explain how to obtain differentially expressed genes under these conditions.

How can I define my reference for the expression analysis when I have three treatments? How can I define my formula?

is it ~ A + B?

when I show over-expressed genes, should I have two lists? e.g.,

A vs Control B vs Control

comparisons DESEQ2 expression multiple differential analysis • 240 views
ADD COMMENT
1
Entering edit mode
7 weeks ago

Hi, you should have a single column in your metadata for condition, and this should contain A, B, and control.

For example:

metadata <- data.frame(
  condition = c(
    'A','A','A',
    'B','B','B',
    'control','control','control'))
metadata$condition <- factor(metadata$condition,
  levels = c('control','A','B'))

So, your design formula would be:

~ condition

Then, you can compare different levels like this: DESeq2 compare all levels

Kevin

ADD COMMENT
1
Entering edit mode

Thank you so much Kevin. This helps a lot.

ADD REPLY

Login before adding your answer.

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