emmeans中所有预测的平均值

Question

我在使用emmeans评估所有预测的平均值(或加权平均值)时遇到了问题。例如，混合模型：

library(emmeans)
library(lme4)
m1 <- lmer(mpg ~ 1 + wt + (1|cyl),data=mtcars)

固定效果wt成功：

emmeans(m1,specs="wt")
   wt emmean   SE   df lower.CL upper.CL
 3.22   20.2 1.71 1.83     12.1     28.3

然而，为了计算预测的平均值，以下方法以前有效(大约12个月前)，但现在失败了：

emmeans(m1,specs="1")
NOTE: Results may be misleading due to involvement in interactions
Error in `[[<-.data.frame`(`*tmp*`, ".wgt.", value = 1) : 
  replacement has 1 row, data has 0

对于简单的线性模型，也会出现同样的错误。非常感谢您的帮助。

Answer 1

当我遇到问题中描述的问题时，我认为我使用的是当前版本的emmeans (1.4.8)。然而，我实际上可能一直在使用emmeans 1.4.6 (请参阅下面Russ been的评论)。我恢复到emmeans v1.4.3，代码正常工作。然后我更新到emmeans的当前版本(1.4.8)，代码继续工作。最有可能的原因是我使用emmeans 1.4.6，它有一个已知的bug。请查看this github entry了解更多信息。