# NOT RUN {
# model with population-level effects only
fit1 <- brm(rating ~ treat + period + carry,
data = inhaler)
waic1 <- waic(fit1)
# model with an additional varying intercept for subjects
fit2 <- brm(rating ~ treat + period + carry + (1|subject),
data = inhaler)
waic2 <- waic(fit2)
# compare both models
compare_ic(waic1, waic2)
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab