Skip to contents

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

Source: R/Ex2delayedEffect.R
Ex2delayedEffect.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(Ex2delayedEffect)

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

Ex1delayedEffect, Ex3curewithph, Ex4belly, Ex5widening, Ex6crossing

Examples

library(survival)

data(Ex2delayedEffect)
km1 <- with(Ex2delayedEffect, survfit(Surv(month, evntd) ~ trt))
km1
#> Call: survfit(formula = Surv(month, evntd) ~ trt)
#> 
#>         n events median 0.95LCL 0.95UCL
#> trt=0 137    123   2.84    2.18    3.50
#> trt=1 135    105   3.45    2.13    5.07
plot(km1)

with(subset(Ex2delayedEffect, trt == 1), survfit(Surv(month, evntd) ~ trt))
#> Call: survfit(formula = Surv(month, evntd) ~ trt)
#> 
#>        n events median 0.95LCL 0.95UCL
#> [1,] 135    105   3.45    2.13    5.07
with(subset(Ex2delayedEffect, trt == 0), survfit(Surv(month, evntd) ~ trt))
#> Call: survfit(formula = Surv(month, evntd) ~ trt)
#> 
#>        n events median 0.95LCL 0.95UCL
#> [1,] 137    123   2.84    2.18     3.5