Reduce interval for root finding

Merck · Dec 19, 2023 · afc5e6a · afc5e6a
1 parent 6d4371a
commit afc5e6a
Show file tree

Hide file tree

Showing 4 changed files with 8 additions and 5 deletions.
diff --git a/DESCRIPTION b/DESCRIPTION
@@ -1,6 +1,6 @@
 Package: gsDesign2
 Title: Group Sequential Design with Non-Constant Effect
-Version: 1.1.0.1
+Version: 1.1.0.2
 Authors@R: c(
     person("Keaven", "Anderson", email = "[email protected]", role = c("aut")),
     person("Yilong", "Zhang", email = "[email protected]", role = c("aut")),

diff --git a/R/expected_time.R b/R/expected_time.R
@@ -95,6 +95,7 @@ expected_time <- function(
     target_event = 150,
     ratio = 1,
     interval = c(.01, 100)) {
+  print(interval)
   # ----------------------------#
   #    check inputs             #
   # ----------------------------#

diff --git a/R/gs_power_ahr.R b/R/gs_power_ahr.R
@@ -52,7 +52,8 @@
 #'   Normally, `r` will not be changed by the user.
 #' @param tol Tolerance parameter for boundary convergence (on Z-scale).
 #' @param interval An interval that is presumed to include the time at which
-#'   expected event count is equal to targeted event.
+#'   expected event count is equal to targeted event. If `analysis_time` is
+#'   defined, the upper limit of the interval is set to its maximum.
 #'
 #' @return A tibble with columns `Analysis`, `Bound`, `Z`, `Probability`,
 #'   `theta`, `Time`, `AHR`, `Events`.
@@ -171,7 +172,7 @@ gs_power_ahr <- function(
     info_scale = c("h0_h1_info", "h0_info", "h1_info"),
     r = 18,
     tol = 1e-6,
-    interval = c(.01, 100)) {
+    interval = c(.01, if(is.null(analysis_time)) 100 else max(analysis_time))) {
   # Get the number of analysis
   n_analysis <- max(length(event), length(analysis_time), na.rm = TRUE)
 

diff --git a/man/gs_power_ahr.Rd b/man/gs_power_ahr.Rd