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