Dear Statalisters,

I have the following easy issue to solve. I would like to plot margins of employ2, which has four levels and a missing category.

svy: mlogit uniontype2 i.durationtype i.agesm2_f##i.mage_gap_cat i.wave5 i.employ2 ib0.evunionpr i.shared i.educ2_f i.homog if marital==3&duratyrs3>0&duratyrs3<6&homo3==0&sampst ==1&ivfioall==1&(minagentry>18&minagentry<37)&sex! =.&educ2_f!=4&educ2_m!=4&agesm2_f!=., level(90) rrr base

margins i.employ2 , level(90) pr(out(1)) pr(out(2))

mplotoffset , xdim(employ2) legend( order(4 "Marriage" 3 "Dissolution" ) ring(5)) ytitle("Probability") title("Employment") offset(0.15) recast(scatter) plotregion(lcol(black)) saving(employ2, replace) xlab(,labsize(small) alternate) xtitle("") plot2opts(msymbol(D))

I would like to plot the four categories, but not the missing. one.

However, I would avoid restrcting with "if" in the following way,

margins i.employ2 if employ2<5 , level(90) pr(out(1)) pr(out(2))
as I have slight different estimates than if computed margins for all the categories.

Thank you and best,
Lydia