I want to use quadratic terms to fit my general linear mixed model with id as a random effect, using the lme4 package. It’s about how the distance to settlements influences the probability of occurrence of an animal. I use the following code (I hope it is correct): glmer_dissettl <- glmer(case ~ poly(dist_settlements,2) + (1|id),… Read More How to interpret the output of R's glmer() with quadratic terms