Runs a series of convergence checks, prints statistics to the console, and makes traceplots of the worst converged parameter per selection.
# S3 method for emc
check(
emc,
selection = c("mu", "sigma2", "alpha"),
digits = 3,
plot_worst = TRUE,
...
)check(emc, ...)
a list with the statistics for the worst converged parameter per selection
An emc object
A Character vector. Indicates which parameter types to check (e.g., alpha, mu, sigma2, correlation).
Integer. How many digits to round the ESS and Rhat to in the plots
Boolean. If TRUE also plots the chain plots for the worst parameter
Optional arguments that can be passed to get_pars or plot.default (see par())
Note that the Rhat is calculated by doubling the number of chains by
first splitting chains into first and second half, so it also a test of
stationarity.
Efficiency of sampling is indicated by the effective
sample size (ESS) (from the coda R package).
Full range of possible samples manipulations described in get_pars.