In metrumresearchgroup/shredder: Rstan Manipulation API

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)

Set up functions

tidybayes

tb <- function(mod){
  mod%>%
    tidybayes::spread_draws(mu_alpha,mu_beta)  
}

tidybayes with shredder select

sh_select <- function(mod){
  mod%>%
    stan_select(mu_alpha,mu_beta)%>%
    tidybayes::spread_draws(mu_alpha,mu_beta)  
}

tidybayes with shredder thin (every 20th sample)

sh_thin <- function(mod){
  mod%>%
    stan_thin_n(20)%>%
    tidybayes::spread_draws(mu_alpha,mu_beta)
}

Benchmark

Run benchmark

bench <- microbenchmark(
  "tidybayes" = {tb(rats)},
  "shredder select" = {sh_select(rats)},
  "shredder thin" = {sh_thin(rats)})

Output

ggplot2::autoplot(bench)

metrumresearchgroup/shredder documentation built on Sept. 1, 2020, 5:36 p.m.