Как включить фиктивные переменные в модель ivreg? - PullRequest
Новая
       18

Как включить фиктивные переменные в модель ivreg?

0 голосов
/ 04 июля 2018

У меня есть следующая модель: 

ivreg ldemand social_housing transport year (lprice = utilities)

Однако я хочу сделать year фиктивной переменной.

Как я могу сделать это в Stata?

Использование i.year не работает для команды ivreg.

Перекрестно опубликовано Статистик .

1 Ответ

0 голосов
/ 04 июля 2018

Команда ivreg не допускает факторные переменные: 

. webuse hsng2, clear

. ivreg rent pcturban i.region (hsngval = faminc)
factor variables not allowed
r(101);

Однако вы можете использовать префикс xi для создания манекенов на лету: 

.  xi: ivreg rent pcturban  i.region (hsngval = faminc)
i.region          _Iregion_1-4        (naturally coded; _Iregion_1 omitted)

Instrumental variables (2SLS) regression

      Source |       SS           df       MS      Number of obs   =        50
-------------+----------------------------------   F(5, 44)        =      9.10
       Model |  12735.4667         5  2547.09334   Prob > F        =    0.0000
    Residual |  48507.6533        44  1102.44667   R-squared       =    0.2079
-------------+----------------------------------   Adj R-squared   =    0.1179
       Total |    61243.12        49  1249.85959   Root MSE        =    33.203

------------------------------------------------------------------------------
        rent |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     hsngval |   .0038683   .0008958     4.32   0.000     .0020629    .0056737
    pcturban |  -.4980121   .5179779    -0.96   0.342    -1.541928    .5459039
  _Iregion_2 |   1.528672   15.14086     0.10   0.920    -28.98572    32.04306
  _Iregion_3 |    7.74279   15.10906     0.51   0.611    -22.70752     38.1931
  _Iregion_4 |  -40.61235   19.60999    -2.07   0.044    -80.13369   -1.091002
       _cons |   88.26681   31.69154     2.79   0.008     24.39671    152.1369
------------------------------------------------------------------------------
Instrumented:  hsngval
Instruments:   pcturban _Iregion_2 _Iregion_3 _Iregion_4 faminc
------------------------------------------------------------------------------

Важно отметить, что в соответствии с файлом справки команды:

Устаревшая команда

ivreg - это устаревшая команда, начиная с Stata 10. ivreg заменена на команду ivregress .

Таким образом, вместо этого лучше переключиться на ivregress: 

. ivregress 2sls rent pcturban i.region (hsngval = faminc), small

Instrumental variables (2SLS) regression

      Source |       SS       df       MS         Number of obs   =         50
-------------+------------------------------      F(  5,    44)   =       9.10
       Model |  12735.4667     5  2547.09334      Prob > F        =     0.0000
    Residual |  48507.6533    44  1102.44667      R-squared       =     0.2079
-------------+------------------------------      Adj R-squared   =     0.1179
       Total |    61243.12    49  1249.85959      Root MSE        =     33.203

------------------------------------------------------------------------------
        rent |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     hsngval |   .0038683   .0008958     4.32   0.000     .0020629    .0056737
    pcturban |  -.4980121   .5179779    -0.96   0.342    -1.541928    .5459039
             |
      region |
    N Cntrl  |   1.528672   15.14086     0.10   0.920    -28.98572    32.04306
      South  |    7.74279   15.10906     0.51   0.611    -22.70752     38.1931
       West  |  -40.61235   19.60999    -2.07   0.044    -80.13369   -1.091002
             |
       _cons |   88.26681   31.69154     2.79   0.008     24.39671    152.1369
------------------------------------------------------------------------------
Instrumented:  hsngval
Instruments:   pcturban 2.region 3.region 4.region faminc

Введите help ivregress из командной строки Stata для получения более подробной информации.

...