RESET $


? ********************************************************************
?
?     Sept 28, 2005
?
?     REStud
?     Au and Henderson, "Are Chinese Cities Too Small?"
?     LIMDEP program for regressions in Table 2
? 
?     Prefecture level cities 1997 on 1990 instruments
?
?     1. Table 2, 1st Col. - IV Estimation, structural model
?     2. Table 2, 2nd Col. - Ordinary NLLS, structural model
?
? ********************************************************************

READ; File = restud.xls
    ; Format = XLS
    ; Names
    ; Nvar = 27
    ; Nobs = 205 $

CREATE; rysna = 1/(1+ryms)
      ; area9015 = area90^(1.5) 
      ; ms90ar90 = ryms90*area90 $

TYPE; -------------------------------------------------------------------- $
TYPE; 1. Table 2, 1st Col. - IV Estimation, structural model
TYPE; -------------------------------------------------------------------- $

NLSQ; Lhs = lnyna
    ; Fcn = constant + a0*lnk + a1*log(lna-a2*(lna^1.5)) 
                + a3*rysna*log(lna-a2*(lna^1.5)) + _pshs90*pshs90 
                + _rfdlnaN*rfdlnaN + _lmp6976*lmp6_976 + _tint6M6*tcint6M6 
    ; Labels = constant, a0, a1, a2, a3, _pshs90, _rfdlnaN, _lmp6976,
               _tint6M6
    ; Start = .2, .40, 1.1, -.15, -.05, .007, 1.0, .5, 4
    ; Inst = ONE,area90,book90,doctor90,tel90,lnroad90,drw,drc,agn90,
             lnrkl90,nofdi90,lnisiv90,ms90ar90,rfdlN90,lmp6_906,tnt6M906,
             pshs90,area9015 
    ; Res = residual $

TYPE; First stage regressions $

REGRESS; Lhs = lnk
       ; Rhs = ONE,area90,book90,doctor90,tel90,lnroad90,drw,drc,agn90,
               lnrkl90,nofdi90,lnisiv90,ms90ar90,rfdlN90,lmp6_906,tnt6M906,
               pshs90,area9015 $

REGRESS; Lhs = lna
       ; Rhs = ONE,area90,book90,doctor90,tel90,lnroad90,drw,drc,agn90,
               lnrkl90,nofdi90,lnisiv90,ms90ar90,rfdlN90,lmp6_906,tnt6M906,
               pshs90,area9015 $

REGRESS; Lhs = rysna
       ; Rhs = ONE,area90,book90,doctor90,tel90,lnroad90,drw,drc,agn90,
               lnrkl90,nofdi90,lnisiv90,ms90ar90,rfdlN90,lmp6_906,tnt6M906,
               pshs90,area9015 $

TYPE; Note: *** pshs90 is a regressor as well as an instrument *** $
TYPE; $

REGRESS; Lhs = rfdlnaN
       ; Rhs = ONE,area90,book90,doctor90,tel90,lnroad90,drw,drc,agn90,
               lnrkl90,nofdi90,lnisiv90,ms90ar90,rfdlN90,lmp6_906,tnt6M906,
               pshs90,area9015 $

REGRESS; Lhs = lmp6_976
       ; Rhs = ONE,area90,book90,doctor90,tel90,lnroad90,drw,drc,agn90,
               lnrkl90,nofdi90,lnisiv90,ms90ar90,rfdlN90,lmp6_906,tnt6M906,
               pshs90,area9015 $

REGRESS; Lhs = tcint6M6
       ; Rhs = ONE,area90,book90,doctor90,tel90,lnroad90,drw,drc,agn90,
               lnrkl90,nofdi90,lnisiv90,ms90ar90,rfdlN90,lmp6_906,tnt6M906,
               pshs90,area9015 $

TYPE; -------------------------------------------------------------------- $
TYPE; Residual regression on instruments $

REGRESS; Lhs = residual
       ; Rhs = ONE,area90,book90,doctor90,tel90,lnroad90,drw,drc,agn90,
               lnrkl90,nofdi90,lnisiv90,ms90ar90,rfdlN90,lmp6_906,tnt6M906,
               pshs90,area9015 $




TYPE; -------------------------------------------------------------------- $
TYPE; 2. Table 2, 2nd Col. - Ordinary NLLS, structural model
TYPE; -------------------------------------------------------------------- $

NLSQ; Lhs = lnyna
    ; Fcn = constant + a0*lnk + a1*log(lna-a2*(lna^1.5)) 
                + a3*rysna*log(lna-a2*(lna^1.5)) + _pshs90*pshs90 
                + _rfdlnaN*rfdlnaN + _lmp6976*lmp6_976 + _tint6M6*tcint6M6 
    ; Labels = constant, a0, a1, a2, a3, _pshs90, _rfdlnaN, _lmp6976,
               _tint6M6
    ; Start = .2, .40, 1.1, -.15, -.05, .007, 1.0, .5, 4 $