I'm trying to obtain the robust standard errors for a multinomial ordered logit
model:
mod6 <- lrm(wdlshea ~ initdesch + concap + capasst + qualrat +
terrain,data=full2)
The model is fine but when I try to get the RSE I get an error.
coeftest(mod6, vcov = vcovHAC(mod6))
Error in match.arg(type) :
'arg' should be one of ordinary, score, score.binary, pearson,
deviance, pseudo.dep, partial, ........etc.
I'm a novice R user and am not sure how to address this problem. I have also
tried to use alternatives (zelig, polr) but have had no luck. Any assistance
on generating RSE for a multinomial order logit model would be appreciated
Thanks,
Patrick
_________________________________________________________________
Hotmail: Trusted email with powerful SPAM protection.
[[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.
