sysdir set PERSONAL C:\users\aushahn\Documents\stataado

set logtype t
log using C:\users\aushahn\Documents\hahn_thesis3.18,replace

*import data + rescale + dummies
import delimited "C:\Users\aushahn\Documents\econ thesis\data\acs_2019\masters\masterdata_2.csv"
tabulate year, generate(yeardum)
replace pop = pop/1000 
xtile inc = fam, n(5)
tab inc, gen(inc)
summ
*birthrate as dependent 
*OLS
regress brate_n bach_p black_p enroll_p inc1-inc4 fert_pe pop   public unemploy shooting_1 shooting_2 deaths_1 deaths_2 year
regress birthrate bach_p black_p enroll_p inc1-inc4 fert_pe pop   public unemploy shooting_1 shooting_2 deaths_1 deaths_2 year

*random, fixed
xtset fips year
*NCHS
*
*
*
*
*NCHS - regular
xtreg brate_n bach_p black_p enroll_p inc1-inc4   fert_pe pop public unemploy shooting_1 shooting_2 deaths_1 deaths_2 yeardum1- yeardum10, re
xtoverid
*tests time invariant = 0 
xttest0
xtreg brate_n bach_p black_p enroll_p inc1-inc4   fert_pe pop public unemploy shooting_1 shooting_2 deaths_1 deaths_2 yeardum1- yeardum10, fe

*NCHS - individual variables 
*shooting 1
xtreg brate_n bach_p black_p enroll_p inc1-inc4   fert_pe pop public unemploy shooting_1 yeardum1- yeardum10, re
xtoverid
*tests time invariant = 0 
xttest0
xtreg brate_n bach_p black_p enroll_p inc1-inc4   fert_pe pop public unemploy shooting_1 yeardum1- yeardum10, fe

*shooting 2
xtreg brate_n bach_p black_p enroll_p inc1-inc4   fert_pe pop public unemploy shooting_2 yeardum1- yeardum10, re
xtoverid
*tests time invariant = 0 
xttest0
xtreg brate_n bach_p black_p enroll_p inc1-inc4   fert_pe pop public unemploy shooting_2 yeardum1- yeardum10, fe

*deaths 1
xtreg brate_n bach_p black_p enroll_p inc1-inc4   fert_pe pop public unemploy deaths_1 yeardum1- yeardum10, re
xtoverid
*tests time invariant = 0 
xttest0
xtreg brate_n bach_p black_p enroll_p inc1-inc4   fert_pe pop public unemploy deaths_1 yeardum1- yeardum10, fe

*deaths 2
xtreg brate_n bach_p black_p enroll_p inc1-inc4   fert_pe pop public unemploy deaths_2 yeardum1- yeardum10, re
xtoverid
*tests time invariant = 0 
xttest0
xtreg brate_n bach_p black_p enroll_p inc1-inc4   fert_pe pop public unemploy deaths_2 yeardum1- yeardum10, fe

*ACS
*
*
*
*
*ACS - Regular
xtreg birthr bach_p black_p enroll_p inc1-inc4   fert_pe pop public unemploy shooting_1 shooting_2 deaths_1 deaths_2 yeardum1- yeardum10, re
xtoverid
*tests time invariant = 0 
xttest0
xtreg birthr bach_p black_p enroll_p inc1-inc4   fert_pe pop public unemploy shooting_1 shooting_2 deaths_1 deaths_2 yeardum1- yeardum10, fe

*ACS - individual variables 
*shooting 1
xtreg birthr bach_p black_p enroll_p inc1-inc4   fert_pe pop public unemploy shooting_1 yeardum1- yeardum10, re
xtoverid
*tests time invariant = 0 
xttest0
xtreg birthr bach_p black_p enroll_p inc1-inc4   fert_pe pop public unemploy shooting_1 yeardum1- yeardum10, fe

*shooting 2
xtreg birthr bach_p black_p enroll_p inc1-inc4   fert_pe pop public unemploy shooting_2 yeardum1- yeardum10, re
xtoverid
*tests time invariant = 0 
xttest0
xtreg birthr bach_p black_p enroll_p inc1-inc4   fert_pe pop public unemploy shooting_2 yeardum1- yeardum10, fe

*deaths 1
xtreg brate_n bach_p black_p enroll_p inc1-inc4   fert_pe pop public unemploy deaths_1 yeardum1- yeardum10, re
xtoverid
*tests time invariant = 0 
xttest0
xtreg birthr bach_p black_p enroll_p inc1-inc4   fert_pe pop public unemploy deaths_1 yeardum1- yeardum10, fe

*deaths 2
xtreg birthr bach_p black_p enroll_p inc1-inc4   fert_pe pop public unemploy deaths_2 yeardum1- yeardum10, re
xtoverid
*tests time invariant = 0 
xttest0
xtreg birthr bach_p black_p enroll_p inc1-inc4   fert_pe pop public unemploy deaths_2 yeardum1- yeardum10, fe

*
*
* interaction models
*
*
*

corr shooting_1 shooting_2 deaths_1 deaths_2

*create interaction terms
gen shoot1_white = (shooting_1yr * white_perc)
gen death1_white = (deaths_1yr * white_perc)
gen shoot1_inc1 = (shooting_1yr * inc1)
gen shoot1_inc2 = (shooting_1yr * inc2)
gen shoot1_inc3 = (shooting_1yr * inc3)
gen shoot1_inc4 = (shooting_1yr * inc4)
gen death1_inc1 = (shooting_1yr * inc1)
gen death1_inc2 = (shooting_1yr * inc2)
gen death1_inc3 = (shooting_1yr * inc3)
gen death1_inc4 = (shooting_1yr * inc4)
gen shoot1_enroll = (shooting_1yr * enroll_perc)
gen death1_enroll = (deaths_1yr * enroll_perc)

