Dear list members,

I am looking for an alternative to get the hat matrix  (or leverages) of a model fitted in lmer(), with fixed and random effects. 

The model looks like: example<-lmer(y~(x1)+(1|x2)+(1|x3)+(1|x4),weights=w)

I'm not sure, but I guess that by using the hat() and model.matrix()  functions I will  take in account only the fixed part of model. Is this true? I thought  about to use the matrices  in 'example@Zt' and 'example@X' to build the hat matrix, but I suppose may exist a more efficient alternative, since the design matrices are large. 
I'm new in the list and I apologise if this issue was discussed previously, but I did not find  anything related in the previous topics.


Thanks in advance,

Haroldo
	[[alternative HTML version deleted]]


