0

emmeans我正在尝试使用包中的功能在 R 中进行事后测试emmeans。但是,我不知道应该在 specs 参数中添加什么。据我了解,这是我放置要对比的变量(我的自变量)的地方。但是,当我放置 IV/IV 时,它会出现错误。我将我的代码和错误放在下面:

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功能做到这一点?
此致

4

1 回答 1

1

来自?emmeans

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

您的spec论点X包含“A”、“B”和“C”(重复 50 次)。但是回归中的变量是 X、Y 及其交互作用。你需要包括类似的东西

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

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

简短版本:您将字符"X"与变量混淆了X

于 2018-05-12T13:19:47.070 回答