[R-sig-ME] response residuals from a logistic mixed regression model using lmer

vahl wouter.vahl at helsinki.fi
Tue Mar 10 16:50:19 CET 2009


Using the lmer function, I have fitted the following logistic mixed 
regression model on an experimental data set containing three fixed 
factors (A, B & C) and two random variables (day & cage):

model1 <- lmer(cbind(k1, k2) ~A*B*C + (1|day) + 

where k1 is the number of successes and k2 the number of failures in a 

Presenting all data in the same figure turns out to result in a graph 
that is hard to read. Therefore, I am now planning to present my data in 
two graphs, one containing information on A and B, the other containing 
information on B & C.

In presenting my data regarding factor A and B visually, I would like to 
take out the (experimentally imposed) variation explained by the factor 
C and the two random block factors. To do so, I need the raw, response 
residuals, that is the residuals on the logit scale. I have studied the 
help-file of the lmer function as well as that of the mer class, and I 
have read the question posted by Andy Fugard on this list as well as the 
reply of Douglas Bates (both posted at 14 march 2008), but am not able 
to find the response residuals. Any advice would be kindly appreciated.

Wouter Vahl

More information about the R-sig-mixed-models mailing list