fit_i <- brm(rating ~ treat + period + carry + (1|subject), data = inhaler, family = "cumulative")
#extract posterior samples of fixed effects
samples1 <- posterior.samples(fit_i, "b_")
head(samples1)
#extract posterior samples of standard deviations of random effects
samples2 <- posterior.samples(fit_i, "sd_")
head(samples2)
Run the code above in your browser using DataLab