deviance function

deviance returns the deviance based on the conditional likelihood associated with the survival part.

deviance(object, M = 1000, conditional = "survival", verbose = TRUE)

Arguments

• object: an object inheriting from class 'Bqrjm'.
• M: an integer indicating the number of draws used for the approximation of the integral with respect to random effects, M=1000 by default.
• conditional: is "survival" by default (only this one is implemented until now).
• verbose: A logical indicating if information about method's progress (included progress bars for each step) must be printed (default to TRUE). Adds a small extra overload.

Returns

An object which is a list with the following elements:

• deviance: Numerical object returning the deviance
• likelihood: (Conditional) likelihood
• sims.list: list of individual quantities as likelihood, draws of random effects, hazard and survival functions
• control: list of arguments giving details about the deviance

Examples

#---- load data
data(dataLong)

#---- Fit quantile regression joint model for the median
qrjm_50 <- qrjm(formFixed = y ~ visit,
               formRandom = ~ visit,
               formGroup = ~ ID,
               formSurv = survival::Surv(time, event) ~ X1 + X2,
               survMod = "weibull",
               param = "value",
               timeVar= "visit",
               data = dataLong,
               save_va = TRUE,
               parallel = FALSE,
               tau = 0.5)

deviance(qrjm_50, M=200)

Author(s)

Antoine Barbieri and Baptiste Courrèges