mm_glmm_3 <- glmer(cbind(MissedMeds_N, TotalAdministrations) ~ month_id + factor(CareHomeSize) + factor(Ratings) + (1|FacilityKey),
                   family = "binomial", data = facs_3mth, na.action = "na.omit", nAGQ = 0)

> summary(mm_glmm_3) # negative month_id, significant
Generalized linear mixed model fit by maximum likelihood (Adaptive Gauss-Hermite Quadrature, nAGQ = 0) ['glmerMod']
 Family: binomial  ( logit )
Formula: cbind(MissedMeds_N, TotalAdministrations) ~ month_id + factor(CareHomeSize) +      factor(Ratings) + (1 | FacilityKey)
   Data: facs_3mth

     AIC      BIC   logLik deviance df.resid
 21058.9  21095.1 -10522.5  21044.9     1288

Scaled residuals:
     Min       1Q   Median       3Q      Max
-16.4066  -1.0167  -0.0825   0.8157  29.3459

Random effects:
 Groups      Name        Variance Std.Dev.
 FacilityKey (Intercept) 2.372    1.54
Number of obs: 1295, groups:  FacilityKey, 171

Fixed effects:
                        Estimate Std. Error z value Pr(>|z|)
(Intercept)           -5.0890719  0.4455409 -11.422   <2e-16 ***
month_id              -0.0062290  0.0006099 -10.213   <2e-16 ***
factor(CareHomeSize)2  0.2843053  0.3637577   0.782    0.434
factor(CareHomeSize)3  0.4440538  0.3985768   1.114    0.265
factor(Ratings)3       0.4166752  0.3246113   1.284    0.199
factor(Ratings)4       0.1053528  0.4212772   0.250    0.803
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) mnth_d f(CHS)2 f(CHS)3 fc(R)3
month_id    -0.011
fctr(CrHS)2 -0.712  0.000
fctr(CrHS)3 -0.749  0.000  0.754
fctr(Rtng)3 -0.689  0.000  0.106   0.243
fctr(Rtng)4 -0.429  0.001 -0.012   0.045   0.583