[R-sig-ME] extracting residuals in lmer glmm

Martin Henry H. Stevens HStevens at muohio.edu
Thu Aug 14 15:32:05 CEST 2008


Hi folks,
I seem unable to figure out how to extract residuals from a Poisson  
lmer model.

And, I am not sure what to do with the working residuals.

Any references or ideas would be appreciated.

# Example:
y <- rpois(100, lambda=10)
xz <- expand.grid(x=gl(2,5), z=gl(10,1))
m1 <- lmer(y ~ x + (1|z), xz, family="poisson")

resid(m1)
Error: 'resid' is not implemented yet

  residuals(m1)
Error: 'residuals' is not implemented yet

  summary(m1 @ wrkres)
Min. 1st Qu.  Median    Mean 3rd Qu.    Max.
1.38    1.98    2.20    2.22    2.50    3.13


Hank Stevens




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