plothist 1.5.0__tar.gz → 1.7.0__tar.gz

{plothist-1.5.0 → plothist-1.7.0}/.github/workflows/cd.yaml

@@ -15,7 +15,7 @@ jobs:
     runs-on: ubuntu-latest
 
     steps:
-      - uses: actions/checkout@v4
+      - uses: actions/checkout@v5
         with:
           persist-credentials: false
       - name: Install uv and set the python version
@@ -46,7 +46,7 @@ jobs:
 
     steps:
       - name: Download all the dists
-        uses: actions/download-artifact@v4
+        uses: actions/download-artifact@v5
         with:
           name: python-package-distributions
           path: dist/
@@ -71,7 +71,7 @@ jobs:
 
     steps:
       - name: Download all the dists
-        uses: actions/download-artifact@v4
+        uses: actions/download-artifact@v5
         with:
           name: python-package-distributions
           path: dist/

{plothist-1.5.0 → plothist-1.7.0}/.github/workflows/ci.yaml

@@ -16,7 +16,7 @@ jobs:
         python-version: ["3.9", "3.10", "3.11", "3.12", "3.13"]
 
     steps:
-      - uses: actions/checkout@v4
+      - uses: actions/checkout@v5
 
       - name: Install uv and set the python version
         uses: astral-sh/setup-uv@v6

{plothist-1.5.0 → plothist-1.7.0}/.github/workflows/coverage.yaml

@@ -11,7 +11,7 @@ jobs:
     runs-on: ubuntu-latest
     steps:
       - name: Checkout
-        uses: actions/checkout@v4
+        uses: actions/checkout@v5
         with:
           fetch-depth: 2
 

{plothist-1.5.0 → plothist-1.7.0}/.github/workflows/pr.yaml

@@ -16,6 +16,6 @@ jobs:
     name: Validate PR title
     runs-on: ubuntu-latest
     steps:
-      - uses: amannn/action-semantic-pull-request@v5
+      - uses: amannn/action-semantic-pull-request@v6
         env:
           GITHUB_TOKEN: ${{ secrets.GITHUB_TOKEN }}

{plothist-1.5.0 → plothist-1.7.0}/.gitignore

