#delimit ;
set more 1 ;

* THIS PROGRAM ESTIMATES THE MODELS PRESENTED IN 
Dirk W. Early and Edgar O. Olsen (2002) 
"Subsidized Housing, Emergency Shelters, and Homelessness: 
An Empirical Investigation Using Data from the 1990 Census", 
Advances in Economic Analysis & Policy: Vol. 2: No. 1, Article 2.
http://www.bepress.com/bejeap/advances/vol2/iss1/art2 ;

* DATA SET IS Early-Olsen.dta ;

* ==================== DESCRIPTIVE STATISTICS ====================== ;
* TABLE 1 TABLE 1 TABLE 1 TABLE 1 TABLE 1 TABLE 1 TABLE 1 TABLE 1  ;

#delimit ;
set more 1 ;
mvdecode _all, mv(-99) ;
sum shelter street poor2 homepov2 shltpov2 strtpov2 fracshlt nbeds units pi
 fbedpov2 fitbdpov fbedhome fitbdhom rshltql hpi_r vr rc 
 rafdc rga rmcapov ur halfpov2 blackpov fempov avetemp temp averain precip
 rmental unitpov2 less5k 
 crime neast south west ada rpland rboeckh
 rmedinc per1000 ngov92 ngovcap ;


* SECOND STAGE REGRESSIONS ;
* ================ linear - bedspov2 endogenous ===================== ;
* TABLE 3 TABLE 3 TABLE 3 TABLE 3 TABLE 3 TABLE 3 TABLE 3 ;


* Table 3 (1) ;
regress homepov2 halfpov2 ur hpi_r vr fbedpov2 rshltql rmcapov zunitpov z01unitp 
  zless5k z01less5 rc rafdc zrga z01rga rmental blackpov fempov avetemp averain zcrime z01crime 
  neast south west
                (halfpov2 ur hpi_r vr          rshltql rmcapov zunitpov z01unitp  
  zless5k z01less5 rc rafdc zrga z01rga rmental blackpov fempov avetemp averain zcrime z01crime 
  neast south west
  zrpland z01rplan rboeckh rmedinc ada ngov92 ngovcap per1000) ;


* Table 3 (2) ;
regress shltpov2 halfpov2 ur hpi_r vr fbedpov2 rshltql rmcapov zunitpov z01unitp 
  zless5k z01less5 rc rafdc zrga z01rga rmental blackpov fempov avetemp averain zcrime z01crime 
  neast south west
                (halfpov2 ur hpi_r vr          rshltql rmcapov zunitpov z01unitp  
  zless5k z01less5 rc rafdc zrga z01rga rmental blackpov fempov avetemp averain zcrime z01crime 
  neast south west
  zrpland z01rplan rboeckh rmedinc ada ngov92 ngovcap per1000)  ;


* Table 3 (3) ;
regress strtpov2 halfpov2 ur hpi_r vr fbedpov2 rshltql rmcapov zunitpov z01unitp 
  zless5k z01less5 rc rafdc zrga z01rga rmental blackpov fempov avetemp averain zcrime z01crime 
  neast south west
                (halfpov2 ur hpi_r vr          rshltql rmcapov zunitpov z01unitp  
  zless5k z01less5 rc rafdc zrga z01rga rmental blackpov fempov avetemp averain zcrime z01crime 
  neast south west
  zrpland z01rplan rboeckh rmedinc ada ngov92 ngovcap per1000) ;



* ================ linear - bedspov2 exogenous ===================== ;
* TABLE 4  TABLE 4  TABLE 4  TABLE 4  TABLE 4  TABLE 4  TABLE 4  TABLE 4 ;	

* Table 4 (1) ;
regress homepov2 halfpov2 ur hpi_r vr fbedpov2 rshltql rmcapov zunitpov z01unitp 
  zless5k z01less5 rc rafdc zrga z01rga rmental blackpov fempov avetemp averain zcrime z01crime 
  neast south west ;

* Table 4 (2) ;
regress shltpov2 halfpov2 ur hpi_r vr fbedpov2 rshltql rmcapov zunitpov z01unitp 
  zless5k z01less5 rc rafdc zrga z01rga rmental blackpov fempov avetemp averain zcrime z01crime 
  neast south west ;

* Table 4 (3) ;
regress strtpov2 halfpov2 ur hpi_r vr fbedpov2 rshltql rmcapov zunitpov z01unitp 
  zless5k z01less5 rc rafdc zrga z01rga rmental blackpov fempov avetemp averain zcrime z01crime 
  neast south west ;




* ====== endog bedspov2 - linear, weights, log, double log  === ;
* TABLE 5  TABLE 5  TABLE 5  TABLE 5  TABLE 5  TABLE 5  TABLE 5  TABLE 5  TABLE 5 ;


* ================ interactive terms zless5k hpi_r z_age62p zfamsize ================ ;

