diff --git a/qhtcp-workflow/apps/r/calculate_interaction_zscores.R b/qhtcp-workflow/apps/r/calculate_interaction_zscores.R
index 1c35e639..e019df67 100644
--- a/qhtcp-workflow/apps/r/calculate_interaction_zscores.R
+++ b/qhtcp-workflow/apps/r/calculate_interaction_zscores.R
@@ -1043,7 +1043,7 @@ main <- function() {
     df_na_stats <- ss$df_with_stats
     write.csv(df_na_ss, file = file.path(out_dir, "summary_stats_all_strains.csv"), row.names = FALSE)
     # For plotting (ggplot warns on NAs)
-    df_na_stats_filtered <- df_na_stats %>% filter(across(all_of(summary_vars), is.finite))
+    df_na_stats_filtered <- df_na_stats %>% filter(if_all(all_of(summary_vars), is.finite))
 
     df_na_stats <- df_na_stats %>%
       mutate(