# 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 <- imFitModel(abc, lst.var);
rst.imp <- imImpAll(abc, rst.fit, deltas=c(-0.25,0,0.25),
normal=TRUE, chains = 4, iter = 2000, warmup = 1000);
imPlotImputed(rst.imp, deltas=c(-0.25,0,0.25), xlim=c(0,100), endp=FALSE);
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab