Skip to contents

Trace and rank plots for the worst-converging structural parameters

Source: R/plot_worst_offender_diagnostics.R
plot_probit_worst_offender_diagnostics.Rd

Auto-selects the union of the top n_worst parameters by Rhat and the top n_worst by lowest bulk ESS (from fit$diagnostics_detail$monitored_parameters), and shows their chain trace and rank-overlay plots with human-readable facet labels (covariate x class, or correlation pair) instead of raw Stan parameter names – enough to tell slow-mixing, chain separation, a sticky region, and multimodality apart, which a single Rhat number cannot.

Usage

plot_probit_worst_offender_diagnostics(
  fit,
  draws,
  title_base = "",
  n_worst = 6L
)

Arguments

fit

List returned by fit_bayesian_multivariate_probit().

draws

A posterior::draws_array (typically fit$draws).

title_base

Character. Prefixed to plot titles.

n_worst

Integer. How many parameters to select by each criterion (Rhat, bulk ESS) before deduplicating. Default 6.

Value

Named list with trace and rank ggplot objects, or NULL if no structural parameter diagnostics are available.