Rizopoulos, D. (2012) Joint Models for Longitudinal and Time-to-Event Data: with
Applications in R. Boca Raton: Chapman and Hall/CRC.
Rizopoulos, D. (2010) JM: An R Package for the Joint Modelling of Longitudinal and Time-to-Event Data.
Journal of Statistical Software35 (9), 1–33. http://www.jstatsoft.org/v35/i09/
See Also
jointModel
Examples
## Not run:
# linear mixed model fit without treatment effect
fitLME.null <- lme(sqrt(CD4) ~ obstime,
random = ~ 1 | patient, data = aids)
# cox model fit without treatment effect
fitCOX.null <- coxph(Surv(Time, death) ~ 1,
data = aids.id, x = TRUE)
# joint model fit without treatment effect
fitJOINT.null <- jointModel(fitLME.null, fitCOX.null,
timeVar = "obstime", method = "weibull-PH-aGH")
# linear mixed model fit with treatment effect
fitLME.alt <- lme(sqrt(CD4) ~ obstime * drug - drug,
random = ~ 1 | patient, data = aids)
# cox model fit with treatment effect
fitCOX.alt <- coxph(Surv(Time, death) ~ drug,
data = aids.id, x = TRUE)
# joint model fit with treatment effect
fitJOINT.alt <- jointModel(fitLME.alt, fitCOX.alt, timeVar = "obstime",
method = "weibull-PH-aGH")
# likelihood ratio test for treatment effect
anova(fitJOINT.null, fitJOINT.alt)
## End(Not run)