Dear R users,
I am trying to fit a mixed effects model using the lmer command.

fit1<-lmer(Rate~T+offset(log(Area))+(1|S/T), family=poisson)

The rate is count data, T is treatment a fixed effect, S is site a random 
effect and the interaction S*T a random effect. 

I have been looking at the posts but I do not understand how you check
the fit of the model. The commands resid, qqnorm etc do not seem to work.
Thanks is advance
Beth Strain


