The per-observation residuals are returned, i.e., the difference of the observation and the fitted value including random effects. With type one can specify whether the weights should be used or not.
## S3 method for class 'rlmerMod'residuals(object, type = c("response","weighted"), scaled =FALSE,...)
Arguments
object: rlmerMod object
type: type of residuals
scaled: scale residuals by residual standard deviation (=scale parameter)?
...: ignored
Examples
## Not run: fm <- rlmer(Yield ~(1|Batch), Dyestuff) stopifnot(all.equal(resid(fm, type="weighted"), resid(fm)* getME(fm,"w_e")))## End(Not run)