Re: [R] Follow-up about ordinal logit with mixtures: how about 'continuation ratio' strategy?

From: Frank E Harrell Jr <>
Date: Thu, 10 May 2007 07:38:54 -0500

Paul Johnson wrote:
> This is a follow up to the message I posted 3 days ago about how to
> estimate mixed ordinal logit models. I hope you don't mind that I am
> just pasting in the code and comments from an R file for your
> feedback. Actual estimates are at the end of the post.

. . .


lrm does not give an incorrect sign on the intercepts. Just look at how it states the model in terms of Prob(Y>=j) so that its coefficients are consistent with the way people state binary models.

I'm not clear on your generation of simulated data. I specify the population logit, anti-logit that, and generate binary responses with those probabilities. I don't use rlogis.

See if using the PO model with lrm with penalization on the factor does what you need.

lrm is not set up to omit an intercept with the -1 notation.

My book goes into details about the continuation ratio model.

Frank Harrell mailing list PLEASE do read the posting guide and provide commented, minimal, self-contained, reproducible code. Received on Thu 10 May 2007 - 12:45:07 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Fri 11 May 2007 - 18:31:14 GMT.

Mailing list information is available at Please read the posting guide before posting to the list.