set scheme sj
use inteff3_example

sjlog using inteff1, replace
probit particip female child uni fem_child fem_uni child_uni fem_chi_uni age age_sq german year2 year3 year4 year5 year6 year7 state1 state2 state3 state4 state6 state7 state8 state9 state10 state11 state12 state13 state14 state15 state16
inteff3
sjlog close, replace

sjlog using inteff2, replace
reg particip female child uni fem_child fem_uni child_uni fem_chi_uni age age_sq german year2 year3 year4 year5 year6 year7 state1 state2 state3 state4 state6 state7 state8 state9 state10 state11 state12 state13 state14 state15 state16
sjlog close, replace

qui probit particip female child uni fem_child fem_uni child_uni fem_chi_uni age age_sq german year2 year3 year4 year5 year6 year7 state1 state2 state3 state4 state6 state7 state8 state9 state10 state11 state12 state13 state14 state15 state16
sjlog using inteff3, replace
inteff3, at(0.5225 0      0     39.2785 1636.299 0.085 0.846 0.8545 0.8595 0.867 0.8305 0.855 0.97 0.9585 0.989 0.9115 0.7915 0.9405 0.9395 0.8825 0.8445 0.9475 0.973 0.9615 0.971 0.9725 0.942 1)
sjlog close, replace

sjlog using inteff4, replace
dprobit particip female child uni fem_child fem_uni child_uni fem_chi_uni age age_sq german year2 year3 year4 year5 year6 year7 state1 state2 state3 state4 state6 state7 state8 state9 state10 state11 state12 state13 state14 state15 state16
sjlog close, replace

qui probit particip female child uni fem_child fem_uni child_uni fem_chi_uni age age_sq german year2 year3 year4 year5 year6 year7 state1 state2 state3 state4 state6 state7 state8 state9 state10 state11 state12 state13 state14 state15 state16

sjlog using inteff5, replace
inteff3, average pex1(pe1) pex1x2x3(pe123) sx1(se1) sx1x2x3(se123)
sjlog close, replace

histogram pe1, saving(pe1, replace)
histogram pe123, saving(pe123, replace)