> table(s_complete$cat_SNB_new) 0 1 2148 1360 > model <- lmer(cat_SNB_new ~ Condition_SNB * Dotcounting_SNB + (1 | subject_SNB), data = na.omit(s_complete)) Error in mkRespMod(fr, REML = REMLpass) : response must be numeric

>  model <- lmer(SNB_new ~ Condition_SNB * Dotcounting_SNB + (1 | subject_SNB), data = na.omit(s_complete))
> summary(model)
Linear mixed model fit by REML ['lmerMod']
Formula: SNB_new ~ Condition_SNB * Dotcounting_SNB + (1 | subject_SNB)
   Data: na.omit(s_complete)

REML criterion at convergence: 628.8

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-3.3291 -0.8359 -0.2435  0.7277  2.9290 

Random effects:
 Groups   Name        Variance Std.Dev.
 subject_SNB (Intercept) 0.2609  0.5107  
 Residual               1.4618  1.2081  
Number of obs: 3506, groups:  subject_SNB, 209

Fixed effects:
                                  Estimate Std. Error