from arviz_base import load_arviz_data
from arviz_stats import loo

import arviz_plots as azp

azp.style.use("arviz-variat")

dt = load_arviz_data("radon")
elpd_data = loo(dt, pointwise=True)

pc = azp.plot_khat(
    elpd_data,
    threshold=0.7,
    show_hlines=True,
    show_bins=True,
    backend="none",  # change to preferred backend
)

pc.show()