# NOT RUN {
fit <- brm(count ~ log_Age_c + log_Base4_c * Trt_c
+ (1|patient) + (1|visit),
data = epilepsy, family = "poisson")
plot(fit)
## plot population-level effects only
plot(fit, pars = "^b_")
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab