# 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))
# }
Run the code above in your browser using DataLab