## Not run:
# model <- brm(count ~ log_Age_c + log_Base4_c * Trt_c
# + (1|patient) + (1|visit),
# data = epilepsy, family = "poisson")
# # plot 95% CIs
# stanplot(model, type = "plot", ci_level = 0.95)
# # equivalent to
# stan_plot(model$fit, ci_level = 0.95)
#
# # only show fixed effects in the plots
# # this will not work when calling stan_plot directly
# stanplot(model, pars = "^b", type = "plot", ci_level = 0.95)
#
# # plot some diagnostics on the sampler
# stanplot(model, type = "diag")
# # equivalent to
# stan_diag(model$fit)
# ## End(Not run)
Run the code above in your browser using DataLab