## Not run: 
# ## poisson regression without random effects. 
# fit1 <- brm(count ~ log_Age_c + log_Base4_c * Trt_c, 
#             data = epilepsy, family = poisson())
# summary(fit1) 
# plot(fit1)             
#     
# ## poisson regression with random intercepts of patients and visits
# ## as well as normal priors for fixed effects parameters.    
# fit2 <- brm(count ~ log_Age_c + log_Base4_c * Trt_c 
#             + (1|patient) + (1|visit), 
#             data = epilepsy, family = poisson(), 
#             prior = set_prior("normal(0,5)"))
# summary(fit2) 
# plot(fit2)
# ## End(Not run)
 
Run the code above in your browser using DataLab