Posterior interval

posterior_interval(model2, prob = 0.80)
d %>% 
  add_fitted_draws(model2, n = 20) %>%
  ggplot(aes(x = ADF, y = ADG)) +
  geom_line(aes(y = .value, group = paste(Studycode, .draw))) + 
  facet_wrap(~ .draw) +theme_bw()
d %>%
  add_predicted_draws(model2, n = 100) %>%
  ggplot(aes(x = ADF)) +
  geom_density(aes(x = .prediction, group = .draw)) +
  xlim(-50,200) +theme_bw()

predict_study <- posterior_predict(
  model2, 
  newdata = data.frame(Studycode = c("V0011273"),
                       ADF = c(83),
                       Included2017=c("Yes")))
mcmc_areas(predict_study, prob = 0.8)+
  ggplot2::scale_y_discrete(labels = c("V0011273"))