sysuse nlsw88, clear recode occupation (11/12=4) /// (9/10=13) gen ln_w = ln(wage) qui xi: reg ln_w i.occupation married never_married sheafcoef, latent(occ:_I* ; marst:married never_married) sheafcoef, latent(occ:_I* ; marst:married never_married) beta