*Shooting interaction

*NCHS
xtreg brate_n bach_p black_p enroll_p inc1-inc4 fert_pe pop public unemploy shooting_1 shoot1_w shoot1_inc1-shoot1_inc4 shoot1_e yeardum1- yeardum10, re
xtoverid
*tests time invariant = 0 
xttest0
xtreg brate_n bach_p black_p enroll_p inc1-inc4 fert_pe pop public unemploy shooting_1 shoot1_w shoot1_inc1-shoot1_inc4 shoot1_e yeardum1- yeardum10, fe
*whiteness
xtreg brate_n bach_p black_p enroll_p inc1-inc4 fert_pe pop public unemploy shooting_1 shoot1_w yeardum1- yeardum10, fe
*income
xtreg brate_n bach_p black_p enroll_p inc1-inc4 fert_pe pop public unemploy shooting_1 shoot1_inc1-shoot1_inc4 yeardum1- yeardum10, fe
*enrollment
xtreg brate_n bach_p black_p enroll_p inc1-inc4 fert_pe pop public unemploy shooting_1 shoot1_e yeardum1- yeardum10, fe

*Deaths interaction

*NCHS
xtreg brate_n bach_p black_p enroll_p inc1-inc4   fert_pe pop public unemploy deaths_1 death1_w death1_inc1-death1_inc4 death1_e yeardum1- yeardum10, re
xtoverid
*tests time invariant = 0 
xttest0
xtreg brate_n bach_p black_p enroll_p inc1-inc4   fert_pe pop public unemploy deaths_1 death1_w death1_inc1-death1_inc4 death1_e yeardum1- yeardum10, fe
*white
xtreg brate_n bach_p black_p enroll_p inc1-inc4   fert_pe pop public unemploy deaths_1 death1_w yeardum1- yeardum10, fe
*income
xtreg brate_n bach_p black_p enroll_p inc1-inc4   fert_pe pop public unemploy deaths_1 death1_inc1-death1_inc4  yeardum1- yeardum10, fe
*enroll
xtreg brate_n bach_p black_p enroll_p inc1-inc4   fert_pe pop public unemploy deaths_1 death1_e yeardum1- yeardum10, fe


*2 year vars
gen shoot2_white = (shooting_2yr * white_perc)
gen death2_white = (deaths_2yr * white_perc)
gen shoot2_inc1 = (shooting_2yr * inc1)
gen shoot2_inc2 = (shooting_2yr * inc2)
gen shoot2_inc3 = (shooting_2yr * inc3)
gen shoot2_inc4 = (shooting_2yr * inc4)
gen death2_inc1 = (shooting_2yr * inc1)
gen death2_inc2 = (shooting_2yr * inc2)
gen death2_inc3 = (shooting_2yr * inc3)
gen death2_inc4 = (shooting_2yr * inc4)
gen shoot2_enroll = (shooting_2yr * enroll_perc)
gen death2_enroll = (deaths_2yr * enroll_perc)

*Shooting interaction2

*NCHS
xtreg brate_n bach_p black_p enroll_p inc1-inc4 fert_pe pop public unemploy shooting_2 shoot2_w shoot2_inc1-shoot2_inc4 shoot2_e yeardum1- yeardum10, re
xtoverid
*tests time invariant = 0 
xttest0
xtreg brate_n bach_p black_p enroll_p inc1-inc4 fert_pe pop public unemploy shooting_2 shoot2_w shoot2_inc1-shoot2_inc4 shoot2_e yeardum1- yeardum10, fe
*whiteness2
xtreg brate_n bach_p black_p enroll_p inc1-inc4 fert_pe pop public unemploy shooting_2 shoot2_w yeardum1- yeardum10, fe
*income2
xtreg brate_n bach_p black_p enroll_p inc1-inc4 fert_pe pop public unemploy shooting_2 shoot2_inc1-shoot2_inc4 yeardum1- yeardum10, fe
*enrollment2
xtreg brate_n bach_p black_p enroll_p inc1-inc4 fert_pe pop public unemploy shooting_2 shoot2_e yeardum1- yeardum10, fe

*Deaths interaction2

*NCHS
xtreg brate_n bach_p black_p enroll_p inc1-inc4   fert_pe pop public unemploy deaths_2 death2_w death2_inc1-death2_inc4 death2_e yeardum1- yeardum10, re
xtoverid
*tests time invariant = 0 
xttest0
xtreg brate_n bach_p black_p enroll_p inc1-inc4   fert_pe pop public unemploy deaths_2 death2_w death2_inc1-death2_inc4 death2_e yeardum1- yeardum10, fe
*white2
xtreg brate_n bach_p black_p enroll_p inc1-inc4   fert_pe pop public unemploy deaths_2 death2_w yeardum1- yeardum10, fe
*income2
xtreg brate_n bach_p black_p enroll_p inc1-inc4   fert_pe pop public unemploy deaths_2 death2_inc1-death2_inc4  yeardum1- yeardum10, fe
*enroll
xtreg brate_n bach_p black_p enroll_p inc1-inc4   fert_pe pop public unemploy deaths_2 death2_e yeardum1- yeardum10, fe





log close