fit <- coxdual(Surv(start, stop, status) ~ cluster(id)
+ trans(from, to) + z, data = dualrc, init.coxph = TRUE)
haz <- step2jump(fit$coxph$basehaz, stratum = 3)
Haz <- jump2step(haz, stratum = 3)
all(Haz == fit$coxph$basehaz)
Run the code above in your browser using DataLab