Аргумент 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
,