fix subset_ccs

cole-trapnell-lab · Jan 5, 2024 · 4c358b8 · 4c358b8
1 parent 0f1b96f
commit 4c358b8
Showing 1 changed file with 4 additions and 2 deletions.
diff --git a/R/methods-cell_count_model.R b/R/methods-cell_count_model.R
@@ -97,15 +97,17 @@ subset_ccs = function(ccs, ...) {
   group_ids = filtered_cds_coldata %>% pull(group_id)
 
   cell_group_assignments = ccs@metadata$cell_group_assignments[ccs@metadata$cell_group_assignments$group_id %in% group_ids,]
-  samples = cell_group_assignments$sample %>% unique()
-  cell_groups = cell_group_assignments$cell_group %>% unique()
+  samples = intersect( colnames(ccs), unique(cell_group_assignments$sample))
+  cell_groups = intersect(rownames(ccs), unique(cell_group_assignments$cell_group))
   sub_ccs = ccs[cell_groups, samples]
   sub_ccs@metadata$cell_group_assignments = cell_group_assignments
 
   return(sub_ccs)
 
 }
 
+
+
 #' #' subset ccs by cell groups
 #' #' @param ccs
 #' #' @param cell_groups