## Not run:
# fit <- brm(count ~ log_Age_c + log_Base4_c * Trt_c + (1+Trt_c|visit),
# data = epilepsy, family = "poisson", chains = 1)
# ## return the means of random effects covariances
# (vc <- VarCorr(fit))
# as.data.frame(vc)
#
# ## return 2.5% and 97.5% quantiles of random effects covariances
# VarCorr(fit, estimate = "quantile", probs = c(0.025, 0.975))
# ## End(Not run)
Run the code above in your browser using DataLab