Making marginal effects appear in a graph after marginsplot

[Luis Ortiz]

Dear members of the list,

I'm generating the marginal effects of gender after the following model

xtreg work_authority age gender_val2 i.work_dimensions2 female##children_16_simp##highest_educ3 i.cycle

Next, the marginal effects that I request:

margins, dydx(female) at(children_16_simp=(0 1 2 3) highest_educ3=(1 2 3))

In the following graph, generated with marginsplot (below), I would like the marginal effects (figures) to appear in the graph by each value in the interaction; that is, for any combination of 'children_16_simp' (categories of motherhood) and 'highest_educ3' (three levels of educational attainment).

marginsplot, xdimension(children_16_simp) plotregion(margin(large))

Do you know of any way of doing that?

Many thanks for your attention

And best wishes

Luis Ortiz