Hi,
I am puzzled about this situation: if we have estimated a ComBat model for batch effects in dataset A, can we use this model to adjust batch effect for dataset B?
In another case, there are two batches, i.e., batch 1 and batch 2, in your dataset, but you just want to use 70% data of high quality to estimate and remove batch effects for batch 2 (if we set batch 1 as the reference batch). After estimating the model, we also need to adjust batch effects for the rest 30% data of batch 2, how can we do this?
Thanks in advance,
BP
1 answer
I would question this approach, as it assumes that the batch effects are the same in the new data - this is a big assumption to make, though. If at all possible, just include batch as a covariate in your design formula(e) while making statistical inferences.
If you must directly remove batch somewhere, then use removeBatchEffect() from limma.
Log in to answer this question.
-
combat_met error: mu_hat_mat[, batches_ind[[jj]]] : subscript out of bounds
written by ramiro.barrantes •
-
ComBat-seq batch correction scRNA-seq
written by luciano.t.simonetta •
-
ComBat-seq adjustment of batch effects
written by chatterjee.arumoy •
-
sva combat removes wanted biological variation
written by exin
-
Batch effects : ComBat or removebatcheffects (limma package) ?
written by lessismore
-
SVA : Setting up the data from an ExpressionSet and applying sva function to extimate batches
written by lessismore
-
How to use Combat in R??
written by madkitty
-
combat correction or include batch effect in the model?
written by Inbalor •
-
ComBat (SVA) Application - Two Technical Effects
written by andrew.j.skelton73
-
How to adjust data for "continuous" covariates ?
written by housebie •
