library(survival)
fit <- survreg(Surv(time, status) ~ ph.ecog + age + sex + strata(sex),
data = lung)
nd = transform(expand.grid(ph.ecog=0:1, sex=1:2), age=60)
simulate(fit, seed=1002, newdata=nd)
simulate(fit, seed=1002, newdata=nd, t0=500)
Run the code above in your browser using DataLab