Greice,

Você pode obter esses resíduos fazendo a diferença entre os valores observados (y) e os preditos com efeitos fixos (y_x) e os preditos com efeitos fixos e aleatórios (y_xz), matricialmente

r_cond  = y - X\beta - Zb
r_marg = y - X\beta

X é a matriz relativa aos efeitos fixos, Z é a matriz relativa aos efeitos aleatórios, \beta é o vetor de estimativas dos efeitos fixos e b é o vetor de efeitos aleatórios preditos. V eja como extrair as matrizes em http://r-br.2285057.n4.nabble.com/R-br-Matrix-Modelo-td3574029.html#a3574331.

À disposição.
Walmes.

==========================================================================
Walmes Marques Zeviani
LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
Departamento de Estatística - Universidade Federal do Paraná
fone: (+55) 41 3361 3573
VoIP: (3361 3600) 1053 1173
e-mail: walmes@ufpr.br
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================