[R-sig-ME] residuals in lme4

[lunina at access.unizh.ch]

Hello

I want to get the residuals of a heteroscedastic linear mixed-effects-model. 
How can I fit such one ond how is it possible to extract the residuals?

I saw that the the usual function "resid(.)" works for lmer, but not for 
lmer2. But no one of these two shows any change in the model values if I 
define some "weights". Is "weights" not yet implementet? And is there 
another possibility to calculate a weighted linear mixed-effects-model?

I already tried it with "lme(.)" of the package "nlme", but the models with 
"lme" fit much worse than those with "lmer", also if I compare the 
unweighted ones...

Thanks a lot

Monika