## Not run:
# # Fit a joint model with bivariate longitudinal outcomes
#
# data(heart.valve)
# hvd <- heart.valve[!is.na(heart.valve$log.grad) & !is.na(heart.valve$log.lvmi), ]
#
# fit2 <- mjoint(
# formLongFixed = list("grad" = log.grad ~ time + sex + hs,
# "lvmi" = log.lvmi ~ time + sex),
# formLongRandom = list("grad" = ~ 1 | num,
# "lvmi" = ~ time | num),
# formSurv = Surv(fuyrs, status) ~ age,
# data = list(hvd, hvd),
# inits = list("gamma" = c(0.11, 1.51, 0.80)),
# timeVar = "time",
# verbose = TRUE)
# sampleData(fit2, size = 10)
# ## End(Not run)
Run the code above in your browser using DataLab