[R] likelyhood maximization problem with polr

Peng, C cpeng.usm at gmail.com
Tue Sep 7 07:19:57 CEST 2010


Since the default initial value is not good enough. You should choose one
based on your experience or luck. I choose start=rep(1,5) since there are
parameters in the model.

> polr(Species~Sepal.Length+Sepal.Width+Petal.Length+Petal.Width,iris,
> start=rep(1,6), method = "logistic") 
Call:
polr(formula = Species ~ Sepal.Length + Sepal.Width + Petal.Length + 
    Petal.Width, data = iris, start = rep(1, 6), method = "logistic")

Coefficients:
Sepal.Length  Sepal.Width Petal.Length  Petal.Width 
   -2.471271    -6.675614     9.432753    18.266280 

Intercepts:
   setosa|versicolor versicolor|virginica 
            5.513223            42.598774 

Residual Deviance: 11.89856 
AIC: 23.89856 


-- 
View this message in context: http://r.789695.n4.nabble.com/likelyhood-maximization-problem-with-polr-tp2528818p2529174.html
Sent from the R help mailing list archive at Nabble.com.



More information about the R-help mailing list