handling no-count-peaks

immunogenomics · Aug 14, 2024 · 2cdf42f · 2cdf42f
1 parent 734f864
commit 2cdf42f
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/R/SCENTfunctions.R b/R/SCENTfunctions.R
@@ -166,7 +166,9 @@ SCENT_algorithm <- function(object, celltype, ncores, regr = "poisson", bin = TR
 
     #binarize peaks:
     if(bin){
-      atac_target[atac_target$atac>0,]$atac<-1
+      if(nrow(atac_target[atac_target$atac>0,])>0){
+        atac_target[atac_target$atac>0,]$atac<-1
+      }
     }
 
     mrna_target <- object@rna[gene,]