# NOT RUN {
fit <- brm(rating ~ treat + period + carry + (1|subject),
data = inhaler, family = "cumulative")
# extract posterior samples of population-level effects
samples1 <- posterior_samples(fit, pars = "^b")
head(samples1)
# extract posterior samples of group-level standard deviations
samples2 <- posterior_samples(fit, pars = "^sd_")
head(samples2)
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab