From 270cc8def5251153c5e72f5ea7a0792af68e587f Mon Sep 17 00:00:00 2001 From: Bryan Roessler Date: Thu, 12 Sep 2024 19:24:53 -0400 Subject: [PATCH] Gracefully handle empty var cols --- qhtcp-workflow/apps/r/calculate_interaction_zscores.R | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/qhtcp-workflow/apps/r/calculate_interaction_zscores.R b/qhtcp-workflow/apps/r/calculate_interaction_zscores.R index fcb4f48c..98ca2cc3 100644 --- a/qhtcp-workflow/apps/r/calculate_interaction_zscores.R +++ b/qhtcp-workflow/apps/r/calculate_interaction_zscores.R @@ -171,7 +171,7 @@ calculate_summary_stats <- function(df, variables, group_vars = c("OrfRep", "con ), .names = "{.fn}_{.col}") ) - print(summary_stats) + # print(summary_stats) # Prevent .x and .y suffix issues by renaming columns df_cleaned <- df %>% @@ -224,10 +224,10 @@ calculate_interaction_scores <- function(df, max_conc, variables, group_vars = c across(all_of(variables), list( mean = ~mean(., na.rm = TRUE), median = ~median(., na.rm = TRUE), - max = ~max(., na.rm = TRUE), - min = ~min(., na.rm = TRUE), + max = ~ifelse(all(is.na(.)), NA, max(., na.rm = TRUE)), + min = ~ifelse(all(is.na(.)), NA, min(., na.rm = TRUE)), sd = ~sd(., na.rm = TRUE), - se = ~sd(., na.rm = TRUE) / sqrt(sum(!is.na(.)) - 1) + se = ~ifelse(sum(!is.na(.)) > 1, sd(., na.rm = TRUE) / sqrt(sum(!is.na(.)) - 1), NA) ), .names = "{.fn}_{.col}") ) %>% ungroup()