Statalist The Stata Listserver


[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]

[no subject]



Random-effects logistic regression              Number of obs      =       675
Group variable (i): e_farm                      Number of groups   =        28

Random effects u_i ~ Gaussian                   Obs per group: min =         1
                                                               avg =      24.1
                                                               max =       124

                                                Wald chi2(16)      =     55.51
Log likelihood  = -363.08026                    Prob > chi2        =    0.0000

------------------------------------------------------------------------------
      status |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
  _Iv_type_1 |   .8567931   .5292698     1.62   0.105    -.1805566    1.894143
  _Iv_type_3 |   .7836975   .5244423     1.49   0.135    -.2441906    1.811586
   floorarea |   .0005784   .0005917     0.98   0.328    -.0005813    .0017382
        bqs2 |    .522315    .322696     1.62   0.106    -.1101576    1.154788
        bqs5 |   .3837457   .2552641     1.50   0.133    -.1165628    .8840542
    _Ibqs7_1 |  -.5865056   .5794704    -1.01   0.311    -1.722247    .5492356
    _Ibqs7_2 |  -.1258125   .3440977    -0.37   0.715    -.8002316    .5486067
    _Ipps1_0 |  -.1305904   .4061407    -0.32   0.748    -.9266116    .6654309
    _Ipps1_1 |   .1964106   .3752239     0.52   0.601    -.5390148    .9318359
    _Ibis1_1 |   1.217941   1.285577     0.95   0.343    -1.301743    3.737625
    _Ibis1_2 |   .4460985   .3017643     1.48   0.139    -.1453485    1.037546
        bis5 |   .1678659    .266148     0.63   0.528    -.3537747    .6895064
        sas1 |   .9023956   .3671265     2.46   0.014     .1828409     1.62195
     _IrdB_2 |  -1.049914   .8792765    -1.19   0.232    -2.773265    .6734358
     _IrdB_3 |  -.4266696   .5833107    -0.73   0.464    -1.569938    .7165984
     _IrdB_4 |   .4683452   .3644601     1.29   0.199    -.2459835    1.182674
       _cons |  -3.388897   .7438128    -4.56   0.000    -4.846743   -1.931051
-------------+----------------------------------------------------------------
    /lnsig2u |        -15   419.6638                     -837.5258    807.5258
-------------+----------------------------------------------------------------
     sigma_u |   .0005531   .1160549                      1.4e-182    2.2e+175
         rho |   9.30e-08    .000039                             0           .
------------------------------------------------------------------------------
Likelihood-ratio test of rho=0: chibar2(01) =     0.00 Prob >= chibar2 = 1.000



**** Logit results *****************************************************

. xi: logit status i.v_type floorarea bqs2 bqs5 i.bqs7 i.pps1 i.bis1 bis5 sas1 i.rd
> B,  nolog
i.v_type          _Iv_type_1-3        (_Iv_type_2 for v_t~e==horizontal omitted)
i.bqs7            _Ibqs7_0-2          (naturally coded; _Ibqs7_0 omitted)
i.pps1            _Ipps1_0-2          (naturally coded; _Ipps1_2 omitted)
i.bis1            _Ibis1_0-2          (naturally coded; _Ibis1_0 omitted)
i.rdB             _IrdB_2-9           (naturally coded; _IrdB_9 omitted)

Logistic regression                               Number of obs   =        675
                                                  LR chi2(16)     =      78.07
                                                  Prob > chi2     =     0.0000
Log likelihood = -363.08026                       Pseudo R2       =     0.0971

------------------------------------------------------------------------------
      status |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
  _Iv_type_1 |   .8567931   .5292696     1.62   0.105    -.1805564    1.894143
  _Iv_type_3 |   .7836975   .5244422     1.49   0.135    -.2441902    1.811585
   floorarea |   .0005784   .0005917     0.98   0.328    -.0005813    .0017382
        bqs2 |    .522315   .3226957     1.62   0.106     -.110157    1.154787
        bqs5 |   .3837457   .2552639     1.50   0.133    -.1165624    .8840538
    _Ibqs7_1 |  -.5865056   .5794702    -1.01   0.311    -1.722246    .5492351
    _Ibqs7_2 |  -.1258125   .3440976    -0.37   0.715    -.8002315    .5486065
    _Ipps1_0 |  -.1305904   .4061405    -0.32   0.748    -.9266111    .6654304
    _Ipps1_1 |   .1964106   .3752237     0.52   0.601    -.5390143    .9318355
    _Ibis1_1 |   1.217941   1.285577     0.95   0.343    -1.301743    3.737625
    _Ibis1_2 |   .4460985    .301764     1.48   0.139     -.145348    1.037545
        bis5 |   .1678659   .2661478     0.63   0.528    -.3537743     .689506
        sas1 |   .9023956   .3671263     2.46   0.014     .1828413     1.62195
     _IrdB_2 |  -1.049914   .8792763    -1.19   0.232    -2.773264    .6734355
     _IrdB_3 |  -.4266696   .5833101    -0.73   0.464    -1.569936    .7165972
     _IrdB_4 |   .4683452   .3644599     1.29   0.199     -.245983    1.182673
       _cons |  -3.388897   .7438126    -4.56   0.000    -4.846743   -1.931051
