График предельного эффекта различных регрессий на одном графике

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 имеет одну линию на манекен каждого года.

Answer 1

Вместо этого вы можете использовать команду сообщества coefplot:

sysuse auto, clear
estimates clear

egen price2 = std(price)
egen rep782 = std(rep78)
replace rep782 = rep78 * -10

regress mpg i.foreign##c.price2 rep78
margins, dydx(foreign) at(price2=(-3(0.5)3)) post
estimates store m1

regress mpg i.foreign##c.rep782 price2
margins, dydx(foreign) at(rep782=(-3(0.5)3)) post
estimates store m2

coefplot m1 m2, at