This function plots a variety of diagnostics from a Variability Model.
These include a histogram of the Rhat values (so-called percent scale reduction
factors). An Rhat value of 1 indicates that no reduction in the variability of
the estimates is possible from running the chain longer. Values below 1.10 or 1.05
are typically considered indicative of convergence, with higher values indicating
the model did not converge and should be changed or run longer.
A histogram of the effective sample size indicates for every parameter estimated how
many effective posterior samples are available for inference. Low values may indicate
high autocorrelation in the samples and may be a sign of failure to converge.
The maximum possible will be the total iterations available.
Histograms of the posterior medians for the latent variability and intercept estimates
are also shown.
Usage
vm_diagnostics(object, plot = TRUE, ...)
Arguments
object
Results from running varian.
plot
Logical whether to plot the results or just return the grob
for the plots. Defaults to TRUE.