Extract log-likelihood from an mjoint object
Extract log-likelihood from an mjoint object.
## S3 method for class 'mjoint' logLik(object, ...)
object |
an object inheriting from class |
... |
additional arguments; currently none are used. |
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.
Graeme L. Hickey (graemeleehickey@gmail.com)
Henderson R, Diggle PJ, Dobson A. Joint modelling of longitudinal measurements and event time data. Biostatistics. 2000; 1(4): 465-480.
logLik for the generic method description.
## 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)Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.