Compute posterior draws of residuals/predictive errors
## S3 method for class 'brmcoda'residuals(object,...)
Arguments
object: An object of class brmcoda.
...: Further arguments passed to residuals.brmsfit.
Returns
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
## fit a modelif(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)}