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

Я пытаюсь выяснить, чтобы сделать 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 функционировать?
С наилучшими пожеланиями

1 ответ

Решение

От ?emmeans:

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

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

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

оценить предельные средние для этих переменных.

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

Другие вопросы по тегам