Entering edit mode
21 months ago
bioinformatics.girl
▴
20
I am integrating 5 different datasets for a scATAC-seq analysis.
I have had the following error persist for weeks; have looked online as well for possible solutions yet nothing has worked. Has anyone run into this issue?
The following code is:
object.list = list(sample1,sample2, sample3, sample4, sample5)
integration.anchors <- FindIntegrationAnchors(
object.list = object.list,
anchor.features = rownames(sample1),
reduction = "rlsi",
dims = 2:30
)
library(future)
plan("sequential")
options(future.globals.maxSize = 16000 * 1024^2)
integrated<-IntegrateData(
integration.anchors,
new.assay.name = "integrated",
normalization.method = c("LogNormalize", "SCT"),
dims = 1:27,
sd.weight = 1,
eps = 0,
verbose = TRUE
)