Dear R-Users: I am trying to find the robust (or sandwich) estimates of the standard error of fixed effects parameter estimates using the package "lmer2". In model-1, I used "robust=TRUE" on the other, in model-2, I used "robust=FALSE". Both models giving me the same estimates. So my question is, does the robust option works in lmer2 to get the robust estimates of the standard error? If anybody could offer me a suggestion I would greatly appreciate it. Thank you. Model-1:> p.mle<-lmer2(ddimer~race+steroid+psi+sofa+apache + (apache|subject), method="ML", data=final, robust=TRUE, cluster="id", weights=final$w) > beta=fixef(p.mle) > Vcov=vcov(p.mle, useScale=FALSE) > se=sqrt(diag(Vcov)) > beta(Intercept) race steroid psi sofa apache 5.826489820 -0.001920670 -0.242040171 0.005293996 0.075468340 0.009245152> se[1] 0.108325229 0.058921371 0.055975547 0.001285687 0.018119089 0.002559902 Model-2:> p.mle<-lmer2(ddimer~race+steroid+psi+sofa+apache + (apache|subject), method="ML", data=final, robust=FALSE, cluster="id", weights=final$w) > beta=fixef(p.mle) > Vcov=vcov(p.mle, useScale=FALSE) > se=sqrt(diag(Vcov)) > beta(Intercept) race steroid psi sofa apache 5.826489820 -0.001920670 -0.242040171 0.005293996 0.075468340 0.009245152> se[1] 0.108325229 0.058921371 0.055975547 0.001285687 0.018119089 0.002559902 Best Regards, Sattar ____________________________________________________________________________________ , and more! [[alternative HTML version deleted]]
Abdus Sattar said the following on 9/19/2007 7:03 AM:> Dear R-Users: > > I am trying to find the robust (or sandwich) estimates of the standard error of fixed effects parameter estimates using the package "lmer2". In model-1, I used "robust=TRUE" on the other, in model-2, I used "robust=FALSE". Both models giving me the same estimates. So my question is, does the robust option works in lmer2 to get the robust estimates of the standard error? If anybody could offer me a suggestion I would greatly appreciate it. Thank you. > > Model-1: > >> p.mle<-lmer2(ddimer~race+steroid+psi+sofa+apache + (apache|subject), method="ML", data=final, robust=TRUE, cluster="id", weights=final$w) >> beta=fixef(p.mle) >> Vcov=vcov(p.mle, useScale=FALSE) >> se=sqrt(diag(Vcov)) >> beta > (Intercept) race steroid psi sofa apache > 5.826489820 -0.001920670 -0.242040171 0.005293996 0.075468340 0.009245152 >> se > [1] 0.108325229 0.058921371 0.055975547 0.001285687 0.018119089 0.002559902 > > Model-2: > >> p.mle<-lmer2(ddimer~race+steroid+psi+sofa+apache + (apache|subject), method="ML", data=final, robust=FALSE, cluster="id", weights=final$w) >> beta=fixef(p.mle) >> Vcov=vcov(p.mle, useScale=FALSE) >> se=sqrt(diag(Vcov)) >> beta > (Intercept) race steroid psi sofa apache > 5.826489820 -0.001920670 -0.242040171 0.005293996 0.075468340 0.009245152 >> se > [1] 0.108325229 0.058921371 0.055975547 0.001285687 0.018119089 0.002559902 > > > Best Regards, > > Sattar > >The help page to ?lmer2 in the lme4 package makes no mention of "cluster" or "robust" arguments. To me, that would mean these arguments are ignored. HTH, --sundar
Hi Doran, My interests to estimate fixed effects parameters in a mixed model (not the random effects parameters) and, in addition, in the model there are some nuisace parameters. For this estimation I am using pseudo-maximum likelihood methods. First I am estimating nuisance parameters then using these estimates in my likelihood function to get the estimates of fixed effects parameters using lmer2. Therefore, to get the correct estimates of the standard errors, I need robust (or sandwich) estiamtes of the SE. Thank you and thanks to Thomas his remarks, Sattar ----- Original Message ---- From: Thomas Lumley <tlumley@u.washington.edu> To: "Doran, Harold" <HDoran@air.org> Cc: Sundar Dorai-Raj <sundar.dorai-raj@pdf.com>; Abdus Sattar <upsattar@yahoo.com>; R-help@stat.math.ethz.ch; gsrwork@yahoo.com Sent: Wednesday, September 19, 2007 1:13:20 PM Subject: Re: [R] Robust or Sandwich estimates in lmer2 On Wed, 19 Sep 2007, Doran, Harold wrote:> This has come up before and I'll again ask the question "why would you > want robust standard errors in lmer"?And I'll again answer: using lmer() does not automatically guarantee correct model specification, either for the correlation structure or for the marginal variance. -thomas Thomas Lumley Assoc. Professor, Biostatistics tlumley@u.washington.edu University of Washington, Seattle ____________________________________________________________________________________ ttp://surveylink.yahoo.com/gmrs/yahoo_panel_invite.asp?a=7 [[alternative HTML version deleted]]
Abdus Sattar <upsattar at yahoo.com> asked:> > I am trying to find the robust (or sandwich) estimates of the standard > error of fixed effects parameter estimates using the package "lmer2". >Others have already pointed out that this is not implemented in lmer2. You could try a delete-n jackknife, which you would have to implement yourself. An ordinary delete-1 jackknife does not work for clustered/correlated data. This would also give you bias-corrected point estimates. I should add that, in the limited simulations I have done (for Gaussian mixed model analysis of pedigree data), the jackknife standard errors seemed a bit conservative (too big). David Duffy. -- | David Duffy (MBBS PhD) ,-_|\ | email: davidD at qimr.edu.au ph: INT+61+7+3362-0217 fax: -0101 / * | Epidemiology Unit, Queensland Institute of Medical Research \_,-._/ | 300 Herston Rd, Brisbane, Queensland 4029, Australia GPG 4D0B994A v
Possibly Parallel Threads
- Extracting df (degree of freedom) & estfun (estimating function) from model built in lmer or lmer2
- estfun & df
- lmer2 error under Mac OS X on PowerPC G5 but not on Dual-Core Intel Xeon
- How to get p-values, seperate vectors of regression coefficients and their s.e. from the "yags" output?
- Library & Package for Tobit regression