## Load data and fit the model
data(bmt)
bmt = transform(bmt, d4=d2+d3)
A = as.numeric(bmt$group>1)
## Model with competing risk data
fit1 = surv.tteICE(A, bmt$t2, bmt$d4, 'composite')
## Plot asymptotic confidence intervals based on explicit formulas
plot_ate(fit1, ylim=c(-0.4,0.4))
# \donttest{
## Plot bootstrap confidence intervals (may take some seconds)
plot_ate(fit1, nboot=200, ylim=c(-0.4,0.4))
# }
## Model with semicompeting risk data
fit2 = scr.tteICE(A, bmt$t1, bmt$d1, bmt$t2, bmt$d2, "composite")
## Plot asymptotic confidence intervals based on explicit formulas
plot_ate(fit2, ylim=c(-0.4,0.4),
plot.configs=list(add.null.line=FALSE))
## Plot bootstrap confidence intervals
plot_ate(fit2, nboot=200, ylim=c(-0.4,0.4),
plot.configs=list(add.null.line=FALSE, lty=2, main=""))
Run the code above in your browser using DataLab