View source: R/logLik.mjoint.R
logLik.mjoint | R Documentation |
mjoint
objectExtract 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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.