arviz 0.20.0__py3-none-any.whl → 0.21.0__py3-none-any.whl

arviz/preview.py

@@ -1,17 +1,48 @@
-# pylint: disable=unused-import,unused-wildcard-import,wildcard-import
+# pylint: disable=unused-import,unused-wildcard-import,wildcard-import,invalid-name
 """Expose features from arviz-xyz refactored packages inside ``arviz.preview`` namespace."""
+import logging
+
+_log = logging.getLogger(__name__)
+
+info = ""
 
 try:
     from arviz_base import *
+
+    status = "arviz_base available, exposing its functions as part of arviz.preview"
+    _log.info(status)
 except ModuleNotFoundError:
-    pass
+    status = "arviz_base not installed"
+    _log.info(status)
+except ImportError:
+    status = "Unable to import arviz_base"
+    _log.info(status, exc_info=True)
+
+info += status + "\n"
 
 try:
-    import arviz_stats
+    from arviz_stats import *
+
+    status = "arviz_stats available, exposing its functions as part of arviz.preview"
+    _log.info(status)
 except ModuleNotFoundError:
-    pass
+    status = "arviz_stats not installed"
+    _log.info(status)
+except ImportError:
+    status = "Unable to import arviz_stats"
+    _log.info(status, exc_info=True)
+info += status + "\n"
 
 try:
     from arviz_plots import *
+
+    status = "arviz_plots available, exposing its functions as part of arviz.preview"
+    _log.info(status)
 except ModuleNotFoundError:
-    pass
+    status = "arviz_plots not installed"
+    _log.info(status)
+except ImportError:
+    status = "Unable to import arviz_plots"
+    _log.info(status, exc_info=True)
+
+info += status + "\n"

arviz/stats/__init__.py

