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

mergeron/examples/safeharbor_boundaries_for_symmetric_firm_mergers.py

@@ -6,17 +6,16 @@ This boundary is limited to share-margin space, given restrictions from symmetry
 
 from pathlib import Path
 
+import matplotlib as mpl
 import numpy as np
 from matplotlib import cm as colormgr
 from matplotlib import colors as mcolors
 from matplotlib.ticker import AutoMinorLocator, MultipleLocator, StrMethodFormatter
 
 import mergeron.core.guidelines_standards as gsf
+from mergeron import DATA_DIR
 
 PROG_PATH = Path(__file__)
-data_path = Path.home() / PROG_PATH.parents[1].stem
-if not data_path.is_dir():
-    data_path.mkdir(parents=True)
 
 
 def _main() -> None:
@@ -43,14 +42,11 @@ def _main() -> None:
 
     print("Generate data for plots")
     # Mgn coords, and div-ratio coords
-    dh_bar, g_bar, r_bar, dr_bar = 0.01, 0.06, 0.80, 0.20
+    g_bar, r_bar, dr_bar = 0.06, 0.80, 0.20
 
     step_size = 10**-5
     sym_shr_vec = np.arange(gsf.shr_from_gbd(g_bar), 0.5 + step_size, step_size)
     mst_vec = (g_bar / r_bar) * (1 - sym_shr_vec) / sym_shr_vec
-    mst_vec_gsf = np.where(
-        sym_shr_vec > gsf.shr_from_gbd(dr_bar, m_star=1.00), np.nan, mst_vec
-    )
 
     print("Setup basic figure and axes for plots of safe harbor boundaries.")
 
@@ -120,7 +116,7 @@ def _main() -> None:
         linestyle="-",
         linewidth=0.75,
         alpha=1,
-        color=colormgr.cividis(dr_bar_mgn[0]),
+        color=mpl.colormaps["cividis"](dr_bar_mgn[0]),
         zorder=5.5,
     )
 
@@ -173,10 +169,7 @@ def _main() -> None:
     _minorLocator = AutoMinorLocator(5)
     _majorLocator = MultipleLocator(0.05)
     for _axs in ax1.xaxis, ax1.yaxis:
-        if _axs == ax1.xaxis:
-            _majorticklabels_rot = 45
-        elif _axs == ax1.yaxis:
-            _majorticklabels_rot = 0
+        _majorticklabels_rot = 45 if _axs == ax1.xaxis else 0
         # x-axis
         _axs.set_major_locator(_majorLocator)
         _axs.set_minor_locator(_minorLocator)
@@ -206,7 +199,7 @@ def _main() -> None:
     cm_plot.ax.set_ylim(0, 1.0)
     cm_plot.outline.set_visible(False)
 
-    plt.savefig(data_path.joinpath(f"{PROG_PATH.stem}.pdf"))
+    plt.savefig(DATA_DIR / PROG_PATH.with_suffix(".pdf").name)
 
 
 if __name__ == "__main__":

mergeron/examples/sound_guppi_safeharbor.py

@@ -13,15 +13,15 @@ import mergeron.gen.guidelines_tests as gtl
 import mergeron.gen.investigations_stats as isl
 from mergeron import DATA_DIR
 
