# \donttest{
fit <- 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)
oldpar <- par(mfrow = c(2, 2))
plot(fit)
par(oldpar)
# }
Run the code above in your browser using DataLab