Using purrr to do bootstrap estimation of the mean

purrr_bootstrap.R

library(tidyverse)

replicates <- (1:100000)%>%
                map(~sample(faithful$waiting, replace = T))%>%
                map(mean)%>%
                simplify()

data_frame(replicates = replicates)%>%
  ggplot(aes(replicates))+
    stat_density()