Re: [R] Problem with ordered probit model in MASS

From: Prof Brian Ripley <ripley_at_stats.ox.ac.uk>
Date: Wed 24 Jan 2007 - 14:55:01 GMT

On Wed, 24 Jan 2007, justin bem wrote:

(Again: this is a duplicate post.)

> Dear all,
>
> I got this message, while using the polr function in MASS
>
>> EQ<-as.formula("dep~fpta+tcdv+cdta+cmta+prcd+patc+lactifs+excta")
>> Estim<-polr(EQ,don, subset=(cote!=0),method="probit",na.action=na.omit)
>
> Error in polr(EQ, don, subset = (cote != 0), method = "probit", na.action = na.omit) :
> attempt for find suitable starting values failed
> In addition: Warning messages:
> 1: algorithm did not converge in: glm.fit(X, y1, wt, family = binomial("probit"), offset = offset)
> 2: fitted probabilities numerically 0 or 1 occurred in: glm.fit(X, y1, wt, family = binomial("probit"), offset = offset)
>
> how can I initialise starting values ?

What do you think the 'start' argument to polr() is for?

If you are asking how you find suitable values, I cannot help you as I know nothing about your problem, and failing to find starting values usually means that the model is very far from appropriate.

> Justin BEM
> Elève Ingénieur Statisticien Economiste
> BP 294 Yaoundé.
> Tél (00237)9597295.

-- 
Brian D. Ripley,                  ripley@stats.ox.ac.uk
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272866 (PA)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595


______________________________________________ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

Received on Thu Jan 25 16:35:16 2007

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Thu 25 Jan 2007 - 06:30:29 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.