library(rstan) library(shredder) library(tidybayes) library(microbenchmark) rstan_options(auto_write = TRUE)
Set seed
set.seed(1234)
Run Model
rats <- shredder::rats_example(nCores = 4)
tb <- function(mod){ mod%>% tidybayes::spread_draws(mu_alpha,mu_beta) }
sh_select <- function(mod){ mod%>% stan_select(mu_alpha,mu_beta)%>% tidybayes::spread_draws(mu_alpha,mu_beta) }
sh_thin <- function(mod){ mod%>% stan_thin_n(20)%>% tidybayes::spread_draws(mu_alpha,mu_beta) }
bench <- microbenchmark( "tidybayes" = {tb(rats)}, "shredder select" = {sh_select(rats)}, "shredder thin" = {sh_thin(rats)})
ggplot2::autoplot(bench)