Bookmark and Share

Notice: On April 23, 2014, Statalist moved from an email list to a forum, based at statalist.org.


[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

st: Boxtid r(430) flat region resulting in a missing likelihood


From   sela <[email protected]>
To   [email protected]
Subject   st: Boxtid r(430) flat region resulting in a missing likelihood
Date   Sun, 27 Mar 2011 06:51:53 -0700 (PDT)

Dear Stata Users,

I'm building a logistic regression modell and therefore am using the
-boxtid- command in order to test for linearity in the logit. I'm using the
following variables:
change: 0 , 1 (dependent variable)
retro: -2 until 2 (5 point scale, independent variable)
polint: 1 until 5 (9 point scale, independent variable)
use: -10 until 10 (21 point scale, independent variable)

when doing the following command
    boxtid logit change retro polint use [pw= cggesow]
the nonlinearity test for both retro and polint is significant.

This would be the output:


Logistic regression                               Number of obs   =       
710
                                                       Wald chi2(4)    =     
15.33
                                                       Prob > chi2     =    
0.0041
Log pseudolikelihood = -275.28524         Pseudo R2       =     0.0597

------------------------------------------------------------------------------
                  |               Robust
    change    |      Coef.   Std. Err.      z    P>|z|     [95% Conf.
Interval]
-------------+----------------------------------------------------------------
    Iretr__1 |   -.045204   .0565584    -0.80   0.424    -.1560563   
.0656484
    Iretr_p1 |   .0188079    .040473     0.46   0.642    -.0605176   
.0981334
    Iuse___1 |   2.302325   75.66987     0.03   0.976    -146.0079   
150.6125
    Iuse__p1 |   -.970768    17.1763    -0.06   0.955    -34.63569   
32.69415
       _cons |   -1.93799    .255246    -7.59   0.000    -2.438263  
-1.437717
------------------------------------------------------------------------------
retro    |  -.5394969   .2083624     -2.589   Nonlin. dev. 4.721   (P =
0.030)
      p1 |   3.535742   2.306073      1.533
------------------------------------------------------------------------------
polint   |  -.1599401   .2681029     -0.597   Nonlin. dev. 10.404  (P =
0.001)
      p1 |  -290.0341          .          .
------------------------------------------------------------------------------
use      | -.1236776   .1175078     -1.053   Nonlin. dev. 1.632   (P =
0.201)
      p1 |  -.2366999   2.137676     -0.111
------------------------------------------------------------------------------
Deviance:  539.095.



I then used exp(x) to change both polint and retro (so it's now polint_new =
exp(polint) and retro_new = exp(retro) ) but this caused the error message
"r(430) flat region resulting in a missing likelihood".

When doing the -boxtid- command only bivariate or also when excluding one
independent variable respectively, it does present results.

Is my idea of using exp(x) to transform the variables wrong? What remedy
would help me avoiding the flat region (and: what exactly does Stata mean by
flat region?)?

I hope someone wants to help me and can give me some advice since I've been
playing with those variables for quite sometime...

Thank you!

--
View this message in context: http://statalist.1588530.n2.nabble.com/Boxtid-r-430-flat-region-resulting-in-a-missing-likelihood-tp6212352p6212352.html
Sent from the Statalist mailing list archive at Nabble.com.
*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/


© Copyright 1996–2018 StataCorp LLC   |   Terms of use   |   Privacy   |   Contact us   |   Site index