-tests_of_interest: tuple[gtl.UPPTestSpec, ...] = (
+tests_of_interest: tuple[gtl.UPPTestRegime, ...] = (
     (isl.PolicySelector.CLRN, gtl.GUPPIWghtngSelector.MAX, gtl.GUPPIWghtngSelector.MAX),
     (isl.PolicySelector.ENFT, gtl.GUPPIWghtngSelector.MIN, gtl.GUPPIWghtngSelector.MIN),
 )
 
-mod_path = Path(__file__)
+PROG_PATH = Path(__file__)
 
 
-def analyze_inv_data(
+def analyze_invres_data(
     _sample_size: int = 10**6,
     _hmg_std_pub_year: Literal[1992, 2010, 2023] = 1992,
     _test_sel: tuple[
@@ -50,16 +50,16 @@ def analyze_inv_data(
         If True, simulated data are save to file (hdf5 format)
 
     """
-    _inv_parm_vec = gsf.GuidelinesStandards(_hmg_std_pub_year).presumption[2:]
+    _invres_parm_vec = gsf.GuidelinesStandards(_hmg_std_pub_year).presumption
 
-    _save_data_to_file: Literal[False] | tuple[Literal[True], ptb.File, str] = False
+    _save_data_to_file: gtl.SaveData = False
     if save_data_to_file_flag:
         _h5_hier_pat = re.compile(r"\W")
         _blosc_filters = ptb.Filters(
             complevel=3, complib="blosc:lz4hc", fletcher32=True
         )
         _h5_datafile = ptb.open_file(
-            str(DATA_DIR / f"{mod_path.stem}_sound.h5"),
+            DATA_DIR / PROG_PATH.with_suffix(".h5").name,
             mode="w",
             title=f"GUPPI Safeharbor {_test_sel[0].capitalize()} Rate Module",
             filters=_blosc_filters,
@@ -115,27 +115,32 @@ def analyze_inv_data(
 
         _ind_sample_spec = dgl.MarketSampleSpec(
             _sample_size,
-            _inv_parm_vec[0],
+            _invres_parm_vec.rec,
             dgl.PRIConstants.SYM,
             share_spec=dgl.ShareSpec(
-                dgl.SHRConstants.UNI, _recapture_spec_test, dgl.EMPTY_ARRAY_DEFAULT
+                _recapture_spec_test,
+                dgl.SHRConstants.UNI,
+                dgl.EMPTY_ARRAY_DEFAULT,
+                dgl.FCOUNT_WTS_DEFAULT,
             ),
             pcm_spec=dgl.PCMSpec(
                 _pcm_dist_type_test, _pcm_dist_firm2_test, _pcm_dist_parms_test
             ),
         )
 
-        _inv_cnts_kwargs = {
-            "sim_inv_sel": _test_sel,
+        _invres_cnts_kwargs = {
+            "sim_test_regime": _test_sel,
             "save_data_to_file": _save_data_to_file,
         }
 
         _start_time = datetime.now()
         (
-            _inv_rate_sim_byfirmcount_array,
-            _inv_rate_sim_bydelta_array,
-            _inv_rate_sim_byconczone_array,
-        ) = gtl.sim_inv_cnts_ll(_inv_parm_vec, _ind_sample_spec, _inv_cnts_kwargs)
+            _invres_rate_sim_byfirmcount_array,
+            _invres_rate_sim_bydelta_array,
+            _invres_rate_sim_byconczone_array,
+        ) = gtl.sim_invres_cnts_ll(
+            _invres_parm_vec, _ind_sample_spec, _invres_cnts_kwargs
+        )
         _run_duration = datetime.now() - _start_time
         print(
             f"Simulation completed in {_run_duration / timedelta(seconds=1):.6f} secs.",
@@ -143,8 +148,8 @@ def analyze_inv_data(
             sep=", ",
         )
 
-        _stats_hdr_list, _stats_dat_list = isl.latex_tbl_inv_stats_1dim(
-            _inv_rate_sim_bydelta_array,
+        _stats_hdr_list, _stats_dat_list = isl.latex_tbl_invres_stats_1dim(
+            _invres_rate_sim_bydelta_array,
             return_type_sel=isl.StatsReturnSelector.RPT,
             sort_order=isl.SortSelector.REV,
         )
@@ -162,8 +167,8 @@ def analyze_inv_data(
         del _pcm_dist_type_test, _pcm_dist_parms_test
 
     if save_data_to_file_flag:
-        _h5_datafile.close()
+        _save_data_to_file[1].close()  # type: ignore
 
 
 if __name__ == "__main__":
-    analyze_inv_data(10**7, 2023, tests_of_interest[1], save_data_to_file_flag=False)
+    analyze_invres_data(10**7, 2023, tests_of_interest[1], save_data_to_file_flag=False)

mergeron/examples/testIntrinsicClearanceRates.py

@@ -1,11 +1,11 @@
 """
 
-    Estimation of intrinsic clearance rates for ΔHHI safeharbor and GUPPI safeharbor
-    Demonstrates MC integration with closed-form integrand and indicator-function
-    integrand (Tested with Python 3.8 and 3.9; NumPy module is required.)
+Estimation of intrinsic clearance rates for ΔHHI safeharbor and GUPPI safeharbor
+Demonstrates MC integration with closed-form integrand and indicator-function
+integrand (Tested with Python 3.8 and 3.9; NumPy module is required.)
 
-    S. Murthy Kambhampaty, © 2019. This work is licensed under a CC BY-NC-SA 4.0 License
-    https://creativecommons.org/licenses/by-nc-sa/4.0/
+S. Murthy Kambhampaty, © 2019. This work is licensed under a CC BY-NC-SA 4.0 License
+https://creativecommons.org/licenses/by-nc-sa/4.0/
 
 """
 

mergeron/examples/visualize_empirical_margin_distribution.py

@@ -97,4 +97,4 @@ if __name__ == "__main__":
     mgn_ax.set_ylabel("Relative Frequency", fontsize=10)
 
     mgn_fig.tight_layout()
-    plt.savefig(DATA_DIR / f"{PROG_PATH.stem}.pdf")
+    plt.savefig(PROG_PATH.with_suffix(".pdf"))

mergeron/examples/{visualize_guidelines_tests_scatterplots.py → visualize_guidelines_tests.py}

@@ -19,7 +19,6 @@ import numpy as np
 import tables as ptb  # type: ignore
 from matplotlib import cm, colors
 from matplotlib.ticker import StrMethodFormatter
-from numpy import arange, argsort, column_stack, einsum, ones_like
 from numpy.typing import NDArray
 
 import mergeron.core.guidelines_standards as gsf
@@ -39,9 +38,9 @@ blosc_filters = ptb.Filters(
 
 def gen_plot_data(
     _market_data: dgl.MarketsSample,
-    _g_bar: float,
+    _std_vec: gsf.GuidelinesSTD,
     _pcm_firm2_star: float = 0.30,
-    _inv_sel: gtl.UPPTestSpec = (
+    _test_regime: gtl.UPPTestRegime = (
         isl.PolicySelector.CLRN,
         gtl.GUPPIWghtngSelector.MAX,
         None,
@@ -54,29 +53,20 @@ def gen_plot_data(
         f"{_pcm_firm2_star * 100:03.1f}".replace(".", "dot")
     )
 
-    _pcm_array = column_stack((
+    _pcm_array = np.column_stack((
         _m1 := _market_data.pcm_array[:, [0]],
-        _pcm_firm2_star * ones_like(_m1),
+        _pcm_firm2_star * np.ones_like(_m1),
     ))
     del _m1
 
-    _guppi_array = einsum("ij,ij->ij", _pcm_array[:, ::-1], _market_data.divr_array)
-
-    _guppi_est = (
-        _guppi_array.max(axis=1, keepdims=True)
-        if _inv_sel[1] == gtl.GUPPIWghtngSelector.MAX
-        else _guppi_array.min(axis=1, keepdims=True)
-    )
-
-    _gbd_test = (
-        (_guppi_est < _g_bar)
-        if _inv_sel[0] == isl.PolicySelector.CLRN
-        else (_guppi_est >= _g_bar)
+    _upp_tests = gtl.gen_upp_arrays(
+        _std_vec,
+        dgl.MarketsSample(_market_data.frmshr_array, _pcm_array, *_market_data[2:]),
+        _test_regime,
     )
 
-    del _guppi_est
-
-    _gbd_test_rows = np.where(_gbd_test)[0]
+    _gbd_test_rows = np.where(_upp_tests.guppi_test_simple)[0]
+    del _upp_tests
 
     _qtyshr_firm1_inv, _qtyshr_firm2_inv = (
         _market_data.frmshr_array[_gbd_test_rows][:, [0]],
@@ -86,7 +76,7 @@ def gen_plot_data(
         _pcm_array[_gbd_test_rows][:, [0]],
         _pcm_array[_gbd_test_rows][:, [1]],
     )
-    del _gbd_test, _gbd_test_rows
+    del _gbd_test_rows
 
     _pcm_plotter = _pcm_firm1_inv
 
@@ -108,8 +98,8 @@ def gen_plot_data(
                 title=f"{_array_name}",
             )
 
-    _pcm_sorter = argsort(_pcm_plotter, axis=0)
-    if _inv_sel[0] != isl.PolicySelector.CLRN:
+    _pcm_sorter = np.argsort(_pcm_plotter, axis=0)
+    if _test_regime[0] != isl.PolicySelector.CLRN:
         _pcm_sorter = _pcm_sorter[::-1, :]
     _qtyshr_firm1_plotter = _qtyshr_firm1_inv[_pcm_sorter]
     _qtyshr_firm2_plotter = _qtyshr_firm2_inv[_pcm_sorter]
@@ -130,13 +120,16 @@ def gen_plot_data(
 def _main(
     _hmg_pub_year: gsf.HMGPubYear,
     _market_sample_spec: dgl.MarketSampleSpec,
-    _inv_sel: gtl.UPPTestSpec,
+    _test_regime: gtl.UPPTestRegime,
     _save_data_to_file: gtl.SaveData,
 ) -> None:
-    _r_bar, _g_bar, _divr_bar, *_ = getattr(
+    guidelins_std_vec = getattr(
         gsf.GuidelinesStandards(_hmg_pub_year),
-        "safeharbor" if _inv_sel[0] == isl.PolicySelector.ENFT else "presumption",
-    )[2:]
+        "safeharbor" if _test_regime[0] == isl.PolicySelector.ENFT else "presumption",
+    )
+
+    _, _r_bar, _g_bar, _divr_bar, *_ = guidelins_std_vec
+
     market_data = dgl.gen_market_sample(_market_sample_spec, seed_seq_list=None)
 
     # Set up a plot grid to fill in the various scatterplots
@@ -173,12 +166,12 @@ def _main(
         _ax_now.text(
             0.81,
             0.72,
-            "\n".join([
-                r"$m_2 = m^* = {0:.{1}f}\%$".format(
+            "\n".join((
+                R"$m_2 = m^* = {0:.{1}f}\%$".format(
                     (_pcmv := _pcm_firm2_star * 100), 1 * (_pcmv % 1 > 0)
                 ),
-                r"$m_1 \neq m_2$",
-            ]),
+                R"$m_1 \neq m_2$",
+            )),
             rotation=0,
             ha="right",
             va="top",
@@ -203,9 +196,9 @@ def _main(
 
         _qtyshr_firm1_plotter, _qtyshr_firm2_plotter, _pcm_plotter = gen_plot_data(
             market_data,
-            _g_bar,
+            guidelins_std_vec,
             _pcm_firm2_star,
-            _inv_sel,
+            _test_regime,
             h5handle=_save_data_to_file[1] if _save_data_to_file else None,
         )
 
@@ -231,7 +224,7 @@ def _main(
         ax=_ax_cm,
         orientation="vertical",
         fraction=3.0,
-        ticks=arange(0, 1.2, 0.2),
+        ticks=np.arange(0, 1.2, 0.2),
         format=StrMethodFormatter("{x:>3.0%}"),
     )
     _cm_plot.set_label(label="Firm 1 Price-Cost Margin, $m_1$", fontsize=10)
@@ -242,26 +235,30 @@ def _main(
 
     _cm_plot.outline.set_visible(False)
 
-    _base_name = DATA_DIR.joinpath(
-        f"{PROG_PATH.stem}_{_hmg_pub_year}gbar{f"{_g_bar * 100:02.0f}"}PCT_{market_sample_spec.share_spec[0]}Recapture_{_inv_sel}"
+    _base_name = DATA_DIR / "{}_{}Rate_{}gbar{}PCT_{}Recapture".format(
+        PROG_PATH.stem,
+        f"{_test_regime[0]}".capitalize(),
+        _hmg_pub_year,
+        f"{_g_bar * 100:02.0f}",
+        market_sample_spec.share_spec[0],
     )
     _my_fig_2dsg_savepath = DATA_DIR / f"{_base_name}_2DScatterGrid.pdf"
-    print(f"Save 2D plot to, {_my_fig_2dsg_savepath!r}")
+    print(f"Save 2D plot to, {f'"{_my_fig_2dsg_savepath}"'}")
     _fig_2dsg.savefig(_my_fig_2dsg_savepath, dpi=600)
 
 
 if __name__ == "__main__":
     # Get Guidelines parameter values
     hmg_pub_year: Final = 2023
-    inv_sel: gtl.UPPTestSpec = (
+    test_regime: gtl.UPPTestRegime = (
         isl.PolicySelector.ENFT,
         gtl.GUPPIWghtngSelector.MIN,
         gtl.GUPPIWghtngSelector.MIN,
     )
-    r_bar, g_bar, divr_bar, *_ = getattr(
+    r_bar = getattr(
         gsf.GuidelinesStandards(hmg_pub_year),
-        "safeharbor" if inv_sel[0] == isl.PolicySelector.ENFT else "presumption",
-    )[2:]
+        "presumption" if test_regime[0] == isl.PolicySelector.ENFT else "safeharbor",
+    ).rec
 
     sample_sz = 10**7
 
@@ -269,16 +266,13 @@ if __name__ == "__main__":
         sample_sz,
         r_bar,
         share_spec=dgl.ShareSpec(
-            dgl.RECConstants.INOUT,
-            dgl.SHRConstants.UNI,
-            DIST_PARMS_DEFAULT,
-            dgl.FCOUNT_WTS_DEFAULT,
+            dgl.RECConstants.INOUT, dgl.SHRConstants.UNI, DIST_PARMS_DEFAULT, None
         ),
     )
 
     save_data_to_file_flag = False
     if save_data_to_file_flag:
-        h5path = DATA_DIR / f"{PROG_PATH.stem}.h5"
+        h5path = DATA_DIR / PROG_PATH.with_suffix(".h5").name
         h5datafile = ptb.open_file(
             h5path,
             mode="w",
@@ -293,7 +287,7 @@ if __name__ == "__main__":
     else:
         save_data_to_file = False
 
-    _main(hmg_pub_year, market_sample_spec, inv_sel, save_data_to_file)
+    _main(hmg_pub_year, market_sample_spec, test_regime, save_data_to_file)
 
     if save_data_to_file_flag:
-        h5datafile.close()
+        save_data_to_file[1].close()  # type: ignore

mergeron/ext/__init__.py

@@ -1,5 +1,5 @@
 from importlib.metadata import version
 
-from .. import _PKG_NAME
+from .. import _PKG_NAME  # noqa: TID252
 
 __version__ = version(_PKG_NAME)

mergeron/gen/__init__.py

@@ -1,5 +1,5 @@
 from importlib.metadata import version
 
-from .. import _PKG_NAME
+from .. import _PKG_NAME  # noqa: TID252
 
 __version__ = version(_PKG_NAME)

mergeron/gen/data_generation.py

@@ -12,12 +12,12 @@ from .. import _PKG_NAME  # noqa: TID252
 __version__ = version(_PKG_NAME)
 
 import enum
-from typing import Literal, NamedTuple
+from typing import Literal, NamedTuple, TypeVar
 
 import attrs
 import numpy as np
 from numpy.random import SeedSequence
-from numpy.typing import NDArray
+from numpy.typing import NBitBase, NDArray
 
 from ..core.damodaran_margin_data import resample_mgn_data  # noqa: TID252
 from ..core.pseudorandom_numbers import (  # noqa: TID252
@@ -29,6 +29,8 @@ from ..core.pseudorandom_numbers import (  # noqa: TID252
 EMPTY_ARRAY_DEFAULT = np.zeros(2)
 FCOUNT_WTS_DEFAULT = ((_nr := np.arange(1, 6)[::-1]) / _nr.sum()).astype(np.float64)
 
+TF = TypeVar("TF", bound=NBitBase)
+
 
 @enum.unique
 class PRIConstants(tuple[bool, str | None], enum.ReprEnum):
@@ -378,7 +380,7 @@ class ShareDataSample(NamedTuple):
 class PriceDataSample(NamedTuple):
     """Container for generated price array, and related."""
 
-    price_array: NDArray[np.float64]
+    price_array: NDArray[np.floating]
     """Merging-firms' prices"""
 
     hsr_filing_test: NDArray[np.bool_]
@@ -660,7 +662,7 @@ def _gen_share_data(
 
 def _gen_market_shares_uniform(
     _s_size: int = 10**6,
-    _dist_parms_mktshr: NDArray[np.floating] | None = DIST_PARMS_DEFAULT,
+    _dist_parms_mktshr: NDArray[np.floating[TF]] | None = DIST_PARMS_DEFAULT,
     _mktshr_rng_seed_seq: SeedSequence | None = None,
     _nthreads: int = 16,
     /,
@@ -685,7 +687,9 @@ def _gen_market_shares_uniform(
     """
 
     _frmshr_array = np.empty((_s_size, 2), dtype=np.float64)
-    _dist_parms_mktshr = _dist_parms_mktshr or DIST_PARMS_DEFAULT
+    _dist_parms_mktshr = (
+        DIST_PARMS_DEFAULT if _dist_parms_mktshr is None else _dist_parms_mktshr  # type: ignore
+    )
     _mrng = MultithreadedRNG(
         _frmshr_array,
         dist_type="Uniform",
@@ -723,8 +727,8 @@ def _gen_market_shares_dirichlet_multisample(
     _s_size: int = 10**6,
     _recapture_spec: RECConstants = RECConstants.INOUT,
     _dist_type_dir: SHRConstants = SHRConstants.DIR_FLAT,
-    _dist_parms_dir: NDArray[np.floating] | None = None,
-    _firm_count_wts: NDArray[np.floating] | None = None,  # type: ignore
+    _dist_parms_dir: NDArray[np.floating[TF]] | None = None,
+    _firm_count_wts: NDArray[np.floating[TF]] | None = None,  # type: ignore
     _fcount_rng_seed_seq: SeedSequence | None = None,
     _mktshr_rng_seed_seq: SeedSequence | None = None,
     _nthreads: int = 16,
@@ -858,7 +862,7 @@ def _gen_market_shares_dirichlet_multisample(
 
 
 def _gen_market_shares_dirichlet(
-    _dir_alphas: NDArray[np.floating],
+    _dir_alphas: NDArray[np.floating[TF]],
     _s_size: int = 10**6,
     _recapture_spec: RECConstants = RECConstants.INOUT,
     _mktshr_rng_seed_seq: SeedSequence | None = None,
@@ -940,8 +944,8 @@ def _gen_market_shares_dirichlet(
 
 
 def _gen_pr_ratio(
-    _frmshr_array: NDArray[np.floating],
-    _nth_firm_share: NDArray[np.floating],
+    _frmshr_array: NDArray[np.floating[TF]],
+    _nth_firm_share: NDArray[np.floating[TF]],
     _mkt_sample_spec: MarketSampleSpec,
     _seed_seq: SeedSequence | None = None,
     /,
@@ -1030,10 +1034,10 @@ def _gen_pr_ratio(
 
 
 def gen_divr_array(
-    _frmshr_array: NDArray[np.floating],
+    _frmshr_array: NDArray[np.floating[TF]],
     _r_bar: float,
     _recapture_spec: RECConstants = RECConstants.INOUT,
-    _aggregate_purchase_prob: NDArray[np.floating] = EMPTY_ARRAY_DEFAULT,
+    _aggregate_purchase_prob: NDArray[np.floating[TF]] = EMPTY_ARRAY_DEFAULT,
     /,
 ) -> NDArray[np.float64]:
     """
@@ -1090,10 +1094,10 @@ def gen_divr_array(
 
 
 def _gen_pcm_data(
-    _frmshr_array: NDArray[np.floating],
+    _frmshr_array: NDArray[np.floating[TF]],
     _mkt_sample_spec: MarketSampleSpec,
-    _price_array: NDArray[np.floating],
-    _aggregate_purchase_prob: NDArray[np.floating],
+    _price_array: NDArray[np.floating[TF]],
+    _aggregate_purchase_prob: NDArray[np.floating[TF]],
     _pcm_rng_seed_seq: SeedSequence,
     _nthreads: int = 16,
     /,
@@ -1174,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
     ----------
@@ -1187,16 +1191,13 @@ 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)
 
 
-def beta_located_bound(_dist_parms: NDArray[np.floating], /) -> NDArray[np.float64]:
+def beta_located_bound(_dist_parms: NDArray[np.floating[TF]], /) -> NDArray[np.float64]:
     R"""
     Return shape parameters for a non-standard beta, given the mean, stddev, range
 
@@ -1204,7 +1205,7 @@ def beta_located_bound(_dist_parms: NDArray[np.floating], /) -> NDArray[np.float
     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
     ----------
@@ -1221,7 +1222,7 @@ def beta_located_bound(_dist_parms: NDArray[np.floating], /) -> NDArray[np.float
 
     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