sysuse nlsw88, clear
recode occupation (11/12=4) ///
                  (9/10=13)

xi: logit union i.occupation married never_married
sheafcoef, latent( occ: _I* ; marst: married never_married)
sheafcoef, latent( occ: _I* ; marst: married never_married) eform