clear set more off cd "$constructeddirectory\" use alldhsirclean, clear cd "$outdirectory\" keep if polyg~=. global outopts dec(3) label nocons excel ctitle(" ") cttop(" ") nonotes *Main controls global controls $giscontrols $indivcontrolsnorel ***Missions and schooling/polygamy cap gen br=. cap gen bf=. cap gen aet=. ***Regression foreach depvar in polyg { foreach rhsvar in catholic_lt20 protestant_lt20 { cap erase tableA5_`depvar'_`rhsvar'_controlsample.xml cap erase tableA5_`depvar'_`rhsvar'_controlsample.txt reg `depvar' `rhsvar' $controls if `rhsvar'_sample==1, cl(clusterid) replace br=_b[`rhsvar'] outreg2 `rhsvar' using tableA5_`depvar'_`rhsvar'_controlsample, $outopts /// addtext(Outcome, `depvar', Other cont., Geo./IndNoRelig., F.E., None, Clustering, Cluster) /// append areg `depvar' `rhsvar' if `rhsvar'_sample==1, a(countryroundid) cl(clusterid) replace bf=_b[`rhsvar'] replace aet=bf/(br-bf) scalar aet=aet[1] outreg2 `rhsvar' using tableA5_`depvar'_`rhsvar'_controlsample if `rhsvar'_sample==1, $outopts /// addtext(Outcome, `depvar', Other cont., None, F.E., CountryRound, Clustering, Cluster) /// addstat(AET,scalar(aet)) append areg `depvar' `rhsvar' $controls if `rhsvar'_sample==1, a(countryroundid) cl(clusterid) replace bf=_b[`rhsvar'] replace aet=bf/(br-bf) scalar aet=aet[1] outreg2 `rhsvar' using tableA5_`depvar'_`rhsvar'_controlsample, $outopts /// addtext(Outcome, `depvar', Other cont., Geo./IndNoRelig., F.E., CountryRound, Clustering, Cluster) /// addstat(AET,scalar(aet)) append areg `depvar' `rhsvar' $controls if `rhsvar'_sample==1, a(regionid) cl(clusterid) replace bf=_b[`rhsvar'] replace aet=bf/(br-bf) scalar aet=aet[1] outreg2 `rhsvar' using tableA5_`depvar'_`rhsvar'_controlsample, $outopts /// addtext(Outcome, `depvar', Other cont., Geo./IndNoRelig., F.E., Region, Clustering, Cluster) /// addstat(AET,scalar(aet)) append cap erase tableA5_`depvar'_`rhsvar'_controlsample.txt } }