Re: [R] bugs and misfeatures in polr(MASS).... fixed!

2014-08-12 Thread tjb
The official maintainers were dismissive when I suggested there were some problems I could fix with the then implementation of polr. I haven't looked at it since, sorry. On Tue, Aug 12, 2014 at 7:44 PM, Guido Biele [via R] < ml-node+s789695n4695392...@n4.nabble.com> wrote: > I modified (where ne

Re: [R] bugs and misfeatures in polr(MASS).... fixed!

2012-10-25 Thread tjb
This could well be out of date because I have not paid any attention to the official POLR code in a year or more. Attached is fixed-polr.R. cheers, Tim On Wed, Oct 24, 2012 at 10:38 PM, ahs [via R] < ml-node+s789695n4647311...@n4.nabble.com> wrote: > Great! > You can skip my question about s0 th

Re: [R] bugs and misfeatures in polr(MASS).... fixed!

2012-10-24 Thread tjb
;t seem to get the change of code right. I redefine the function polr by > the lines > > tjb wrote > if(missing(start)) { > # try something that should always work -tjb > u <- as.integer(table(y)) > u <- (cumsum(u)/sum(u))[1:q] > zetas &

Re: [R] likelyhood maximization problem with polr

2010-11-09 Thread tjb
blackscorpio wrote: > > Thank you for your answer. I have already tried lrm and it's true that it > works better than polr in such a case. Nevertheless lrm does not work with > the addterm and dropterm functions (to my knowledge) and I need to use > them. Maybe do you know alternate functions th

Re: [R] bugs and misfeatures in polr(MASS).... fixed!

2010-11-06 Thread tjb
Note that that the enhancements in my original post solve the unresolved problem of Chaehyung Ahn (22 Mar 2005) whose data I reproduce: y,x,lx 0,3.2e-02,-1.49485 0,3.2e-02,-1.49485 0,1.0e-01,-1.0 0,1.0e-01,-1.0 0,3.2e-01,-0.49485 0,3.2e-01,-0.49485 1,1.0e+00,0.0 0,1.0e+00,0.0 1,3.

Re: [R] bugs and misfeatures in polr(MASS).... fixed!

2010-11-06 Thread tjb
hat should always work -tjb u <- as.integer(table(y)) u <- (cumsum(u)/sum(u))[1:q] zetas <- switch(method, "logistic"= qlogis(u), "probit"= qnorm(u), "cauchit"= qcauchy(u),

Re: [R] likelyhood maximization problem with polr

2010-11-06 Thread tjb
blackscorpio wrote: > > Dear community, > > I am currently trying to fit an ordinal logistic regression model with the > polr function. I often get the same error message : > > "attempt to find suitable starting values failed", for example with : > ... > Does anyone have a clue ? > Yes. The