joineRML (version 0.1.1)

logLik.mjoint: Extract log-likelihood from an mjoint object

Description

Extract log-likelihood from an mjoint object.

Usage

"logLik"(object, ...)

Arguments

object

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

...

additional arguments; currently none are used.

Value

Returns an object of class logLik. This is a number with two attributes: "df" (degrees of freedom), giving the number of parameters in the model, and "nobs", the number of observations used in estimation.

References

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

See Also

logLik for the generic method description.

Examples

## 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)
# 
# logLik(fit2)
# ## End(Not run)