I am looking for a command that will display estimates of a two-part ML estimation. This model predicts the effects of variables age, female, family, educ_uni on parameters phi and sigma.

Code:
#delimit ;
            ml model lf eut_crra
                (phi:    choice 
                        p_ah x_ah p_al x_al 
                        p_bh x_bh p_bl x_bl 
                        = age female family educ_uni) 
                (sigma: = age female family educ_uni), 
                    cluster(pid) 
                    technique(nr) 
                    maximize
                    iterate(300)
                    difficult;
            #delimit cr
            estimates store est1
            
    * tabulate results
    estimates table est1 , b(%9.3f) stats(ll N N_clust) stfmt(%9.0f)
The results are:
Array

While I can call up the estimates of variables age ... _cons on phi as _b[age] ... _b[_cons],

Code:
di _b[age]
global _b_age = _b[age]
what would be the command to call up estimates of age ... _cons on sigma?

Thank you in advance,
Michael