# \donttest{
# Fit a joint model
fit <- jmcs(ydata = ydata, cdata = cdata,
long.formula = response ~ time + x1,
surv.formula = Surv(surv, failure_type) ~ x1 + x2,
random = ~ time| ID)
fit2 <- jmcs(ydata = ydata, cdata = cdata,
long.formula = response ~ time + gender + x1 + race,
surv.formula = Surv(surv, failure_type) ~ x1 + gender + x2 + race,
random = ~ time| ID)
anova(fit, fit2)
# }
Run the code above in your browser using DataLab