fix bug about reduction dimensions

stemangiola · Oct 28, 2023 · 078bef9 · 078bef9
1 parent 67b61d5
commit 078bef9
Show file tree

Hide file tree

Showing 4 changed files with 24 additions and 2 deletions.
diff --git a/R/methods.R b/R/methods.R
@@ -1140,6 +1140,17 @@ setGeneric("reduce_dimensions", function(.data,
 	col_names = get_abundance_norm_if_exists(.data, .abundance)
 	.abundance = col_names$.abundance
 
+	# adjust top for the max number of features I have
+	if(top > .data |> distinct(!!.feature) |> nrow()){
+	  warning(sprintf(
+	    "tidybulk says: the \"top\" argument %s is higher than the number of features %s", 
+	    top, 
+	    .data |> distinct(!!.feature) |> nrow()
+	  ))
+
+	  top = min(top, .data |> distinct(!!.feature) |> nrow())
+	}
+
 	# Validate data frame
 	if(do_validate()) {
 	validation(.data, !!.element, !!.feature, !!.abundance)

diff --git a/R/methods_SE.R b/R/methods_SE.R
@@ -467,6 +467,17 @@ setMethod("cluster_elements",
   if(.abundance |> quo_is_symbolic()) my_assay = quo_name(.abundance)
   else my_assay = get_assay_scaled_if_exists_SE(.data)
 
+  # adjust top for the max number of features I have
+  if(top > nrow(.data)){
+    warning(sprintf(
+      "tidybulk says: the \"top\" argument %s is higher than the number of features %s", 
+      top, 
+      nrow(.data)
+    ))
+
+    top = min(top, nrow(.data))
+  }
+
 	my_assay =
 		.data %>%
 

diff --git a/man/get_reduced_dimensions_UMAP_bulk.Rd b/man/get_reduced_dimensions_UMAP_bulk.Rd
diff --git a/man/get_reduced_dimensions_UMAP_bulk_SE.Rd b/man/get_reduced_dimensions_UMAP_bulk_SE.Rd