glmer predict：关于结果变量不是一个因素的警告消息

我有一个混合效应模型，二项结果与

glmer

拟合。出于绘图目的，我想预测一个小数据集的总体水平值

下面是一个示例，说明了我的方法：

silly <- glmer(Sex ~ distance +age + (1|Subject), data=Orthodont, family=binomial)

sillypred <- expand.grid(distance=c(20, 25), age=unique(Orthodont$age))
sillypred$fitted <- predict(silly, sillypred, re.form=NA, type="response")

但是，当我检查时，它看起来是：

str(Orthodont["Sex"])

变量

fitted

仍在创建中，且值有意义，但我对这条消息很好奇。有什么我应该关心的吗？否则，此消息的目的是什么

---

这似乎是一个微不足道的问题（毕竟，这一切似乎都很有效），在这种情况下，我表示歉意，但我想确保我没有忽略一些重要的事情。

这似乎是一个无害的错误，当从具有因子响应的模型进行预测时总会出现（问题是我们使用

xlev

参数来

model.frame

，包括响应变量的级别）。我已将此发布在。感谢您的报告

您可以忽略警告，也可以将响应变量强制为二进制值（这将给出相同的结果）

str(Orthodont["Sex"])