21.4 reg_intervals from {rsample}

Similar purpose (?) as {infer} package, few supported models and interval types

reg_intervals(
  my_formula,
  race_top_results,
  model_fn = "glm",
  times = 2000,
  type = "percentile",
  keep_reps = TRUE
)

## # A tibble: 2 × 7
##   term            .lower .estimate   .upper .alpha .method           .replicates
##   <chr>            <dbl>     <dbl>    <dbl>  <dbl> <chr>      <list<tibble[,2]>>
## 1 aid_stations -0.0286   -0.0128   0.00200    0.05 percentile        [2,000 × 2]
## 2 participants  0.000122  0.000440 0.000904   0.05 percentile        [2,000 × 2]