Some residual methods for tidybayes

tidy_residuals.R

tidy_residuals <- function(model){
  UseMethod("tidy_residuals")
}

tidy_residuals.brmsfit <- function(model){
  res <- residuals(model, summary=FALSE)
  props <- summary(model)
  nchains <- props$chains
  iter <- props$iter - props$warmup
  
  start <- seq(1, nrow(res), nrow(res)/nchains)
  
  chains <- map(start, ~as.mcmc(res[.x:(.x+iter-1),])) %>%
    coda::as.mcmc.list()
  
  tidy_draws(chains)
  
}

gather_residuals <- function(model){
  tidy_residuals(model) %>% 
    gather(.variable, .value, -.chain, -.iteration, -.draw)
}

spread_residuals <- function(model){
  tidy_residuals(model)
}