model1 <- lmer(categorical_FOD_FODs~Condition_FODs*Gender_FODs + (1|subject_FODs), data=indvar_FODs, control = lmerControl( optimizer ='bobyqa', optCtrl=list(method='nlminb', maxfun=6000000), calc.derivs = FALSE)) Error in mkRespMod(fr, REML = REMLpass) : response must be numeric Please edit the code according to the categorical dependent variable which is categorical_FOD_FODs

function that: replace 0 with NA when there are missing data points