Update utility_scripts.R

biobakery · Feb 24, 2025 · 598fd78 · 598fd78
1 parent f0a3caf
commit 598fd78
Showing 1 changed file with 4 additions and 2 deletions.
diff --git a/R/utility_scripts.R b/R/utility_scripts.R
@@ -315,8 +315,10 @@ write_results <- function(output,
     ###########################################
 
     significant_results <-
-        fit_data[fit_data$qval_joint <= max_significance | 
-                    fit_data$qval_individual <= max_significance &
+        fit_data[(fit_data$qval_joint <= max_significance &
+                      !is.na(fit_data$qval_joint)) | 
+                    (fit_data$qval_individual <= max_significance &
+                         !is.na(fit_data$qval_individual)) &
                     is.na(fit_data$error),]
     significant_results$error <- NULL
     significant_results_file <-