人均预测的均值

Question

我很难使用平均数来评估所有预测的均值（或加权均值）。例如，混合模型：

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 Lenth的评论）。我回到了emmeans v1.4.3，并且代码正常工作。然后，我更新到了emmeans的当前版本（1.4.8），并且代码继续起作用。最可能的原因是我使用emmeans 1.4.6，它有一个已知的bug。有关更多信息，请参见this github entry。