mergeron 2024.738953.0__py3-none-any.whl → 2024.738953.1__py3-none-any.whl

mergeron/examples/example_parameterizations.py

@@ -1,143 +0,0 @@
-"""
-
-A few parameters/relations employed in policy analysis using this pacakage
-are demonstrated here
-
-"""
-
-from typing import NamedTuple
-
-import numpy as np
-from numpy.typing import NDArray
-
-from mergeron.core.pseudorandom_numbers import prng
-from mergeron.gen.data_generation import (
-    RECConstants,
-    _gen_market_shares_dirichlet,
-    _gen_market_shares_uniform,
-)
-
-
-def gen_rval_mnl(
-    _ssz: int = 10**8, _r_bar: float = 0.80, *, pcm_dist_type: str = "Uniform"
-) -> NamedTuple:
-    """
-
-    :param _ssz: sample size
-    :param _r_bar: recapture rate
-    :param pcm_dist_type: margin distribution name
-    :return:
-    """
-
-    # Define output type
-    class RvalMNL(NamedTuple):
-        qtyshr_array: NDArray[np.floating]
-        pcm_array: NDArray[np.floating]
-        mnl_test_rows: NDArray[np.integer]
-
-    _qtyshr_array = _gen_market_shares_uniform(_ssz).mktshr_array[:, :2]
-
-    _qtyshr_min = _qtyshr_array.min(axis=1, keepdims=True, initial=None)
-    _cprob = np.divide(_r_bar, 1 - (1 - _r_bar) * _qtyshr_min)
-    del _qtyshr_min
-
-    _purchprob_array = _cprob * _qtyshr_array
-
-    _pcm0: NDArray[np.floating]
-    if pcm_dist_type == "Uniform":
-        _pcm0 = prng().uniform(size=(_ssz, 1))
-    elif pcm_dist_type == "Beta":
-        _pcm0 = prng().beta(10, 10, size=(_ssz, 1))
-    else:
-        raise ValueError("Invalid type for distribution of margins")
-
-    _pcm1 = np.divide(
-        _pcm0 * (1 - _purchprob_array[:, [0]]), (1 - _purchprob_array[:, [1]])
-    )
-
-    _mnl_test_rows = _pcm1.__gt__(0) & _pcm1.__lt__(1)
-    _pcm_array = np.column_stack((_pcm0, _pcm1))
-
-    return RvalMNL(_qtyshr_array, _pcm_array, _mnl_test_rows)
-
-
-def gen_rval_ssp(_ssz: int = 10**6, _r_bar: float = 0.8) -> NamedTuple:
-    R"""Given r_1 (_r_bar), generates r_2 and diversion ratios under share proportionality
-
-    More testing with::
-
-        from sympy import solve, symbols
-        s1, s2, m1, mstar, rbar, dstar = symbols('s_1, s_2, m_1, m^{*}, \overline{r}, \delta^{*}')
-        infl2_pts = solve(
-            (
-                s2 - dstar + dstar * rbar * s1 + dstar * (1 - rbar) * s2,
-                s1 - (1 - s2) * dstar * mstar / m1
-            ),
-            s1, s2)
-
-    Substitute various values into the solution to test against plots generated elsewhere
-    """
-
-    # Define output type
-    class RvalSSP(NamedTuple):
-        """Container for share, recapture ratio, and diversion rate arrays"""
-
-        qtyshr_array: NDArray[np.floating]
-        r_val: NDArray[np.floating]
-        divr_array: NDArray[np.floating]
-
-    _qtyshr_array = _gen_market_shares_uniform(_ssz).mktshr_array[:, :2]
-
-    _qtyshr_min = _qtyshr_array.min(axis=1, keepdims=True, initial=None)
-    _cprob = np.divide(_r_bar, 1 - (1 - _r_bar) * _qtyshr_min)
-    del _qtyshr_min
-
-    _r_val = np.divide(_cprob * (1 - _qtyshr_array), 1 - _cprob * _qtyshr_array)
-
-    _divr_array = _r_val * _qtyshr_array[:, ::-1] / (1 - _qtyshr_array)
-
-    return RvalSSP(_qtyshr_array, _r_val, _divr_array)
-
-
-def gen_implied_mkt_shr_1(_fcount: int = 5) -> NDArray[np.floating]:
-    """
-    Generate implied market shares for firm 1 with all pairs of products
-    in a putative market of `_fcount` firms.
-
-    Parameters
-    ----------
-    _fcount: Firm count for the given market
-
-    Returns
-    -------
-    Vector of implied market shares
-
-    """
-
-    _mkt_sample = _gen_market_shares_dirichlet([1] * _fcount, 10, RECConstants.OUTIN)
-
-    _mktshr_array = _mkt_sample.mktshr_array
-    _chprob_outside_good = _mkt_sample.choice_prob_outgd
-
-    _recapture_array = np.divide(
-        (1 - _chprob_outside_good) * (1 - _mktshr_array),
-        1 - (1 - _chprob_outside_good) * _mktshr_array,
-    )
-
-    _divratio_1j = np.divide(
-        _recapture_array[:, :1] * _mktshr_array[:, 1:], 1 - _mktshr_array[:, :1]
-    )
-    _divratio_j1 = np.divide(
-        _recapture_array[:, 1:] * _mktshr_array[:, :1], 1 - _mktshr_array[:, 1:]
-    )
-
-    _implied_mkt_shr_1 = np.divide(
-        _divratio_j1 * (_recapture_array[:, :1] - _divratio_1j),
-        _recapture_array[:, :1] * _recapture_array[:, 1:] - _divratio_1j * _divratio_j1,
-    )
-
-    return _implied_mkt_shr_1
-
-
-if __name__ == "__main__":
-    print(gen_implied_mkt_shr_1(10))

