from arviz_base import load_arviz_data

import arviz_plots as azp

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

data = load_arviz_data("centered_eight")
pc = azp.plot_pair_focus(
    data,
    var_names=["theta","tau"],
    focus_var="mu",
    backend="none", # change to preferred backend
)
pc.show()