clear
cd directory
set more off
set mem 2000000

use birth_data, clear

*state specific quadratic time trend
gen time = year - 1980
gen time2 = time^2
xi i.state i.state*time i.state*time2 i.year 

* Panel A
reg married bldtest _I* hs_grad black dmage [w=count], cluster(state)
est store all
reg married bldtest _I* hs_grad dmage [w=count] if black==1, cluster(state)
est store black
reg married bldtest _I* black dmage [w=count] if hs_grad==0, cluster(state)
est store noHS
reg married bldtest _I* hs_grad black dmage [w=count] if dmage<25, cluster(state)
est store young
est table all black noHS young, b(%9.4f) se(%9.4f) keep(bldtest hs_grad black dmage)

* Panel B
replace distance=distance/100
reg married distance _I* hs_grad black dmage [w=count], cluster(state)
est store all
reg married distance _I* hs_grad dmage [w=count] if black==1, cluster(state)
est store black
reg married distance _I* black dmage [w=count] if hs_grad==0, cluster(state)
est store noHS
reg married distance _I* hs_grad black dmage [w=count] if dmage<25, cluster(state)
est store young
est table all black noHS young, b(%9.4f) se(%9.4f) keep(distance hs_grad black dmage)

sum married if hs_grad~=. [w=count]
sum married [w=count] if black==1 & hs_grad~=. 
sum married [w=count] if hs_grad==0 & hs_grad~=.
sum married [w=count] if dmage<25 & hs_grad~=.


log close