Question

我试图用emmeans包中的emmeans函数来计算R中的posthoc测试。但是，我无法找出我应该在规格参数中加入什么。据我所知，它是我想要对比的变量（我的自变量）。但是，当我把我的IV / IVs它给出错误。我将我的代码和错误放在下面：

X= as.factor(rep(c("A", "B", "C"), each= 50))
Y= as.factor(c("K", "L", "M"), times= 50)
Z= rnorm(150)
Model= lm(Z ~ X+Y+X:Y)
emmeans(Model, X)

`错误没有在参考网格中命名为a的变量。

我认为规格论证与我的想法完全不同。我想对比主要效果和互动效果。如何使用emmeans功能执行此操作？最好的问候

Answer 1

来自?emmeans：

spec一个字符向量，指定需要EMM的预测变量的名称。

您的spec参数为X，其中包含＆＃34; A＆＃34;，＆＃34; B＆＃34;和＆＃34; C＆＃34; （重复50次）。但是回归中的变量是X，Y和它们的相互作用。你需要包括像

这样的东西

emmeans(Model, c("X", "Y"))

估算这些变量的边际均值。

简短版：您将字符"X"与变量X混淆。