summary.mjoint: Summary of an 'mjoint' object

View source: R/summary.mjoint.R

summary.mjointR Documentation

Summary of an mjoint object


This function provides a summary of an mjoint object.


## S3 method for class 'mjoint'
summary(object, bootSE = NULL, ...)



an object inheriting from class mjoint for a joint model of time-to-event and multivariate longitudinal data.


an object inheriting from class bootSE for the corresponding model. If bootSE=NULL, the function will attempt to utilize approximate standard error estimates (if available) calculated from the empirical information matrix.


additional arguments; currently none are used.


A list containing the coefficient matrices for the longitudinal and time-to-event sub-models; variance-covariance matrix for the random effects; residual error variances; log-likelihood of joint model; AIC and BIC statistics; and model fit objects.


Graeme L. Hickey (


Wulfsohn MS, Tsiatis AA. A joint model for survival and longitudinal data measured with error. Biometrics. 1997; 53(1): 330-339.

Henderson R, Diggle PJ, Dobson A. Joint modelling of longitudinal measurements and event time data. Biostatistics. 2000; 1(4): 465-480.

Lin H, McCulloch CE, Mayne ST. Maximum likelihood estimation in the joint analysis of time-to-event and multiple longitudinal variables. Stat Med. 2002; 21: 2369-2382.

See Also

mjoint, mjoint.object, and summary for the generic method description.


## Not run: 
# Fit a joint model with bivariate longitudinal outcomes

hvd <- heart.valve[!$log.grad) & !$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)

## End(Not run)

joineRML documentation built on Jan. 22, 2023, 1:18 a.m.