Dear Stata users,
I try to improve my study which use probit and marginal partial effect in education. However, the command "estsimp" does not work longer. I tried to see if I need to update it or which other command replaced it. I use Stata13. I compute the probability to get high score in 2 regions. I write below my command hoping someone can help me to replace the command or the lines or code
I also attach the figure I created almost 2 years ago when estsimp worked. My variable is High_Score and others are dummies
I hope some of you can help me.
Thank you very much
With kind regards,
Simona
estsimp logit High_Score Gender Attend_ISCED0 FamSouth SouthItaly HISCED Lang_Home FAMSTRUC IMMIG CULTPOSS HOMSCH, vce (robust)
gen p_mean=.
gen p_low=.
gen p_high=.
forvalues x=1(1)6 {
setx HISCED `x' Gender 1 FamSouth 0 SouthItaly 0 Attend_ISCED0 1 Lang_Home 1 FAMSTRUC 1 IMMIG 1
simqi, genpr(High_Score0 High_Score1)
summarize High_Score1, meanonly
replace p_mean=r(mean) in `x'
centile High_Score1, c(2.5 97.5)
replace p_low=r(c_1) in `x'
replace p_high=r(c_2) in `x'
drop High_Score0 High_Score1
}
gen x_axis=_n if _n>0&_n<6
graph twoway (line p_mean x_axis) (line p_low x_axis, lpattern(dash)) (line p_high x_axis, lpattern(dash)), yscale(range(0 1.0)) xlabel(0(1)6) ylabel(0(.1) 1.0) ytitle("Pr(High_Score)")xtitle("HISCED") legend(order(1 2) label(1 "Pr(High_Score)") label(2 "95% CI"))graph1.pdf
I try to improve my study which use probit and marginal partial effect in education. However, the command "estsimp" does not work longer. I tried to see if I need to update it or which other command replaced it. I use Stata13. I compute the probability to get high score in 2 regions. I write below my command hoping someone can help me to replace the command or the lines or code
I also attach the figure I created almost 2 years ago when estsimp worked. My variable is High_Score and others are dummies
I hope some of you can help me.
Thank you very much
With kind regards,
Simona
estsimp logit High_Score Gender Attend_ISCED0 FamSouth SouthItaly HISCED Lang_Home FAMSTRUC IMMIG CULTPOSS HOMSCH, vce (robust)
gen p_mean=.
gen p_low=.
gen p_high=.
forvalues x=1(1)6 {
setx HISCED `x' Gender 1 FamSouth 0 SouthItaly 0 Attend_ISCED0 1 Lang_Home 1 FAMSTRUC 1 IMMIG 1
simqi, genpr(High_Score0 High_Score1)
summarize High_Score1, meanonly
replace p_mean=r(mean) in `x'
centile High_Score1, c(2.5 97.5)
replace p_low=r(c_1) in `x'
replace p_high=r(c_2) in `x'
drop High_Score0 High_Score1
}
gen x_axis=_n if _n>0&_n<6
graph twoway (line p_mean x_axis) (line p_low x_axis, lpattern(dash)) (line p_high x_axis, lpattern(dash)), yscale(range(0 1.0)) xlabel(0(1)6) ylabel(0(.1) 1.0) ytitle("Pr(High_Score)")xtitle("HISCED") legend(order(1 2) label(1 "Pr(High_Score)") label(2 "95% CI"))graph1.pdf
