fix the bug on p_obs

Merck · Jan 31, 2024 · b1bc45f · b1bc45f
1 parent ef046d3
commit b1bc45f
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/R/calc_seq_p.R b/R/calc_seq_p.R
@@ -129,7 +129,7 @@ calc_seq_p <- function(
 
     p_diff <- NULL
     for (hhh in all_hypothesis) {
-      p_diff_new <- p_obs[[hhh]] - ans[[hhh]]
+      p_diff_new <- (p_obs %>% filter(analysis <= test_analysis))[[hhh]] - ans[[hhh]]
       p_diff <- c(p_diff, p_diff_new)
     }