Dear all,

I have constructed a generalized linear mixed model using "lmer" with continuous covariates, a binary response of zero or one (i.e. logistic) and have included a single random effect for each individual (IND). I would like to plot the residuals but need to distinguish the residuals based their associated random effect factor. In other words I need to relate each residual to each random effect factor for display purposes. 

>model<-lmer(response~cov1+cov2.....+(1|IND),data,family=binomial(link="logit"))


I have extracted the residuals using: 
> resids<-residuals(model)

but I now need a way identify which residuals are linked to each factor of the random effect.

I tried to using "cbind" to join the residuals back the original data set but the for some reason the vectors are of different length.

Any guidance would be greatly appreciated.

Thanks,
 
  
Ben Wilson
	[[alternative HTML version deleted]]

