Controlling for batch and estrus cycle using DESeq2 in RNA-seq analysis

Hi,

I was wondering if I could get some advice on controlling for random effects using DESeq2. For example, given simulated data:

Treatment   Batch   Estrus
Control      1          1
Control      2          2
Control      3          1
Control      4          1
Control      1          2
Control      2          2
Control      3          2
Control      4          1
Control      1          3
Control      2      1
Treated      3          3
Treated      4          3
Treated      1          2
Treated      2          3
Treated      3          2
Treated      4          3
Treated      1          2
Treated      2          1
Treated      3          1
Treated      4          2

What is the best way to control for batch and estrus? I understand that DESeq2 does not support random effect analysis. As such, I was wondering if I could be advised on a work around for this type of analysis?

Thanks

random-effects rna-seq batch-effect

• 1,531 views • link

updated 2.3 years ago by Ram

• written 9.0 years ago by Dipro •

Hello Dipro!

It appears that your post has been cross-posted to another site: https://support.bioconductor.org/p/98425/ (and answered there)

This is typically not recommended as it runs the risk of annoying people in both communities.

• 0 views • link

written 9.0 years ago by WouterDeCoster

0 answers

No answers yet.

More posts like this

Tutorial: removeBatchEffect explained using base R linear models
written 2.3 years ago by nhaus •

Hi all, I always had a some difficulties to understand what removing the effect of a covariates actually means and what happens "under the hood". …
Batch Effect Removal on non-Linearly independent samples
written 3.1 years ago by James •

The wet lab that I work with did a bulk RNA-Seq experiment. In the experiment, they had wild type and diseased cells. They then treated …
Two rounds of batch effect removal
written 4.3 years ago by Kazuna •

Hi, we are currently doing RNAseq analysis and would like to remove batch effects. I would like to compare my data with previously published data, …
rna-seq/DEseq2 analysis - is it a linear combination problem?
written 5.4 years ago by aabbcc1 •

Hi, I need some advice concerning my data analysis. It's RNA-seq data looking at gene expression differences with a treatment. Also recorded is the sampling …
limma::removeBatchEffect and using DESeq2 for RNA-seq differential expression
written 6.5 years ago by 2405592M

Hi everyone, I have an RNA-seq experiment that I am analysing. I received cell culture replicates in batches. I used the `limma::removeBatchEffect` to remove any …
Compare different library
written 6.8 years ago by inminin •

Hello, I have total RNA sequencing data (treatment 4, control 5), but we used different library sequencing for one treatment sample. Can I do DESeq2 …
DESeq2 design with multiple conditions
written 6.9 years ago by maria2019

Hi, I have 16 samples and multiple conditions (treatment and sex): (4 Untreated_male, 4 Untreated_female, 4 Treated_male, 4 Treated_female) I want to do 3 types …
Using DESeq2 to analyze similarities in treatment effect - 2 conditions and 4 genotypes
written 7.1 years ago by martin.strand •

I have a mRNASec dataset with 2 conditions and 4 genotypes and I'm interested in finding genes that have a similar change in expression across …
RNA-seq batch effect correction
written 9.2 years ago by endikavarela •

Hi, I'm new to RNA-seq data analysis and i need help to remove some batch effect from my RNA-seq data. First, my data structure is …
Differential Expression with Factorial Design and adjust for Lane+Batch Effect
written 10.7 years ago by Sam

Hi, so recently we have got a data back with the following design ``` Condition Treatment Lane Batch Case Treated 1 1 Case Treated 2 …