# 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) > mod_FODs_optimizer <- glmer(y ~ x1 * x2 + (1|subject_FODs), family = binomial, data = indvar_FODs, + control = glmerControl(optimizer = "bobyqa")) boundary (singular) fit: see help('isSingular') add optCtrl according to the data. improve this test write the edited version of the code

# 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)