Even if the batches were meaningfully distributed, OP doesn't need to use ComBat-seq for batch especially when they're trying to do DE - they need to use batch as covariate then. But like you say, the sample type and dataset number are confounded - there is no way to attribute any change observed to just biology or just batch.
How to adjust for multiple batches using Combat-Seq
Dear all,
I have RNA-seq data that I've collected from 4 datasets, 3 of which correspond to the same type of sample (control, green). The other dataset (4) contains only cancer samples (red). All these samples are integrated in the same matrix, and the goal is to perform DE analysis between control and cancer patients.
Therefore, I want to perform ComBat to adjust for batch effects. However, I have two batch effects: one is Dataset, and the other is Hospital (for dataset 1, samples were collected from two different hospitals adding another layer of variability).
How should I proceed? I've tried to find the answer looking in different forums, but nothing helps.
Thanks in advance!
• 1,699 views
•
link
updated
by
Ram
•
written
by
sabriledesma13 •
•
written
by
sabriledesma13 •
1 answer
Log in to answer this question.
More posts like this
-
combat_met error: mu_hat_mat[, batches_ind[[jj]]] : subscript out of bounds
written by ramiro.barrantes •
-
How to handle Single-Level Batch Variable for Batch Correction in RNA-Seq Data
written by Nelo •
-
Tutorial: removeBatchEffect explained using base R linear models
written by nhaus •
-
creating batch colum for batch correction
written by Expert •
-
ComBat-seq batch correction scRNA-seq
written by luciano.t.simonetta •
-
ComBat-seq problem
written by SKY •
-
Batch Adjusted counts for Downstream Analysis
written by saicharanp18 •
-
Can we use previous estimated ComBat model to correct batch effects for new data?
written by Min Dai
-
Help with multiple batch effects
written by fp89 •
-
EdgeR-GLMs correcting for batch effects
written by Matina
