glmer预测：关于结果变量的警告信息不是一个因素

Question

我有一个混合效果模型，二项式结果符合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")

我收到以下警告信息：

Warning message:
In model.frame.default(delete.response(Terms), newdata, na.action = na.action,  :
  variable 'Sex' is not a factor

然而，当我检查时，它看起来像是：

str(Orthodont["Sex"])

仍然会创建变量fitted并且值很有意义，但我对此消息感到好奇。有什么我应该关注的吗？否则，此消息的目的是什么。

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

Answer 1

这似乎是一个无害的错误，当从具有因子响应的模型进行预测时总会发生（问题是我们使用xlev参数到model.frame包括响应变量的级别）。我在https://github.com/lme4/lme4/issues/205发布了这个帖子。谢谢你的报告！

您可以忽略警告或将响应变量强制转换为二进制值（这将产生相同的结果）。