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

2007-01-25 Thread David Duffy
Justin Bem wrote:

> I got this message, while using the polr function in MASS
> 
> Error in polr(EQ, don, subset = (cote != 0), method = "probit", na.action = 
> na.omit) :
> attempt for find suitable starting values failed
> 
> how can I initialise starting values ?

You might try running an ordered logistic in lrm() and using those estimates
as starting values (if it runs OK).

David Duffy.

__
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.


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

2007-01-24 Thread Prof Brian Ripley

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,  [EMAIL PROTECTED]
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, UKFax:  +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.


[R] Problem with ordered probit model in MASS

2007-01-24 Thread justin bem
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 ?
 
Justin BEM
Elève Ingénieur Statisticien Economiste
BP 294 Yaoundé.
Tél (00237)9597295.










___ 
Découvrez une nouvelle façon d'obtenir des réponses à toutes vos questions ! 
Profitez des connaissances, des opinions et des expériences des internautes sur 
Yahoo! Questions/Réponses 

[[alternative HTML version deleted]]

__
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.


[R] Problem with ordered probit model in MASS

2007-01-23 Thread justin bem
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 ?
 
Justin BEM
Elève Ingénieur Statisticien Economiste
BP 294 Yaoundé.
Tél (00237)9597295.








___ 
Découvrez une nouvelle façon d'obtenir des réponses à toutes vos questions ! 
Profitez des connaissances, des opinions et des expériences des internautes sur 
Yahoo! Questions/Réponses 

[[alternative HTML version deleted]]

__
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.