set more off set memory 2g set matsize 800 version 9 cd c:\shin\project\nhi_health\JDE_revision\2nd_revision\data_program capture log close log using Table6&8.log, replace capture program drop action program action args tech samp2 dep cd c:\shin\project\nhi_health\JDE_revision\2nd_revision\data_program\Table6_8 use `samp2', clear gen Q=`tech' gen NHIPRI_Q=NHIPRI_`tech' gen NHIFRM_Q=NHIFRM_`tech' gen PRI_Q=PRI`tech' gen FRM_Q=FRM`tech' gen NHIQ=NHI`tech' summarize Q if PRI==1 & NHI==1 [aw=number] scalar am_pri=r(mean) summarize Q if FRM==1 & NHI==1 [aw=number] scalar am_frm=r(mean) summarize Q if (PRI==0 & FRM==0) & NHI==1 [aw=number] scalar am_gov=r(mean) summarize Q if PRI==1 & NHI==0 [aw=number] scalar bm_pri=r(mean) summarize Q if FRM==1 & NHI==0 [aw=number] scalar bm_frm=r(mean) summarize Q if (PRI==0 & FRM==0) & NHI==0 [aw=number] scalar bm_gov=r(mean) egen NHI_county=group(county1 NHI) xi i.yearn i.county1 cd c:\shin\project\nhi_health\JDE_revision\2nd_revision\data_program local outputfile "output\Table6&8_`tech'" local ind_vars = `"NHIPRI_Q NHIFRM_Q PRI_Q FRM_Q NHIQ NHI_PRI NHI_FRM PRI FRM NHI Q"' /***************************/ /* reg `dep' `ind_vars' _Iyearn_1991-_Iyearn_1993 _Iyearn_1996-_Iyearn_2001 _Icounty1_* [aw=number] lincom NHI_PRI+Q*((am_pri-bm_pri)-(am_gov-bm_gov))+PRI_Q*(am_pri-bm_pri)+NHIQ*(am_pri-am_gov)+NHIPRI_Q*am_pri scalar PRI_NHI_Q=r(estimate) scalar se_p=r(se) gen t_PRI_NHI_Q= PRI_NHI_Q/se_p sum t_PRI_NHI_Q lincom NHI_FRM+Q*((am_frm-bm_frm)-(am_gov-bm_gov))+FRM_Q*(am_frm-bm_frm)+NHIQ*(am_frm-am_gov)+NHIFRM_Q*am_frm scalar FRM_NHI_Q=r(estimate) scalar se_f=r(se) gen t_FRM_NHI_Q= FRM_NHI_Q/se_f sum t_FRM_NHI_Q outreg2 `ind_vars' using `outputfile',excel stat(coef se tstat) nocons bd(3) td(2) /* */ ctitle("`samp2'") /* */ addstat(PRI_NHI_Q, PRI_NHI_Q, std-p, se_p, t-stat-p, t_PRI_NHI_Q, FRM_NHI_Q, FRM_NHI_Q, std-f, se_f, t-stat-f, t_FRM_NHI_Q) /* */ addtext(Model, Grouped, Cluster, No) drop t_PRI_NHI_Q t_FRM_NHI_Q scalar drop PRI_NHI_Q FRM_NHI_Q se_p se_f */ /***************************/ reg `dep' `ind_vars' _Iyearn_1991-_Iyearn_1993 _Iyearn_1996-_Iyearn_2001 _Icounty1_* [aw=number], cluster(NHI_county) lincom NHI_PRI+Q*((am_pri-bm_pri)-(am_gov-bm_gov))+PRI_Q*(am_pri-bm_pri)+NHIQ*(am_pri-am_gov)+NHIPRI_Q*am_pri scalar PRI_NHI_Q=r(estimate) scalar se_p=r(se) gen t_PRI_NHI_Q= PRI_NHI_Q/se_p sum t_PRI_NHI_Q lincom NHI_FRM+Q*((am_frm-bm_frm)-(am_gov-bm_gov))+FRM_Q*(am_frm-bm_frm)+NHIQ*(am_frm-am_gov)+NHIFRM_Q*am_frm scalar FRM_NHI_Q=r(estimate) scalar se_f=r(se) gen t_FRM_NHI_Q= FRM_NHI_Q/se_f sum t_FRM_NHI_Q outreg2 `ind_vars' using `outputfile',excel stat(coef se tstat) nocons bd(3) td(2) /* */ ctitle("`samp2'") /* */ addstat(PRI_NHI_Q, PRI_NHI_Q, std-p, se_p, t-stat-p, t_PRI_NHI_Q, FRM_NHI_Q, FRM_NHI_Q, std-f, se_f, t-stat-f, t_FRM_NHI_Q) /* */ addtext(Model, Grouped, Cluster, County_NHI) drop t_PRI_NHI_Q t_FRM_NHI_Q scalar drop PRI_NHI_Q FRM_NHI_Q se_p se_f /***************************/ /* reg `dep' `ind_vars' _Iyearn_1991-_Iyearn_1993 _Iyearn_1996-_Iyearn_2001 _Icounty1_* [aw=number], cluster(ins_county) lincom NHI_PRI+Q*((am_pri-bm_pri)-(am_gov-bm_gov))+PRI_Q*(am_pri-bm_pri)+NHIQ*(am_pri-am_gov)+NHIPRI_Q*am_pri scalar PRI_NHI_Q=r(estimate) scalar se_p=r(se) gen t_PRI_NHI_Q= PRI_NHI_Q/se_p sum t_PRI_NHI_Q lincom NHI_FRM+Q*((am_frm-bm_frm)-(am_gov-bm_gov))+FRM_Q*(am_frm-bm_frm)+NHIQ*(am_frm-am_gov)+NHIFRM_Q*am_frm scalar FRM_NHI_Q=r(estimate) scalar se_f=r(se) gen t_FRM_NHI_Q= FRM_NHI_Q/se_f sum t_FRM_NHI_Q outreg2 `ind_vars' using `outputfile',excel stat(coef se tstat) nocons bd(3) td(2) /* */ ctitle("`samp2'") /* */ addstat(PRI_NHI_Q, PRI_NHI_Q, std-p, se_p, t-stat-p, t_PRI_NHI_Q, FRM_NHI_Q, FRM_NHI_Q, std-f, se_f, t-stat-f, t_FRM_NHI_Q) /* */ addtext(Model, Grouped, Cluster, Group_County) drop t_PRI_NHI_Q t_FRM_NHI_Q scalar drop PRI_NHI_Q FRM_NHI_Q se_p se_f /***************************/ reg `dep' `ind_vars' _Iyearn_1991-_Iyearn_1993 _Iyearn_1996-_Iyearn_2001 _Icounty1_* [aw=number], cluster(county1) lincom NHI_PRI+Q*((am_pri-bm_pri)-(am_gov-bm_gov))+PRI_Q*(am_pri-bm_pri)+NHIQ*(am_pri-am_gov)+NHIPRI_Q*am_pri scalar PRI_NHI_Q=r(estimate) scalar se_p=r(se) gen t_PRI_NHI_Q= PRI_NHI_Q/se_p sum t_PRI_NHI_Q lincom NHI_FRM+Q*((am_frm-bm_frm)-(am_gov-bm_gov))+FRM_Q*(am_frm-bm_frm)+NHIQ*(am_frm-am_gov)+NHIFRM_Q*am_frm scalar FRM_NHI_Q=r(estimate) scalar se_f=r(se) gen t_FRM_NHI_Q= FRM_NHI_Q/se_f sum t_FRM_NHI_Q outreg2 `ind_vars' using `outputfile',excel stat(coef se tstat) nocons bd(3) td(2) /* */ ctitle("`samp2'") /* */ addstat(PRI_NHI_Q, PRI_NHI_Q, std-p, se_p, t-stat-p, t_PRI_NHI_Q, FRM_NHI_Q, FRM_NHI_Q, std-f, se_f, t-stat-f, t_FRM_NHI_Q) /* */ addtext(Model, Grouped, Cluster, County) drop t_PRI_NHI_Q t_FRM_NHI_Q scalar drop PRI_NHI_Q FRM_NHI_Q se_p se_f */ /***************************/ xtset ins_county yearn newey2 `dep' `ind_vars' _Iyearn_1991-_Iyearn_1993 _Iyearn_1996-_Iyearn_2001 _Icounty1_* [aw=number], lag(2) force lincom NHI_PRI+Q*((am_pri-bm_pri)-(am_gov-bm_gov))+PRI_Q*(am_pri-bm_pri)+NHIQ*(am_pri-am_gov)+NHIPRI_Q*am_pri scalar PRI_NHI_Q=r(estimate) scalar se_p=r(se) gen t_PRI_NHI_Q= PRI_NHI_Q/se_p sum t_PRI_NHI_Q lincom NHI_FRM+Q*((am_frm-bm_frm)-(am_gov-bm_gov))+FRM_Q*(am_frm-bm_frm)+NHIQ*(am_frm-am_gov)+NHIFRM_Q*am_frm scalar FRM_NHI_Q=r(estimate) scalar se_f=r(se) gen t_FRM_NHI_Q= FRM_NHI_Q/se_f sum t_FRM_NHI_Q outreg2 `ind_vars' using `outputfile',excel stat(coef se tstat) nocons bd(3) td(2) /* */ ctitle("`samp2'") /* */ addstat(PRI_NHI_Q, PRI_NHI_Q, std-p, se_p, t-stat-p, t_PRI_NHI_Q, FRM_NHI_Q, FRM_NHI_Q, std-f, se_f, t-stat-f, t_FRM_NHI_Q) /* */ addtext(Model, Newey2) drop t_PRI_NHI_Q t_FRM_NHI_Q scalar drop PRI_NHI_Q FRM_NHI_Q se_p se_f drop Q NHIPRI_Q NHIFRM_Q PRI_Q FRM_Q NHIQ end foreach i in "p_pedia" "p_obst" "p_bbaby" "p_pedia_lag1" "p_obst_lag1" "p_bbaby_lag1"{ action `i' whole postdeath1 foreach j in 0 1{ action `i' hs_grad`j' postdeath1 } foreach j in 0 1{ action `i' fullterm`j' postdeath1 } foreach j in 1 0{ action `i' lbw`j' postdeath1 } } exit