emmeans 中所有预测的平均值

Mean of all predictions in emmeans

我在使用 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

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

当我遇到问题中描述的问题时,我以为我使用的是当前版本的 emmeans (1.4.8)。但是,我实际上可能一直在使用 emmeans 1.4.6(请参阅下面 Russ Lenth 的评论)。我恢复到 emmeans v1.4.3 并且代码有效。然后我更新到当前版本的 emmeans (1.4.8),代码继续工作。最有可能的原因是我使用了 emmeans 1.4.6,它有一个已知的错误。请参阅 this github entry 了解更多信息。