Compute posterior draws of residuals/predictive errors
Usage
# S3 method for class 'brmcoda'
residuals(object, ...)
Arguments
- object
An object of class
brmcoda
.- ...
Further arguments passed to
residuals.brmsfit
.
Value
An array
of predictive error/residual draws. If
summary = FALSE
the output resembles those of
predictive_error.brmsfit
. If summary = TRUE
the output
is an N x E matrix, where N is the number of observations and E denotes
the summary statistics computed from the draws.
Examples
# \donttest{
## fit a model
if(requireNamespace("cmdstanr")){
m <- brmcoda(complr = complr(data = mcompd, sbp = sbp,
parts = c("TST", "WAKE", "MVPA", "LPA", "SB"),
idvar = "ID", total = 1440),
formula = Stress ~ bilr1 + bilr2 + bilr3 + bilr4 +
wilr1 + wilr2 + wilr3 + wilr4 + (1 | ID),
chain = 1, iter = 500,
backend = "cmdstanr")
## extract residuals
res <- residuals(m)
head(res)
}# }
#> Error: CmdStan path has not been set yet. See ?set_cmdstan_path.