mergeron/examples/guidelines_enforcement_patterns.py

@@ -1,74 +0,0 @@
-from __future__ import annotations
-
-import numpy as np
-
-import mergeron.core.proportions_tests as pcl
-
-enf_counts = np.array(
-    (
-        (24, 25, 33, 33),
-        (30, 33, 27, 32),
-        (14, 15, 9, 11),
-        (13, 22, 6, 11),
-        (6, 8, 7, 9),
-        (1, 3, 2, 4),
-        (1, 3, 2, 4),
-        (1e-15, 1e-15, 1, 1),
-    ),
-    np.int16,
-)
-
-if __name__ == "__main__":
-    adj_n = len(enf_counts)
-
-    print()
-    print("C.I.s for relative frequency of II Requests")
-    cis_1996_2003, cis_2004_2011 = (
-        pcl.propn_ci_multinomial(_c) for _c in (enf_counts[:, _i] for _i in (1, 3))
-    )
-    for cdx, ci in enumerate(cis_1996_2003):
-        print(ci, ";", cis_2004_2011[cdx, :])
-
-    print()
-    print(
-        "C.I.s for relative frequency of II Requests, 1996-2003",
-        "(Bonnferoni-adjusted Wilson C.I.s)",
-    )
-    for count in (ic_1996_2003 := enf_counts[:, 1]):
-        print(
-            pcl.propn_ci(
-                count, ic_1996_2003.sum(), alpha=0.05 / adj_n, method="Wilson"
-            )[2:]
-        )
-
-    print("Conf. intervals for relative frequency of II Requests (Goodman, 1965)")
-    for goodman_alternative in "default", "simplified":
-        for goodman_method in "goodman", "quesenberry-hurst":
-            print(f"Method, {goodman_method!r}; alternative, {goodman_alternative!r}")
-            for _cis in np.column_stack([
-                pcl.propn_ci_multinomial(
-                    enf_counts[:, _cidx],
-                    method=goodman_method,
-                    alternative=goodman_alternative,
-                )
-                for _cidx in (1, 3)
-            ]):
-                print(_cis)
-            print()
-
-    print()
-    print("Conf. intervals for differences in proportions enforced (Goodman, 1964)")
-    print(repr(pcl.propn_diff_ci_multinomial(enf_counts[:, [1, 3]])))
-
-    print()
-    print("Goodman's chi-squared test for homogeneity of enforcement patterns")
-    print(repr(pcl.propn_test_multinomial(enf_counts[:, [1, 3]])))
-
-    print()
-    print(
-        "C.I.s for differences in proportions enforced",
-        "(Bonnferoni-adjusted Newcombe C.I.s)",
-    )
-    for counts in enf_counts:
-        print(pcl.propn_diff_ci(*counts, alpha=0.05 / adj_n, method="Newcombe"))
-    print()

