I am trying to use the predict function to estimate a dependent variable in my mixed model given value of the two independent variables that are in the model. The estimable function in lmer for some reasons doesn't work with glmer. Can someone help with how I could estimate this? see my code. I got error that my object is not a matrix.

model1=glmer(response~treatment+month+treatment*month +

(1|id),

data=data1,

family=binomial(link = "logit"),

nAGQ=100)

summary(model1)

nedff = data.frame(treatment =1, month =12)

predict(model1, nedff, se.fit=T, interval ="confidence")