Skip to content
Prev 4843 / 398506 Next

polr question

On Sat, 11 Mar 2000, Troels Ring wrote:

            
No, but there are at least three conventions here.  The simplest way to
interpret a polr model is as a discretized linear model with a logistic
error distribution, and that is what polr does.  So think of the
coefficients as being those of the variation with x of the (mean, mode) of
the underlying logistic, and the intercepts as the breakpoints for the
discretization.
(Um, it's MASS, and you are only getting away with this on Windows.
Don't trust what is displayed in Explorer by default.)
lower categories is higher ?

I don't see any mention of a carrier here, nor of x?  But if you had
the table

y=    1   2   3
x=0   19  29  24
x=1  497 560 269

I can reproduce your results, and the larger x does tend to give a lower
response.  (That does seem the natural sign for the x coefficient to me.)

Note that glm will give the same sign (apart from the intercepts) as polr:
Call:  glm(formula = yy ~ x, family = binomial, weights = freq) 

Coefficients:
(Intercept)            x  
     1.0258      -0.5142
Call:  glm(formula = (yy > 2) ~ x, family = binomial, weights = freq) 

Coefficients:
(Intercept)            x  
    -0.6931      -0.6753  

The sign of the intercepts are different, as you need to subtract the
breakpoint to get a zero-breakpoint model.
[...]

I think that has one of the other conventions, which are logistic models
for P(Y <= k | x)  or P(Y > k | x).