Create Analyse function for piecewise exponential model

Create Analyse function for piecewise exponential model

Usage

analyse_piecewise_exponential(cuts, testing_only = FALSE)

Arguments

cuts

interval boundaries for the piecewise exponential model

testing_only

if set to TRUE omits all statistics in the intervals and just returns the p value of the global test.

Value

an analyse function that can be used in runSimulation

Details

If there's any time interval no patients ever enter, NA is returned for all time intervals. This behavior will likely change in future package versions.

Examples

condition <- merge(
   assumptions_delayed_effect(),
   design_fixed_followup(),
   by=NULL
 ) |>
 head(1)
dat <- generate_delayed_effect(condition)
analyse_piecewise_exponential(cuts=c(90, 360))(condition, dat)
#> $p
#> [1] 0.0003369932
#> 
#> $hr_s
#> # A tibble: 1 × 3
#>   `trt:intervalI1` `trt:intervalI2` `trt:intervalI3`
#>              <dbl>            <dbl>            <dbl>
#> 1            0.367            0.486            0.680
#> 
#> $p_s
#> # A tibble: 1 × 3
#>   `trt:intervalI1` `trt:intervalI2` `trt:intervalI3`
#>              <dbl>            <dbl>            <dbl>
#> 1           0.0565           0.0169          0.00330
#> 
#> $hr_s_lower
#> # A tibble: 1 × 3
#>   `trt:intervalI1` `trt:intervalI2` `trt:intervalI3`
#>              <dbl>            <dbl>            <dbl>
#> 1            -2.14            -1.33           -0.643
#> 
#> $hr_s_upper
#> # A tibble: 1 × 3
#>   `trt:intervalI1` `trt:intervalI2` `trt:intervalI3`
#>              <dbl>            <dbl>            <dbl>
#> 1          -0.0286           -0.142           -0.127
#> 
#> $N_pat
#> [1] 300
#> 
#> $N_evt
#> [1] 300
#> 
#> $interval_table
#>     
#>        0   1
#>   I1 282  18
#>   I2 234  48
#>   I3   0 234
#>