Skip to content

Commit ae89008

Browse files
committed
Add check to MergeSeurat for single-cell objects
1 parent 6c9b31a commit ae89008

File tree

1 file changed

+10
-0
lines changed

1 file changed

+10
-0
lines changed

singlecell/resources/chunks/MergeSeurat.R

Lines changed: 10 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -15,6 +15,11 @@ mergeBatch <- function(dat) {
1515
} else {
1616
toMerge[[datasetId]] <- readSeuratRDS(dat[[datasetId]])
1717
}
18+
19+
if (ncol(toMerge[[datasetId]]) == 1) {
20+
logger::log_info(paste0('Dataset has single cell, skipping: ', datasetId))
21+
toMerge[[datasetId]] <- NULL
22+
}
1823
}
1924

2025
if (!is.null(assaysToDrop)) {
@@ -28,6 +33,10 @@ mergeBatch <- function(dat) {
2833
}
2934
}
3035

36+
if (length(toMerge) == 0) {
37+
stop('There were no passing seurat objects!')
38+
}
39+
3140
seuratObj <- CellMembrane::MergeSeuratObjs(toMerge, projectName = projectName, doGC = doDiet, errorOnBarcodeSuffix = errorOnBarcodeSuffix)
3241
return(seuratObj)
3342
}
@@ -61,6 +70,7 @@ if (length(seuratObjects) == 1) {
6170
unlink(mergedObjectFiles[[1]])
6271
} else {
6372
logger::log_info('performing final merge')
73+
# TODO: check for single cell in object
6474
seuratObj <- readRDS(mergedObjectFiles[[1]])
6575
unlink(mergedObjectFiles[[1]])
6676

0 commit comments

Comments
 (0)