@@ -9,6 +9,7 @@ from .stats_utils import *
 
 __all__ = [
     "apply_test_function",
+    "bayes_factor",
     "bfmi",
     "compare",
     "hdi",

arviz/stats/diagnostics.py

@@ -744,8 +744,8 @@ def _ess_sd(ary, relative=False):
     ary = np.asarray(ary)
     if _not_valid(ary, shape_kwargs=dict(min_draws=4, min_chains=1)):
         return np.nan
-    ary = _split_chains(ary)
-    return min(_ess(ary, relative=relative), _ess(ary**2, relative=relative))
+    ary = (ary - ary.mean()) ** 2
+    return _ess(_split_chains(ary), relative=relative)
 
 
 def _ess_quantile(ary, prob, relative=False):
@@ -838,13 +838,15 @@ def _mcse_sd(ary):
     ary = np.asarray(ary)
     if _not_valid(ary, shape_kwargs=dict(min_draws=4, min_chains=1)):
         return np.nan
-    ess = _ess_sd(ary)
+    sims_c2 = (ary - ary.mean()) ** 2
+    ess = _ess_mean(sims_c2)
+    evar = (sims_c2).mean()
+    varvar = ((sims_c2**2).mean() - evar**2) / ess
+    varsd = varvar / evar / 4
     if _numba_flag:
-        sd = float(_sqrt(svar(np.ravel(ary), ddof=1), np.zeros(1)).item())
+        mcse_sd_value = float(_sqrt(np.ravel(varsd), np.zeros(1)))
     else:
-        sd = np.std(ary, ddof=1)
-    fac_mcse_sd = np.sqrt(np.exp(1) * (1 - 1 / ess) ** (ess - 1) - 1)
-    mcse_sd_value = sd * fac_mcse_sd
+        mcse_sd_value = np.sqrt(varsd)
     return mcse_sd_value
 
 
@@ -973,19 +975,21 @@ def _multichain_statistics(ary, focus="mean"):
         # ess mean
         ess_mean_value = _ess_mean(ary)
 
-        # ess sd
-        ess_sd_value = _ess_sd(ary)
-
         # mcse_mean
-        sd = np.std(ary, ddof=1)
-        mcse_mean_value = sd / np.sqrt(ess_mean_value)
+        sims_c2 = (ary - ary.mean()) ** 2
+        sims_c2_sum = sims_c2.sum()
+        var = sims_c2_sum / (sims_c2.size - 1)
+        mcse_mean_value = np.sqrt(var / ess_mean_value)
 
         # ess bulk
         ess_bulk_value = _ess(z_split)
 
         # mcse_sd
-        fac_mcse_sd = np.sqrt(np.exp(1) * (1 - 1 / ess_sd_value) ** (ess_sd_value - 1) - 1)
-        mcse_sd_value = sd * fac_mcse_sd
+        evar = sims_c2_sum / sims_c2.size
+        ess_mean_sims = _ess_mean(sims_c2)
+        varvar = ((sims_c2**2).mean() - evar**2) / ess_mean_sims
+        varsd = varvar / evar / 4
+        mcse_sd_value = np.sqrt(varsd)
 
         return (
             mcse_mean_value,

arviz/stats/stats.py

@@ -27,6 +27,7 @@ from ..utils import Numba, _numba_var, _var_names, get_coords
 from .density_utils import get_bins as _get_bins
 from .density_utils import histogram as _histogram
 from .density_utils import kde as _kde
+from .density_utils import _kde_linear
 from .diagnostics import _mc_error, _multichain_statistics, ess
 from .stats_utils import ELPDData, _circular_standard_deviation, smooth_data
 from .stats_utils import get_log_likelihood as _get_log_likelihood
@@ -41,6 +42,7 @@ from ..labels import BaseLabeller
 
 __all__ = [
     "apply_test_function",
+    "bayes_factor",
     "compare",
     "hdi",
     "loo",
@@ -2337,3 +2339,87 @@ def _cjs_dist(draws, weights):
     bound = cdf_p_int + cdf_q_int
 
     return np.sqrt((cjs_pq + cjs_qp) / bound)
+
+
+def bayes_factor(idata, var_name, ref_val=0, prior=None, return_ref_vals=False):
+    r"""Approximated Bayes Factor for comparing hypothesis of two nested models.
+
+    The Bayes factor is estimated by comparing a model (H1) against a model in which the
+    parameter of interest has been restricted to be a point-null (H0). This computation
+    assumes the models are nested and thus H0 is a special case of H1.
+
+    Notes
+    -----
+    The bayes Factor is approximated as the Savage-Dickey density ratio
+    algorithm presented in [1]_.
+
+    Parameters
+    ----------
+    idata : InferenceData
+        Any object that can be converted to an :class:`arviz.InferenceData` object
+        Refer to documentation of :func:`arviz.convert_to_dataset` for details.
+    var_name : str, optional
+        Name of variable we want to test.
+    ref_val : int, default 0
+        Point-null for Bayes factor estimation.
+    prior : numpy.array, optional
+        In case we want to use different prior, for example for sensitivity analysis.
+    return_ref_vals : bool, optional
+        Whether to return the values of the prior and posterior at the reference value.
+        Used by :func:`arviz.plot_bf` to display the distribution comparison.
+
+
+    Returns
+    -------
+    dict : A dictionary with BF10 (Bayes Factor 10 (H1/H0 ratio), and BF01 (H0/H1 ratio).
+
+    References
+    ----------
+    .. [1] Heck, D., 2019. A caveat on the Savage-Dickey density ratio:
+       The case of computing Bayes factors for regression parameters.
+
+    Examples
+    --------
+    Moderate evidence indicating that the parameter "a" is different from zero.
+
+    .. ipython::
+
+        In [1]: import numpy as np
+           ...: import arviz as az
+           ...: idata = az.from_dict(posterior={"a":np.random.normal(1, 0.5, 5000)},
+           ...: prior={"a":np.random.normal(0, 1, 5000)})
+           ...: az.bayes_factor(idata, var_name="a", ref_val=0)
+
+    """
+
+    posterior = extract(idata, var_names=var_name).values
+
+    if ref_val > posterior.max() or ref_val < posterior.min():
+        _log.warning(
+            "The reference value is outside of the posterior. "
+            "This translate into infinite support for H1, which is most likely an overstatement."
+        )
+
+    if posterior.ndim > 1:
+        _log.warning("Posterior distribution has {posterior.ndim} dimensions")
+
+    if prior is None:
+        prior = extract(idata, var_names=var_name, group="prior").values
+
+    if posterior.dtype.kind == "f":
+        posterior_grid, posterior_pdf, *_ = _kde_linear(posterior)
+        prior_grid, prior_pdf, *_ = _kde_linear(prior)
+        posterior_at_ref_val = np.interp(ref_val, posterior_grid, posterior_pdf)
+        prior_at_ref_val = np.interp(ref_val, prior_grid, prior_pdf)
+
+    elif posterior.dtype.kind == "i":
+        posterior_at_ref_val = (posterior == ref_val).mean()
+        prior_at_ref_val = (prior == ref_val).mean()
+
+    bf_10 = prior_at_ref_val / posterior_at_ref_val
+    bf = {"BF10": bf_10, "BF01": 1 / bf_10}
+
+    if return_ref_vals:
+        return (bf, {"prior": prior_at_ref_val, "posterior": posterior_at_ref_val})
+    else:
+        return bf

arviz/tests/base_tests/test_data.py

@@ -3,6 +3,7 @@
 
 import importlib
 import os
+import warnings
 from collections import namedtuple
 from copy import deepcopy
 from html import escape
@@ -32,7 +33,13 @@ from ... import (
     extract,
 )
 
-from ...data.base import dict_to_dataset, generate_dims_coords, infer_stan_dtypes, make_attrs
+from ...data.base import (
+    dict_to_dataset,
+    generate_dims_coords,
+    infer_stan_dtypes,
+    make_attrs,
+    numpy_to_data_array,
+)
 from ...data.datasets import LOCAL_DATASETS, REMOTE_DATASETS, RemoteFileMetadata
 from ..helpers import (  # pylint: disable=unused-import
     chains,
@@ -231,6 +238,17 @@ def test_dims_coords_skip_event_dims(shape):
     assert "z" not in coords
 
 
+@pytest.mark.parametrize("dims", [None, ["chain", "draw"], ["chain", "draw", None]])
+def test_numpy_to_data_array_with_dims(dims):
+    da = numpy_to_data_array(
+        np.empty((4, 500, 7)),
+        var_name="a",
+        dims=dims,
+        default_dims=["chain", "draw"],
+    )
+    assert list(da.dims) == ["chain", "draw", "a_dim_0"]
+
+
 def test_make_attrs():
     extra_attrs = {"key": "Value"}
     attrs = make_attrs(attrs=extra_attrs)
@@ -921,7 +939,7 @@ class TestInferenceData:  # pylint: disable=too-many-public-methods
         data = np.random.normal(size=(4, 500, 8))
         idata = data_random
         with pytest.warns(UserWarning, match="The group.+not defined in the InferenceData scheme"):
-            idata.add_groups({"new_group": idata.posterior})
+            idata.add_groups({"new_group": idata.posterior}, warn_on_custom_groups=True)
         with pytest.warns(UserWarning, match="the default dims.+will be added automatically"):
             idata.add_groups(constant_data={"a": data[..., 0], "b": data})
         assert idata.new_group.equals(idata.posterior)
@@ -962,8 +980,8 @@ class TestInferenceData:  # pylint: disable=too-many-public-methods
         with pytest.raises(ValueError, match="join must be either"):
             idata.extend(idata2, join="outer")
         idata2.add_groups(new_group=idata2.prior)
-        with pytest.warns(UserWarning):
-            idata.extend(idata2)
+        with pytest.warns(UserWarning, match="new_group"):
+            idata.extend(idata2, warn_on_custom_groups=True)
 
 
 class TestNumpyToDataArray:
@@ -1156,11 +1174,17 @@ def test_bad_inference_data():
         InferenceData(posterior=[1, 2, 3])
 
 
-def test_inference_data_other_groups():
+@pytest.mark.parametrize("warn", [True, False])
+def test_inference_data_other_groups(warn):
     datadict = {"a": np.random.randn(100), "b": np.random.randn(1, 100, 10)}
     dataset = convert_to_dataset(datadict, coords={"c": np.arange(10)}, dims={"b": ["c"]})
-    with pytest.warns(UserWarning, match="not.+in.+InferenceData scheme"):
-        idata = InferenceData(other_group=dataset)
+    if warn:
+        with pytest.warns(UserWarning, match="not.+in.+InferenceData scheme"):
+            idata = InferenceData(other_group=dataset, warn_on_custom_groups=True)
+    else:
+        with warnings.catch_warnings():
+            warnings.simplefilter("error")
+            idata = InferenceData(other_group=dataset, warn_on_custom_groups=False)
     fails = check_multiple_attrs({"other_group": ["a", "b"]}, idata)
     assert not fails
 

arviz/tests/base_tests/test_diagnostics.py

@@ -120,10 +120,11 @@ class TestDiagnostics:
         ```
         Reference file:
 
-        Created: 2020-08-31
-        System: Ubuntu 18.04.5 LTS
-        R version 4.0.2 (2020-06-22)
-        posterior 0.1.2
+        Created: 2024-12-20
+        System: Ubuntu 24.04.1 LTS
+        R version 4.4.2 (2024-10-31)
+        posterior version from https://github.com/stan-dev/posterior/pull/388
+            (after release 1.6.0 but before the fixes in the PR were released).
         """
         # download input files
         here = os.path.dirname(os.path.abspath(__file__))

arviz/tests/base_tests/test_plots_matplotlib.py

@@ -2,21 +2,22 @@
 
 # pylint: disable=redefined-outer-name,too-many-lines
 import os
+import re
 from copy import deepcopy
 
 import matplotlib.pyplot as plt
 import numpy as np
 import pytest
+import xarray as xr
 from matplotlib import animation
 from pandas import DataFrame
 from scipy.stats import gaussian_kde, norm
-import xarray as xr
 
 from ...data import from_dict, load_arviz_data
 from ...plots import (
     plot_autocorr,
-    plot_bpv,
     plot_bf,
+    plot_bpv,
     plot_compare,
     plot_density,
     plot_dist,
@@ -43,20 +44,20 @@ from ...plots import (
     plot_ts,
     plot_violin,
 )
+from ...plots.dotplot import wilkinson_algorithm
+from ...plots.plot_utils import plot_point_interval
 from ...rcparams import rc_context, rcParams
 from ...stats import compare, hdi, loo, waic
 from ...stats.density_utils import kde as _kde
-from ...utils import _cov, BehaviourChangeWarning
-from ...plots.plot_utils import plot_point_interval
-from ...plots.dotplot import wilkinson_algorithm
+from ...utils import BehaviourChangeWarning, _cov
 from ..helpers import (  # pylint: disable=unused-import
+    RandomVariableTestClass,
     create_model,
     create_multidimensional_model,
     does_not_warn,
     eight_schools_params,
     models,
     multidim_models,
-    RandomVariableTestClass,
 )
 
 rcParams["data.load"] = "eager"
@@ -1236,6 +1237,23 @@ def test_plot_hdi_dataset_error(models):
         plot_hdi(np.arange(8), hdi_data=hdi_data)
 
 
+def test_plot_hdi_string_error():
+    """Check x as type string raises an error."""
+    x_data = ["a", "b", "c", "d"]
+    y_data = np.random.normal(0, 5, (1, 200, len(x_data)))
+    hdi_data = hdi(y_data)
+    with pytest.raises(
+        NotImplementedError,
+        match=re.escape(
+            (
+                "The `arviz.plot_hdi()` function does not support categorical data. "
+                "Consider using `arviz.plot_forest()`."
+            )
+        ),
+    ):
+        plot_hdi(x=x_data, y=y_data, hdi_data=hdi_data)
+
+
 def test_plot_hdi_datetime_error():
     """Check x as datetime raises an error."""
     x_data = np.arange(start="2022-01-01", stop="2022-03-01", dtype=np.datetime64)
@@ -2082,7 +2100,5 @@ def test_plot_bf():
     idata = from_dict(
         posterior={"a": np.random.normal(1, 0.5, 5000)}, prior={"a": np.random.normal(0, 1, 5000)}
     )
-    bf_dict0, _ = plot_bf(idata, var_name="a", ref_val=0)
-    bf_dict1, _ = plot_bf(idata, prior=np.random.normal(0, 10, 5000), var_name="a", ref_val=0)
-    assert bf_dict0["BF10"] > bf_dict0["BF01"]
-    assert bf_dict1["BF10"] < bf_dict1["BF01"]
+    _, bf_plot = plot_bf(idata, var_name="a", ref_val=0)
+    assert bf_plot is not None

arviz/tests/base_tests/test_stats.py

@@ -18,6 +18,7 @@ from ...data import concat, convert_to_inference_data, from_dict, load_arviz_dat
 from ...rcparams import rcParams
 from ...stats import (
     apply_test_function,
+    bayes_factor,
     compare,
     ess,
     hdi,
@@ -871,3 +872,13 @@ def test_priorsens_coords(psens_data):
     assert "mu" in result
     assert "theta" in result
     assert "school" in result.theta_t.dims
+
+
+def test_bayes_factor():
+    idata = from_dict(
+        posterior={"a": np.random.normal(1, 0.5, 5000)}, prior={"a": np.random.normal(0, 1, 5000)}
+    )
+    bf_dict0 = bayes_factor(idata, var_name="a", ref_val=0)
+    bf_dict1 = bayes_factor(idata, prior=np.random.normal(0, 10, 5000), var_name="a", ref_val=0)
+    assert bf_dict0["BF10"] > bf_dict0["BF01"]
+    assert bf_dict1["BF10"] < bf_dict1["BF01"]

arviz/utils.py

@@ -330,6 +330,7 @@ class Numba:
     """A class to toggle numba states."""
 
     numba_flag = numba_check()
+    """bool: Indicates whether Numba optimizations are enabled. Defaults to False."""
 
     @classmethod
     def disable_numba(cls):
@@ -732,7 +733,10 @@ class Dask:
     """
 
     dask_flag = False
+    """bool: Enables Dask parallelization when set to True. Defaults to False."""
     dask_kwargs = None
+    """dict: Additional keyword arguments for Dask configuration.
+    Defaults to an empty dictionary."""
 
     @classmethod
     def enable_dask(cls, dask_kwargs=None):

{arviz-0.20.0.dist-info → arviz-0.21.0.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: arviz
-Version: 0.20.0
+Version: 0.21.0
 Summary: Exploratory analysis of Bayesian models
 Home-page: http://github.com/arviz-devs/arviz
 Author: ArviZ Developers

{arviz-0.20.0.dist-info → arviz-0.21.0.dist-info}/RECORD

@@ -1,15 +1,15 @@
-arviz/__init__.py,sha256=cZ-_ahEB5206wiaQaaTe3LjM8VIC04GAkV-noExc2kw,10397
+arviz/__init__.py,sha256=1PWUe3oeouP_C8JmAXhMrHarplymbqLyWqm950C-aAI,10397
 arviz/labels.py,sha256=w4-t0qdJzjKrqRyhzbtk6ucoMIAxle1HpHYlH7up06Q,6828
-arviz/preview.py,sha256=5HwHycvbSAkWsAEYHBfd3Crwmeq8T06ysSgB528iQnA,375
+arviz/preview.py,sha256=Fmff8j9Zlvgi5w2PRwnbkEOioJY8fK9p1SWQWjTl4N8,1314
 arviz/py.typed,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 arviz/rcparams.py,sha256=lQnBprbDihEcsP5FujoJetGna4T8nc20JnuVf1Wn-1k,21031
 arviz/sel_utils.py,sha256=xvAYENhGXDTrhaT4Itlk1SJQUUGZ6BGcR04fPFgvzdM,6951
-arviz/utils.py,sha256=sOUIkqmWtwY3PM-rQpN6arW9MsPt-MfvzdwAe78rT8U,26425
+arviz/utils.py,sha256=-q3eAeficZcLakOt6UXN0eyVx6uIVVpkqA_ky8FRSjE,26699
 arviz/data/__init__.py,sha256=SG2umdZ8uVNYKVBoVYPy5tNxJnzLdyw0spPMjMTae5k,1558
-arviz/data/base.py,sha256=ptyUVyp57mVq0QVrPRYzRT1EAisSYZqZucq5_V_ZdX8,21369
+arviz/data/base.py,sha256=PALdVidyCxJqje_za4XPwXH010qAMJt4VzzeOflUCC8,21365
 arviz/data/converters.py,sha256=4QU5_0OU84HV3gSXXY_N5bd6S9ft81lfpnCT-8E8qyc,7933
 arviz/data/datasets.py,sha256=wPi23KZI-w4HrhjmY1LUg0Rj0aJobSJ4WO3LBrjfZQc,5392
-arviz/data/inference_data.py,sha256=3hVaL26YmWeblrZhVKmnaKON3Itn_EsGNssEgSQiEuc,92745
+arviz/data/inference_data.py,sha256=GYmQzz26c614IX_XVvmoTA0ycaC2pZs3B48hd2HYhrA,93619
 arviz/data/io_beanmachine.py,sha256=QQVBD6rftvs6_kLIb4Vm1QzQ6BsS0J9DTrzw2Jj4ob8,3745
 arviz/data/io_cmdstan.py,sha256=8YX9Nfkx4LqjfKms4s4GTOkOjZNelb6SukvRZRHY6iM,38994
 arviz/data/io_cmdstanpy.py,sha256=iSr8ciKBFoIa1tJGHEX-2JKkUJRyaTXzRXf-5mu8q5U,42991
@@ -31,8 +31,8 @@ arviz/data/example_data/data/centered_eight.nc,sha256=jvw6uv4MeW65rqe2lJDU4kAKM8
 arviz/data/example_data/data/non_centered_eight.nc,sha256=r7kyd10HyJTTRQs4OlSCXPVt3T-nLsPd3g-bcPYnPmA,836647
 arviz/plots/__init__.py,sha256=atWhfLXHAD6zaTjbdWdTaJrPTNBl3XpcSCY4etgw_cY,1514
 arviz/plots/autocorrplot.py,sha256=f-rYtnZtxgOns5sIT6pdw-EdU74lAVmAeV9b6rKLsgw,5666
-arviz/plots/bfplot.py,sha256=s8bizLCkF5uZnz66UKmuUJs_AXuMVrgGObTa2ofUF-Q,5124
-arviz/plots/bpvplot.py,sha256=QR4znezjiSLMyMagvnJLb9Y1SIlT99gJ52o4FnX5tZg,12098
+arviz/plots/bfplot.py,sha256=TKCkk60dgIk70CNWz9pHDXb1HGwHT_aLfNgFH6jDz9c,4367
+arviz/plots/bpvplot.py,sha256=Pf0ME21_teSQd_7CcuYE2MxfDJHcy5eW7Ms7yIDZSCA,12505
 arviz/plots/compareplot.py,sha256=Z8usSMEeQKs4GmkziDR4dVzSh3Ocd4ySfiNDZVaFOUc,6078
 arviz/plots/densityplot.py,sha256=6477ZljpBCcZRw0SUwcTO4FYjxqw_qYsJupWNo-jCok,10895
 arviz/plots/distcomparisonplot.py,sha256=gVNQUN0VX7hC527fcUk1oxtQRdIl5mrltU95c0Nra9k,7184
@@ -43,7 +43,7 @@ arviz/plots/elpdplot.py,sha256=NKqPkwTj9BWDzwMnG0cPeLmYBGMX_meP9F6bqvTwLKY,6433
 arviz/plots/energyplot.py,sha256=znEDPYpWaTIX0XpdVoyhXOITJ4A8BYkI9t1TVhJq4Qo,4797
 arviz/plots/essplot.py,sha256=ch0DjUQDILk4ohpSUR-9VHejGFb6Xmelly-qa-fKb9E,11741
 arviz/plots/forestplot.py,sha256=37Wj4wFGjydZS4SpdqZ9qLxhZBo4rFxV_MWQnZAS7DA,11897
-arviz/plots/hdiplot.py,sha256=cNj2r0dPxtquZNoChgjFO-wmqgU-W-dq_ed1kWJ2vHI,7589
+arviz/plots/hdiplot.py,sha256=Pii9ZsuejEM-I24dn39muUL-yYKTfe2RWzAuU0W-3SI,7798
 arviz/plots/kdeplot.py,sha256=55mQPyIrvPKeJbjTsqWTExcff6QPvXFXG47XtDfPs_Y,11924
 arviz/plots/khatplot.py,sha256=u03gmBG1xwxG07ASLZEJB-GsRhRHtUCtbglpC1N0aIg,8086
 arviz/plots/lmplot.py,sha256=LxR7RXkaAi5J8076isebVrtdk6UwbcTRekEymM9S6cY,11726
@@ -51,7 +51,7 @@ arviz/plots/loopitplot.py,sha256=bFUO_Fy4z6u6E39GdaF4rIvc-COWNwF4A0n2kcmZBfA,832
 arviz/plots/mcseplot.py,sha256=rsiz4E9M9p58YetAaF75gbenGIj4M0hapWnh9UJOXzY,6829
 arviz/plots/pairplot.py,sha256=yVxyLC7ms0SKpTunifacUyOblH8BW9CsQqE2Hv0ARP4,10407
 arviz/plots/parallelplot.py,sha256=ZBEsHvnlmSXLRpSwP-KUwzgWBC2S4siKXFGJnLf7uAE,7125
-arviz/plots/plot_utils.py,sha256=mawuYnwb9AfdLQ9HAerTOGT9iKNjXrQraXvtDd9tFco,18243
+arviz/plots/plot_utils.py,sha256=VyVR50HrZegdkWa6ZxtRnC_WJstooYvaB-xsDHf6kaQ,18337
 arviz/plots/posteriorplot.py,sha256=pC-5SQZOGq1F4opM_sQLxn4ZG2we4Y9ULV8yhxjGVdo,10952
 arviz/plots/ppcplot.py,sha256=UPTtXDWHf3wFAb-apNPGcz8qw9CQwINGml_2YkYI-iM,13967
 arviz/plots/rankplot.py,sha256=lz0swHs6EBe-gXn4udP1Um3RS-EatsOAmguYqGMlIjU,8648
@@ -59,39 +59,39 @@ arviz/plots/separationplot.py,sha256=Fx_QVeFUcF45fm7nn06pt0qubOzvH8QMU1cw5RLyaik
 arviz/plots/traceplot.py,sha256=dwcF7rsjMAIxZ_LPv7Z8os01uQZHXTkDFWEBtsbzI9k,10216
 arviz/plots/tsplot.py,sha256=haTyvfGX5fA8Zle9bzllybG5n307BUJIxGywNAnOsU0,15925
 arviz/plots/violinplot.py,sha256=yxoEMGTIt4CDinZaNHPYI5MqFvXB2J2gyKdKJ47PKdk,7129
-arviz/plots/backends/__init__.py,sha256=LZxXo7ogt7ZVrdUea0xkxpx5YV0wpV20MHwdItgk3oU,7752
+arviz/plots/backends/__init__.py,sha256=LQbEHNuGuD2CRNmO-Djnfuk4kMF58guFGjEYr68mAxI,7867
 arviz/plots/backends/bokeh/__init__.py,sha256=e2wfZNdGTFU5GjsLokCLpknweaNgjZ5v3k7NB0gry6g,4877
 arviz/plots/backends/bokeh/autocorrplot.py,sha256=9CruVndfMnqsi4waav7Gcr7pn0nmUB0mkJxkzOjtGtw,2463
 arviz/plots/backends/bokeh/bfplot.py,sha256=ydjomuA5iTw9LE2_eq9_u6Ox3MCy0ulE_DU4qgH0MO8,406
-arviz/plots/backends/bokeh/bpvplot.py,sha256=btwA4qbO8iCtZ-vriscqNlwn0zEooObh7DD6-fsyM44,6642
+arviz/plots/backends/bokeh/bpvplot.py,sha256=FItG_h39UeoZX9nNE7GfM95cvT2gstlf374b2XuCEsM,6673
 arviz/plots/backends/bokeh/compareplot.py,sha256=52jOjhrRqB-Dgdw8odKfM4aPXk9med6MopCiHxKNgJM,4736
-arviz/plots/backends/bokeh/densityplot.py,sha256=TQ5OkHHnihT1NBa4vtt6wJUsll_K3xaeEFzDtqrmcVU,6286
+arviz/plots/backends/bokeh/densityplot.py,sha256=LHkRpz7QJDobpBgw6Wuy23nFr8Ruc2CUBUDVwbj5A9c,6356
 arviz/plots/backends/bokeh/distcomparisonplot.py,sha256=o8FHMb1ZzKPpt7fXhwBr6HGhqpclO1Qk9o6aTGypgv0,431
 arviz/plots/backends/bokeh/distplot.py,sha256=a2yY4waIPdwGhDpUtYi87Ra-TJiAA67oQumIiU-nXiA,4851
-arviz/plots/backends/bokeh/dotplot.py,sha256=qy8B4QL6B_ZfoVtDDA6c7HcsQofQLlE1tTnkcJO1mPA,2828
+arviz/plots/backends/bokeh/dotplot.py,sha256=xnbaCSyGjCSfdiIIRf5zuMeogfNdNPK-1LenasjTAa8,2923
 arviz/plots/backends/bokeh/ecdfplot.py,sha256=Zrinhu6ViG3BbClxnRFZ_zIZm7S1nARiydHi-l1_i-g,1680
 arviz/plots/backends/bokeh/elpdplot.py,sha256=5bn_rH1Aixm0--BArAP4m4kuZMjWxR8ox-8T-xK_GMY,6487
 arviz/plots/backends/bokeh/energyplot.py,sha256=FKPYRaWwM32Vw0AHE0MIWBn9wZQeoYdP_YpYZGijIq0,4571
-arviz/plots/backends/bokeh/essplot.py,sha256=zghMEYUH3bAljqNt72-dRa8IFom27F7dMklHNQOrzBM,5494
-arviz/plots/backends/bokeh/forestplot.py,sha256=DIoGzSIAMA9vmKwpoNxVoNKFpDBwpfjdGd6KTl9u_OY,27219
+arviz/plots/backends/bokeh/essplot.py,sha256=0i2E1TvlzzP9mTavlGlJ5hPoFW5tGWCkAPAKrVHbgYQ,5560
+arviz/plots/backends/bokeh/forestplot.py,sha256=TAaXMqF2F1G7BIeDdWas14C6pCKZ3RTNhLIjUm-_IQs,27466
 arviz/plots/backends/bokeh/hdiplot.py,sha256=bAhTPi9D7cw2ytSvjSIRD3g-PqpB9OFwusrQAnmKCYY,1538
 arviz/plots/backends/bokeh/kdeplot.py,sha256=nKEgJfnP7NK2Y2cipF_RglEC6GpowaagH5wZP3iZq-U,9329
-arviz/plots/backends/bokeh/khatplot.py,sha256=J5m9vugSzpi9sE-T_2PRzI8pA100Hf8rUMhpmkTP8n4,4599
-arviz/plots/backends/bokeh/lmplot.py,sha256=kqfkQ1HAahi_XwkWxrub73zPM_4dXEFvKAd0TFKsoXQ,5392
+arviz/plots/backends/bokeh/khatplot.py,sha256=Iz2C6YaQOhnucAFbjSWTha2HNLFmctXc4tt_rLdS6Ko,4663
+arviz/plots/backends/bokeh/lmplot.py,sha256=0VIu1CPngKFQZayOqASAYHfXU8s_SRzysjdmxkELRJo,5600
 arviz/plots/backends/bokeh/loopitplot.py,sha256=FWjcsSWGJNy4wM63_N2qpg5oECb3Cq-uLy3xDS5x6j4,7172
-arviz/plots/backends/bokeh/mcseplot.py,sha256=teyAG-A4vIHor_Qr4iB3l5nLJ--EzvPvw1x44hUOVSI,5924
-arviz/plots/backends/bokeh/pairplot.py,sha256=2sByGaDuDT_ebyhPLIP3jRO1i43DLVDTPKWCpjd9MXc,13146
+arviz/plots/backends/bokeh/mcseplot.py,sha256=QHyXeANcVSeLRdPJGh5RihZmaYq8cAPjSAn6FGJXQRk,5960
+arviz/plots/backends/bokeh/pairplot.py,sha256=kNvB3pVcpWXz7GMjBKc0pNyJdbx7VYgyEQOZeFbVgKo,13416
 arviz/plots/backends/bokeh/parallelplot.py,sha256=SNChOLWvcKxXuuJsIfWs9CNj7qDuVb95UZyp4CP1BQE,2230
 arviz/plots/backends/bokeh/posteriorplot.py,sha256=yBAlGo3lQpFBVro724eZtKnOHgVkY7LXupa-8_VtE7c,9380
-arviz/plots/backends/bokeh/ppcplot.py,sha256=mfDzTXxJU26CcK8ZKno3Y_yK_y3VYeUT3vaJqUoGEq4,13269
-arviz/plots/backends/bokeh/rankplot.py,sha256=UvzzzhOEnHTX4BkO30JO7iKyU0dwbq7kqo0rV6LlXEU,4474
+arviz/plots/backends/bokeh/ppcplot.py,sha256=fDeXXn1WtjjSs-285nDo76NKJTp2m5hIZ7V9S6iqA30,13307
+arviz/plots/backends/bokeh/rankplot.py,sha256=cRTAU-cWt__0J4CzsHHSf47SA_vB1n4mQ16XIfve7aE,4531
 arviz/plots/backends/bokeh/separationplot.py,sha256=_VhJUnPW03xVpy2y4h8npLyPPD49H76WQMmj3hOq2Hc,2416
-arviz/plots/backends/bokeh/traceplot.py,sha256=mHyJQqj0IzwKzSsK0eJ7u05V-WObZCJKuLbJ9EpdPtM,14208
-arviz/plots/backends/bokeh/violinplot.py,sha256=_GaLDsq6muqnDr1BQCFClCGqy7hyEk7moYpMd0MNwNA,4359
+arviz/plots/backends/bokeh/traceplot.py,sha256=mnw2rBoLNVrWYDR-m0HyECHFD1-dlsS18gj6D5qGGQ4,14246
+arviz/plots/backends/bokeh/violinplot.py,sha256=lfAzkv6qABLSz5R4uzFUotVExnG_c8HP3vkT_7Yc1dQ,4389
 arviz/plots/backends/matplotlib/__init__.py,sha256=LBEWakXN4QFoIXp_aPXPMTzXnA8VJt24k5RaowPCQoY,3629
 arviz/plots/backends/matplotlib/autocorrplot.py,sha256=ahyNnwyNrLubZTsxNvTi5hAZ5gV9dmo_wecrTYyMMDI,1807
 arviz/plots/backends/matplotlib/bfplot.py,sha256=00-xGO_VpmTxCkYiC1cGGsAW0HNO3bhQJkFxA6ssMh0,1828
-arviz/plots/backends/matplotlib/bpvplot.py,sha256=SeAGSOxhUFMIlQcluVnGEF5SnW2qjg1HGx4St-aiFz8,6262
+arviz/plots/backends/matplotlib/bpvplot.py,sha256=EJ79XptBAPrERh59wg0PCkKliIS1oPgcdzMkLmSFcno,6301
 arviz/plots/backends/matplotlib/compareplot.py,sha256=qvjSkNDBTZh3vxIw4pNEyRnaAM8qFG9vWw0X23hscsI,3695
 arviz/plots/backends/matplotlib/densityplot.py,sha256=zlqzYvH3VXWPiiIvL--fxtMp3hIGrcttKkQxKPV1U2s,5479
 arviz/plots/backends/matplotlib/distcomparisonplot.py,sha256=XZY2jITNKtcIMsg5tl_lzuwI-2DcdUdCWqrK7bsdvWE,3568
@@ -136,28 +136,28 @@ arviz/plots/styles/arviz-white.mplstyle,sha256=p3dbvWzOKhA-u8r3BmTF-bR5bhKh87iIk
 arviz/plots/styles/arviz-whitegrid.mplstyle,sha256=IMjjlfG3wg7heUjcVrkez1SNoiMI6BLztGPmsUp1iws,1072
 arviz/static/css/style.css,sha256=wcC7rvCT4E6TycEiw7YqxwyaaZ4tTRDGqMkYYjAqrao,5910
 arviz/static/html/icons-svg-inline.html,sha256=t-ChbtS1Gv8uZxc31DCJS8SuXDsLGUHoKgwv8zu6j2M,1343
-arviz/stats/__init__.py,sha256=jWXBXngHmYFy8m_3QJKYRvLszI4L5Q1aIBS79PC9Gms,700
+arviz/stats/__init__.py,sha256=kvrANzMkqyHMTdry7N5w836E2OP0tJM6bm5-G8OZaA0,721
 arviz/stats/density_utils.py,sha256=WCEkXCQsaycY_usA6xkorhIqr3_5ru7wgwE9eJJgdH8,32216
-arviz/stats/diagnostics.py,sha256=qGIrq258E0vl_tJwZSCvmp0dbQxvysrdJ8zzIsaK50A,32501
+arviz/stats/diagnostics.py,sha256=COTy2c5ROAirCAK_UNo7kQnggN71maBRPwy54ZdabKE,32656
 arviz/stats/ecdf_utils.py,sha256=Wy38wL-bsHDlpyU9qnDjedYBvbP_6ZrzJuWo6NzD2Xg,11835
-arviz/stats/stats.py,sha256=_4_L4_dOmyeWploQMQuUPjPFzjC3DMEetn_RIl3IP8c,86678
+arviz/stats/stats.py,sha256=0cpdhmTwEVNY6VedeYhl2D3zlXj2hn2Mp2yqVkufP94,89884
 arviz/stats/stats_refitting.py,sha256=trbPC7LCnsb-n5D6g7J0bzXJCPfcDukJDniB4ud1z9E,5415
 arviz/stats/stats_utils.py,sha256=XG8ILPVs8Jbh_v7jzLfwMkm2HraT2j2-Hxe_kEYlLjQ,20076
 arviz/tests/__init__.py,sha256=TiS6C1IzwAXmNa8u36Y2xzL1CTTZm2PwzAtmZgoqepE,18
 arviz/tests/conftest.py,sha256=6U9WpKmYf38EVRoFZNBpV0CunQvESBFJG2SJ8IBEkL4,1270
 arviz/tests/helpers.py,sha256=qhsOhLtfyz-dC2yuT6ug0frYZlbims06BljJuEVDP6E,23593
 arviz/tests/base_tests/__init__.py,sha256=zg7V5_0DZrCz7ozETqET6bUpAvUUmmkSoLpJLwaIj2E,23
-arviz/tests/base_tests/test_data.py,sha256=QTKeQtTflU5TSnGk-rFLV9OQhmvoQ7RT962MLpu2YXU,63543
+arviz/tests/base_tests/test_data.py,sha256=L3vfscLm8aMVjd3MswJqmJ_J7Y2sNjwjCeJQb4TX9_Q,64296
 arviz/tests/base_tests/test_data_zarr.py,sha256=sPWnIQ7vPhN5Ql3Dq4JENuSwQV5IeignQjy9BAYe1_A,5441
-arviz/tests/base_tests/test_diagnostics.py,sha256=w8yT2WxelnxH-ynAN9lvspQTRR_UlDmnXIRMe7fus1c,20219
+arviz/tests/base_tests/test_diagnostics.py,sha256=pbuy1-nvTKWSHv0nnhXOhpG4e2uy-4GGZb4lxAdoMpw,20353
 arviz/tests/base_tests/test_diagnostics_numba.py,sha256=2G5O-7Hz66DSaHIZtjs2XL45RezYnXQZH6Dg2Ow-p4Q,2847
 arviz/tests/base_tests/test_helpers.py,sha256=PogHpWCMBEtkuzKt9jGQ8uIPg0cLDwztXxOnPSPNyVE,669
 arviz/tests/base_tests/test_labels.py,sha256=X08vTMmcgXkYGbE2Qnu_UUDSTAIvSNKdnyqoWwmj008,1686
 arviz/tests/base_tests/test_plot_utils.py,sha256=lwDZYDNrlEOKP-asJv6qu3sH_4y-FiHcFlqnMTpZyhw,11771
 arviz/tests/base_tests/test_plots_bokeh.py,sha256=1JqUqLKUb1g4c4w41K3j_LCT4eqb3u1qGnYbUuJMHPE,39148
-arviz/tests/base_tests/test_plots_matplotlib.py,sha256=l9usdmjdI8RV954rknfJvcC-wwnoaRtqJ8wN2SukH1I,65818
+arviz/tests/base_tests/test_plots_matplotlib.py,sha256=sRn-jWzBeCkk49i8DO_a1l9WXoia269sANOQfPpqYjQ,66197
 arviz/tests/base_tests/test_rcparams.py,sha256=b9ueOXd9C0xiUIqgS0qnzvalHFgTFK7sUqL8UAzgJNs,10851
-arviz/tests/base_tests/test_stats.py,sha256=uBrlHdwLQ9w9bP4Cq454qoZ754Xll1kYOpxEr2AwZF4,32385
+arviz/tests/base_tests/test_stats.py,sha256=zEK8uKvXBFpnLQPxoHkM5L9-KWqHQXvUcWAzKGMGgFQ,32809
 arviz/tests/base_tests/test_stats_ecdf_utils.py,sha256=p1FnQzlC0fjjKDFfhbHIrrbwAVhLiygH4J0aarx89A0,6038
 arviz/tests/base_tests/test_stats_numba.py,sha256=wGXgNuSO_gwJajoYtXSgpIe88PcBRyIkRihxC8paR-o,1582
 arviz/tests/base_tests/test_stats_utils.py,sha256=Udkw8tODs8mLt3_hO3HgNczrU0n09IJrML2agXF-upQ,13864
@@ -176,8 +176,8 @@ arviz/wrappers/__init__.py,sha256=d8GTUuBW_30LyDyk6qn2MAnvg-GZCeUw_i5SUPqaa1w,35
 arviz/wrappers/base.py,sha256=Vvh330pdzIvBEaikHsDP1ej6L2jCZZ0Dqj5TvUbYesI,9134
 arviz/wrappers/wrap_pymc.py,sha256=ltKv55aG0WTWXVDJuff5TXkgJJ_ESLvlT-JPlh3yHAg,1143
 arviz/wrappers/wrap_stan.py,sha256=c40brlajoPcc3xk00xI9Hqc-y0xcbAmFAIZOtfXWeqo,5525
-arviz-0.20.0.dist-info/LICENSE,sha256=xllut76FgcGL5zbIRvuRc7aezPbvlMUTWJPsVr2Sugg,11358
-arviz-0.20.0.dist-info/METADATA,sha256=xxMeqVgWrA6YkcGKDSJhbCyn-FKRuIBqlUIl8t1Qess,8848
-arviz-0.20.0.dist-info/WHEEL,sha256=eOLhNAGa2EW3wWl_TU484h7q1UNgy0JXjjoqKoxAAQc,92
-arviz-0.20.0.dist-info/top_level.txt,sha256=5MFvqrTtYRWsIx-SjKuFIUHtrnVJq0Ngd0Nc2_etQhE,6
-arviz-0.20.0.dist-info/RECORD,,
+arviz-0.21.0.dist-info/LICENSE,sha256=xllut76FgcGL5zbIRvuRc7aezPbvlMUTWJPsVr2Sugg,11358
+arviz-0.21.0.dist-info/METADATA,sha256=3zSEB4npy44RXZtTxQVEkIvjvCqbDXQ9s6vsAyxvXmg,8848
+arviz-0.21.0.dist-info/WHEEL,sha256=tZoeGjtWxWRfdplE7E3d45VPlLNQnvbKiYnx7gwAy8A,92
+arviz-0.21.0.dist-info/top_level.txt,sha256=5MFvqrTtYRWsIx-SjKuFIUHtrnVJq0Ngd0Nc2_etQhE,6
+arviz-0.21.0.dist-info/RECORD,,

{arviz-0.20.0.dist-info → arviz-0.21.0.dist-info}/WHEEL

@@ -1,5 +1,5 @@
 Wheel-Version: 1.0
-Generator: bdist_wheel (0.44.0)
+Generator: bdist_wheel (0.45.1)
 Root-Is-Purelib: true
 Tag: py3-none-any