·foreach var of varlist DD_chn DD_f_chn poll_emp_so2_chn poll_emp_pt_chn poll_emp_no2_chn poll_emp_pca_chn {
        sort prefect_adj age_group year
        bysort prefect_adj age_group: gen f_`var'=`var'[_n+1]
        }
·ivreg2 d_All_mortality (f_shock_pca_chn=f_shock_pca_f_chn) i.region_year if age_group==1 & year>=1990 & year<=2000 [aweight=weight], endog(f_shock_pca_chn) first ffirst savefirst partial(i.region_year) cluster(province)
·matrix stat=e(first)
·local temp1=stat[15,1]
·outreg2 using "$dir/result/result_age0_future_f.tex", keep(f_shock_pca_chn f_DD_chn) dec(3) addstat("Fstat PollExpShock", `temp1') tex replace
·ivreg2 d_All_mortality (f_DD_chn=f_DD_f_chn) i.region_year if age_group==1 & year>=1990 & year<=2000 [aweight=weight], endog(f_DD_chn) first ffirst savefirst partial(i.region_year) cluster(province)
·matrix stat=e(first)
·local temp1=stat[15,1]
·outreg2 using "$dir/result/result_age0_future_f.tex", keep(f_shock_pca_chn f_DD_chn) dec(3) addstat("Fstat ExpShock", `temp1') tex append
·ivreg2 d_All_mortality (f_shock_pca_chn f_DD_chn=f_shock_pca_f_chn f_DD_f_chn) i.region_year if age_group==1 & year>=1990 & year<=2000 [aweight=weight], endog(f_shock_pca_chn f_DD_chn) first ffirst savefirst partial(i.region_year) cluster(province)
·matrix stat=e(first)
·local temp1=stat[15,1]
·local temp2=stat[15,2]
·outreg2 using "$dir/result/result_age0_future_f.tex", keep(f_shock_pca_chn f_DD_chn) dec(3) addstat("Fstat PollExpShock", `temp1',"Fstat ExpShock", `temp2') tex append