# NOT RUN {
lst.var <- list(trt="TRT", surv="SURV", outcome=c("Y1","Y2"), y0=NULL,
endp=c("Y2"), unitTime="days",
trt.label = c("UC+SBT", "SAT+SBT"),
cov=c("AGE"), endfml="Y2", duration=365, bounds=c(0,100));
rst.fit <- imFit(abc, lst.var);
rst.imp <- imImpAll(abc, rst.fit, lst.var, deltas=c(-0.25,0,0.25),
normal=TRUE, iter=300, n.imp=2, thin=10, p.scale=10);
rst.est <- imEstimate(rst.imp, quantiles=c(0.25,0.5,0.75));
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab