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