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.137   0.137   0.137   0.137   0.137   0.137 
summary(residuals(mod, type = 'working')[[1]])
#>    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
#>   1.159   1.159   1.159   1.159   1.159   1.159