Initial approach to bootstrapping a mean value, but is much slower than just using base::sample() so used that...

bootstrap-mean-rsample.R

library(tidymodels)

### MUCH faster computationally to use base R `sample()` for this step... so did not use this approach
resamples <- rsample::bootstraps(preds, 5000)

avg_diff_sample <- function(split){
  analysis(split) %>% 
    summarise(diff = mean(diff_abs_resids)) %>% 
    pull(diff)
}

resamples <- resamples %>% 
  transmute(diff = map_dbl(splits, avg_diff_sample))