reduce code

n.R

@@ -25,47 +25,20 @@ binding = FALSE
 info_frac = NULL
 analysis_time = c(24, 36)
 
-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)
+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
+)
 
-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
+event_new = c(floor(event[1]), ceiling(event[2]))
+sample_size_new = (ceiling(x$analysis$n[2] / 2) * 2)  # 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_new$analysis$n
-n - round(n)
+with(enroll_rate_new, sum(rate * duration)) - 454