[R] weird (?) resid(.) from 2-level lme()

Dieter Menne dieter.menne at menne-biomed.de
Thu Oct 11 08:38:51 CEST 2007


Irene Mantzouni <ima <at> difres.dk> writes:

> I am trying to fit a 2-level hierarchical lme(). 
> I can extract ranef() and coef() without problems for levels=1:2. 
> However, when it comes to resid() at level=2, the resulting list has some
unnamed entries (label="NA"). I
> checked with my data (NAs have been omitted) and I found out that it gives 2-3
residuals less in each
> subgroup. 

Please post a reproducible example

Dieter



More information about the R-help mailing list