fit <- brm(count ~ log_Age_c + log_Base4_c * Trt_c + (1+Trt_c|visit),
data = epilepsy, family = "poisson", chains = 1)
## group-level means with corresponding covariances
rf <- ranef(fit, var = TRUE)
attr(rf, "var")
## group-level medians
ranef(fit, estimate = "median")
Run the code above in your browser using DataLab