
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]]

