/* DO file for analysis at TTWA level, as reported in: Hilber, C.A.L. and W. Vermeulen (2014), The impact of supply constraints on house prices in England, Economic Journal. Tables and figures are produced in order of appearance in the paper. Christian Hilber and Wouter Vermeulen, 22 July 2014 */ clear clear matrix set memory 700m set matsize 2000 set more 1 adopath + g:\stata // path of folder that contains ado files (xtivreg2, outreg2) cd "M:\p_wvm\SERC Projecten\Supply constraints in the UK (LUCK)\Published data and code" // path of folder that contains data capture log close log using "analysis FUR.log", replace use "data FUR.dta", clear ******************************************************************************************** * * * DATA PREPARATIONS * * * ******************************************************************************************** * generate year dummies tab year, gen(year_) * standardize variables and create interaction effect with log earnings * note that variables are standardized to their standard deviation at LPA level!!! gen sd_refusal_maj_7908 = 0.08728701 gen sd_delchange_maj1 = 0.2200945 gen sd_pdevel90_m2 = 0.2331092 gen sd_pop_density_1911_imp = 2561.628 gen sd_range_meters = 171.1572 gen sd_labourvotes1983 = 0.0905755 #delimit ; foreach vari in refusal_maj_7908 delchange_maj1 pdevel90_m2 pop_density_1911_imp range_meters labourvotes1983 {; egen m_`vari' = mean(`vari'); gen st_`vari' = (`vari' - m_`vari') / sd_`vari'; gen `vari'_learn = st_`vari' * lmale_earn_real; }; #delimit cr ******************************************************************************************** * * * TABLE 4: BASELINE RESULTS AT FUR LEVEL USING BOTH INSTRUMENTS * * * ******************************************************************************************** #delimit ; xtivreg2 lrindex2 (refusal_maj_7908_learn pdevel90_m2_learn = delchange_maj1_learn labourvotes1983_learn pop_density_1911_imp_learn) lmale_earn_real range_meters_learn year_2-year_35 [aweight = furhou], fe robust cluster(furcode) ffirst; outreg2 using "Table 4 FUR", drop(year_*) nocon nolabel aster replace word nor2 title("Robustness to geographical scale") ctitle("FUR") addstat("Kleibergen-Paap F", e(widstat)); #delimit cr log close