Extract residuals from an object of class JointAI

# S3 method for JointAI
residuals(object, type = c("working", "pearson", "response"), warn = TRUE, ...)

## Arguments

object

object inheriting from class 'JointAI'

type

type of residuals: "deviance", "response", "working"

warn

logical; should warnings be given? Default is TRUE.

...

currently not used

## Note

• For mixed models residuals are currently calculated using the fixed effects only.

• For ordinal (mixed) models and parametric survival models only type = "response" is available.

• For Cox proportional hazards models residuals are not yet implemented.

## Examples

mod <- glm_imp(B1 ~ C1 + C2 + O1, data = wideDF, n.iter = 100,
family = binomial(), mess = FALSE)
summary(residuals(mod, type = 'response')[[1]])
#>    Min. 1st Qu.  Median    Mean 3rd Qu.    Max.
#>  0.1471  0.1471  0.1471  0.1471  0.1471  0.1471
summary(residuals(mod, type = 'working')[[1]])
#>    Min. 1st Qu.  Median    Mean 3rd Qu.    Max.
#>   1.172   1.172   1.172   1.172   1.172   1.172