## Not run:
# fit <- brm(time ~ age * sex, data = kidney)
#
# lp <- log_posterior(fit)
# head(lp)
#
# np <- nuts_params(fit)
# str(np)
# # extract the number of divergence transitions
# sum(subset(np, Parameter == "divergent__")$Value)
#
# head(rhat(fit))
# head(neff_ratio(fit))
# ## End(Not run)
Run the code above in your browser using DataLab