[R] help with summary(polr_model)
Luca Braglia
lbraglia at gmail.com
Thu Mar 13 23:52:43 CET 2008
hello everybody
I'm a newbie with ordered probit and with polr too.
The problem is that I have a dependent variable I need to explain
with an ordered probit that is
> head(dfscale$sod.sit.ec.fam,100)
> [1] 5 7 5 6 5 5 6 8 6 8 8 8 6 6 6 5 0 5 NA 6
> [21] 7 NA NA 0 0 2 5 3 6 6 7 6 NA 8 6 6 7 NA NA NA
> [41] 4 5 NA 10 3 9 10 10 7 5 5 5 NA NA 5 7 1 1 7 3
> [61] 6 7 NA 4 6 NA 3 4 5 3 4 4 4 5 NA 2 4 4 8 5
> [81] 5 4 6 6 5 6 7 7 7 7 7 6 8 7 6 5 5 5 6 7
so when I estimate the model
> base.op <- polr(
> as.ordered(sod.sit.ec.fam) ~
> log(y) +
> log(ncomp) +
> eta +
> I(eta^2) +
> uomo +
> elem +
> laurea +
> saluteno +
> extra
> ,
> weights = dfscale$w,
> data = dfscale,
> method = "probit"
> )
I get
> Warning in eval(expr, envir, enclos) :
> non-integer #successes in a binomial glm
Then when I try to summarize the model, in order to use xtable
> summary(base.op)
> Re-fitting to get Hessian
> Errore in optim(start, fmin, gmin, method = "BFGS", hessian =
> Hess, ...) :
> valore iniziale in 'vmmin' non finito
that freely translated to english means "starting value in
'vmmin' not finished (or finite, more likely)"
what does this mean? How can I have my summary?
Thank you
Luca
--
Luca Braglia, aka Bragliozzo
http://bragliozzo.altervista.org
More information about the R-help
mailing list