mergeron 2024.738936.0__py3-none-any.whl → 2024.738936.2__py3-none-any.whl

mergeron/core/guidelines_standards.py

@@ -17,6 +17,7 @@ from typing import Any, Literal, TypeAlias
 import numpy as np
 from mpmath import mp, mpf  # type: ignore
 from numpy.typing import NDArray
+from scipy.spatial.distance import minkowski as distance_function
 
 mp.prec = 80
 mp.trap_complex = True
@@ -168,7 +169,7 @@ def lerp(
     /,
 ) -> float | mpf | NDArray[np.float64]:
     """
-    From the C++ standard function of the same name.
+    From the function of the same name in the C++ standard [2]_
 
     Constructs the weighted average, :math:`w_1 x_1 + w_2 x_2`, where
     :math:`w_1 = 1 - r` and :math:`w_2 = r`.
@@ -178,21 +179,35 @@ def lerp(
     _x1, _x2
         bounds :math:`x_1, x_2` to interpolate between.
     _r
-        interpolation weight :math:`r` assigned to :math:`x_2`.
+        interpolation weight :math:`r` assigned to :math:`x_2`
 
     Returns
     -------
-        The linear interpolation, or weighted average, :math:`x_1 (1 - r) + x_2 r`.
+        The linear interpolation, or weighted average,
+        :math:`x_1 + r \\cdot (x_1 - x_2) \\equiv (1 - r) \\cdot x_1 + r \\cdot x_2`.
 
     Raises
     ------
     ValueError
         If the interpolation weight is not in the interval, :math:`[0, 1]`.
 
+    References
+    ----------
+
+    .. [2] C++ Reference, https://en.cppreference.com/w/cpp/numeric/lerp
+
     """
+
     if not 0 <= _r <= 1:
         raise ValueError("Specified interpolation weight must lie in [0, 1].")
