def METHOD_NAME(cls):
    formula_outcome = 'bweight ~ prenatal1_ + mmarried_ + mage + fbaby_'
    mod = OLS.from_formula(formula_outcome, dta_cat)
    tind = np.asarray(dta_cat['mbsmoke_'])
    cls.teff = TreatmentEffect(mod, tind, results_select=res_probit)