DESeq2 design formula for time - condition

DESeq2 design formula for time - condition - donor

Hi, I need advice on how to set up the formula for the experimental design in DESeq2 in a complex analysis in which I have three DONORS and for each I have three TIMES and for each time I have an untreated control and a treated: coldata structure

The goal of the study is to find genes whose expression varies according to condition and time and independently of the donor. Please I would really appreciate your help. THANK YOU!

Ale

deseq2

• 866 views • link

updated 2.9 years ago by jared.andrews07

• written 2.9 years ago by asilvestris •

1 answer

I believe your design would be:

dds <- DESeqDataSet(dds, design = ~donor + condition + time + condition:time)
dds <- DESeq(dds, test = "LRT", reduced = ~Culture + H3_Status + DaysDiff)

Alternatively, you could collapse condition & time into a single column do pair-wise comparisons at each timepoint.

dds$timecondition <- paste0(dds$time, dds$condition)
dds <- DESeqDataSet(dds, design = ~donor + timecondition)
dds <- DESeq(dds)

• 0 views • link

written 2.9 years ago by jared.andrews07

More posts like this

Perform DESeq2 with One factor with three levels vs Two factor with two levels
written 1.8 years ago by وفاء •

Hi everyone I have patients RNA-seq data before and after treatment for one type of cancer. My goal is to use machine learning to predict …
design formula DEG paired design
written 3.3 years ago by Caro •

``` Hi all I have a single cell data experimental design which looks like this: cluster_metadata Sample_N Condition Treatment Sample_ID Sample_Name 1 healthy untreated S1 …
Multi factor analysis in DESeq2 and design formula
written 3.8 years ago by nik.kraemer •

Hi all, I am trying to do a multifactor analysis in DESeq2 with RNA-Seq data from mammalian cells. I have 12 different samples: name,cell,treatment Sample …
Alternate DESeq2 designs to model donor variance
written 3.9 years ago by swbarnes2

I have a set of RNASeq data; 3 treated replicates, 3 control replicates, from each of 4 donors. I observe in the PCA plot that …
Interaction between the treatment and the drug
written 4.5 years ago by wmsalsah •

My main factor is treatment (before and after treatment). I would like to know the gene differentially expressed because of the treatment, and to know …
DESeq2 design
written 5.0 years ago by Rajesh •

Hi Experts, i have this feature counts matrix for 23 samples (each containing multiple replicates ranging 3-16). Among these 23 samples some are untreated for …
DESeq2 design question
written 5.0 years ago by Lepomis_8 •

I have a count matrix from an RNA-seq experiment that I'd like to normalize using DESeq2 and perform DE analysis on. My code is below: …
DESeq2 and integers values
written 7.4 years ago by Morris_Chair

Dear all, I have an experiment with two different conditions, two treated and two untreated samples from where I performed RNA sequencing. here is the …
DESeq2 model matrix with interaction term
written 10.3 years ago by george.bask •

Hi all, I am trying to figure out the best way to make a model matrix for the following RNA-seq experiment using DESeq2. The experimental …
DESeq2 model design: dose and time effect
written 11.0 years ago by CuriusScientist

I have three samples each with 3 biological replicates i.e. one control and other at dose 12ug/ml and 25ug/ml at 24hr ``` > library('DESeq2') > …