Allow either q-values in significant results

biobakery · Nov 15, 2024 · 3093a65 · 3093a65
1 parent e7644c3
commit 3093a65
Show file tree

Hide file tree

Showing 3 changed files with 1,381 additions and 1,152 deletions.
diff --git a/R/utility_scripts.R b/R/utility_scripts.R
@@ -333,7 +333,8 @@ write_results <- function(output,
     ###########################################
 
     significant_results <-
-        fit_data[fit_data$qval_joint <= max_significance &
+        fit_data[fit_data$qval_joint <= max_significance | 
+                    fit_data$qval_individual <= max_significance &
                     is.na(fit_data$error),]
     significant_results$error <- NULL
     significant_results_file <-