Merge pull request #470 from Merck/469-info-scale-when-updating-design-from-initial-alpha-to-full-alpha

Correct `info_scale` in `gs_update_ahr()` when updating design from initial alpha to full alpha

Author: keaven

R/gs_design_ahr.R

@@ -360,6 +360,7 @@ gs_design_ahr <- function(
     enroll_rate = enroll_rate, fail_rate = fail_rate,
     alpha = alpha, beta = beta, ratio = ratio,
     info_frac = info_frac, analysis_time = analysis_time,
+    info_scale = info_scale,
     upper = upper, upar = upar,
     lower = lower, lpar = lpar,
     test_upper = test_upper, test_lower = test_lower,

R/gs_power_ahr.R

@@ -261,6 +261,7 @@ gs_power_ahr <- function(
   input <- list(
     enroll_rate = enroll_rate, fail_rate = fail_rate,
     event = event, analysis_time = analysis_time,
+    info_scale = info_scale,
     upper = upper, upar = upar,
     lower = lower, lpar = lpar,
     test_lower = test_lower, test_upper = test_upper,

R/gs_update_ahr.R

@@ -296,8 +296,16 @@ gs_update_ahr <- function(
     x_updated_h0 <- gs_power_npe(theta = 0,
                                  theta0 = 0,
                                  theta1 = x$analysis$theta,
+                                 # statistical information at all analyses for input `theta`
                                  info = x$analysis$info0,
-                                 info_scale = "h0_info",
+                                 # statistical information under null hypothesis,
+                                 # impacts null hypothesis bound calculation.
+                                 info0 = x$analysis$info0,
+                                 # statistical information under hypothesis used for
+                                 # futility bound calculation if different from `info`
+                                 # impacts futility hypothesis bound calculation.
+                                 info1 = x$analysis$info,
+                                 info_scale = x$input$info_scale,
                                  upper = x$input$upper, upar = upar_update,
                                  test_upper = x$input$test_upper,
                                  lower = x$input$lower, lpar = x$input$lpar,
@@ -308,9 +316,16 @@ gs_update_ahr <- function(
     x_updated_h1 <- gs_power_npe(theta = x$analysis$theta,
                                  theta0 = 0,
                                  theta1 = x$analysis$theta,
+                                 # statistical information under null hypothesis,
+                                 # impacts null hypothesis bound calculation.
                                  info0 = x$analysis$info0,
+                                 # statistical information at all analyses for input `theta`
                                  info = x$analysis$info,
-                                 info_scale = "h0_info",
+                                 # statistical information under hypothesis used for
+                                 # futility bound calculation if different from `info`
+                                 # impacts futility hypothesis bound calculation.
+                                 info1 = NULL,
+                                 info_scale = x$input$info_scale,
                                  upper = x$input$upper, upar = upar_update,
                                  test_upper = x$input$test_upper,
                                  lower = x$input$lower, lpar = x$input$lpar,