This do-file creates Table 2 (p. 63), Table 3 (p. 64-65), and Table 4 (p. 67) in the SMH paper.
// Task: models // Project: insecurity in context ****************************************** *** 0. Program setup ****************************************** version 11.2 set linesize 80 set maxvar 10000 clear all macro drop _all ****************************************** *** 1. Open the data set ****************************************** use "2013-05-11 vars generated.dta" drop if age > 65 drop if selfemployed == 1 ****************************************** *** 2. Analytic sample and weight ****************************************** misschk generalhappy mntlhlth jobinsecurity logincome female age black divorced widow nvmarried babies preteen teens selfemployed PT year, gen(m) // weight set-up svyset sampcode [pw = wtssnr], strata(year) ****************************************** *** 3. Prepare vars ****************************************** sum incomerec if mnumber == 0, d gen income1 = (incomerec <= 18750) if mnumber == 0 & incomerec < . gen income2 = (incomerec > 18750 & incomerec <= 32500) if mnumber == 0 & incomerec < . gen income3 = (incomerec > 32500 & incomerec <= 55000) if mnumber == 0 & incomerec < . gen income4 = (incomerec > 55000) if mnumber == 0 & incomerec < . gen incomer = 1 if income1 == 1 replace incomer = 2 if income2 == 1 replace incomer = 3 if income3 == 1 replace incomer = 4 if income4 == 1 gen jobinsecurityyear2010 = jobinsecurity * year2010 label var jobinsecurityyear2010 "jobinsecurity * year2010" global covar "female age agesq black divorced widow nvmarried babies preteen teens selfemployed PT" // global output global outreg1 "outreg2 using "gss1", excel bdec(3) alpha(0.001, 0.01, 0.05, 0.1) symbol(***, **, *, +)" global outreg2 "outreg2 using "gss2", excel bdec(3) alpha(0.001, 0.01, 0.05, 0.1) symbol(***, **, *, +)" global outreg3 "outreg2 using "gss3", excel bdec(3) alpha(0.001, 0.01, 0.05, 0.1) symbol(***, **, *, +)" ****************************************** *** 4. Tables ****************************************** ************************** *** 4.1. Table 2 ************************** foreach var of varlist generalhappy { svy, subpop(if mnumber == 0): reg `var' year2010 jobinsecurity $covar logincome if mnumber == 0 $outreg1 svy, subpop(if mnumber == 0 & year == 2006): reg `var' year2010 jobinsecurity $covar logincome if mnumber == 0 & year == 2006 est store m1 $outreg1 svy, subpop(if mnumber == 0 & year == 2010): reg `var' year2010 jobinsecurity $covar logincome if mnumber == 0 & year == 2010 est store m2 $outreg1 } suest m1 m2 test [m1]jobinsecurity = [m2]jobinsecurity foreach var of varlist mntlhlth { svy, subpop(if mnumber == 0): nbreg `var' year2010 jobinsecurity $covar logincome if mnumber == 0 $outreg1 svy, subpop(if mnumber == 0 & year == 2006): nbreg `var' year2010 jobinsecurity $covar logincome if mnumber == 0 & year == 2006 est store m3 $outreg1 svy, subpop(if mnumber == 0 & year == 2010): nbreg `var' year2010 jobinsecurity $covar logincome if mnumber == 0 & year == 2010 est store m4 $outreg1 } suest m3 m4 test [m3_mntlhlth]jobinsecurity = [m4_mntlhlth]jobinsecurity est clear ************************** *** 4.2. Table 3 ************************** foreach var of varlist generalhappy { svy, subpop(if mnumber == 0 & income1 == 1): reg `var' year2010 jobinsecurity $covar if mnumber == 0 & income1 == 1 est store m1 $outreg2 svy, subpop(if mnumber == 0 & income2 == 1): reg `var' year2010 jobinsecurity $covar if mnumber == 0 & income2 == 1 est store m2 $outreg2 svy, subpop(if mnumber == 0 & income3 == 1): reg `var' year2010 jobinsecurity $covar if mnumber == 0 & income3 == 1 est store m3 $outreg2 svy, subpop(if mnumber == 0 & income4 == 1): reg `var' year2010 jobinsecurity $covar if mnumber == 0 & income4 == 1 est store m4 $outreg2 svy, subpop(if mnumber == 0 & income1 == 1): reg `var' year2010 jobinsecurity jobinsecurityyear2010 $covar if mnumber == 0 & income1 == 1 est store m5 $outreg2 svy, subpop(if mnumber == 0 & income2 == 1): reg `var' year2010 jobinsecurity jobinsecurityyear2010 $covar if mnumber == 0 & income2 == 1 est store m6 $outreg2 svy, subpop(if mnumber == 0 & income3 == 1): reg `var' year2010 jobinsecurity jobinsecurityyear2010 $covar if mnumber == 0 & income3 == 1 est store m7 $outreg2 svy, subpop(if mnumber == 0 & income4 == 1): reg `var' year2010 jobinsecurity jobinsecurityyear2010 $covar if mnumber == 0 & income4 == 1 est store m8 $outreg2 } suest m1 m2 m3 m4 test [m1]jobinsecurity = [m2]jobinsecurity = [m3]jobinsecurity = [m4]jobinsecurity suest m5 m6 m7 m8 test [m5]jobinsecurityyear2010 = [m6]jobinsecurityyear2010 = [m7]jobinsecurityyear2010 = [m8]jobinsecurityyear2010 est clear foreach var of varlist mntlhlth { svy, subpop(if mnumber == 0 & income1 == 1): nbreg `var' year2010 jobinsecurity $covar if mnumber == 0 & income1 == 1 est store m1 $outreg2 svy, subpop(if mnumber == 0 & income2 == 1): nbreg `var' year2010 jobinsecurity $covar if mnumber == 0 & income2 == 1 est store m2 $outreg2 svy, subpop(if mnumber == 0 & income3 == 1): nbreg `var' year2010 jobinsecurity $covar if mnumber == 0 & income3 == 1 est store m3 $outreg2 svy, subpop(if mnumber == 0 & income4 == 1): nbreg `var' year2010 jobinsecurity $covar if mnumber == 0 & income4 == 1 est store m4 $outreg2 svy, subpop(if mnumber == 0 & income1 == 1): nbreg `var' year2010 jobinsecurity jobinsecurityyear2010 $covar if mnumber == 0 & income1 == 1 est store m5 $outreg2 svy, subpop(if mnumber == 0 & income2 == 1): nbreg `var' year2010 jobinsecurity jobinsecurityyear2010 $covar if mnumber == 0 & income2 == 1 est store m6 $outreg2 svy, subpop(if mnumber == 0 & income3 == 1): nbreg `var' year2010 jobinsecurity jobinsecurityyear2010 $covar if mnumber == 0 & income3 == 1 est store m7 $outreg2 svy, subpop(if mnumber == 0 & income4 == 1): nbreg `var' year2010 jobinsecurity jobinsecurityyear2010 $covar if mnumber == 0 & income4 == 1 est store m8 $outreg2 } suest m1 m2 m3 m4 test [m1_mntlhlth]jobinsecurity = [m2_mntlhlth]jobinsecurity = [m3_mntlhlth]jobinsecurity = [m4_mntlhlth]jobinsecurity suest m5 m6 m7 m8 test [m5_mntlhlth]jobinsecurityyear2010 = [m6_mntlhlth]jobinsecurityyear2010 = [m7_mntlhlth]jobinsecurityyear2010 = [m8_mntlhlth]jobinsecurityyear2010 est clear ************************** *** 4.3. Table 4 ************************** foreach var of varlist generalhappy { svy, subpop(if mnumber == 0 & female == 0): reg `var' year2010 jobinsecurity $covar logincome if mnumber == 0 & female == 0 est store m1 $outreg3 svy, subpop(if mnumber == 0 & female == 1): reg `var' year2010 jobinsecurity $covar logincome if mnumber == 0 & female == 1 est store m2 $outreg3 svy, subpop(if mnumber == 0 & female == 0): reg `var' year2010 jobinsecurity jobinsecurityyear2010 $covar logincome if mnumber == 0 & female == 0 est store m3 $outreg3 svy, subpop(if mnumber == 0 & female == 1): reg `var' year2010 jobinsecurity jobinsecurityyear2010 $covar logincome if mnumber == 0 & female == 1 est store m4 $outreg3 } suest m1 m2 test [m1]jobinsecurity = [m2]jobinsecurity suest m3 m4 test [m3]jobinsecurityyear2010 = [m4]jobinsecurityyear2010 est clear foreach var of varlist mntlhlth { svy, subpop(if mnumber == 0 & female == 0): nbreg `var' year2010 jobinsecurity $covar logincome if mnumber == 0 & female == 0 est store m1 $outreg3 svy, subpop(if mnumber == 0 & female == 1): nbreg `var' year2010 jobinsecurity $covar logincome if mnumber == 0 & female == 1 est store m2 $outreg3 svy, subpop(if mnumber == 0 & female == 0): nbreg `var' year2010 jobinsecurity jobinsecurityyear2010 $covar logincome if mnumber == 0 & female == 0 est store m3 $outreg3 svy, subpop(if mnumber == 0 & female == 1): nbreg `var' year2010 jobinsecurity jobinsecurityyear2010 $covar logincome if mnumber == 0 & female == 1 est store m4 $outreg3 } suest m1 m2 test [m1_mntlhlth]jobinsecurity = [m2_mntlhlth]jobinsecurity suest m3 m4 test [m3_mntlhlth]jobinsecurityyear2010 = [m4_mntlhlth]jobinsecurityyear2010 est clear