Added a fix to do_GroupwiseDEHeatmap where pvalues of 0 would cause an error.

Author: enblacar

NEWS.md

@@ -111,6 +111,7 @@ The reason of these modification is to allow for a much clearer and concise outp
 
 ## do_GroupwiseDEHeatmap()
 -   Reduced the legend title texts to allow for more room in the plot.
+-   Fixed a bug in which providing a DE genes object with `p_val_adj = 0` would cause the function to crash.
 
 ## do_LigandReceptorPlot()
 -   Added a new parameter `top_interactions_by_group` which when set to `TRUE` will report for each pair of `source` and `target`, as many interactions as stated in `top_interactions`.

R/do_GroupwiseDEHeatmap.R

@@ -210,6 +210,7 @@ do_GroupwiseDEHeatmap <- function(sample,
               dplyr::group_by(.data$cluster) %>%
               dplyr::filter(.data[[specificity]] <= p.cutoff) %>% 
               dplyr::slice_head(n = top_genes)
+
 
   gene.order <- data.use$gene %>% unique()
 
@@ -224,6 +225,10 @@ do_GroupwiseDEHeatmap <- function(sample,
               dplyr::arrange(dplyr::desc(.data[[magnitude]]), dplyr::desc(.data$p.adj.log10.minus)) %>% 
               dplyr::select(-dplyr::all_of(c(specificity)))
 
+  data.use$p.adj.log10.minus[data.use$p.adj.log10.minus == Inf] <- .Machine$double.xmax
+  
+  
+  
   # Workaround parameter deprecation.
   if (base::isTRUE(utils::packageVersion("Seurat") < "4.9.9")){
     data <- Seurat::GetAssayData(object = sample,