Re: [R] Robust SE for lrm object

2010-03-06 Thread Frank E Harrell Jr
Achim Zeileis wrote: On Sat, 6 Mar 2010, David Winsemius wrote: 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) T

Re: [R] Robust SE for lrm object

2010-03-05 Thread Achim Zeileis
On Sat, 6 Mar 2010, David Winsemius wrote: 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 w

Re: [R] Robust SE for lrm object

2010-03-05 Thread David Winsemius
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.

[R] Robust SE for lrm object

2010-03-05 Thread Patrick Shea
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(