[R-sig-ME] BLUP in lmer

[Roberts, J. Kyle]

Olivier,

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,
Kyle

***************************************
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 ?

Thanks,
olivier

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