qqnorm() is an S3 generic without a method for lmer. It doesn't even have methods for other models based on S3:
methods("qqnorm")
[1] qqnorm.default So why do you think qqnorm(fittedModel) will work? You can do qqnorm(residuals(crop.lme)) Andy
From: Jacob Michaelson Does anyone know how to extract residuals in lmer? Here's the error I get:
>
crop.lme=lmer(response~variety*irrigation*pesticide+(1|rep)+(1|rep: pesticide)+(1|rep:pesticide:irrigation), crop.data)
> qqnorm(crop.lme)
Error in qqnorm.default(crop.lme) : y is empty or has only NAs
> resid(crop.lme)
NULL Thanks! --Jake
______________________________________________ R-help at stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html