more debugging

n.R

@@ -1,38 +1,71 @@
 install.packages('gsDesign2', quiet = TRUE)
 
-ratio <- 1
+library(gsDesign2)
 
-enroll_rate <- gsDesign2::define_enroll_rate(duration = c(2, 2, 8),
-                                             rate = c(1, 2, 3))
+ratio = 1
 
-fail_rate <- gsDesign2::define_fail_rate(duration = c(4, Inf), 
-                                         fail_rate = log(2) / 12,
-                                         hr = c(1, .6),
-                                         dropout_rate = .001)
+enroll_rate = define_enroll_rate(duration = c(2, 2, 8), rate = c(1, 2, 3))
 
-alpha <- 0.025
-beta <- 0.1
+fail_rate = define_fail_rate(
+  duration = c(4, Inf), fail_rate = log(2) / 12, hr = c(1, .6), dropout_rate = .001
+)
 
-upper <- gsDesign2::gs_spending_bound
-upar <- list(sf = gsDesign::sfLDOF, total_spend = alpha)
-test_upper <- rep(TRUE, 2)
+alpha = 0.025
+beta = 0.1
 
-lower <- gsDesign2::gs_spending_bound
-lpar <- list(sf = gsDesign::sfLDOF, total_spend = beta)
-test_lower <- c(TRUE, FALSE)
-binding <- FALSE
+upper = gs_spending_bound
+upar = list(sf = gsDesign::sfLDOF, total_spend = alpha)
+test_upper = rep(TRUE, 2)
+
+lower = gs_spending_bound
+lpar = list(sf = gsDesign::sfLDOF, total_spend = beta)
+test_lower = c(TRUE, FALSE)
+binding = FALSE
 
 info_frac = NULL
 analysis_time = c(24, 36)
 
-x <- gsDesign2::gs_design_ahr(enroll_rate = enroll_rate, fail_rate = fail_rate, 
+x = gs_design_ahr(enroll_rate = enroll_rate, fail_rate = fail_rate, 
                               alpha = alpha, beta = beta, ratio = ratio,
                               info_frac = info_frac, analysis_time = analysis_time, 
                               upper = upper, upar = upar, test_upper = test_upper,
                               lower = lower, lpar = lpar, test_lower = test_lower,
-                              binding = binding) |> gsDesign2::to_integer()
+                              binding = binding)
+
+n_analysis <- length(x$analysis$analysis)
+multiply_factor <- x$input$ratio + 1
+event = x$analysis$event
+event_new = c(floor(event[1:(n_analysis - 1)]), ceiling(event[n_analysis]))
+sample_size_new = (ceiling(x$analysis$n[n_analysis] / multiply_factor) * multiply_factor)  # 454L
+enroll_rate = x$enroll_rate
+enroll_rate_new = within(enroll_rate, {
+  rate = rate * sample_size_new / x$analysis$n[2]
+})
+upar_new = x$input$upar
+info_with_new_event = gs_info_ahr(
+  enroll_rate = enroll_rate_new,
+  fail_rate = x$input$fail_rate,
+  ratio = x$input$ratio,
+  event = event_new,
+  analysis_time = NULL
+)
+upar_new$timing = info_with_new_event$info / max(info_with_new_event$info)
+lpar_new = x$input$lpar
+lpar_new$timing = upar_new$timing
+x_new = gs_power_ahr(
+  enroll_rate = enroll_rate_new,
+  fail_rate = x$input$fail_rate,
+  event = event_new,
+  analysis_time = NULL,
+  ratio = x$input$ratio,
+  upper = x$input$upper, upar = upar_new,
+  lower = x$input$lower, lpar = lpar_new,
+  test_upper = x$input$test_upper,
+  test_lower = x$input$test_lower,
+  binding = x$input$binding,
+  info_scale = x$input$info_scale, r = x$input$r, tol = x$input$tol,
+  interval = c(0.01, max(x$analysis$time) + 100)
+)
 
-n <- x$analysis$n
+n = x_new$analysis$n
 n - round(n)
-as.integer(n)
-length(sample('All', n[2], TRUE))