Re: [R] Problem with ordered probit model in MASS
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
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
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
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.