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 $indivcontrols ***Missions and schooling/polygamy cap gen br=. cap gen bf=. cap gen aet=. foreach depvar in polyg { foreach rhsvar in lncatholic_dist lnprotestant_dist { cap erase tableA4_`depvar'_`rhsvar'_rdum.xml cap erase tableA4_`depvar'_`rhsvar'_rdum.txt reg `depvar' `rhsvar' $controls, cl(clusterid) replace br=_b[`rhsvar'] outreg2 `rhsvar' using tableA4_`depvar'_`rhsvar'_rdum, $outopts /// addtext(Outcome, `depvar', Other cont., Geo./Ind., F.E., None, Clustering, Cluster) /// append areg `depvar' `rhsvar', a(countryroundid) cl(clusterid) replace bf=_b[`rhsvar'] replace aet=bf/(br-bf) scalar aet=aet[1] outreg2 `rhsvar' using tableA4_`depvar'_`rhsvar'_rdum, $outopts /// addtext(Outcome, `depvar', Other cont., None, F.E., CountryRound, Clustering, Cluster) /// addstat(AET,scalar(aet)) append areg `depvar' `rhsvar' $controls, a(countryroundid) cl(clusterid) replace bf=_b[`rhsvar'] replace aet=bf/(br-bf) scalar aet=aet[1] outreg2 `rhsvar' using tableA4_`depvar'_`rhsvar'_rdum, $outopts /// addtext(Outcome, `depvar', Other cont., Geo./Ind., F.E., CountryRound, Clustering, Cluster) /// addstat(AET,scalar(aet)) append areg `depvar' `rhsvar' $controls, a(regionid) cl(clusterid) replace bf=_b[`rhsvar'] replace aet=bf/(br-bf) scalar aet=aet[1] outreg2 `rhsvar' using tableA4_`depvar'_`rhsvar'_rdum, $outopts /// addtext(Outcome, `depvar', Other cont., Geo./Ind., F.E., Region, Clustering, Cluster) /// addstat(AET,scalar(aet)) append cap erase tableA4_`depvar'_`rhsvar'_rdum.txt } }