predictions2

predictions <- posterior_predict(model2, newdata = d)
ppc_intervals(d$ADG, 
              yrep = predictions, 
              x = d$ADF, 
              prob = 0.5, prob_outer = 0.95)
prediction_summary(model2, data = d)
cv_procedure <- prediction_summary_cv(model = model2,
                                      group = "Studycode",
                                      data  = d, 
                                      k = 10)
cv_procedure$cv
model_elpd <- loo(model2, save_psis = TRUE)
print(model_elpd)
model_elpd$estimates
plot(model_elpd)
yrep <- posterior_predict(model2)
ppc_loo_pit_overlay(
  y = d$ADG,
  yrep = yrep,
  lw = weights(model_elpd$psis_object)
)
ppc_loo_pit_qq(
  y = d$ADG,
  yrep = yrep,
  lw = weights(model_elpd$psis_object)
)