Skip to contents

Time-to-event data example 3 for non-proportional hazards working group

Source: R/ex3_cure_with_ph.R
ex3_cure_with_ph.Rd

Survival objects reverse-engineered datasets from published Kaplan-Meier curves. Individual trials are de-identified since the data are only approximations of the actual data. Data are intended to evaluate methods and designs for trials where non-proportional hazards may be anticipated for outcome data.

Usage

data(ex3_cure_with_ph)

Format

Data frame with 4 variables:

  • id: Sequential numbering of unique identifiers.

  • month: Time-to-event.

  • event: 1 for event, 0 for censored.

  • trt: 1 for experimental, 0 for control.

References

Lin, Ray S., Ji Lin, Satrajit Roychoudhury, Keaven M. Anderson, Tianle Hu, Bo Huang, Larry F Leon, Jason J.Z. Liao, Rong Liu, Xiaodong Luo, Pralay Mukhopadhyay, Rui Qin, Kay Tatsuoka, Xuejing Wang, Yang Wang, Jian Zhu, Tai-Tsang Chen, Renee Iacona & Cross-Pharma Non-proportional Hazards Working Group. 2020. Alternative analysis methods for time to event endpoints under nonproportional hazards: A comparative analysis. Statistics in Biopharmaceutical Research 12(2): 187–198.

See also

ex1_delayed_effect, ex2_delayed_effect, ex4_belly, ex5_widening, ex6_crossing

Examples

library(survival)

data(ex3_cure_with_ph)
km1 <- with(ex3_cure_with_ph, survfit(Surv(month, evntd) ~ trt))
km1
#> Call: survfit(formula = Surv(month, evntd) ~ trt)
#> 
#>         n events median 0.95LCL 0.95UCL
#> trt=0 137    101   1.05   0.523    1.74
#> trt=1 143     86   1.74   1.158    3.13
plot(km1)