Hello, I am looking for a way to extract fitted values of a 2-level generalized linear mixed model using lme4-package. As it is a Poisson-model, I have to use lmer(.), e.g.: (y ~ x1 + x2 + (1|PtID), data = mydata, family = "poisson", method = "ML") If it was a linear mixed model with normally distributed response, that would be very easy, e.g.: ### 2-level LME fm1 <- lme(pixel ~ day + I(day^2), data = Pixel, random = list(Dog = ~ day, Side = ~ 1)) fitted(fm1, level = 0:2) plot(augPred(fm1, level = c(0,1,2))) How can I extract the fitted values for different levels so that I can plot it afterwards ? Is there something similar to augPred in lmer() ??? Thanks in advance for your answers, Christine Adrion. -- View this message in context: http://www.nabble.com/lmer%28.%29%3A-extract-fitted-values-tf4653387.html#a13295162 Sent from the R help mailing list archive at Nabble.com.