def METHOD_NAME(model_data: AttrDict) -> None:
    res = model_data.res
    stata_stats = model_data.stata_result.stats
    for i, key in enumerate(res.equations):
        eq = res.equations[key]
        stat = eq.f_statistic.stat
        stata_stat = stata_stats.loc[f"F_{i + 1}"].squeeze()
        if np.isnan(stata_stat):
            stata_stat = stata_stats.loc[f"chi2_{i + 1}"].squeeze()
        assert_allclose(stat, stata_stat)
        pval = eq.f_statistic.pval
        stata_pval = stata_stats.loc[f"p_{i + 1}"]
        assert_allclose(pval, stata_pval, atol=1e-6)