Skip to contents

gs_power_ahr: computes power using average hazard ratio under non-proportional hazards

Source: vignettes/usage_gs_power_ahr.Rmd
usage_gs_power_ahr.Rmd
library(tibble)
library(dplyr)
#library(gsDesign2)
devtools::load_all()

Introduction of gs_power_ahr()

Use cases of gs_power_ahr()

Example 1 

x <- gs_power_ahr(enrollRates = tibble(Stratum = "All", duration = c(2, 2, 10), rate = c(3, 6, 9)),
                  failRates = tibble(Stratum = "All", duration = c(3, 100), failRate = log(2) / c(9, 18), 
                                     hr = c(.9, .6), dropoutRate = rep(.001, 2)),
                  analysisTimes = c(12, 24, 36), events = NULL,
                  binding = TRUE,
                  upper = gs_spending_bound,
                  upar = list(sf = gsDesign::sfLDOF, total_spend = 0.025, param = NULL, timing = NULL),
                  lower = gs_spending_bound,
                  lpar = list(sf = gsDesign::sfLDOF, total_spend = 0.025, param = NULL, timing = NULL))
x
## $enrollRates
## # A tibble: 3 × 3
##   Stratum duration  rate
##   <chr>      <dbl> <dbl>
## 1 All            2     3
## 2 All            2     6
## 3 All           10     9
## 
## $failRates
## # A tibble: 2 × 5
##   Stratum duration failRate    hr dropoutRate
##   <chr>      <dbl>    <dbl> <dbl>       <dbl>
## 1 All            3   0.0770   0.9       0.001
## 2 All          100   0.0385   0.6       0.001
## 
## $bounds
## # A tibble: 6 × 7
##   Analysis Bound Probability Probability0      Z `~HR at bound` `Nominal p`
##      <int> <chr>       <dbl>        <dbl>  <dbl>          <dbl>       <dbl>
## 1        1 Upper   0.000370     0.0000538  3.87           0.178   0.0000538
## 2        1 Lower   0.0000612    0.0000538 -3.40           4.55    1.00     
## 3        2 Upper   0.116        0.00921    2.36           0.506   0.00919  
## 4        2 Lower   0.00907      0.00921   -1.20           1.42    0.885    
## 5        3 Upper   0.324        0.0250     2.01           0.608   0.0222   
## 6        3 Lower   0.0250       0.0250    -0.473          1.12    0.682    
## 
## $analysis
## # A tibble: 3 × 10
##   Analysis  Time     N Events   AHR theta  info info0    IF   IF0
##      <int> <dbl> <dbl>  <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1        1    12    90   20.4 0.811 0.210  5.03  5.10 0.309 0.308
## 2        2    24   108   49.1 0.715 0.335 12.0  12.3  0.738 0.741
## 3        3    36   108   66.2 0.683 0.381 16.3  16.6  1     1    
## 
## attr(,"class")
## [1] "ahr"       "gs_design" "list"

Inner Logic of gs_power_ahr()

TODO