## Not run:
# model <- brm(count ~ log_Age_c + log_Base4_c * Trt_c
# + (1|patient) + (1|visit),
# data = epilepsy, family = "poisson")
#
# # plot posterior intervals
# stanplot(model)
#
# # only show population-level effects in the plots
# stanplot(model, pars = "^b_")
#
# # show histograms of the posterior distributions
# stanplot(model, type = "hist")
#
# # plot some diagnostics of the sampler
# stanplot(model, type = "neff")
# stanplot(model, type = "rhat")
#
# # plot some diagnostics specific to the NUTS sampler
# stanplot(model, type = "nuts_acceptance")
# stanplot(model, type = "nuts_divergence")
# ## End(Not run)
Run the code above in your browser using DataLab