аргумент specs в функции emmeans с R

Question

Я пытаюсь выяснить, как выполнить posthoc-тест в R с помощью функции emmeans из пакета emmeans. Тем не менее, я не мог выяснить, что я должен поставить в аргумент спецификации. насколько я понимаю, это где я помещаю переменные, которые я хочу контрастировать (мои независимые переменные). Однако, когда я ставлю свои 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)

`ошибка - переменная не названа в сетке ссылок.

Полагаю, аргумент о спецификациях совершенно отличается от того, что я думаю. Я хочу сравнить как основные эффекты, так и эффекты взаимодействия. Как я могу сделать это с функцией emmeans?
С наилучшими пожеланиями

Answer 1

С ?emmeans:

spec Вектор символов, указывающий имена предикторов, для которых требуются EMM.

Ваш аргумент spec равен X который содержит «A», «B» и «C» (повторяется 50 раз).Но ваши переменные в регрессии - это X, Y и их взаимодействие.Вам необходимо добавить что-то вроде

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

для оценки среднего значения этих переменных.

Краткая версия: вы путаете символ "X" с переменной X.