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 了解更多信息。
我在使用 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 了解更多信息。