deviance {BeQut}

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

Description

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

Usage

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.

Value

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

Author(s)

Antoine Barbieri and Baptiste Courrèges

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)