Thanks! Yes I’m using the latest versions of the packages installed from the CRAN. I posted in the emmeans github and received the following response:
I think perhaps that what is happening is that
epred() doesn’t produce numerical predictions. The brms package has its own emmeans support, and you may need to ask the developer of that package for details, and/or ask if the
epred option should work.
I’m attaching some data and a reproducible example:
Example.csv (6.5 KB)
model<-brm(formula=Presence ~ Region + Sex + Age + DD + PMI + (1 | ID), data=example, family=“bernoulli”,warmup=2000,iter=4000,seed=34)
emmeans(model, ~ Region, epred=TRUE, re_formula=NULL)
The warning only appears with a bernoulli model (I’ve run negative binomial and zero-inflated beta without the warning).