Non-significant variables (P > 0.1) were discarded
using a manual stepwise backward selection procedure
until reaching the Minimal Adequate Model.
Calculations were performed using the glmmPQL function in the MASS package for
R (Venables and Ripley, 2002).