On Mar 5, 2010, at 11:54 PM, Patrick Shea wrote:
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
Have you loaded the library that contains the vcovHAC function?
(And do you know whether coeftest works with Design/rms objects?)
--
David Winsemius, MD
West Hartford, CT
______________________________________________
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.
