R/percentile-data-frame.R

In DavisVaughan/rboot: What the Package Does (One Line, Title Case)

boot_pctl_df <- function(data, ..., values = NULL, alpha = 0.05, times = 1000) {

  summarise_exprs <- rlang::enquos(...)

  data_bootstrapped <- strapgod::bootstrapify(data, times = times)

  estimates <- dplyr::summarise(data_bootstrapped, !!! summarise_exprs)

  syms_groups <- dplyr::groups(data)
  vars_estimates <- tidyselect::vars_select(names(estimates), !!rlang::enquo(values))
  vars_estimates <- setdiff(vars_estimates, ".bootstrap")

  # keep all
  if (length(vars_estimates) == 0L) {
    vars_estimates <- names(summarise_exprs)
  }

  syms_estimates <- rlang::syms(vars_estimates)

  estimates <- dplyr::select(estimates, !!!syms_groups, !!!syms_estimates)

  estimates <- tidyr::gather(estimates, ".statistic", ".estimate", !!!syms_estimates)

  estimates <- dplyr::group_by(estimates, .statistic, add = TRUE)

  out <- dplyr::group_map(estimates, ~calc_boot_pctl(.x$.estimate, alpha = alpha))
  out <- dplyr::ungroup(out)

  out
}