Batch correction of ADT data

[santiiiiiii]

Hi, I have used DSB (https://github.com/niaid/dsb) to normalise my ADT data. I however have a batch issue and would like to use rpca + integration anchor to correct ADT batch effects. I tried running the following:

subsample_list <- SplitObject(combined, split.by = "Batch")
subsample_list <- lapply(X = subsample_list, FUN = function(x) {
x <- FindVariableFeatures(x, selection.method = "vst")
})

I end up with the following error:

Error in simpleLoess(y, x, w, span, degree = degree, parametric = parametric, :
invalid 'x'