@@ -15,6 +15,7 @@ temp_img/
 html/
 sg_execution_times.rst
 docs/examples/*/*.svg
+src/plothist/examples/*/*.svg
 example_gallery/
 .svg_metadata.yaml
 tests/*.yaml

{plothist-1.5.0 → plothist-1.7.0}/.pre-commit-config.yaml

@@ -5,7 +5,7 @@ ci:
 
 repos:
   - repo: https://github.com/adamchainz/blacken-docs
-    rev: "1.19.1"
+    rev: "1.20.0"
     hooks:
       - id: blacken-docs
         args: ["--line-length", "120"]
@@ -13,7 +13,7 @@ repos:
           - black==24.2.0
 
   - repo: https://github.com/pre-commit/pre-commit-hooks
-    rev: v5.0.0
+    rev: v6.0.0
     hooks:
       - id: check-added-large-files
       - id: check-case-conflict
@@ -30,7 +30,7 @@ repos:
         args: ["--pytest-test-first"]
 
   - repo: https://github.com/astral-sh/ruff-pre-commit
-    rev: v0.12.2
+    rev: v0.13.3
     hooks:
       - id: ruff
         args: ["--fix", "--show-fixes"]
@@ -44,13 +44,13 @@ repos:
           - tomli
 
   - repo: https://github.com/python-jsonschema/check-jsonschema
-    rev: 0.33.2
+    rev: 0.34.0
     hooks:
       - id: check-github-workflows
         args: ["--verbose"]
 
   - repo: https://github.com/pre-commit/mirrors-mypy
-    rev: v1.16.1
+    rev: v1.18.2
     hooks:
       - id: mypy
         files: src
@@ -66,13 +66,13 @@ repos:
       - id: rst-inline-touching-normal
 
   - repo: https://github.com/python-jsonschema/check-jsonschema
-    rev: 0.33.2
+    rev: 0.34.0
     hooks:
       - id: check-readthedocs
       - id: check-dependabot
       - id: check-github-workflows
 
   - repo: https://github.com/henryiii/validate-pyproject-schema-store
-    rev: 2025.06.23
+    rev: 2025.10.03
     hooks:
       - id: validate-pyproject

plothist-1.7.0/AUTHORS.md

@@ -0,0 +1,2 @@
+* Cyrille Praz ([@cyrraz](https://github.com/cyrraz))
+* Tristan Fillinger ([@0ctagon](https://github.com/0ctagon))

plothist-1.7.0/CONTRIBUTING.md

@@ -0,0 +1 @@
+For contributing, refer to the chapter [Contributing](https://plothist.readthedocs.io/en/latest/CONTRIBUTING.html) in the documentation.

{plothist-1.5.0 → plothist-1.7.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: plothist
-Version: 1.5.0
+Version: 1.7.0
 Summary: Plot histograms in a scalable way and a beautiful style.
 Project-URL: Homepage, https://github.com/cyrraz/plothist
 Project-URL: Documentation, https://plothist.readthedocs.io/

{plothist-1.5.0 → plothist-1.7.0}/codecov.yml

@@ -1,6 +1,3 @@
 comment: false
 
 github_checks: false
-
-ignore:
-  - "src/plothist/scripts"

{plothist-1.5.0 → plothist-1.7.0}/docs/basics/variable_registry.rst

@@ -25,7 +25,7 @@ For each variable in the list, the following information is stored by default in
 
     variable_0:
         name: variable_0
-        bins: 50
+        bins: auto
         range:
         - min
         - max
@@ -63,24 +63,24 @@ Update the registry
 
 Multiple functions are available to modify the plotting information of the variables in the registry, add or remove some parameters.
 
-Ranges
-------
+Binning and ranges
+------------------
 
-The :func:`update_variable_registry_ranges() <plothist.variable_registry.update_variable_registry_ranges>` function automatically updates the range parameter in the ``yaml`` file to the ``min`` and ``max`` values of the variable in the dataset:
+The :func:`update_variable_registry_binning() <plothist.variable_registry.update_variable_registry_binning>` function automatically updates the number of bins parameter in the ``yaml`` file to the length of [``numpy.histogram_bin_edges``](https://numpy.org/doc/2.1/reference/generated/numpy.histogram_bin_edges.html#numpy-histogram-bin-edges) minus one (the bins are regular) and automatically updates the range parameter in the ``yaml`` file to the ``min`` and ``max`` values of the variable in the dataset:
 
 .. code-block:: python
 
-    from plothist import update_variable_registry_ranges
+    from plothist import update_variable_registry_binning
 
-    update_variable_registry_ranges(df, variable_keys)
+    update_variable_registry_binning(df, variable_keys)
 
-The range has been updated for all the variables in ``variables_keys``. The ``yaml`` file is now:
+The number of bins and the range has been updated for all the variables in ``variables_keys``. The ``yaml`` file is now:
 
 .. code-block:: yaml
 
     variable_0:
         name: variable_0
-        bins: 50
+        bins: 121 # = len(numpy.histogram_bin_edges(df["variable_0"], bins="auto")) - 1
         range:
         - -10.55227774892869    # min(df["variable_0"])
         - 10.04658448558009     # max(df["variable_0"])
@@ -94,10 +94,9 @@ The range has been updated for all the variables in ``variables_keys``. The ``ya
     variable_1:
         ...
 
-Then, you may manually modify the ``yaml`` to get a more suitable range to display in the plot.
-
-Calling this function again on the same variable keys will not overwrite their ``range`` parameter, unless the ``overwrite`` parameter is set to ``True``.
+Then, you may manually modify the ``yaml`` to get a more suitable binning and range to display in the plot.
 
+Calling this function again on the same variable keys will not overwrite their ``bins`` or ``range`` parameter, unless the ``overwrite`` parameter is set to ``True``.
 
 Add or modify variable properties
 ---------------------------------
@@ -124,7 +123,7 @@ This will add the new properties to the ``yaml`` file to all the variables in ``
 
     variable_0:
         name: variable_0
-        bins: 50
+        bins: 121
         range:
         - -10.55227774892869
         - 10.04658448558009
@@ -160,7 +159,7 @@ To remove a parameter from the plotting information, you can use the :func:`remo
 
     from plothist import remove_variable_registry_parameters
 
-    remove_variable_registry_parameters(["range", "log", "legend_ncols", "new_property"], variable_keys)
+    remove_variable_registry_parameters(["bins", "range", "log", "legend_ncols", "new_property"], variable_keys)
 
 The ``yaml`` file is updated:
 
@@ -168,7 +167,6 @@ The ``yaml`` file is updated:
 
     variable_0:
         name: variable_0
-        bins: 50
         label: variable_0
         legend_location: best
         docstring: ''
@@ -197,7 +195,7 @@ Here is an example of how to create, update, and use the variable registry to pl
         plot_hist,
         create_variable_registry,
         update_variable_registry,
-        update_variable_registry_ranges,
+        update_variable_registry_binning,
         get_variable_from_registry,
         add_text,
     )
@@ -208,8 +206,8 @@ Here is an example of how to create, update, and use the variable registry to pl
     # Create the registry
     create_variable_registry(variable_keys)
 
-    # Update the ranges
-    update_variable_registry_ranges(df, variable_keys)
+    # Update the number of bins and range
+    update_variable_registry_binning(df, variable_keys)
 
     # Add custom info
     update_variable_registry({"text": "my analysis"}, variable_keys)
@@ -245,7 +243,7 @@ Example: to plot a zoom on a variable but still keep the original one, you can c
 
     variable_0:
         name: variable_0
-        bins: 50
+        bins: 121
         range:
         - -10
         - 10
@@ -257,7 +255,7 @@ Example: to plot a zoom on a variable but still keep the original one, you can c
 
     variable_0_zoom:
         name: variable_0
-        bins: 50
+        bins: 121
         range:
         - -1
         - 1

{plothist-1.5.0 → plothist-1.7.0}/docs/conf.py

@@ -1,5 +1,8 @@
 import sys
 
+from sphinx_gallery.scrapers import matplotlib_scraper
+from sphinx_gallery.sorting import ExplicitOrder
+
 import plothist
 
 #
@@ -175,9 +178,6 @@ def reset_mpl(gallery_conf, fname):
     set_style("default")
 
 
-from sphinx_gallery.scrapers import matplotlib_scraper
-
-
 class matplotlib_svg_scraper:
     def __repr__(self):
         return self.__class__.__name__
@@ -186,8 +186,6 @@ class matplotlib_svg_scraper:
         return matplotlib_scraper(*args, format="svg", bbox_inches="tight", **kwargs)
 
 
-from sphinx_gallery.sorting import ExplicitOrder
-
 sphinx_gallery_conf = {
     # path to your example scripts
     "examples_dirs": ["../src/plothist/examples"],

plothist-1.7.0/docs/documentation/statistics.rst

@@ -0,0 +1,21 @@
+.. _documentation-statistics-label:
+
+===================
+Notes on statistics
+===================
+
+In the :func:`plot_error_hist() <plothist.plotters.plot_error_hist>` function, the represented uncertainties can be:
+
+- ``symmetrical``: use the square root of the variance of the histogram
+- ``asymmetrical``: frequentist approach based on a Poisson confidence interval (see `page 16 on this presentation <https://www.nikhef.nl/~ivov/Statistics/PoissonError/2017_05_15_PoissonError_LHCb_IvovanVulpen.pdf>`_)
+
+    - the upper limit on bins with no entries can be interpreted differently. By default, a one-sided 68% confidence bound is used. The behavior can be changed using the ``uncertainty_type = asymmetrical_double_sided_zeros``, which will use a two-sided 68% confidence bound.
+    - this method doesn't support weighted histograms. For more details, see the `note about uncertainties when the histogram entries are weighted <https://www.pp.rhul.ac.uk/~cowan/stat/notes/errors_with_weights.pdf>`_.
+
+.. image:: ../img/uncertainty_types.png
+   :alt: Simple error hist
+   :width: 500
+
+The uncertainty on the efficiency (ratio of two histograms h1 and h2 when the entries of h1 are a subset of the entries of h2) in :func:`get_efficiency() <plothist.comparison.get_efficiency>` is calculated following equation (19) of:
+
+- `T. Ullrich, Z. Xu, Treatment of Errors in Efficiency Calculations, arxiv: 0701199v1. <https://arxiv.org/pdf/physics/0701199v1.pdf>`_

{plothist-1.5.0 → plothist-1.7.0}/docs/usage/installation.rst

@@ -64,11 +64,11 @@ To update ``plothist`` to its latest stable version, follow the instructions abo
 Install the development version
 ===============================
 
-1. Clone the plothist repository and go to the package folder:
+1. `Fork <https://docs.github.com/en/pull-requests/collaborating-with-pull-requests/working-with-forks/fork-a-repo?tool=webui>`_ and clone locally the plothist repository, then go to the package folder:
 
     .. code-block:: bash
 
-       git clone git@github.com:cyrraz/plothist.git
+       git clone git@github.com:USERNAME/plothist.git
 
        cd plothist
 

{plothist-1.5.0 → plothist-1.7.0}/pyproject.toml

@@ -83,7 +83,9 @@ extend-select = [
     "PD",          # pandas-vet
 ]
 ignore = [
-    "E402",  # module level import not at top of file
     "NPY002",  # Replace legacy `np.random` call with `np.random.Generator`
-    "PD901",  # Avoid using the generic variable name `df` for DataFrames
+]
+[tool.ruff.lint.per-file-ignores]
+"src/plothist/examples/**/*.py" = [
+    "E402",  # module level import not at top of file
 ]

