# dependent variable y <- as.factor(indvar_FODs$categorical_FOD_FODs) # independent variables x1 <- indvar_FODs$Condition_FODs x2 <- indvar_FODs$Language_used_FODs # model with optimizer argument mod_FODs <- glmer(y ~ x1 * x2 + (1|subject_FODs), family = binomial, data = indvar_FODs, control = glmerControl(optimizer = "bobyqa")) # print model summary(mod_FODs) include optimizer because of convergence error
options(contrasts = c("contr.sum", "contr.poly"))