Stacked bar graphs of margins

[Pina Valle]

Hello,

Is it possible for Stata to generate a stacked bar graph of predictive margins after mlogit? Example below:

use https://stats.idre.ucla.edu/stat/data/hsbdemo, clear

qui: mlogit prog i.ses i.female
margins ses, predict(outcome(1))

Predictive margins                              Number of obs     =        200
Model VCE    : OIM

Expression   : Pr(prog==general), predict(outcome(1))

------------------------------------------------------------------------------
             |            Delta-method
             |     Margin   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
         ses |
        low  |   .3459232   .0704257     4.91   0.000     .2078913     .483955
     middle  |   .2093635   .0416951     5.02   0.000     .1276426    .2910843
       high  |     .15396   .0472522     3.26   0.001     .0613474    .2465726
------------------------------------------------------------------------------

qui: mlogit prog i.ses i.female
margins ses, predict(outcome(2))

Predictive margins                              Number of obs     =        200
Model VCE    : OIM

Expression   : Pr(prog==academic), predict(outcome(2))

------------------------------------------------------------------------------
             |            Delta-method
             |     Margin   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
         ses |
        low  |   .3985931    .071842     5.55   0.000     .2577855    .5394008
     middle  |   .4646755   .0512312     9.07   0.000     .3642643    .5650867
       high  |   .7256886   .0585014    12.40   0.000      .611028    .8403492
------------------------------------------------------------------------------

qui: mlogit prog i.ses i.female
margins ses, predict(outcome(3))

Predictive margins                              Number of obs     =        200
Model VCE    : OIM

Expression   : Pr(prog==vocation), predict(outcome(3))

------------------------------------------------------------------------------
             |            Delta-method
             |     Margin   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
         ses |
        low  |   .2554837   .0642489     3.98   0.000     .1295582    .3814092
     middle  |    .325961   .0481665     6.77   0.000     .2315564    .4203656
       high  |   .1203514   .0427004     2.82   0.005     .0366601    .2040427
------------------------------------------------------------------------------

After getting the results, I manually created the figure below using Excel but I'm wondering if Stata can produce something similar. Thank you!