fit_e <- brm(count ~ log_Age_c + log_Base4_c * Trt_c + (1+Trt_c|visit),
data = epilepsy, family = "poisson", n.chains = 1)
## random effects means centered around zero with corresponding covariances
rf <- ranef(fit_e, var = TRUE)
attr(rf, "var")
## random effects medians centered around the corresponding fixed effect
ranef(fit_e, estimate = "median", center.zero = FALSE)
Run the code above in your browser using DataLab