12.6 Posterior Distribution

# shortcut instead of running stan_glm() again
equality_model <- update(equality_model_prior, prior_PD = FALSE)

12.6.1 Checks

mcmc_trace(equality_model)

mcmc_dens_overlay(equality_model)

mcmc_acf(equality_model)

12.6.2 Posterior Predictive Check

set.seed(1)
pp_check(equality_model, plotfun = "hist", nreps = 5) + 
  xlab("laws")

pp_check(equality_model) + 
  xlab("laws")