------------------------------------------------------------------------------



************************************************************************
And from Stata 9.2, where I have used the -iter(3)- option to force
-xtlogit- to cease and desist, we obtain the same results with lnsig2u
further towards negative infinity (but not convergence):
************************************************************************

. xi: xtlogit status i.v_type floorarea bqs2 bqs5 i.bqs7 i.pps1 i.bis1 bis5 sas1 i.
> rdB,  nolog iter(3)
i.v_type          _Iv_type_1-3        (_Iv_type_2 for v_t~e==horizontal omitted)
i.bqs7            _Ibqs7_0-2          (naturally coded; _Ibqs7_0 omitted)
i.pps1            _Ipps1_0-2          (naturally coded; _Ipps1_2 omitted)
i.bis1            _Ibis1_0-2          (naturally coded; _Ibis1_0 omitted)
i.rdB             _IrdB_2-9           (naturally coded; _IrdB_9 omitted)

Random-effects logistic regression              Number of obs      =       675
Group variable (i): e_farm                      Number of groups   =        28

Random effects u_i ~ Gaussian                   Obs per group: min =         1
                                                               avg =      24.1
                                                               max =       124

                                                Wald chi2(16)      =     55.51
Log likelihood  = -363.08026                    Prob > chi2        =    0.0000

------------------------------------------------------------------------------
      status |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
  _Iv_type_1 |   .8567931   .5292696     1.62   0.105    -.1805564    1.894143
  _Iv_type_3 |   .7836975   .5244422     1.49   0.135    -.2441902    1.811585
   floorarea |   .0005784   .0005917     0.98   0.328    -.0005813    .0017382
        bqs2 |    .522315   .3226957     1.62   0.106     -.110157    1.154787
        bqs5 |   .3837457   .2552639     1.50   0.133    -.1165624    .8840538
    _Ibqs7_1 |  -.5865056   .5794702    -1.01   0.311    -1.722246    .5492351
    _Ibqs7_2 |  -.1258125   .3440976    -0.37   0.715    -.8002315    .5486065
    _Ipps1_0 |  -.1305904   .4061405    -0.32   0.748    -.9266111    .6654304
    _Ipps1_1 |   .1964106   .3752237     0.52   0.601    -.5390143    .9318355
    _Ibis1_1 |   1.217941   1.285577     0.95   0.343    -1.301743    3.737625
    _Ibis1_2 |   .4460985    .301764     1.48   0.139     -.145348    1.037545
        bis5 |   .1678659   .2661478     0.63   0.528    -.3537743     .689506
        sas1 |   .9023956   .3671263     2.46   0.014     .1828413     1.62195
     _IrdB_2 |  -1.049914   .8792763    -1.19   0.232    -2.773264    .6734355
     _IrdB_3 |  -.4266696   .5833101    -0.73   0.464    -1.569936    .7165972
     _IrdB_4 |   .4683452   .3644599     1.29   0.199     -.245983    1.182673
       _cons |  -3.388897   .7438126    -4.56   0.000    -4.846743   -1.931051
-------------+----------------------------------------------------------------
    /lnsig2u |  -52.57779          .                             .           .
-------------+----------------------------------------------------------------
     sigma_u |   3.83e-12          .                             .           .
         rho |   4.45e-24          .                             .           .
------------------------------------------------------------------------------
Likelihood-ratio test of rho=0: chibar2(01) =     0.00 Prob >= chibar2 = 1.000

*****************************************************************************

*
*   For searches and help try:
*   http://www.stata.com/support/faqs/res/findit.html
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/



© Copyright 1996–2014 StataCorp LP   |   Terms of use   |   Privacy   |   Contact us   |   What's new   |   Site index