We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
1 parent c26427e commit a9236c1Copy full SHA for a9236c1
1 file changed
R/Filter_QC.R
@@ -881,12 +881,15 @@ filterQC <- function(object,
881
## Annotate Doublets: ####
882
## Gene filter does not effect doublet ident and so not recalculated
883
884
+
885
if( do.doublets.fitler==T){
- sce <- as.SingleCellExperiment(so)
886
+ sce <- as.SingleCellExperiment(so)
887
888
set.seed(123)
- sce.dbl <- scDblFinder(sce)%>%suppressWarnings()
- sce.class <- sce.dbl$scDblFinder.class
889
- so <- AddMetaData(so,sce.class,"Doublet")
+ RNGkind("default", "default", "default")
890
+ sce.dbl <- scDblFinder(sce,BPPARAM=SerialParam(RNGseed = 123))%>%suppressWarnings()
891
+ sce.class <- sce.dbl$scDblFinder.class
892
+ so <- AddMetaData(so,sce.class,"Doublet")
893
}else{ print('doublets Identification not Run')}
894
return(so)
895
})
0 commit comments