mergeron/examples/investigations_stats_obs_tables.py

@@ -1,481 +0,0 @@
-"""
-
-Construct observed clearance rates and enforcement rates, by specification,
-from FTC merger investigations data
-
-"""
-
-import sys
-import warnings
-from collections.abc import Mapping, Sequence
-
-from numpy import einsum, row_stack, unique
-
-import mergeron.core.ftc_merger_investigations_data as fid
-import mergeron.gen.investigations_stats as isl
-from mergeron import DATA_DIR
-from mergeron.gen import INVResolution
-
-if not sys.warnoptions:
-    warnings.simplefilter("ignore")
-
-INVRES_RATIO_FORMAT_STR = "{: >3.0f}/{:<3.0f}"
-INVDATA_DOTTEX_FORMAT_STR = "{}.tex".format(
-    "_".join(("FTCMergerInvestigationsDataTables", "{}", "OBS"))
-)
-
-
-def invres_stats_odds_ratio_byhhianddelta(
-    _data_array_dict: Mapping[str, Mapping[str, Mapping[str, fid.INVData]]],
-    _data_periods: tuple[str, str],
-    _merger_classes: Sequence[isl.INDGRPConstants | isl.EVIDENConstants],
-    /,
-) -> tuple[str, ...]:
-    """
-    Reconstruct tables by HHI and Delta.
-
-    Source tables as well as tables from constructed periods.
-    """
-    if not all(_dpd in _data_array_dict for _dpd in _data_periods):
-        raise ValueError(
-            f"All given data periods, {_data_periods!r} must be contained "
-            f"in {tuple(_data_array_dict.keys())!r}"
-        )
-
-    print("Odds ratios by HHI and Delta:")
-    _stats_group = isl.StatsGrpSelector.HD
-    _invres_rate_table_content = isl.StatsContainer()
-    _invres_rate_table_design = isl.latex_jinja_env.get_template(
-        "ftcinvdata_byhhianddelta_table_template.tex.jinja2"
-    )
-    _invres_rate_table_content.obs_summary_type = f"{_stats_group}"
-
-    _output_dottex_pathlist: tuple[str, ...] = ()
-    for _merger_class in _merger_classes:
-        _table_ind_group = (
-            _merger_class
-            if isinstance(_merger_class, isl.INDGRPConstants)
-            else isl.INDGRPConstants.ALL
-        )
-        _table_evid_cond = (
-            _merger_class
-            if isinstance(_merger_class, isl.EVIDENConstants)
-            else isl.EVIDENConstants.UR
-        )
-        _invres_rate_table_content.obs_merger_class = f"{_merger_class}"
-
-        for _data_period in _data_periods:
-            _data_array_dict_sub = _data_array_dict[_data_period][f"{_stats_group}"]
-            _table_no = isl.table_no_lku(
-                _data_array_dict_sub,  # type: ignore
-                _table_ind_group,
-                _table_evid_cond,
-            )
-
-            _invres_rate_table_content.table_ref = _table_no
-
-            _data_period_0, _data_period_1 = (int(f) for f in _data_period.split("-"))
-            if _data_period_0 != 1996:
-                _invres_rate_table_content.invdata_notestr = " ".join((
-                    "NOTES:",
-                    isl.LTX_ARRAY_LINEEND,
-                    R"\(\cdot\) Data for period, {}".format(
-                        _data_periods[1].replace("-", "--")
-                    ),
-                    f"calculated by subtracting reported figures for 1996--{_data_period_0 - 1}",
-                    R"from reported figures for 1996--2011",
-                    isl.LTX_ARRAY_LINEEND,
-                    isl.LTX_ARRAY_LINEEND,
-                ))
-
-                _invres_rate_table_content.invdata_sourcestr = " ".join((
-                    "\\(\\cdot\\) Fed. Trade Comm'n ({}), at~\\cref{{fn:{}}},".format(
-                        _data_period_0, f"FTCInvData1996to{_data_period_0}"
-                    ),
-                    isl.LTX_ARRAY_LINEEND,
-                ))
-                _invres_rate_table_content.invdata_sourcestr += " ".join((
-                    "\\(\\cdot\\) Fed. Trade Comm'n ({}), at~\\cref{{fn:{}}},".format(
-                        _data_period_1, f"FTCInvData1996to{_data_period_1}"
-                    ),
-                    isl.LTX_ARRAY_LINEEND,
-                ))
-            else:
-                _invres_rate_table_content.invdata_sourcestr = " ".join((
-                    "\\(\\cdot\\) Fed. Trade Comm'n ({}), at~\\cref{{fn:{}}},".format(
-                        _data_period_1, f"FTCInvData1996to{_data_period_1}"
-                    ),
-                    isl.LTX_ARRAY_LINEEND,
-                ))
-
-            _invres_rate_table_content.obs_merger_class = f"{_merger_class}"
-            _invres_rate_table_content.obs_period = _data_period.split("-")
-
-            _invres_cnts_array = _data_array_dict_sub[_table_no][-1]  # type: ignore
-            _odds_ratio_data_str = ""
-            for _hhi_range_it in unique(_invres_cnts_array[:, 0]):  # type: ignore
-                _invres_cnts_row_for_hhi_range = _invres_cnts_array[
-                    _invres_cnts_array[:, 0] == _hhi_range_it
-                ][:, 2:]  # type: ignore
-                _odds_ratio_data_str += " & ".join([
-                    INVRES_RATIO_FORMAT_STR.format(*g)  # type: ignore
-                    for g in _invres_cnts_row_for_hhi_range
-                ])
-                _odds_ratio_data_str += " & {}".format(
-                    INVRES_RATIO_FORMAT_STR.format(
-                        *einsum("ij->j", _invres_cnts_row_for_hhi_range)
-                    )
-                )
-                _odds_ratio_data_str += isl.LTX_ARRAY_LINEEND
-
-            _invres_cnts_row_for_hhi_tots = row_stack([
-                einsum(
-                    "ij->j", _invres_cnts_array[_invres_cnts_array[:, 1] == f][:, 2:]
-                )
-                for f in unique(_invres_cnts_array[:, 1])
-            ])
-            _odds_ratio_data_str += " & ".join([
-                INVRES_RATIO_FORMAT_STR.format(*f)
-                for f in _invres_cnts_row_for_hhi_tots
-            ])
-            _odds_ratio_data_str += " & {}".format(
-                INVRES_RATIO_FORMAT_STR.format(
-                    *einsum("ij->j", _invres_cnts_row_for_hhi_tots)
-                )
-            )
-            _odds_ratio_data_str += isl.LTX_ARRAY_LINEEND
-            _invres_rate_table_content.invdata_byhhianddelta = _odds_ratio_data_str
-
-            _output_dottex_path = DATA_DIR / INVDATA_DOTTEX_FORMAT_STR.format(
-                f"{_stats_group}_{_data_period_1}_{_merger_class.replace(" ", "")}"
-            )
-            with _output_dottex_path.open(
-                "w", encoding="utf8"
-            ) as _invres_rate_table_dottex:
-                _invres_rate_table_dottex.write(
-                    _invres_rate_table_design.render(
-                        tmpl_data=_invres_rate_table_content
-                    )
-                )
-                print("\n", file=_invres_rate_table_dottex)
-
-            _output_dottex_pathlist += (_output_dottex_path.name,)
-            print(_odds_ratio_data_str)
-            print()
-            del _odds_ratio_data_str
-
-    return _output_dottex_pathlist
-
-
-def invres_stats_obs_setup(
-    _data_array_dict: Mapping,
-    _data_periods: tuple[str, str],
-    _merger_classes: Sequence[isl.INDGRPConstants | isl.EVIDENConstants],
-    _invres_spec: INVResolution = INVResolution.CLRN,
-    /,
-) -> tuple[str, ...]:
-    _notes_str_base = " ".join((
-        "NOTES:",
-        isl.LTX_ARRAY_LINEEND,
-        Rf"\(\cdot\) Data for period, {_data_periods[1].replace("-", "--")}",
-        "calculated by subtracting reported figures for 1996--{}".format(
-            int(_data_periods[1].split("-")[0]) - 1
-        ),
-        "from reported figures for 1996--{}".format(_data_periods[1].split("-")[1]),
-    ))
-    _stats_group_dict = {
-        isl.StatsGrpSelector.FC: {
-            "desc": f"{_invres_spec.capitalize()} rates by Firm Count",
-            "title_str": "By Number of Significant Competitors",
-            "hval": "Firm Count",
-            "hcol_width": 54,
-            "notewidth": 0.67,
-            "notestr": " ".join((
-                _notes_str_base,
-                isl.LTX_ARRAY_LINEEND,
-                isl.LTX_ARRAY_LINEEND,
-                isl.LTX_ARRAY_LINEEND,
-            )),
-        },
-        isl.StatsGrpSelector.DL: {
-            "desc": Rf"{_invres_spec.capitalize()} rates by range of \(\Delta HHI\)",
-            "title_str": R"By Change in Concentration (\Deltah{})",
-            "hval": R"$\Delta HHI$",
-            "hval_plus": R"{ $[\Delta_L, \Delta_H)$ pts.}",
-            "hcol_width": 54,
-            "notewidth": 0.67,
-            "notestr": " ".join((
-                _notes_str_base,
-                isl.LTX_ARRAY_LINEEND,
-                isl.LTX_ARRAY_LINEEND,
-            )),
-            "notestr_plus": " ".join((
-                R"\(\cdot\) Ranges of $\Delta HHI$ are defined as",
-                "half-open intervals with",
-                R"$\Delta_L \leqslant \Delta HHI < \Delta_H$, except that",
-                R"$2500 \text{ pts.} \leqslant \Delta HHI \leqslant 5000 \text{ pts.}$",
-                R"in the closed interval [2500, 5000] pts.",
-                isl.LTX_ARRAY_LINEEND,
-                isl.LTX_ARRAY_LINEEND,
-            )),
-        },
-        isl.StatsGrpSelector.ZN: {
-            "desc": f"{_invres_spec.capitalize()} rates by Approximate Presumption Zone",
-            "title_str": R"By Approximate \textit{2010 Guidelines} Concentration-Based Standards",
-            "hval": "Approximate Standard",
-            "hcol_width": 190,
-            "notewidth": 0.96,
-            "notestr": " ".join((
-                _notes_str_base,
-                isl.LTX_ARRAY_LINEEND,
-                isl.LTX_ARRAY_LINEEND,
-            )),
-        },
-    }
-
-    _output_dottex_pathlist = ()
-    for _stats_group_key in _stats_group_dict:
-        _output_dottex_path = _invres_stats_obs_render(
-            _data_array_dict,
-            _data_periods,
-            _merger_classes,
-            _stats_group_key,
-            _stats_group_dict[_stats_group_key],
-            _invres_spec,
-        )
-        _output_dottex_pathlist += (_output_dottex_path,)  # type: ignore
-
-    return _output_dottex_pathlist
-
-
-def _invres_stats_obs_render(
-    _data_array_dict: Mapping,
-    _data_periods: tuple[str, str],
-    _merger_classes: Sequence[isl.INDGRPConstants | isl.EVIDENConstants],
-    _stats_group: isl.StatsGrpSelector,
-    _stats_group_dict: Mapping,
-    _invres_spec: isl.INVResolution = INVResolution.CLRN,
-    /,
-) -> str:
-    _invres_rate_table_content = isl.StatsContainer()
-    _invres_rate_table_design = isl.latex_jinja_env.get_template(
-        "ftcinvdata_summarypaired_table_template.tex.jinja2"
-    )
-
-    print(
-        f'{_stats_group_dict["desc"]}:', ", ".join([f'"{g}"' for g in _merger_classes])
-    )
-    _invres_rate_table_content.test_regime = _invres_spec.capitalize()
-    _invres_rate_table_content.obs_summary_type = f"{_stats_group}"
-    _invres_rate_table_content.obs_summary_type_title = _stats_group_dict.get(
-        "title_str"
-    )
-    _invres_rate_table_content.hdrcol_raw_width = f'{_stats_group_dict["hcol_width"]}pt'
-
-    _hs1 = _stats_group_dict["hval"]
-    _hs2 = _h if (_h := _stats_group_dict.get("hval_plus", "")) else _hs1
-    _invres_rate_table_content.invdata_hdrcoldescstr = (
-        isl.latex_hrdcoldesc_format_str.format(
-            "hdrcol_raw",
-            f'{_stats_group_dict["hcol_width"]}pt',
-            "hdrcoldesc_raw",
-            "center",
-            " ".join((
-                _hs1 if _hs1 != _hs2 else Rf"{{ \phantom{{{_hs1}}} }}",
-                isl.LTX_ARRAY_LINEEND,
-                _hs2,
-                isl.LTX_ARRAY_LINEEND,
-            )),
-        )
-    )
-    del _hs1, _hs2
-
-    _invres_rate_table_content.obs_merger_class_0 = f"{_merger_classes[0]}"
-    _invres_rate_table_content.obs_merger_class_1 = f"{_merger_classes[1]}"
-    _invres_rate_table_content.obs_periods_str = (
-        Rf"{" & ".join(_data_periods)} \\".replace("-", "--")
-    )
-
-    _invres_rate_table_content.invdata_notewidth = _stats_group_dict["notewidth"]
-    _invres_rate_table_content.invdata_notestr = _stats_group_dict["notestr"]
-    if _n2 := _stats_group_dict.get("notestr_plus", ""):
-        _invres_rate_table_content.invdata_notestr += _n2
-    del _n2
-
-    _invdata_sourcestr_format_str = "{} {}".format(
-        R"\(\cdot\) Fed. Trade Comm'n ({}), at note~\cref{{fn:{}}}",
-        isl.LTX_ARRAY_LINEEND,
-    )
-
-    _table_nos = get_table_nos(
-        _data_array_dict, _merger_classes, _stats_group, _data_periods[0]
-    )
-    _src_table_nos_str = (
-        f'{", ".join(_table_nos[:-1])} and {_table_nos[-1]}'
-        if _merger_classes[0] != isl.EVIDENConstants.NE
-        else "Table~{0}.1, Table~{1}.1, and Table~{1}.2".format(
-            *(4, 10) if _stats_group == "ByFirmCount" else (3, 9)
-        )
-    )
-    _invres_rate_table_content.invdata_sourcestr = _invdata_sourcestr_format_str.format(
-        "2003", "FTCInvData1996to2003"
-    )
-    _invres_rate_table_content.invdata_sourcestr += (
-        _invdata_sourcestr_format_str.format("2011", "FTCInvData1996to2011")
-    )
-
-    _invdata_hdr_list: list[str] = []
-    _invdata_dat_list: list[list[str]] = []
-    _invres_cnt_totals: list[str] = []
-    _sort_order = (
-        isl.SortSelector.UCH
-        if _stats_group == isl.StatsGrpSelector.FC
-        else isl.SortSelector.REV
-    )
-
-    for _merger_class in _merger_classes:
-        _table_ind_group = (
-            _merger_class
-            if isinstance(_merger_class, isl.INDGRPConstants)
-            else isl.INDGRPConstants.ALL
-        )
-        _table_evid_cond = (
-            _merger_class
-            if isinstance(_merger_class, isl.EVIDENConstants)
-            else isl.EVIDENConstants.UR
-        )
-        for _data_period in _data_periods:
-            _invres_cnt_totals += [
-                isl.invres_stats_output(
-                    _data_array_dict,
-                    _data_period,
-                    _table_ind_group,
-                    _table_evid_cond,
-                    _stats_group,
-                    _invres_spec,
-                    return_type_sel=isl.StatsReturnSelector.CNT,
-                    print_to_screen=False,
-                )[1][-1][0]
-            ]
-
-            _invdata_hdr_list_it, _invdata_dat_list_it = isl.invres_stats_output(
-                _data_array_dict,
-                _data_period,
-                _table_ind_group,
-                _table_evid_cond,
-                _stats_group,
-                _invres_spec,
-                return_type_sel=isl.StatsReturnSelector.RPT,
-                sort_order=_sort_order,
-                print_to_screen=False,
-            )
-            _invdata_hdr_list = _invdata_hdr_list_it
-            _invdata_dat_list = (
-                _invdata_dat_list_it[:]
-                if not _invdata_dat_list
-                else [
-                    _invdata_dat_list[_r][:] + _invdata_dat_list_it[_r][:]
-                    for _r in range(len(_invdata_dat_list))
-                ]
-            )
-
-    isl.stats_print_rows(_invdata_hdr_list, _invdata_dat_list)
-
-    _invdata_hdrstr = "".join([
-        f"{_invdata_hdr_list[g]} {isl.LTX_ARRAY_LINEEND}"
-        for g in range(len(_invdata_hdr_list))
-    ])
-
-    _invdata_datstr = "".join([
-        f"{" & ".join(_invdata_dat_list[g])} {isl.LTX_ARRAY_LINEEND}"
-        for g in range(len(_invdata_dat_list))
-    ])
-
-    (
-        _invres_rate_table_content.mkt_counts_str_class_0,
-        _invres_rate_table_content.mkt_counts_str_class_1,
-    ) = (
-        R"{} \\".format(" & ".join([f"Obs. = {f}" for f in g]))
-        for g in [
-            _invres_cnt_totals[: len(_data_periods)],
-            _invres_cnt_totals[len(_data_periods) :],
-        ]
-    )
-
-    _invres_rate_table_content.invdata_numrows = len(_invdata_hdr_list)
-    _invres_rate_table_content.invdata_hdrstr = _invdata_hdrstr
-    _invres_rate_table_content.invdata_datstr = _invdata_datstr
-
-    _output_dottex_path = DATA_DIR / INVDATA_DOTTEX_FORMAT_STR.format(_stats_group)
-    with _output_dottex_path.open("w", encoding="UTF-8") as _output_dottex_file:
-        _output_dottex_file.write(
-            _invres_rate_table_design.render(tmpl_data=_invres_rate_table_content)
-        )
-        print("\n", file=_output_dottex_file)
-    del _invdata_hdrstr, _invdata_datstr
-
-    return _output_dottex_path.name
-
-
-def get_table_nos(
-    _data_array_dict: Mapping[str, fid.INVData],
-    _merger_classes: Sequence[isl.INDGRPConstants | isl.EVIDENConstants],
-    _stats_group: isl.StatsGrpSelector,
-    _data_period: str,
-    /,
-) -> list[str]:
-    _stats_group_major = (
-        "ByFirmCount" if _stats_group == isl.StatsGrpSelector.FC else "ByHHIandDelta"
-    )
-
-    _table_ind_groups = tuple(
-        (_m if isinstance(_m, isl.INDGRPConstants) else isl.INDGRPConstants.ALL)
-        for _m in _merger_classes
-    )
-    _table_evid_conds = tuple(
-        (_m if isinstance(_m, isl.EVIDENConstants) else isl.EVIDENConstants.UR)
-        for _m in _merger_classes
-    )
-
-    return list(
-        dict.fromkeys(
-            isl.table_no_lku(
-                _data_array_dict[_data_period][_stats_group_major],
-                _table_ind_group,
-                _table_evid_cond,
-            )
-            for _table_ind_group in _table_ind_groups
-            for _table_evid_cond in _table_evid_conds
-        )
-    )
-
-
-if __name__ == "__main__":
-    invdata_array_dict = fid.construct_data(
-        fid.INVDATA_ARCHIVE_PATH,
-        flag_backward_compatibility=False,
-        flag_pharma_for_exclusion=True,
-    )
-
-    merger_classes = (
-        isl.EVIDENConstants.NE,
-        isl.EVIDENConstants.ED,
-    )  # clstl.INDGRPConstants.IID)
-    data_periods = ("1996-2003", "2004-2011")
-    test_regime = INVResolution.ENFT
-
-    # Now generate the various tables summarizing merger investigations data
-    invres_stats_byhhianddelta_pathlist = invres_stats_odds_ratio_byhhianddelta(
-        invdata_array_dict, data_periods, merger_classes
-    )
-    print(invres_stats_byhhianddelta_pathlist)
-    invres_stats_allothers_pathlist = invres_stats_obs_setup(
-        invdata_array_dict, data_periods, merger_classes, test_regime
-    )
-    isl.render_table_pdf(
-        invres_stats_byhhianddelta_pathlist,
-        INVDATA_DOTTEX_FORMAT_STR.format("ByHHIandDelta"),
-    )
-    isl.render_table_pdf(
-        invres_stats_allothers_pathlist, INVDATA_DOTTEX_FORMAT_STR.format("AllOthers")
-    )