clear
set mem 32m
set matsize 800
use SB_Paper

gen lfrsvnoc  = log(farsnocc/(vmtrural+vmturban))
gen lfrsvocc  = log(farsocc /(vmtrural+vmturban))
gen lvmtrura  = log(vmtrural)
gen lvmturba  = log(vmturban)
gen ldensrur  = log(densrur)
gen ldensurb  = log(densurb)
gen lusage    = log(usage)
gen lincome   = log(percapin)
gen lage      = log(meanage)
gen lpblack   = log(precentb) 
gen lphisp    = log(precenth) 
gen lviopcap  = log(viopcap)
gen lpropcap  = log(proppcap)
gen lunemp    = log(unemp)
gen lfueltax  = log(fueltax)

* full IV regression

xi: ivreg lfrsvocc (lusage=ds dp dsp) lincome lage lpblack lphisp ldensrur ldensurb lviopcap lpropcap lunemp lvmtrura lvmturba lfueltax lim65 lim70p mlda21 bac08 i.fips i.year, robust
xi: ivreg lfrsvnoc (lusage=ds dp dsp) lincome lage lpblack lphisp ldensrur ldensurb lviopcap lpropcap lunemp lvmtrura lvmturba lfueltax lim65 lim70p mlda21 bac08 i.fips i.year, robust

* standard regression with FE

xi: reg   lfrsvocc lusage lincome lage lpblack lphisp ldensrur ldensurb lviopcap lpropcap lunemp lvmtrura lvmturba lfueltax lim65 lim70p mlda21 bac08 i.fips i.year, robust
xi: reg   lfrsvnoc lusage lincome lage lpblack lphisp ldensrur ldensurb lviopcap lpropcap lunemp lvmtrura lvmturba lfueltax lim65 lim70p mlda21 bac08 i.fips i.year, robust

* standard regression with no FE

xi: reg   lfrsvocc lusage lincome lage lpblack lphisp ldensrur ldensurb lviopcap lpropcap lunemp lvmtrura lvmturba lfueltax lim65 lim70p mlda21 bac08 i.year, robust
xi: reg   lfrsvnoc lusage lincome lage lpblack lphisp ldensrur ldensurb lviopcap lpropcap lunemp lvmtrura lvmturba lfueltax lim65 lim70p mlda21 bac08 i.year, robust

* usage regression (can use boxcox as well)

xi: reg  usage ds dp dsp lincome lage lpblack lphisp ldensrur ldensurb lviopcap lpropcap lunemp lvmtrura lvmturba lfueltax lim65 lim70p mlda21 bac08 i.fips i.year, robust
xi: reg  usage ds dp dsp lincome lage lpblack lphisp ldensrur ldensurb lviopcap lpropcap lunemp lvmtrura lvmturba lfueltax lim65 lim70p mlda21 bac08 i.year, robust