# NOT RUN {
prothro$t0 <- as.numeric(prothro$time == 0)
lmeFit <- lme(pro ~ treat * (time + t0), random = ~ time | id, data = prothro)
survFit <- coxph(Surv(Time, death) ~ treat, data = prothros, x = TRUE)
jointFit <- jointModelBayes(lmeFit, survFit, timeVar = "time")
if (require("xtable")) {
xtable:::xtable(jointFit, math.style.negative = TRUE)
}
# }
Run the code above in your browser using DataLab