if (FALSE) {
require(xtable)
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 <- jointModel(lmeFit, survFit, timeVar = "time",
method = "weibull-PH-aGH")
xtable(jointFit, math.style.negative = TRUE)
}
Run the code above in your browser using DataLab