open "/Users/jcaceres/Google Drive/jcd/carlos3/2015_2016/EA_2016/Web_EconomiaAplicada/materiales/fertility.gdt"

summary weeksm1

genr employed=weeksm1>0

#LINEAR PROB MODEL
ols employed const morekids agem1 black hispan othrace --robust

#PROBIT
probit employed const morekids agem1 black hispan othrace --robust

genr bp=$coeff

#IMPACT OF FERTILITY
genr z0=bp[1]+bp[2]*0+bp[3]*mean(agem1)+bp[4]*mean(black)+bp[5]*mean(hispan)+bp[6]*mean(othrace)
genr z1=bp[1]+bp[2]*1+bp[3]*mean(agem1)+bp[4]*mean(black)+bp[5]*mean(hispan)+bp[6]*mean(othrace)
genr z_mean=bp[1]+bp[2]*mean(morekids)+bp[3]*mean(agem1)+bp[4]*mean(black)+bp[5]*mean(hispan)+bp[6]*mean(othrace)

genr impact_fertility=cdf(N,z1)-cdf(N,z0)

#IMPACT OF AGE

genr impact_age=pdf(N,z_mean)*bp[3]

#LOGIT
logit employed const morekids agem1 black hispan othrace --robust

genr bl=$coeff

genr w0=bl[1]+bl[2]*0+bl[3]*mean(agem1)+bl[4]*mean(black)+bl[5]*mean(hispan)+bl[6]*mean(othrace)
genr w1=bl[1]+bl[2]*1+bl[3]*mean(agem1)+bl[4]*mean(black)+bl[5]*mean(hispan)+bl[6]*mean(othrace)
genr w_mean=bl[1]+bl[2]*mean(morekids)+bl[3]*mean(agem1)+bl[4]*mean(black)+bl[5]*mean(hispan)+bl[6]*mean(othrace)

#IMPACT OF FERTILITY

genr impact_fertilityl=(1/(1+exp(-w1)))-(1/(1+exp(-w0)))

#IMPACT OF AGE

genr impact_agel=(exp(-w_mean)/(1+exp(-w_mean))^2)*bl[3]