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

[Dieter Menne]

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