Sorry, I have been away. I'll have a look at this and get back to you
tomorrow.

Tim

On Thu, Oct 18, 2012 at 10:35 PM, ahs [via R] <
ml-node+s789695n4646600h91@n4.nabble.com> wrote:

> Hello,
> I am trying to use this fix for the convergence problem in polr, but I
> don'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 <-
>          switch(method,
>                 "logistic"= qlogis(u),
>                 "probit"=   qnorm(u),
>                 "cauchit"=  qcauchy(u),
>                 "cloglog"=  -log(-log(u)) )
>       s0 <- c(rep(0,pc),zetas[1],log(diff(zetas)))
>
> but - which part should I erase? The whole old if(missing(start)) {...} or
> something else or nothing? I have tried all ways I could think of.
>
> I also can't see where the new object s0 is used in the old code?
>
> Any input on my problem would be very much appreciated!
>
>
> ------------------------------
>  If you reply to this email, your message will be added to the discussion
> below:
>
> http://r.789695.n4.nabble.com/bugs-and-misfeatures-in-polr-MASS-fixed-tp3024677p4646600.html
>  To unsubscribe from bugs and misfeatures in polr(MASS).... fixed!, click
> here<http://r.789695.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=3024677&code=dGltb3RoeS5iZW5oYW1AdXFjb25uZWN0LmVkdS5hdXwzMDI0Njc3fDE5NTE2NDMxMjk=>
> .
> NAML<http://r.789695.n4.nabble.com/template/NamlServlet.jtp?macro=macro_viewer&id=instant_html%21nabble%3Aemail.naml&base=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespace&breadcrumbs=notify_subscribers%21nabble%3Aemail.naml-instant_emails%21nabble%3Aemail.naml-send_instant_email%21nabble%3Aemail.naml>
>




-----
Tim J. Benham
--
View this message in context: http://r.789695.n4.nabble.com/bugs-and-misfeatures-in-polr-MASS-fixed-tp3024677p4647307.html
Sent from the R help mailing list archive at Nabble.com.
	[[alternative HTML version deleted]]


