***************** *** Table A.2 *** ***************** * To run these files, store the ado files cgmreg.ado and _makestar.ado in your folder that contains .ado files *cgmreg.ado" *_makestar.ado" *helppoor2 use "data\gss_estimation.dta", clear xi: cgmwildboot helppoor2 shock $controls if newsample == 0, cluster(region) bootcluster(region) reps($reps) seed(999) sum helppoor2 if e(N) scalar mean_outc = r(mean) matrix helppoor2 = e(bootresults) foreach var of varlist shock educ married female black unem{ scalar b_`var' = _b[`var'] scalar se_`var' = _se[`var'] } clear matsave helppoor2, p("$dir\data") dropall replace g nobs = e(N) g r2 = e(r2) g nclust = e(N_region) g mean_outc = mean_outc g b_shock = b_shock g b_educ = b_educ g b_married = b_married g b_female = b_female g b_black = b_black g b_unem = b_unem g se_shock = se_shock g se_educ = se_educ g se_married = se_married g se_female = se_female g se_black = se_black g se_unem = se_unem save, replace *welfare use "data\gss_estimation.dta", clear xi: cgmwildboot welfare shock $controls if newsample == 0, cluster(region) bootcluster(region) reps($reps) seed(999) sum welfare if e(N) scalar mean_outc = r(mean) matrix welfare = e(bootresults) foreach var of varlist shock educ married female black unem{ scalar b_`var' = _b[`var'] scalar se_`var' = _se[`var'] } clear matsave welfare, p("$dir\data") dropall replace g nobs = e(N) g r2 = e(r2) g nclust = e(N_region) g mean_outc = mean_outc g b_shock = b_shock g b_educ = b_educ g b_married = b_married g b_female = b_female g b_black = b_black g b_unem = b_unem g se_shock = se_shock g se_educ = se_educ g se_married = se_married g se_female = se_female g se_black = se_black g se_unem = se_unem save, replace *getahead use "data\gss_estimation.dta", clear xi: cgmwildboot getahead shock $controls if newsample == 0, cluster(region) bootcluster(region) reps($reps) seed(999) sum getahead if e(N) scalar mean_outc = r(mean) matrix getahead = e(bootresults) foreach var of varlist shock educ married female black unem{ scalar b_`var' = _b[`var'] scalar se_`var' = _se[`var'] } clear matsave getahead, p("$dir\data") dropall replace g nobs = e(N) g r2 = e(r2) g nclust = e(N_region) g mean_outc = mean_outc g b_shock = b_shock g b_educ = b_educ g b_married = b_married g b_female = b_female g b_black = b_black g b_unem = b_unem g se_shock = se_shock g se_educ = se_educ g se_married = se_married g se_female = se_female g se_black = se_black g se_unem = se_unem save, replace *partyidnew use "data\gss_estimation.dta", clear xi: cgmwildboot partyidnew shock $controls if newsample == 0, cluster(region) bootcluster(region) reps($reps) seed(999) sum partyidnew if e(N) scalar mean_outc = r(mean) matrix partyidnew = e(bootresults) foreach var of varlist shock educ married female black unem{ scalar b_`var' = _b[`var'] scalar se_`var' = _se[`var'] } clear matsave partyidnew, p("$dir\data") dropall replace g nobs = e(N) g r2 = e(r2) g nclust = e(N_region) g mean_outc = mean_outc g b_shock = b_shock g b_educ = b_educ g b_married = b_married g b_female = b_female g b_black = b_black g b_unem = b_unem g se_shock = se_shock g se_educ = se_educ g se_married = se_married g se_female = se_female g se_black = se_black g se_unem = se_unem save, replace *polviews2 use "data\gss_estimation.dta", clear xi: cgmwildboot polviews2 shock $controls if newsample == 0, cluster(region) bootcluster(region) reps($reps) seed(999) sum polviews2 if e(N) scalar mean_outc = r(mean) matrix polviews2 = e(bootresults) foreach var of varlist shock educ married female black unem{ scalar b_`var' = _b[`var'] scalar se_`var' = _se[`var'] } clear matsave polviews2, p("$dir\data") dropall replace g nobs = e(N) g r2 = e(r2) g nclust = e(N_region) g mean_outc = mean_outc g b_shock = b_shock g b_educ = b_educ g b_married = b_married g b_female = b_female g b_black = b_black g b_unem = b_unem g se_shock = se_shock g se_educ = se_educ g se_married = se_married g se_female = se_female g se_black = se_black g se_unem = se_unem save, replace capture log close *presidential2 use "data\gss_estimation.dta", clear xi: cgmwildboot presidential2 shock $controls if newsample == 0, cluster(region) bootcluster(region) reps($reps) seed(999) sum presidential2 if e(N) scalar mean_outc = r(mean) matrix presidential2 = e(bootresults) foreach var of varlist shock educ married female black unem{ scalar b_`var' = _b[`var'] scalar se_`var' = _se[`var'] } clear matsave presidential2, p("$dir\data") dropall replace g nobs = e(N) g r2 = e(r2) g nclust = e(N_region) g mean_outc = mean_outc g b_shock = b_shock g b_educ = b_educ g b_married = b_married g b_female = b_female g b_black = b_black g b_unem = b_unem g se_shock = se_shock g se_educ = se_educ g se_married = se_married g se_female = se_female g se_black = se_black g se_unem = se_unem save, replace capture log close * analytic p-values global outcomes "helppoor2 welfare getahead partyidnew polviews2 presidential2" global coeffs "shock educ married female black unem" local i = 1 foreach out of global outcomes{ use "data\gss_estimation.dta", clear reg `out' shock $controls if newsample == 0, cluster(region) local r`i' = string(e(r2), "%15.2fc") parmest, saving("data\analytic.dta", replace) keep if e(sample) count local n`i' = string(r(N), "%15.0fc") sum `out' local out`i' = string(r(mean), "%15.3fc") keep yearint duplicates drop count local waves`i' = string(r(N), "%15.0fc") use "data\analytic.dta", clear foreach item of global coeffs{ _makestar "`item'" sum estimate if parm == "`item'" local nspec = string(r(mean), "%15.3fc") local b_`item'`i' "`nspec'$star" sum p if parm == "`item'" local p1_`item'`i' = string(r(mean), "%15.3fc") sum stderr if parm == "`item'" local se_`item'`i' = string(r(mean), "%15.3fc") } local i = `i' + 1 } * bootstrap p-values global outcomes "helppoor2 welfare getahead partyidnew polviews2 presidential2" local i = 1 foreach out of global outcomes{ use data/`out'.dta, clear foreach var of varlist shock educ married female black unem{ keep if type == 1 g t_orig = `var' if iter == 1 egen t_orig_max = max(t_orig) drop t_orig ren t_orig_max t_orig g outside = `var' > t_orig if b_`var' > 0 replace outside = `var' < t_orig if b_`var' < 0 sum outside g pval = 2*r(mean) sum pval local p_`var'`i' = string(r(mean), "%15.3fc") drop outside t_orig pval } local i = `i' + 1 } texdoc init "tableA2.tex", replace force tex \begin{tabular}{lccccccccccccccccc} tex \toprule tex & \multicolumn{2}{c}{Help poor} & & \multicolumn{2}{c}{Assistance poor} & & \multicolumn{2}{c}{Work-luck} & & \multicolumn{2}{c}{Party affiliation} & & \multicolumn{2}{c}{Political views} & & \multicolumn{2}{c}{Voting Democrat} \\ tex \cmidrule{2-3}\cmidrule{5-6}\cmidrule{8-9}\cmidrule{11-12}\cmidrule{14-15}\cmidrule{17-18} & orig. & repl. & & orig. & repl. & & orig. & repl. & & orig. & repl. & & orig. & repl. & & orig. & repl. \\ tex & (1) & (2) & & (3) & (4) & & (5) & (6) & & (7) & (8) & & (9) & (10) & & (11) & (12) \\ tex \midrule tex Economic shock & 0.033** & `b_shock1' & & 0.021** & `b_shock2' & & 0.017** & `b_shock3' & & 0.177*** & `b_shock4' & & 0.133*** & `b_shock5' & & 0.043*** & `b_shock6' \\ tex \quad SE & (0.016) & (`se_shock1') & & (0.010) & (`se_shock2') & & (0.008) & (`se_shock3') & & (0.029) & (`se_shock4') & & (0.022) & (`se_shock5') & & (0.009) & (`se_shock6') \\ tex \quad p-value & & [`p1_shock1'] & & & [`p1_shock2'] & & & [`p1_shock3'] & & & [`p1_shock4'] & & & [`p1_shock5'] & & & [`p1_shock6'] \\ tex \quad p-value (WBS) & & [`p_shock1'] & & & [`p_shock2'] & & & [`p_shock3'] & & & [`p_shock4'] & & & [`p_shock5'] & & & [`p_shock6'] \\ tex Years of education & --0.051*** & `b_educ1' & & --0.017*** & `b_educ2' & & 0.006*** & `b_educ3' & & --0.033*** & `b_educ4' & & 0.020*** & `b_educ5' & & 0.001 & `b_educ6' \\ tex \quad SE & (0.006) & (`se_educ1') & & (0.002) & (`se_educ2') & & (0.002) & (`se_educ3') & & (0.008) & (`se_educ4') & & (0.007) & (`se_educ5') & & (0.002) & (`se_educ6') \\ tex \quad p-value & & [`p1_educ1'] & & & [`p1_educ2'] & & & [`p1_educ3'] & & & [`p1_educ4'] & & & [`p1_educ5'] & & & [`p1_educ6'] \\ tex \quad p-value (WBS) & & [`p_educ1'] & & & [`p_educ2'] & & & [`p_educ3'] & & & [`p_educ4'] & & & [`p_educ5'] & & & [`p_educ6'] \\ tex Married & --0.139*** & `b_married1' & & --0.036*** & `b_married2' & & --0.059*** & `b_married3' & & --0.220*** & `b_married4' & & --0.264*** & `b_married5' & & --0.058*** & `b_married6' \\ tex \quad SE & (0.016) & (`se_married1') & & (0.009) & (`se_married2') & & (0.008) & (`se_married3') & & (0.036) & (`se_married4') & & (0.014) & (`se_married5') & & (0.008) & (`se_married6') \\ tex \quad p-value & & [`p1_married1'] & & & [`p1_married2'] & & & [`p1_married3'] & & & [`p1_married4'] & & & [`p1_married5'] & & & [`p1_married6'] \\ tex \quad p-value (WBS) & & [`p_married1'] & & & [`p_married2'] & & & [`p_married3'] & & & [`p_married4'] & & & [`p_married5'] & & & [`p_married6'] \\ tex Female & 0.164*** & `b_female1' & & 0.060*** & `b_female2' & & --0.075*** & `b_female3' & & 0.178*** & `b_female4' & & 0.088*** & `b_female5' & & 0.037*** & `b_female6' \\ tex \quad SE & (0.017) & (`se_female1') & & (0.011) & (`se_female2') & & (0.012) & (`se_female3') & & (0.035) & (`se_female4') & & (0.019) & (`se_female5') & & (0.008) & (`se_female6') \\ tex \quad p-value & & [`p1_female1'] & & & [`p1_female2'] & & & [`p1_female3'] & & & [`p1_female4'] & & & [`p1_female5'] & & & [`p1_female6'] \\ tex \quad p-value (WBS) & & [`p_female1'] & & & [`p_female2'] & & & [`p_female3'] & & & [`p_female4'] & & & [`p_female5'] & & & [`p_female6'] \\ tex Black & 0.635*** & `b_black1' & & 0.276*** & `b_black2' & & 0.101*** & `b_black3' & & 1.468*** & `b_black4' & & 0.296*** & `b_black5' & & 0.449*** & `b_black6' \\ tex \quad SE & (0.036) & (`se_black1') & & (0.012) & (`se_black2') & & (0.010) & (`se_black3') & & (0.064) & (`se_black4') & & (0.021) & (`se_black5') & & (0.016) & (`se_black6') \\ tex \quad p-value & & [`p1_black1'] & & & [`p1_black2'] & & & [`p1_black3'] & & & [`p1_black4'] & & & [`p1_black5'] & & & [`p1_black6'] \\ tex \quad p-value (WBS) & & [`p_black1'] & & & [`p_black2'] & & & [`p_black3'] & & & [`p_black4'] & & & [`p_black5'] & & & [`p_black6'] \\ tex Unemployed & 0.118*** & `b_unem1' & & 0.066*** & `b_unem2' & & 0.058*** & `b_unem3' & & 0.113*** & `b_unem4' & & 0.100** & `b_unem5' & & 0.050*** & `b_unem6' \\ tex \quad SE & (0.036) & (`se_unem1') & & (0.016) & (`se_unem2') & & (0.005) & (`se_unem3') & & (0.040) & (`se_unem4') & & (0.045) & (`se_unem5') & & (0.014) & (`se_unem6') \\ tex \quad p-value & & [`p1_unem1'] & & & [`p1_unem2'] & & & [`p1_unem3'] & & & [`p1_unem4'] & & & [`p1_unem5'] & & & [`p1_unem6'] \\ tex \quad p-value (WBS) & & [`p_unem1'] & & & [`p_unem2'] & & & [`p_unem3'] & & & [`p_unem4'] & & & [`p_unem5'] & & & [`p_unem6'] \\ tex Observations & 24,287 & `n1' & & 15,416 & `n2' & & 30,694 & `n3' & & 43,443 & `n4' & & 38,525 & `n5' & & 27,267 & `n6' \\ tex R2 & 0.09 & `r1' & & 0.07 & `r2' & & 0.02 & `r3' & & 0.11 & `r4' & & 0.05 & `r5' & & 0.15 & `r6' \\ tex Mean of outcome & `out1' & 3.103 & & 2.557 & `out2' & & 1.459 & `out3' & & 3.318 & `out4' & & 3.905 & `out5' & & 0.498 & `out6' \\ tex Number of waves & - & `waves1' & & - & `waves2' & & - & `waves3' & & - & `waves4' & & - & `waves5' & & - & `waves6'\\ tex \bottomrule tex \end{tabular}% texdoc close /* End of file */