Hello Everyone
I have Bulk-RNASeq sequencing data and processing the data in the DESeq2 package. I am concerned about the "batch" effect due to the different sequencing runs performed at different location. Samples of group1 and group2 were sequenced separably at different locations.
How to handle
condition Batch
Sample1 G1 1
Sample2 G1 1
Sample3 G2 2
Sample4 G2 2
dds<-DESeqDataSetFromMatrix(countData = count_f, colData = col_data, design = ~Batch+condition)
Should I include "Batch" or not while performing the differential analysis among two groups (G1 vs G2)?
I would appreciate all the suggestions.
•
written
by
Omics data mining
1 answer
It's perfectly confounded (batch is the exact same as condition), so nothing you can do about it. It's a (sorry to be blunt) poor experimental design. All you can do is to hope that the batch effect (if it exists) does not overshadow the true biological variation. Not nesting batch and condition is really the basis for any experimental design, you should try to avoid that in the future.
Log in to answer this question.
-
pseudobulk DESeq2 design matrix
written by HK •
-
Statistical analysisi in DESeq2 with different condition
written by michelafrancesconi9 •
-
Time-series RNA-seq analysis using DESeq2
written by Debebe •
-
Metagenomic differential abundance analysis
written by imaparna27 •
-
How to create shiny app for deseq2 object creation
written by majeedmj.ict •
-
Batch correction in DESeq2---Normalized Counts
written by mtsompana •
-
How to perform batch correction when we have only single condition?
written by svp
-
Deseq2 Differentially Expressed Genes between conditions
written by bryce.kirby •
-
FeatureCounts from Galaxy to DESeq2 in RStudio
written by Newbieish •
-
Can't group samples in StringTie for differential expression
written by CandiceChuDVM
