[R] residuals in lmer

[Liaw, Andy]

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