{plothist-1.5.0 → plothist-1.7.0}/src/plothist/__init__.py

@@ -9,12 +9,7 @@ from .comparison import (
     get_ratio,
     get_ratio_variances,
 )
-from .histogramming import (
-    create_axis,
-    flatten_2d_hist,
-    make_2d_hist,
-    make_hist,
-)
+from .histogramming import create_axis, flatten_2d_hist, make_2d_hist, make_hist
 from .plothist_style import (
     add_luminosity,
     add_text,
@@ -43,6 +38,7 @@ from .variable_registry import (
     get_variable_from_registry,
     remove_variable_registry_parameters,
     update_variable_registry,
+    update_variable_registry_binning,
     update_variable_registry_ranges,
 )
 
@@ -81,22 +77,24 @@ __all__ = [
     "set_fitting_ylabel_fontsize",
     "set_style",
     "update_variable_registry",
+    "update_variable_registry_binning",
     "update_variable_registry_ranges",
 ]
 
 
-# Get style file and use it
+from importlib import resources
 from importlib.resources import files
 
+import boost_histogram as bh
+import matplotlib.font_manager as fm
 import matplotlib.pyplot as plt
 
+# Get style file and use it
+
 style_file = files("plothist").joinpath("default_style.mplstyle")
 plt.style.use(style_file)
 
 # Install fonts
-from importlib import resources
-
-import matplotlib.font_manager as fm
 
 with resources.as_file(resources.files("plothist_utils") / "fonts") as font_path:
     font_files = fm.findSystemFonts(fontpaths=[str(font_path)])
@@ -104,7 +102,6 @@ with resources.as_file(resources.files("plothist_utils") / "fonts") as font_path
         fm.fontManager.addfont(font)
 
 # Check version of boost_histogram
-import boost_histogram as bh
 
 if tuple(int(part) for part in bh.__version__.split(".")) < (1, 4, 0):
     raise ImportError(

{plothist-1.5.0 → plothist-1.7.0}/src/plothist/_version.py

@@ -1,7 +1,14 @@
 # file generated by setuptools-scm
 # don't change, don't track in version control
 
-__all__ = ["__version__", "__version_tuple__", "version", "version_tuple"]
+__all__ = [
+    "__version__",
+    "__version_tuple__",
+    "version",
+    "version_tuple",
+    "__commit_id__",
+    "commit_id",
+]
 
 TYPE_CHECKING = False
 if TYPE_CHECKING:
@@ -9,13 +16,19 @@ if TYPE_CHECKING:
     from typing import Union
 
     VERSION_TUPLE = Tuple[Union[int, str], ...]
+    COMMIT_ID = Union[str, None]
 else:
     VERSION_TUPLE = object
+    COMMIT_ID = object
 
 version: str
 __version__: str
 __version_tuple__: VERSION_TUPLE
 version_tuple: VERSION_TUPLE
+commit_id: COMMIT_ID
+__commit_id__: COMMIT_ID
 
-__version__ = version = '1.5.0'
-__version_tuple__ = version_tuple = (1, 5, 0)
+__version__ = version = '1.7.0'
+__version_tuple__ = version_tuple = (1, 7, 0)
+
+__commit_id__ = commit_id = None

{plothist-1.5.0 → plothist-1.7.0}/src/plothist/comparison.py

@@ -22,9 +22,16 @@ def _check_uncertainty_type(uncertainty_type: str) -> None:
         If the uncertainty type is not valid.
 
     """
-    if uncertainty_type not in ["symmetrical", "asymmetrical"]:
+    _valid_uncertainty_types = [
+        "symmetrical",
+        "asymmetrical",
+        "asymmetrical_double_sided_zeros",
+        "asymmetrical_one_sided_zeros",
+    ]
+
+    if uncertainty_type not in _valid_uncertainty_types:
         raise ValueError(
-            f"Uncertainty type {uncertainty_type} not valid. Must be 'symmetrical' or 'asymmetrical'."
+            f"Uncertainty type {uncertainty_type} not valid. Must be in {_valid_uncertainty_types}."
         )
 
 
@@ -42,11 +49,12 @@ def _is_unweighted(hist: bh.Histogram) -> bool:
     bool
         True if the histogram is unweighted, False otherwise.
     """
-    return np.allclose(hist.variances(), hist.values())
+    return np.allclose(hist.variances(), hist.values(), equal_nan=True)
 
 
 def get_asymmetrical_uncertainties(
     hist: bh.Histogram,
+    uncertainty_type: str = "asymmetrical",
 ) -> tuple[np.ndarray, np.ndarray]:
     """
     Get Poisson asymmetrical uncertainties for a histogram via a frequentist approach based on a confidence-interval computation.
@@ -57,6 +65,8 @@ def get_asymmetrical_uncertainties(
     ----------
     hist : bh.Histogram
         The histogram.
+    uncertainty_type : str, optional
+        The type of uncertainty to compute for bins with 0 entry. Default is "asymmetrical" (= "asymmetrical_one_sided_zeros"). Use "asymmetrical_double_sided_zeros" to have the double-sided definition. More information in :ref:`documentation-statistics-label`.
 
     Returns
     -------
@@ -72,16 +82,44 @@ def get_asymmetrical_uncertainties(
 
     """
     _check_counting_histogram(hist)
+    _check_uncertainty_type(uncertainty_type)
 
     if not _is_unweighted(hist):
         raise ValueError(
             "Asymmetrical uncertainties can only be computed for an unweighted histogram."
         )
-    conf_level = 0.682689492
-    alpha = 1.0 - conf_level
+
+    alpha = 1.0 - 0.682689492
+    tail_probability = alpha / 2
+
     n = hist.values()
-    uncertainties_low = n - stats.gamma.ppf(alpha / 2, n, scale=1)
-    uncertainties_high = stats.gamma.ppf(1 - alpha / 2, n + 1, scale=1) - n
+
+    lower_bound = np.zeros_like(n, dtype=float)
+    upper_bound = np.zeros_like(n, dtype=float)
+
+    # Two-sided Garwood intervals for n > 0
+    lower_bound[n > 0] = stats.gamma.ppf(q=tail_probability, a=n[n > 0], scale=1)
+    upper_bound[n > 0] = stats.gamma.ppf(
+        q=1 - tail_probability, a=n[n > 0] + 1, scale=1
+    )
+
+    if uncertainty_type == "asymmetrical_double_sided_zeros":
+        # Two-sided Garwood intervals for n == 0
+        upper_bound[n == 0] = stats.gamma.ppf(q=1 - tail_probability, a=1, scale=1)
+    elif uncertainty_type in ["asymmetrical_one_sided_zeros", "asymmetrical"]:
+        # One-sided upper limit for n == 0
+        upper_bound[n == 0] = stats.gamma.ppf(q=1 - 2 * tail_probability, a=1, scale=1)
+    else:
+        raise ValueError(
+            f"Invalid uncertainty type '{uncertainty_type}' for asymmetrical uncertainties."
+        )
+
+    # Compute asymmetric uncertainties
+    uncertainties_low = n - lower_bound
+    uncertainties_high = upper_bound - n
+
+    uncertainties_low = np.nan_to_num(uncertainties_low, nan=0.0)
+    uncertainties_high = np.nan_to_num(uncertainties_high, nan=0.0)
 
     return uncertainties_low, uncertainties_high
 
@@ -174,8 +212,10 @@ def get_pull(
     _check_counting_histogram(h1)
     _check_counting_histogram(h2)
 
-    if h1_uncertainty_type == "asymmetrical":
-        uncertainties_low, uncertainties_high = get_asymmetrical_uncertainties(h1)
+    if "asymmetrical" in h1_uncertainty_type:
+        uncertainties_low, uncertainties_high = get_asymmetrical_uncertainties(
+            h1, h1_uncertainty_type
+        )
         h1_variances = np.where(
             h1.values() >= h2.values(),
             uncertainties_low**2,
@@ -232,8 +272,10 @@ def get_difference(
 
     difference_values = h1.values() - h2.values()
 
-    if h1_uncertainty_type == "asymmetrical":
-        uncertainties_low, uncertainties_high = get_asymmetrical_uncertainties(h1)
+    if "asymmetrical" in h1_uncertainty_type:
+        uncertainties_low, uncertainties_high = get_asymmetrical_uncertainties(
+            h1, h1_uncertainty_type
+        )
 
         difference_uncertainties_low = np.sqrt(uncertainties_low**2 + h2.variances())
         difference_uncertainties_high = np.sqrt(uncertainties_high**2 + h2.variances())
@@ -390,11 +432,13 @@ def get_ratio(
 
     ratio_values = np.where(h2.values() != 0, h1.values() / h2.values(), np.nan)
 
-    if h1_uncertainty_type == "asymmetrical":
-        uncertainties_low, uncertainties_high = get_asymmetrical_uncertainties(h1)
+    if "asymmetrical" in h1_uncertainty_type:
+        uncertainties_low, uncertainties_high = get_asymmetrical_uncertainties(
+            h1, h1_uncertainty_type
+        )
 
     if ratio_uncertainty_type == "uncorrelated":
-        if h1_uncertainty_type == "asymmetrical":
+        if "asymmetrical" in h1_uncertainty_type:
             h1_high = h1.copy()
             h1_high[:] = np.c_[h1_high.values(), uncertainties_high**2]
             h1_low = h1.copy()
@@ -405,7 +449,7 @@ def get_ratio(
             ratio_uncertainties_low = np.sqrt(get_ratio_variances(h1, h2))
             ratio_uncertainties_high = ratio_uncertainties_low
     elif ratio_uncertainty_type == "split":
-        if h1_uncertainty_type == "asymmetrical":
+        if "asymmetrical" in h1_uncertainty_type:
             ratio_uncertainties_low = uncertainties_low / h2.values()
             ratio_uncertainties_high = uncertainties_high / h2.values()
         else:
@@ -493,7 +537,7 @@ def get_comparison(
                 h1, h2, h1_uncertainty_type
             )
         elif comparison == "asymmetry":
-            if h1_uncertainty_type == "asymmetrical":
+            if "asymmetrical" in h1_uncertainty_type:
                 raise ValueError(
                     "Asymmetrical uncertainties are not supported for the asymmetry comparison."
                 )
@@ -501,7 +545,7 @@ def get_comparison(
             lower_uncertainties = uncertainties
             upper_uncertainties = uncertainties
         elif comparison == "efficiency":
-            if h1_uncertainty_type == "asymmetrical":
+            if "asymmetrical" in h1_uncertainty_type:
                 raise ValueError(
                     "Asymmetrical uncertainties are not supported in an efficiency computation."
                 )

{plothist-1.5.0 → plothist-1.7.0}/src/plothist/examples/2d_hist/2d_hist_correlations.py

@@ -20,7 +20,7 @@ from plothist import (
     get_variable_from_registry,
     make_2d_hist,
     plot_2d_hist,
-    update_variable_registry_ranges,
+    update_variable_registry_binning,
 )
 
 # No need to redo this step if the registry was already created before
@@ -28,7 +28,7 @@ variable_keys = ["variable_0", "variable_1", "variable_2"]
 unique_id = str(int(time.time() * 1000))[-8:]  # unique ID based on current time
 temporary_registry_path = f"./_temporary_variable_registry_{unique_id}.yaml"
 create_variable_registry(variable_keys, path=temporary_registry_path)
-update_variable_registry_ranges(df, variable_keys, path=temporary_registry_path)
+update_variable_registry_binning(df, variable_keys, path=temporary_registry_path)
 
 # Get all the correlation plot between the variables
 variable_keys_combinations = list(combinations(variable_keys, 2))