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

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

1
Entering edit mode

Thank you so much Kevin. This helps a lot.