Pass in rstanfit model object, and a threshold Rhat value for convergence. Returns boolean.

is_converged(fitted_model, threshold = 1.05, parameters = c("sigma", "x", "Z"))

Arguments

fitted_model

Samples extracted (with permuted = FALSE) from a Stan model. E.g. output from invert_chains().

threshold

Threshold for maximum Rhat.

parameters

Vector of parameters to be included in convergence determination. Defaults = c("sigma","x","Z"). Other elements can be added including "pred", "log_lik", or "lp__"