* Table 5 (1);
regress lnhpov2 halfpov2 ur hpi_r vr fbedpov2 rshltql rmcapov zunitpov z01unitp 
  zless5k z01less5 rc rafdc zrga z01rga rmental blackpov fempov avetemp averain zcrime z01crime 
  neast south west
                (halfpov2 ur hpi_r vr          rshltql rmcapov zunitpov z01unitp  
  zless5k z01less5 rc rafdc zrga z01rga rmental blackpov fempov avetemp averain zcrime z01crime 
  neast south west
  zrpland z01rplan rboeckh rmedinc ada ngov92 ngovcap per1000)     if lnhpov > -99 ;

* Table 5 (2) ;
regress lnhpov2 lnhalfpo lnur lnhpi_r lnvr lnfbedpo lnrshltq lnrmcapo zunitpov z01unitp 
  zless5k z01less5 rc lnrafdc zrga z01rga lnrmenta lnblackp lnfempov lnavetem lnaverai
  zcrime z01crime neast south west
                (lnhalfpo lnur lnhpi_r lnvr lnfbedpo lnrshltq lnrmcapo zunitpov z01unitp 
  zless5k z01less5 rc lnrafdc zrga z01rga lnrmenta lnblackp lnfempov lnavetem lnaverai
  zcrime z01crime neast south west
  zrpland z01rplan lnrboeck lnrmedin ada lnngov92 lnngovca lnper100)     if lnhpov > -99  ;

* Table 5 (3) ;
regress homepov2 halfpov2 ur hpi_r vr fbedpov2 rshltql rmcapov zunitpov z01unitp 
  zless5k z01less5 int5k_pi int5k_fs z015k_fs rc rafdc zrga z01rga rmental blackpov
  fempov avetemp averain zcrime z01crime neast south west
                (halfpov2 ur hpi_r vr          rshltql rmcapov zunitpov z01unitp  
  zless5k z01less5 int5k_pi int5k_fs z015k_fs rc rafdc zrga z01rga rmental blackpov
  fempov avetemp averain zcrime z01crime neast south west
  zrpland z01rplan rboeckh rmedinc ada ngov92 ngovcap per1000) ;


* ================ fraction of homeless in shelters w/ most compelling variables  =============== ;
* TABLE 7  TABLE 7  TABLE 7  TABLE 7  TABLE 7  TABLE 7  TABLE 7  TABLE 7  TABLE 7 ;

* Table 7 (1) ;
regress fracshlt fbedhome rshltql rmcapov blackpov fempov temp precip zcrime z01crime
  neast south west
 (                                    rshltql rmcapov blackpov fempov temp precip zcrime z01crime  
  neast south west
  zrpland z01rplan rboeckh rmedinc ada ngov92 ngovcap per1000)   if fracshlt > -99 ;

* Table 7 (2) ;
regress fracshlt halfpov2 ur hpi_r vr fbedhome rshltql rmcapov zunitpov z01unitp 
 zless5k z01less5 rc rafdc zrga z01rga rmental blackpov fempov temp precip zcrime z01crime
  neast south west
 (                     halfpov2 ur hpi_r vr                rshltql rmcapov zunitpov z01unitp 
 zless5k z01less5 rc rafdc zrga z01rga rmental blackpov fempov temp precip zcrime z01crime
  neast south west
  zrpland z01rplan rboeckh rmedinc ada ngov92 ngovcap per1000)   if fracshlt > -99 ;


* ===================== FIRST STAGE REGRESSIONS  ====================== ;
* SHELTER BEDS PER POOR PERSON ;
* Appendix A (1) ;
capture drop fitbeds  ;
regress fbedpov2 
  halfpov2 ur hpi_r vr          rshltql rmcapov zunitpov z01unitp  
  zless5k z01less5 rc rafdc zrga z01rga rmental blackpov fempov avetemp averain zcrime z01crime 
  neast south west
  zrpland z01rplan rboeckh rmedinc ada ngov92 ngovcap per1000 ;


* SHELTER BEDS PER HOMELESS PERSON W/ MOST THEORETICALLY COMPELLING VARIABLES ;
* Appendix A (2) ;
capture drop fitbeds2  ;
regress fbedhome 
  rshltql rmcapov blackpov fempov temp precip zcrime z01crime  
  neast south west
  zrpland z01rplan rboeckh rmedinc ada ngov92 ngovcap per1000  if fracshlt ~= . ;


* SHELTER BEDS PER HOMELESS PERSON W/  EXHAUSTIVE LIST OF VARIABLES ;
* Appendix A (3) ;
capture drop fitbeds3  ;
regress fbedhome halfpov2 ur hpi_r vr                rshltql rmcapov zunitpov z01unitp 
 zless5k z01less5 rc rafdc zrga z01rga rmental blackpov fempov temp precip zcrime z01crime
  neast south west
  zrpland z01rplan rboeckh rmedinc ada ngov92 ngovcap per1000   if fracshlt ~= . ;