在一张图中绘制不同回归的边际效应

Question

我正在运行相同的probit回归五次，但使用的是不同的数据集。

我如何在同一张图上绘制每个回归的边际效应曲线，以一种简单的方式比较它们？

下面是我尝试过的代码：

probit one_y_unemp year* MR* AG* if Qualifica2 == 1
margins, dydx(year*) saving(me_intern, replace)

probit one_y_unemp year* MR* AG* if Qualifica3 == 1
margins, dydx(year*) saving(me_seniormanager, replace)

probit one_y_unemp year* MR* AG* if Qualifica4 == 1
margins, dydx(year*) saving(me_whitecollar, replace)

probit one_y_unemp year* MR* AG* if Qualifica5 == 1
margins, dydx(year*) saving(me_bluecollar, replace)

probit one_y_unemp year* MR* AG* if Qualifica6 == 1
margins, dydx(year*) saving(me_juniormanager, replace)

combomarginsplot me_intern me_seniormanager me_whitecollar me_bluecollar me_juniormanager, ///
labels("Intern" "Manager" "White Collar" "Blue Collar" "Junior Manager") xtitle("Years")

此代码的问题在于，我获得的图表没有为每个专业角色绘制一条曲线，而在x轴上具有年份虚拟对象，而是每年有一条虚拟线。