[R-sig-ME] BLUP in lmer

Roberts, J. Kyle jkrobert at bcm.tmc.edu
Fri Mar 23 17:30:15 CET 2007


Do you mean something like this?
fm1 <- lmer(Reaction ~ Days + (Days|Subject), sleepstudy)
with(fm1, xyplot(resid(.) ~ fitted(.))) 

This gives you a plot of the residual versus the fitted.  "sleepstudy" is included in the package.

Hope this helps,

J. Kyle Roberts, Ph.D.
Baylor College of Medicine
Center for Educational Outreach
One Baylor Plaza, MS:  BCM411
Houston, TX   77030-3411
713-798-6672 - 713-798-8201 Fax
jkrobert at bcm.edu

-----Original Message-----
From: r-sig-mixed-models-bounces at r-project.org [mailto:r-sig-mixed-models-bounces at r-project.org] On Behalf Of Olivier MARTIN
Sent: Friday, March 23, 2007 10:30 AM
To: r-sig-mixed-models at r-project.org
Subject: [R-sig-ME] BLUP in lmer

Hi all,

I am using 'lmer' to fit generalized linear mixed-effects models. 
I would like to know if there is a function to estimate the random effects.
And, is there a way to compare the observed values vs. fitted values or fitted values vs. residuals ?


R-sig-mixed-models at r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models

More information about the R-sig-mixed-models mailing list