Hi all,
I am trying to perform an ordinal regression using proportional odds
model.
First I would like to test whether the assumption of proportional odds is
supported by the data.
I got the following error message when I try to fit one of the relevant
models. Can anyone give
some suggestions on how to fix this. Thank you.
Hanna
fit2 <- vglm(PLANHORIZON ~ factor(TotalCompanion) +
factor(NumTrips)+factor(Numnight)+factor(FMWEBSITE)+factor(OTHERWEB)+factor(BROCHU)
+ +factor(TRAAGENT) + factor(WOM) + factor(AAA) + factor(GUIDEBOOK) +
factor(ADVERTISE) + factor(NEWSPAPER) + factor(INFOOTHER) +
+ ExpTotal + factor(AGE) + factor(CHILD) + factor(INCOME), +
family=cumulative(parallel=F), data=two)
Error in if (any(smallmu <- (mu * (1 - mu) < double.eps))) { :
missing value where TRUE/FALSE needed
In addition: Warning messages:
1: In Deviance.categorical.data.vgam(mu = mu, y = y, w = w, residuals =
residuals, :
fitted values close to 0 or 1
2: In Deviance.categorical.data.vgam(mu = mu, y = y, w = w, residuals =
residuals, :
fitted values close to 0 or 1
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org 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.