-    return _x2 * _r + _x1 * (1 - _r)
+    elif _r == 0:
+        return _x1
+    elif _r == 1:
+        return _x2
+    elif _r == 0.5:
+        return 1 / 2 * (_x1 + _x2)
+    else:
+        return _r * _x2 + (1 - _r) * _x1
 
 
 def gbd_from_dsf(
@@ -552,7 +567,7 @@ def delta_hhi_boundary(
 
 
 def combined_share_boundary(
-    _s_incpt: float = 0.0625, /, *, bdry_dps: int = 10
+    _s_intcpt: float = 0.0625, /, *, bdry_dps: int = 10
 ) -> tuple[NDArray[np.float64], float]:
     """
     Share combinations on the merging-firms' combined share boundary.
@@ -563,7 +578,7 @@ def combined_share_boundary(
 
     Parameters
     ----------
-    _s_incpt:
+    _s_intcpt:
         Merging-firms' combined share.
     bdry_dps
         Number of decimal places for rounding reported shares.
@@ -573,16 +588,16 @@ def combined_share_boundary(
         Array of share-pairs, area under boundary.
 
     """
-    _s_incpt = mpf(f"{_s_incpt}")
-    _s_mid = _s_incpt / 2
+    _s_intcpt = mpf(f"{_s_intcpt}")
+    _s_mid = _s_intcpt / 2
 
-    _s1_pts = (0, _s_mid, _s_incpt)
+    _s1_pts = (0, _s_mid, _s_intcpt)
     return (
         np.column_stack((
             np.array(_s1_pts, np.float64),
             np.array(_s1_pts[::-1], np.float64),
         )),
-        round(float(_s_incpt * _s_mid), bdry_dps),
+        round(float(_s_intcpt * _s_mid), bdry_dps),
     )
 
 
@@ -650,17 +665,17 @@ def shrratio_mgnsym_boundary_max(
     # of function call with other shrratio_mgnsym_boundary functions
     del _r_val
     _delta_star = mpf(f"{_delta_star}")
-    _s_incpt = _delta_star
+    _s_intcpt = _delta_star
     _s_mid = _delta_star / (1 + _delta_star)
 
-    _s1_pts = (0, _s_mid, _s_incpt)
+    _s1_pts = (0, _s_mid, _s_intcpt)
 
     return (
         np.column_stack((
             np.array(_s1_pts, np.float64),
             np.array(_s1_pts[::-1], np.float64),
         )),
-        round(float(_s_incpt * _s_mid), gbd_dps),  # simplified calculation
+        round(float(_s_intcpt * _s_mid), gbd_dps),  # simplified calculation
     )
 
 
@@ -710,7 +725,7 @@ def shrratio_mgnsym_boundary_min(
         )
 
     _delta_star = mpf(f"{_delta_star}")
-    _s_incpt = mpf("1.00")
+    _s_intcpt = mpf("1.00")
     _s_mid = _delta_star / (1 + _delta_star)
 
     if recapture_spec == "inside-out":
@@ -725,14 +740,14 @@ def shrratio_mgnsym_boundary_min(
                 _smin_nr / _guppi_bdry_env_dr,
                 _s_mid,
                 _smax_nr / _guppi_bdry_env_dr,
-                _s_incpt,
+                _s_intcpt,
             ),
             np.float64,
         )
 
         _gbd_area = (_smin_nr + (_smax_nr - _smin_nr) * _s_mid) / _guppi_bdry_env_dr
     else:
-        _s1_pts, _gbd_area = np.array((0, _s_mid, _s_incpt), np.float64), _s_mid
+        _s1_pts, _gbd_area = np.array((0, _s_mid, _s_intcpt), np.float64), _s_mid
 
     return np.column_stack((_s1_pts, _s1_pts[::-1])), round(float(_gbd_area), gbd_dps)
 
@@ -842,10 +857,12 @@ def shrratio_mgnsym_boundary_wtd_avg(
     _delta_star = mpf(f"{_delta_star}")
     _s_mid = _delta_star / (1 + _delta_star)
 
+    # initial conditions
     _gbdry_points = [(_s_mid, _s_mid)]
     _s_1_pre, _s_2_pre = _s_mid, _s_mid
     _s_2_oddval, _s_2_oddsum, _s_2_evnsum = True, 0, 0
 
+    # parameters for iteration
     _gbd_step_sz = mp.power(10, -gbd_dps)
     _theta = _gbd_step_sz * (10 if wgtng_policy == "cross-product-share" else 1)
     for _s_1 in mp.arange(_s_mid - _gbd_step_sz, 0, -_gbd_step_sz):
@@ -883,12 +900,12 @@ def shrratio_mgnsym_boundary_wtd_avg(
                     _delta_test = lerp(_de_1, _de_2, _r)
 
             if wgtng_policy == "cross-product-share":
-                test_flag, incr_decr = (_delta_test > _delta_star, -1)
+                _test_flag, _incr_decr = (_delta_test > _delta_star, -1)
             else:
-                test_flag, incr_decr = (_delta_test < _delta_star, 1)
+                _test_flag, _incr_decr = (_delta_test < _delta_star, 1)
 
-            if test_flag:
-                _s_2 += incr_decr * _gbd_step_sz
+            if _test_flag:
+                _s_2 += _incr_decr * _gbd_step_sz
             else:
                 break
 
@@ -916,13 +933,19 @@ def shrratio_mgnsym_boundary_wtd_avg(
     # Area under boundary
     _gbdry_area_total = 2 * _gbd_prtlarea - mp.power(_s_mid, 2)
 
-    _gbdry_points = np.row_stack((
-        _gbdry_points,
-        (
-            mpf("0.0"),
-            _delta_star if wgtng_policy == "cross-product-share" else mpf("1.0"),
-        ),
-    )).astype(np.float64)
+    match wgtng_policy:
+        case "cross-product-share":
+            _s_intcpt = _delta_star
+        case "own-product-share":
+            _s_intcpt = mpf("1.0")
+        case None if avg_method == "distance":
+            _s_intcpt = _delta_star * mp.sqrt("2")
+        case _:
+            _s_intcpt = _s_2_pre
+
+    _gbdry_points = np.row_stack((_gbdry_points, (mpf("0.0"), _s_intcpt))).astype(
+        np.float64
+    )
     # Points defining boundary to point-of-symmetry
     return (
         np.row_stack((np.flip(_gbdry_points, 0), np.flip(_gbdry_points[1:], 1))),
@@ -1009,13 +1032,13 @@ def shrratio_mgnsym_boundary_xact_avg(
     _gbdry_points_start = np.array([(_s_mid, _s_mid)])
     _s_1 = np.array(mp.arange(_s_mid - _gbd_step_sz, 0, -_gbd_step_sz), np.float64)
     if recapture_spec == "inside-out":
-        _s_incpt = mp.fdiv(
+        _s_intcpt = mp.fdiv(
             mp.fsub(
                 2 * _delta_star * _r_val + 1, mp.fabs(2 * _delta_star * _r_val - 1)
             ),
             2 * mpf(f"{_r_val}"),
         )
-        _nr_t1 = 1 + 2 * _delta_star * _r_val * (1 - _s_1) - _s_1 * (1 - _r_val)
+        _nr_t1 = 1 + 2 * _delta_star * _r_val * (1 - _s_1) - _s_1 * (1 - _r_val)  # type: ignore
 
         _nr_sqrt_mdr = 4 * _delta_star * _r_val
         _nr_sqrt_mdr2 = _nr_sqrt_mdr * _r_val
@@ -1040,9 +1063,9 @@ def shrratio_mgnsym_boundary_xact_avg(
 
         _nr_t2_s1 = _nr_sqrt_s1sq + _nr_sqrt_s1 + _nr_sqrt_nos1
 
-        if not np.isclose(
-            np.einsum("i->", _nr_t2_mdr.astype(np.float64)),  # from mpf to float64
-            np.einsum("i->", _nr_t2_s1.astype(np.float64)),
+        if not np.isclose(  # type: ignore
+            np.einsum("i->", _nr_t2_mdr.astype(np.float64)),  # type: ignore from mpf to float64
+            np.einsum("i->", _nr_t2_s1.astype(np.float64)),  # type: ignore from mpf to float64
             rtol=0,
             atol=0.5 * gbd_dps,
         ):
@@ -1054,7 +1077,7 @@ def shrratio_mgnsym_boundary_xact_avg(
         _s_2 = (_nr_t1 - np.sqrt(_nr_t2_s1)) / (2 * _r_val)
 
     else:
-        _s_incpt = mp.fsub(_delta_star + 1 / 2, mp.fabs(_delta_star - 1 / 2))
+        _s_intcpt = mp.fsub(_delta_star + 1 / 2, mp.fabs(_delta_star - 1 / 2))
         _s_2 = (
             (1 / 2)
             + _delta_star
@@ -1069,7 +1092,7 @@ def shrratio_mgnsym_boundary_xact_avg(
         )
 
     _gbdry_points_inner = np.column_stack((_s_1, _s_2))
-    _gbdry_points_end = np.array([(mpf("0.0"), _s_incpt)], np.float64)
+    _gbdry_points_end = np.array([(mpf("0.0"), _s_intcpt)], np.float64)
 
     _gbdry_points = np.row_stack((
         _gbdry_points_end,
@@ -1111,6 +1134,8 @@ def shrratio_mgnsym_boundary_avg(
     Share combinations along the average GUPPI boundary, with
     symmetric merging-firm margins.
 
+    Faster than calling shrratio_mgnsym_boundary_avg(wgtng_policy=None).
+
     Parameters
     ----------
     _delta_star
@@ -1118,7 +1143,7 @@ def shrratio_mgnsym_boundary_avg(
     _r_val
         Recapture ratio.
     avg_method
-        Whether "arithmetic", "geometric", or "root-mean-square".
+        Whether "arithmetic", "geometric", or "distance".
     recapture_spec
         Whether recapture-ratio is MNL-consistent ("inside-out") or has fixed
         value for both merging firms ("proportional").
@@ -1137,7 +1162,7 @@ def shrratio_mgnsym_boundary_avg(
             "Margin-adjusted benchmark share ratio cannot exceed 1."
         )
 
-    if avg_method not in (_avgmthds := ("arithmetic", "geometric", "root-mean-square")):
+    if avg_method not in (_avgmthds := ("arithmetic", "geometric", "distance")):
         raise ValueError(
             f"Averarging method, {f'"{avg_method}"'} is invalid. "
             f"Must be one of, {_avgmthds!r}."
@@ -1151,18 +1176,21 @@ def shrratio_mgnsym_boundary_avg(
 
     _delta_star = mpf(f"{_delta_star}")
     _s_mid = _delta_star / (1 + _delta_star)
-    _gbd_step_sz = mp.power(10, -gbd_dps)
 
+    # initial conditions
     _s_2 = _s_mid
     _s_2_oddval = True
     _s_2_oddsum = 0
     _s_2_evnsum = 0
     _gbdry_points = [(_s_mid, _s_mid)]
+
+    # parameters for iteration
+    _gbd_step_sz = mp.power(10, -gbd_dps)
     for _s_1 in mp.arange(_s_mid, 0, -_gbd_step_sz):
         _s_1 -= _gbd_step_sz
         while True:
-            _shratio_12 = _s_2 / (1 - _s_1)
-            _shratio_21 = (
+            _delta_12 = _s_2 / (1 - _s_1)
+            _delta_21 = (
                 _s_1 / (1 - _s_2)
                 if recapture_spec == "proportional"
                 else _s_1 / (1 - lerp(_s_1, _s_2, _r_val))
@@ -1170,13 +1198,21 @@ def shrratio_mgnsym_boundary_avg(
 
             match avg_method:
                 case "geometric":
-                    _small_d = mp.sqrt(_shratio_12 * _shratio_21)
+                    _delta_test = mp.sqrt(_delta_12 * _delta_21)
                 case "distance":
-                    _small_d = mp.sqrt(mp.fdiv((_shratio_12**2 + _shratio_21**2), "2"))
+                    # _delta_test = mp.sqrt(mp.fdiv((_delta_12**2 + _delta_21**2), "2"))
+                    _delta_test = mp.sqrt(
+                        mp.fdiv(
+                            mp.fsum(
+                                mp.power(f"{_g}", "2") for _g in (_delta_12, _delta_21)
+                            ),
+                            "2",
+                        )
+                    )
                 case _:
-                    _small_d = mp.fdiv(_shratio_12 + _shratio_21, "2")
+                    _delta_test = mp.fdiv(_delta_12 + _delta_21, "2")
 
-            if _small_d < _delta_star:
+            if _delta_test < _delta_star:
                 _s_2 += _gbd_step_sz
             else:
                 break
@@ -1189,12 +1225,12 @@ def shrratio_mgnsym_boundary_avg(
 
     # Starting at _s_id - _gbd_step_sz means _s_1 is not always
     # an even multiple of _gbd_step_sz
-    _s_incpt = _s_2
+    _s_intcpt = _s_2
 
     _gbd_prtlarea = 2 * _gbd_step_sz * (
         mp.fmul(4 / 3, _s_2_oddsum)
         + mp.fmul(2 / 3, _s_2_evnsum)
-        + mp.fmul(1 / 3, _s_mid + _s_incpt)
+        + mp.fmul(1 / 3, _s_mid + _s_intcpt)
     ) - mp.power(_s_mid, 2)
 
     _gbdry_points = np.array(_gbdry_points, np.float64)
@@ -1202,3 +1238,158 @@ def shrratio_mgnsym_boundary_avg(
         np.row_stack((np.flip(_gbdry_points, 0), np.flip(_gbdry_points[1:], 1))),
         round(float(_gbd_prtlarea), gbd_dps),
     )
+
+
+def shrratio_mgnsym_boundary_distance(
+    _delta_star: float = 0.075,
+    _r_val: float = 0.80,
+    /,
+    *,
+    avg_method: Literal["arithmetic", "distance"] = "arithmetic",
+    wgtng_policy: Literal["own-share", "cross-product-share"] | None = "own-share",
+    recapture_spec: Literal["inside-out", "proportional"] = "inside-out",
+    gbd_dps: int = 5,
+) -> tuple[NDArray[np.float64], float]:
+    """
+    Share combinations for the GUPPI boundaries using various aggregators with symmetric merging-firm margins.
+
+    Parameters
+    ----------
+    _delta_star
+        corollary to GUPPI bound (:math:`\\overline{g} / (m^* \\cdot \\overline{r})`)
+    _r_val
+        recapture ratio
+    avg_method
+        Whether "arithmetic", "geometric", or "distance".
+    wgtng_policy
+        Whether "own-share" or "cross-product-share".
+    recapture_spec
+        Whether recapture-ratio is MNL-consistent ("inside-out") or has fixed
+        value for both merging firms ("proportional").
+    gbd_dps
+        Number of decimal places for rounding returned shares and area.
+
+    Returns
+    -------
+        Array of share-pairs, area under boundary.
+
+    """
+
+    if _delta_star > 1:
+        raise ValueError(
+            "Margin-adjusted benchmark share ratio, `_delta_star` cannot exceed 1."
+        )
+
+    _delta_star = mpf(f"{_delta_star}")
+    _s_mid = _delta_star / (1 + _delta_star)
+
+    # initial conditions
+    _gbdry_points = [(_s_mid, _s_mid)]
+    _s_1_pre, _s_2_pre = _s_mid, _s_mid
+    _s_2_oddval, _s_2_oddsum, _s_2_evnsum = True, 0, 0
+
+    # parameters for iteration
+    _weights_base = (mpf("0.5"),) * 2
+    _gbd_step_sz = mp.power(10, -gbd_dps)
+    _theta = _gbd_step_sz * (10 if wgtng_policy == "cross-product-share" else 1)
+    for _s_1 in mp.arange(_s_mid - _gbd_step_sz, 0, -_gbd_step_sz):
+        # The wtd. avg. GUPPI is not always convex to the origin, so we
+        #   increment _s_2 after each iteration in which our algorithm
+        #   finds (s1, s2) on the boundary
+        _s_2 = _s_2_pre * (1 + _theta)
+
+        if (_s_1 + _s_2) > mpf("0.99875"):
+            # 1: # We lose accuracy at 3-9s and up
+            break
+
+        while True:
+            _de_1 = _s_2 / (1 - _s_1)
+            _de_2 = (
+                _s_1 / (1 - lerp(_s_1, _s_2, _r_val))
+                if recapture_spec == "inside-out"
+                else _s_1 / (1 - _s_2)
+            )
+
+            _weights_i = (
+                (
+                    _w1 := mp.fdiv(
+                        _s_2 if wgtng_policy == "cross-product-share" else _s_1,
+                        _s_1 + _s_2,
+                    ),
+                    1 - _w1,
+                )
+                if wgtng_policy
+                else _weights_base
+            )
+
+            match avg_method:
+                case "arithmetic":
+                    _delta_test = distance_function(
+                        (_de_1, _de_2), (0.0, 0.0), p=1, w=_weights_i
+                    )
+                case "distance":
+                    _delta_test = distance_function(
+                        (_de_1, _de_2), (0.0, 0.0), p=2, w=_weights_i
+                    )
+
+            if wgtng_policy == "cross-product-share":
+                _test_flag, _incr_decr = (_delta_test > _delta_star, -1)
+            else:
+                _test_flag, _incr_decr = (_delta_test < _delta_star, 1)
+
+            if _test_flag:
+                _s_2 += _incr_decr * _gbd_step_sz
+            else:
+                break
+
+        # Build-up boundary points
+        _gbdry_points.append((_s_1, _s_2))
+
+        # Build up area terms
+        _s_2_oddsum += _s_2 if _s_2_oddval else 0
+        _s_2_evnsum += _s_2 if not _s_2_oddval else 0
+        _s_2_oddval = not _s_2_oddval
+
+        # Hold share points
+        _s_2_pre = _s_2
+        _s_1_pre = _s_1
+
+    _gbd_prtlarea = _gbd_step_sz * (
+        (4 * _s_2_oddsum + 2 * _s_2_evnsum + _s_mid + _delta_star) / 3
+        if wgtng_policy == "cross-product-share"
+        else (
+            (4 * _s_2_oddsum + 2 * _s_2_evnsum + _s_mid + _s_2_pre) / 3
+            + _s_1_pre * (1 + _s_2_pre) / 2
+        )
+    )
+
+    # Area under boundary
+    _gbdry_area_total = 2 * _gbd_prtlarea - mp.power(_s_mid, 2)
+
+    match wgtng_policy:
+        case "cross-product-share":
+            _s_intcpt = _delta_star
+        case "own-product-share":
+            _s_intcpt = mpf("1.0")
+        case None if avg_method == "distance":
+            _s_intcpt = _delta_star * mp.sqrt("2")
+        case None if avg_method == "arithmetic" and recapture_spec == "inside-out":
+            _s_intcpt = mp.fdiv(
+                mp.fsub(
+                    2 * _delta_star * _r_val + 1, mp.fabs(2 * _delta_star * _r_val - 1)
+                ),
+                2 * mpf(f"{_r_val}"),
+            )
+        case None if avg_method == "arithmetic":
+            _s_intcpt = mp.fsub(_delta_star + 1 / 2, mp.fabs(_delta_star - 1 / 2))
+        case _:
+            _s_intcpt = _s_2_pre
+
+    _gbdry_points = np.row_stack((_gbdry_points, (mpf("0.0"), _s_intcpt))).astype(
+        np.float64
+    )
+    # Points defining boundary to point-of-symmetry
+    return (
+        np.row_stack((np.flip(_gbdry_points, 0), np.flip(_gbdry_points[1:], 1))),
+        round(float(_gbdry_area_total), gbd_dps),
+    )

mergeron/core/proportions_tests.py

@@ -41,7 +41,7 @@ def propn_ci(
 ]:
     """Returns point estimates and confidence interval for a proportion
 
-    Methods "Clopper-Pearson" and "Exact" are synoymous [2]_.  Similarly,
+    Methods "Clopper-Pearson" and "Exact" are synoymous [3]_.  Similarly,
     "Wilson" and "Score" are synonyms here.
 
     Parameters
@@ -66,7 +66,7 @@ def propn_ci(
     References
     ----------
 
-    .. [2] Alan Agresti & Brent A. Coull (1998) Approximate is Better
+    .. [3] Alan Agresti & Brent A. Coull (1998) Approximate is Better
        than “Exact” for Interval Estimation of Binomial Proportions,
        The American Statistician, 52:2, 119-126,
        https://doi.org/10.1080/00031305.1998.10480550
@@ -200,9 +200,9 @@ def propn_diff_ci(
 ) -> tuple[float, float, float, float]:
     R"""Confidence intervals for differences in binomial proportions.
 
-    Methods available are Agresti-Caffo [3]_, Mee [4]_, Meitinen-Nurminen [4]_ [5]_
-    and Newcombe (aka, Score method) [4]_. See also, source code for the
-    R-language function BinomDiffCI, in the module StatsAndCIs [6]_.
+    Methods available are Agresti-Caffo [4]_, Mee [5]_, Meitinen-Nurminen [5]_ [6]_
+    and Newcombe (aka, Score method) [5]_. See also, source code for the
+    R-language function BinomDiffCI, in the module StatsAndCIs [7]_.
 
     Parameters
     ----------
@@ -222,19 +222,19 @@ def propn_diff_ci(
     References
     ----------
 
-    .. [3] Agresti, A., & Caffo, T. (2000). Simple and Effective
+    .. [4] Agresti, A., & Caffo, T. (2000). Simple and Effective
        Confidence Intervals for Proportions and Differences of Proportions
        Result from Adding Two Successes and Two Failures.
        The American Statistician, 54(4), 280--288. https://doi.org/10.2307/2685779
 
-    .. [4] Newcombe, R.G. (1998). Two-sided confidence intervals for
+    .. [5] Newcombe, R.G. (1998). Two-sided confidence intervals for
        the single proportion: comparison of seven methods. Statist. Med., 17: 857-872.
        https://doi.org/10.1002/(SICI)1097-0258(19980430)17:8%3C857::AID-SIM777%3E3.0.CO;2-E
 
-    .. [5] Miettinen, O. and Nurminen, M. (1985). Comparative analysis of two rates.
+    .. [6] Miettinen, O. and Nurminen, M. (1985). Comparative analysis of two rates.
         Statist. Med., 4: 213-226. https://doi.org/10.1002/sim.4780040211; Appendix I
 
-    .. [6] StatsAndCIs.r, function BinomDiffCI, method, "mn"
+    .. [7] StatsAndCIs.r, function BinomDiffCI, method, "mn"
        https://github.com/cran/DescTools/blob/master/R/StatsAndCIs.r
        (R source code is distributed under the CC-BY license.)
 

mergeron/core/pseudorandom_numbers.py

@@ -78,8 +78,7 @@ def gen_seed_seq_list_default(
 
     References
     ----------
-    *See*,
-    https://numpy.org/doc/stable/reference/random/parallel.html
+    *See*, https://numpy.org/doc/stable/reference/random/parallel.html
 
 
     """

mergeron/examples/concentration_as_diversion.py

@@ -20,7 +20,6 @@ systems.
 
 from __future__ import annotations
 
-from collections import OrderedDict
 from collections.abc import Mapping, Sequence
 from datetime import datetime
 from math import sqrt
@@ -39,97 +38,95 @@ from mergeron import DATA_DIR
 
 PROG_PATH = Path(__file__)
 
-recapture_spec = "inside-out"
+RECAPTURE_SPEC = "inside-out"
 # Map boundary forms to titles and generating-function names, with
 #   additional parameters as relevant
-bdry_specs_dict: Mapping[str, Mapping[str, Any]] = OrderedDict({
+BDRY_SPECS_DICT: Mapping[str, Mapping[str, Any]] = {
     "ΔHHI": {
-        "title_str": "ΔHHI boundary boundary",
+        "title_str": "ΔHHI boundary",
+        "sheet_name": "ΔHHI",
         "func_str": R"\Delta HHI",
         "func": gsf.delta_hhi_boundary,
     },
     "OSWAG Own-shr-wtd Div Ratio Index": {
         "title_str": "Aggregated-diversion-ratio boundary, own-share wtd. avg.",
+        "sheet_name": "OSWAG, wtd avg",
         "func_str": R"(s_1 d_{12} + s_2 d_{21}) / s_M",
         "func": gsf.shrratio_mgnsym_boundary_wtd_avg,
-        "func_kwargs": {"wgtng_policy": "own-share", "recapture_spec": recapture_spec},
+        "func_kwargs": {"wgtng_policy": "own-share", "recapture_spec": RECAPTURE_SPEC},
     },
     "OSWAG Own-shr-wtd Div Ratio Distance": {
         "title_str": "Aggregated-diversion-ratio boundary, own-shr. wtd. distance",
+        "sheet_name": "OSWAG, distance",
         "func_str": R"\surd (s_1 d_{12}^2 / s_M + s_2 d_{21}^2 / s_M)",
         "func": gsf.shrratio_mgnsym_boundary_wtd_avg,
         "func_kwargs": {
             "wgtng_policy": "own-share",
-            "recapture_spec": recapture_spec,
-            "avg_method": "root-mean-square",
-        },
-    },
-    "OSWAG Own-shr-wtd Div Ratio Real Distance": {
-        "title_str": "Aggregated-diversion-ratio boundary, own-shr. wtd. (really) distance ",
-        "func_str": R"\surd (d_{12}^{s_1 / s_M} + d_{21}^{s_2 / s_M})",
-        "func": gsf.shrratio_mgnsym_boundary_wtd_avg,
-        "func_kwargs": {
-            "wgtng_policy": "own-share",
-            "recapture_spec": recapture_spec,
+            "recapture_spec": RECAPTURE_SPEC,
             "avg_method": "distance",
         },
     },
     "OSWAG Min Div Ratio": {
         "title_str": "Aggregated-diversion-ratio boundary, minimum",
+        "sheet_name": "OSWAG, minimum",
         "func_str": R"\min (d_{12}, d_{21})",
         "func": gsf.shrratio_mgnsym_boundary_min,
-        "func_kwargs": {"recapture_spec": recapture_spec},
+        "func_kwargs": {"recapture_spec": RECAPTURE_SPEC},
     },
     "SAG Combined Share": {
         "title_str": "Combined Share boundary",
+        "sheet_name": "SAG, combined-share",
         "func_str": R"s_M",
         "func": gsf.combined_share_boundary,
     },
     "SAG Div Ratio Distance": {
         "title_str": "Aggregated-diversion-ratio boundary, distance",
+        "sheet_name": "SAG, distance",
         "func_str": R"\surd (d_{12}^2 / 2 + d_{21}^2 / 2)",
         "func": gsf.shrratio_mgnsym_boundary_avg,
-        "func_kwargs": {
-            "recapture_spec": recapture_spec,
-            "avg_method": "root-mean-square",
-        },
+        "func_kwargs": {"recapture_spec": RECAPTURE_SPEC, "avg_method": "distance"},
     },
     "SAG Average Div Ratio": {
         "title_str": "Aggregated-diversion-ratio boundary, simple average",
+        "sheet_name": "SAG, average",
         "func_str": R"(d_{12} + d_{21}) / 2",
         "func": gsf.shrratio_mgnsym_boundary_xact_avg,
-        "func_kwargs": {"recapture_spec": recapture_spec},
+        "func_kwargs": {"recapture_spec": RECAPTURE_SPEC},
     },
     "CPSWAG Premerger HHI-contribution": {
         "title_str": "Premerger HHI-contribution boundary",
+        "sheet_name": "CPSWAG, HHI-contrib-pre",
         "func_str": R"HHI_M^{pre}",
         "func": gsf.hhi_pre_contrib_boundary,
     },
     "CPSWAG Cross-product-shr-wtd Div Ratio Index": {
         "title_str": "Aggregated-diversion-ratio boundary, cross-product-share wtd. avg.",
+        "sheet_name": "CPSWAG, wtd avg",
         "func_str": R"(s_2 d_{12} / s_M  + s_1 d_{21} / s_M)",
         "func": gsf.shrratio_mgnsym_boundary_wtd_avg,
         "func_kwargs": {
             "wgtng_policy": "cross-product-share",
-            "recapture_spec": recapture_spec,
+            "recapture_spec": RECAPTURE_SPEC,
         },
     },
     "CPSWAG Cross-product-shr-wtd Div Ratio Distance": {
         "title_str": "Aggregated-diversion-ratio boundary, cross-prod-shr. wtd. distance",
+        "sheet_name": "CPSWAG, distance",
         "func_str": R"\surd (s_2 d_{12}^2 / s_M + s_1 d_{21}^2 / s_M)",
         "func": gsf.shrratio_mgnsym_boundary_wtd_avg,
         "func_kwargs": {
             "wgtng_policy": "cross-product-share",
-            "recapture_spec": recapture_spec,
-            "avg_method": "root-mean-square",
+            "recapture_spec": RECAPTURE_SPEC,
+            "avg_method": "distance",
         },
     },
     "CPSWAG Max Div Ratio": {
         "title_str": "Aggregated-diversion-ratio boundary, maximum",
+        "sheet_name": "CPSWAG, maximum",
         "func_str": R"\max (d_{12}, d_{21})",
         "func": gsf.shrratio_mgnsym_boundary_max,
     },
-})
+}
 
 
 def tabulate_boundary_stats(_gpubyr: Literal[1992, 2010, 2023], /) -> None:
@@ -142,7 +139,9 @@ def tabulate_boundary_stats(_gpubyr: Literal[1992, 2010, 2023], /) -> None:
 
     """
     gso = gsf.GuidelinesStandards(_gpubyr)
-    _dhhi_val, _, _r_val, _g_val, *_ = gso.presumption
+    _dhhi_val, _r_val, _g_val = (
+        getattr(gso.presumption, _f) for _f in ("delta", "rec", "guppi")
+    )
 
     _dhhi_seq = (
         (0.01, 0.02, 0.03125, 0.05, _dhhi_val)
@@ -153,7 +152,7 @@ def tabulate_boundary_stats(_gpubyr: Literal[1992, 2010, 2023], /) -> None:
     _bdry_approx_data_dict = {
         "Criterion": {
             _k: R"\({} < {}\)".format(
-                bdry_specs_dict[_k]["func_str"],
+                BDRY_SPECS_DICT[_k]["func_str"],
                 R"\safeharb{d}"
                 if "Div Ratio" in _k
                 else (
@@ -162,7 +161,7 @@ def tabulate_boundary_stats(_gpubyr: Literal[1992, 2010, 2023], /) -> None:
                     else R"\safeharb{H}"
                 ),
             )
-            for _k in bdry_specs_dict
+            for _k in BDRY_SPECS_DICT
             if not _k.endswith("Distance")
         }
     }
@@ -209,7 +208,7 @@ def _dhhi_stats(
 
     _bdry_stats = Parallel(n_jobs=cpu_count() // 2)(
         delayed(_bdry_stats_col)(_bdry_spec, _dhhi_val, _delta_val, _r_val)
-        for _bdry_spec in bdry_specs_dict
+        for _bdry_spec in BDRY_SPECS_DICT
     )
 
     _bounds_string = R"{{ {} \\ {} }}".format(
@@ -234,9 +233,9 @@ def _bdry_stats_col(
         case "CPSWAG Premerger HHI-contribution":
             return _bdry_spec, f"{_hhi_m_pre_prob:6.5f}"
         case _ if "Div Ratio" in _bdry_spec:
-            _gbd_func = bdry_specs_dict[_bdry_spec]["func"]
+            _gbd_func = BDRY_SPECS_DICT[_bdry_spec]["func"]
             _, _within_bdry_area = _gbd_func(
-                _delta_val, _r_val, **bdry_specs_dict[_bdry_spec].get("func_kwargs", {})
+                _delta_val, _r_val, **BDRY_SPECS_DICT[_bdry_spec].get("func_kwargs", {})
             )
             _within_bdry_prob = 2 * _within_bdry_area
             _within_conc_bdry_prob = (
@@ -261,10 +260,11 @@ def plot_and_save_boundary_coords(
 ) -> None:
     gso = gsf.GuidelinesStandards(_gpubyr)
 
-    _hmg_standards_strings = {
+    _hmg_standards_strings_dict = {
         "distributed": ("presumption", "inferred presumption", "safeharbor"),
         "collected": ("safeharbor", "inferred_presumption", "presumption"),
-    }.get(layout, ())
+    }
+    _hmg_standards_strings = _hmg_standards_strings_dict.get(layout, ())
     if not _hmg_standards_strings:
         raise ValueError(
             f"Layout parameter value, {layout!r} is invalid.  "
@@ -279,7 +279,7 @@ def plot_and_save_boundary_coords(
     for _divr_agg_method, _hmg_standards_str in zip(
         _divr_agg_methods, _hmg_standards_strings, strict=True
     ):
-        _r_bar, _g_bar = gso.presumption[2:4]
+        _r_bar, _g_bar = (getattr(gso.presumption, _f) for _f in ("rec", "guppi"))
         _dhhi_val = getattr(gso, _hmg_standards_str)[0]
         _divr_val = (
             _g_bar
@@ -301,36 +301,22 @@ def plot_and_save_boundary_coords(
             str, Sequence[tuple[float]]
         ] = {}  #: Container for boundary coordinates data, by boundary
 
-        for _bdry_spec_key in bdry_specs_dict:
-            _bdry_spec = (_bdry_spec_key, bdry_specs_dict[_bdry_spec_key])
+        for _bdry_spec_key in BDRY_SPECS_DICT:
+            _bdry_spec = (_bdry_spec_key, BDRY_SPECS_DICT[_bdry_spec_key])
 
             if _bdry_spec_key == "ΔHHI":
-                if any((
-                    _hmg_standards_str == "inferred_presumption",
-                    _hmg_standards_str == "safeharbor" and layout == "collected",
-                    _hmg_standards_str == "presumption" and _divr_agg_method != "OSWAG",
-                )):
+                if _hmg_standards_str != _hmg_standards_strings_dict[layout][0]:
                     continue
 
                 _dh_s1, _dh_s2 = gen_plot_boundary(
                     _bndry_data_dict, gso, _hmg_standards_str, _bdry_spec, _ax1
                 )
 
-                if _hmg_standards_str == "safeharbor":
-                    _ax1.fill_between(
-                        (0, *tuple(_dh_s1), 1),
-                        (1, *tuple(_dh_s2), 0),
-                        edgecolor=None,
-                        facecolor="#64bb64",  # "#64ff64"
-                        alpha=0.7,
-                        rasterized=True,
-                    )
-
                 del _dh_s1, _dh_s2
 
             elif _bdry_spec_key.startswith(
                 _divr_agg_method
-            ) and not _bdry_spec_key.endswith("Distance"):
+            ):  # and not _bdry_spec_key.endswith("Distance"):
                 gen_plot_boundary(
                     _bndry_data_dict, gso, _hmg_standards_str, _bdry_spec, _ax1
                 )
@@ -346,8 +332,8 @@ def plot_and_save_boundary_coords(
         _fig_leg.set_in_layout(False)
 
         for _bndry_name in _bndry_data_dict:
-            if _bndry_name == "ΔHHI":  # _divr_agg_method != "OSWAG" and
-                continue
+            # if _bndry_name == "ΔHHI":  # and _divr_agg_method != "OSWAG" and
+            #     continue
 
             boundary_data_to_worksheet(
                 _bndry_name,
@@ -406,10 +392,14 @@ def gen_plot_boundary(
     _pt_mdco: ptcolor.Mcset = ptcolor.tol_cset("medium-contrast")  # type: ignore
 
     _plot_line_width = 1.0
+    _plot_line_alpha = 0.8
     _plot_line_color = _pt_mdco.black
     _plot_line_style = {"OSWAG": "-", "SAG": "-.", "CPSWAG": "--"}.get(
         _bdry_spec_str.split(" ")[0], "-"
     )
+    if _bdry_spec_str.startswith("ΔHHI"):
+        _plot_line_width = 0.5
+        _plot_line_alpha = 1.0
     _zrdr = 5
 
     if not _bdry_spec_str.startswith("ΔHHI"):
@@ -425,12 +415,12 @@ def gen_plot_boundary(
             case _:
                 _plot_line_color = _pt_mdco.dark_red
 
-    _dh_bar = _gso.safeharbor[0]
-    _g_val = _gso.safeharbor[3]
+    _dh_bar = _gso.safeharbor.divr
+    _g_val = _gso.safeharbor.guppi
 
-    _r_bar = _gso.presumption[2]
+    _r_bar = _gso.presumption.rec
 
-    _dhhi_val = getattr(_gso, _gs_str)[0]  # _gso.presumption[0]
+    _dhhi_val = getattr(_gso, _gs_str).delta
     _s_mid = sqrt(_dhhi_val / 2)
     _delta_val = _g_val / _r_bar if _gs_str == "safeharbor" else _s_mid / (1 - _s_mid)
 
@@ -456,7 +446,9 @@ def gen_plot_boundary(
         _plot_label_uom,
     )
 
-    _bndry_data_dict |= {_bdry_spec_str: (_bdry_data, _bdry_area)}  # type: ignore
+    _bndry_data_dict |= {
+        _bdry_spec_str: (_bdry_spec_dict["sheet_name"], _bdry_data, _bdry_area)
+    }  # type: ignore
     _bdry_s1, _bdry_s2 = zip(*_bdry_data, strict=True)
 
     _ax1.plot(
@@ -466,30 +458,26 @@ def gen_plot_boundary(
         color=_plot_line_color,
         linestyle=_plot_line_style,
         linewidth=_plot_line_width,
+        alpha=_plot_line_alpha,
         zorder=_zrdr,
     )
 
     print("\t", _bdry_spec_str, f"{_bdry_s2[0]:.1%}")
-    match _bdry_spec_str:
-        case "ΔHHI":
-            _plot_annotator(
-                _ax1,
-                f"({_bdry_s1[0]:.1%}, {_bdry_s2[0]:.1%})",
-                (_bdry_s1[0], _bdry_s2[0]),
-                (0.005, 0),
-                "left",
-            )
-        case _ if _bdry_spec_str.startswith("SAG") or _bdry_spec_str in (
-            "CPSWAG Premerger HHI-contribution",
-            "CPSWAG Max Div Ratio",
-        ):
-            _plot_annotator(
-                _ax1,
-                f"{_bdry_s2[0]:.1%}",
-                (_bdry_s1[0], _bdry_s2[0]),
-                (-0.005, 0),
-                "right",
-            )
+    if _bdry_spec_str.startswith(("ΔHHI", "OSWAG Min")):
+        _plot_annotator(
+            _ax1,
+            f"({_bdry_s1[1]:.1%}, {_bdry_s2[1]:.1%})",
+            (_bdry_s1[1], _bdry_s2[1]),
+            (0.005, 0),
+            "left",
+        )
+    elif _bdry_spec_str.startswith("SAG") or _bdry_spec_str in (
+        "CPSWAG Premerger HHI-contribution",
+        "CPSWAG Max Div Ratio",
+    ):
+        _plot_annotator(
+            _ax1, f"{_bdry_s2[0]:.1%}", (_bdry_s1[0], _bdry_s2[0]), (-0.005, 0), "right"
+        )
 
     return _bdry_s1, _bdry_s2
 
@@ -546,9 +534,9 @@ def boundary_data_to_worksheet(
         Specified Excel Workbook
 
     """
-    _xl_sheet = _xl_book.add_worksheet(
-        f"{_bndry_name.replace('Boundary', '').strip()}"[:31]
-    )
+    _sheet_name, _bndry_points, _bndry_area = _bndry_data_dict[_bndry_name]
+
+    _xl_sheet = _xl_book.add_worksheet(_sheet_name)
 
     _xl_sheet.write("A1", "Sound GUPPI Safeharbor")
     _xl_sheet.write("A2", "Merger Screens for Unilateral Effects")
@@ -574,7 +562,6 @@ def boundary_data_to_worksheet(
     )
     _xl_sheet.set_footer(f"&L{_left_footer}")
 
-    _bndry_points, _bndry_area = _bndry_data_dict[_bndry_name]
     xlh.scalar_to_sheet(_xl_book, _xl_sheet, "B7", _bndry_area, xlh.CFmt.AREA_NUM)
 
     _results_header_row = 9

mergeron/gen/data_generation.py

@@ -1178,7 +1178,7 @@ def _beta_located(
     """
     Given mean and stddev, return shape parameters for corresponding Beta distribution
 
-    Solve the first two moments of the standard Beta to get the shape parameters. [1]_
+    Solve the first two moments of the standard Beta to get the shape parameters.
 
     Parameters
     ----------
@@ -1191,12 +1191,9 @@ def _beta_located(
     -------
         shape parameters for Beta distribution
 
-    References
-    ----------
-    .. [1] NIST. https://www.itl.nist.gov/div898/handbook/eda/section3/eda366h.htm
-
     """
-    _mul = (_mu - _mu**2 - _sigma**2) / _sigma**2
+
+    _mul = -1 + _mu * (1 - _mu) / _sigma**2
     return np.array([_mu * _mul, (1 - _mu) * _mul], dtype=np.float64)
 
 
@@ -1208,7 +1205,7 @@ def beta_located_bound(_dist_parms: NDArray[np.floating[TF]], /) -> NDArray[np.f
     Recover the r.v.s as
     :math:`\min + (\max - \min) \cdot \symup{Β}(a, b)`,
     with `a` and `b` calculated from the specified mean (:math:`\mu`) and
-    variance (:math:`\sigma`). [7]_
+    variance (:math:`\sigma`). [8]_
 
     Parameters
     ----------
@@ -1225,7 +1222,7 @@ def beta_located_bound(_dist_parms: NDArray[np.floating[TF]], /) -> NDArray[np.f
 
     References
     ----------
-    .. [7] NIST. https://www.itl.nist.gov/div898/handbook/eda/section3/eda366h.htm
+    .. [8] NIST, Beta Distribution. https://www.itl.nist.gov/div898/handbook/eda/section3/eda366h.htm
     """  # noqa: RUF002
 
     _bmu, _bsigma, _bmin, _bmax = _dist_parms

{mergeron-2024.738936.0.dist-info → mergeron-2024.738936.2.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: mergeron
-Version: 2024.738936.0
+Version: 2024.738936.2
 Summary: Analysis of standards defined in Horizontal Merger Guidelines
 License: MIT
 Keywords: merger policy analysis,merger guidelines,merger screening,policy presumptions,concentration standards,upward pricing pressure,GUPPI

{mergeron-2024.738936.0.dist-info → mergeron-2024.738936.2.dist-info}/RECORD

@@ -5,11 +5,11 @@ mergeron/core/__init__.py,sha256=iyfxkX3-SoMS4ZQZKHKPn8JEMN536vpty9oSZf0LHv8,115
 mergeron/core/damodaran_margin_data.py,sha256=N1SU_PDjFMoCzCqone-If-gR4PpcG-MK5TfayusOJLs,8166
 mergeron/core/excel_helper.py,sha256=Zx2pZVjrCV3y0IlpX-_18K3-gbvXxJgnJqPnST7Fg8E,7471
 mergeron/core/ftc_merger_investigations_data.py,sha256=C37gFwSSgUFu98FW8ZDlUTgEaz4qYa1aKPC3SPtXChU,26553
-mergeron/core/guidelines_standards.py,sha256=nPcoKQDVQi_aVYAuxXmOiSCX9Yg8DMQOtYowui2xIC4,36899
-mergeron/core/proportions_tests.py,sha256=PHcaHuWzZC__u4PjN1U11dghL-5G16P6QQhxKJkR77o,15275
-mergeron/core/pseudorandom_numbers.py,sha256=fJQPY9DDqMlQejLyVfjrSylfjpHHn6qoHV_SHGClADQ,9434
+mergeron/core/guidelines_standards.py,sha256=VrPH-WApl7o0AolomnJwx4UqZO1rEzEUiIRT8TRb4_c,43316
+mergeron/core/proportions_tests.py,sha256=zALaWNrGOltsCp2ZSP4pcUfrOCaOoR5YDavqxQHADKA,15275
+mergeron/core/pseudorandom_numbers.py,sha256=ZnIE1ixp3JcgJjlulgMsIcWFVAxeyjC9lsiMIk7OHlM,9430
 mergeron/examples/__init__.py,sha256=iyfxkX3-SoMS4ZQZKHKPn8JEMN536vpty9oSZf0LHv8,115
-mergeron/examples/concentration_as_diversion.py,sha256=5FXpQXSCmKxXiahCy-ga3BCMAa7WdOu-cUZaX2V1N7M,21183
+mergeron/examples/concentration_as_diversion.py,sha256=3re3GKMbV0_g2DFb9_1DXmQZpJeoedAbo2s39GFld_o,20692
 mergeron/examples/example_parameterizations.py,sha256=VP-hi7L0j30ffcEzmJ3P8mOj1VjwEWKCTZSx_CaVQxA,4197
 mergeron/examples/guidelines_enforcement_patterns.py,sha256=gpAW9jLLg3bbf4ueCu71QUKIca_eb3tTtAuTwY14m3c,2235
 mergeron/examples/investigations_stats_obs_tables.py,sha256=IDUZqZHiNvyO05d-yaYUJ_3-G_tchA4Z4MTaEICV9g0,17801
@@ -25,7 +25,7 @@ mergeron/examples/visualize_guidelines_tests.py,sha256=CDnDT8fIu2Ibed4NsUamzamLt
 mergeron/ext/__init__.py,sha256=iyfxkX3-SoMS4ZQZKHKPn8JEMN536vpty9oSZf0LHv8,115
 mergeron/ext/tol_colors.py,sha256=wFOHZXWZonbp9mhmSGu9mVujBYhdTsvx9_WikMpoCmo,22229
 mergeron/gen/__init__.py,sha256=iyfxkX3-SoMS4ZQZKHKPn8JEMN536vpty9oSZf0LHv8,115
-mergeron/gen/data_generation.py,sha256=E0KWhvqU9UKcwqJT1IzsQwmXou-HIgRSHMkhVlyelYI,41453
+mergeron/gen/data_generation.py,sha256=9D-ieJU6oYXf6vMckRRqjz-qdoA7w_3Y4BamlfLCUUU,41348
 mergeron/gen/guidelines_tests.py,sha256=eoDQFIi6CZzSIyxZsaK5pdtwuIXHt8QwBQSV4VGpKs8,15164
 mergeron/gen/investigations_stats.py,sha256=evm8bwaNUPRe9V-RR_DgsumoUYNLGKmkrhSpaBNXXfA,23005
 mergeron/jinja_LaTex_templates/clrrate_cis_summary_table_template.tex.jinja2,sha256=m1fVYQC4Ra3OLk2W_EFeCzwYZ_Ihy9wZbAHs8rBfs00,4836
@@ -36,6 +36,6 @@ mergeron/jinja_LaTex_templates/mergeron.cls,sha256=AV2mk4-uERvAuMkE95Ka7el6LZsb0
 mergeron/jinja_LaTex_templates/mergeron_table_collection_template.tex.jinja2,sha256=nr6xUI0_2KHG4Sz9k1JFVQjs2h9qS9BGt1MeE6Tygs8,2429
 mergeron/jinja_LaTex_templates/setup_tikz_tables.tex.jinja2,sha256=WKVxtp3eoMchfGliQAJMj4w2FtBkWG5z2V3-hBYUYUQ,3292
 mergeron/py.typed,sha256=frcCV1k9oG9oKj3dpUqdJg1PxRT2RSN_XKdLCPjaYaY,2
-mergeron-2024.738936.0.dist-info/METADATA,sha256=iasUgfZLsctlK777sCoJLBYIx02nHamNWBxrA53XVNo,6452
-mergeron-2024.738936.0.dist-info/WHEEL,sha256=FMvqSimYX_P7y0a7UY-_Mc83r5zkBZsCYPm7Lr0Bsq4,88
-mergeron-2024.738936.0.dist-info/RECORD,,
+mergeron-2024.738936.2.dist-info/METADATA,sha256=Ths3nxJp21K9_EB460S3yb88-mF7-VZtTEmITNze86s,6452
+mergeron-2024.738936.2.dist-info/WHEEL,sha256=FMvqSimYX_P7y0a7UY-_Mc83r5zkBZsCYPm7Lr0Bsq4,88
+mergeron-2024.738936.2.dist-info/RECORD,,