Re: [R] standardized random effects with ranef.lme()

From: Spencer Graves <spencer.graves_at_pdf.com>
Date: Mon 31 Jul 2006 - 00:31:55 EST

          Have you tried RSiteSearch("MLWin")? I just got 29 hits. I wonder if any one of these might relate to your question?

          If you would like more help on this issue from this listserve, please submit another post, preferably illustrating your question with the simplest possible self-contained example that illustrates your question, perhaps like the following:

	  fm1.16 <- lme(distance~age, data=Orthodont[1:16,],
            random=~age|Subject)

	  Hope this helps.
	  Spencer Graves

p.s. PLEASE do read the posting guide
"www.R-project.org/posting-guide.html" and provide commented, minimal, self-contained, reproducible code.

Dirk Enzmann wrote:
> Using ranef() (package nlme, version 3.1-75) with an 'lme' object I can
> obtain random effects for intercept and slope of a certain level (say:
> 1) - this corresponds to (say level 1) "residuals" in MLWin. Maybe I'm
> mistaken here, but the results are identical.
>
> However, if I try to get the standardized random effects adding the
> paramter "standard=T" to the specification of ranef(), the results
> differ considerably from the results of MLWin (although MLWin defines
> "standardized" in the same way as "divided by its estimated (diagnostic)
> standard error").
>
> Why do the results differ although the estimates (random effects and
> thus their variances) are almost identical? I noticed that lme() does
> not compute the standard errors of the variances of the random effects -
> for several reasons, but if this is true, how does ranef() calculate the
> standardized random effects (the help says: '"standardized" (i.e.
> divided by the corresponding estimated standard error)').
>
> Is there a way to obtain similar results as in MLWin (or: should I
> prefer the results of ranef() for certain reasons)?
>
> Dirk
>
> -----------------------------
> R version: 2.3.1 Patched (2006-06-21 r38367)
>
>
> *************************************************
> Dr. Dirk Enzmann
> Institute of Criminal Sciences
> Dept. of Criminology
> Edmund-Siemers-Allee 1
> D-20146 Hamburg
> Germany
>
> phone: +49-(0)40-42838.7498 (office)
> +49-(0)40-42838.4591 (Billon)
> fax: +49-(0)40-42838.2344
> email: dirk.enzmann@uni-hamburg.de
> www:
> http://www2.jura.uni-hamburg.de/instkrim/kriminologie/Mitarbeiter/Enzmann/Enzmann.html
>
> ______________________________________________
> R-help@stat.math.ethz.ch 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.



R-help@stat.math.ethz.ch 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. Received on Mon Jul 31 00:43:16 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Tue 01 Aug 2006 - 04:17:06 EST.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.