R help - Mar 2010 - Robust SE for lrm object

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

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