Переменный порядок в терминах взаимодействия
Я пытаюсь подогнать ряд линейных моделей, как показано ниже. Важно, чтобы все термины взаимодействия были отсортированы лексикографически. Обратите внимание, что во второй модели отсутствует основной эффект для х.
x = rnorm(100)
y = rnorm(100)
z = x + y + rnorm(100)
m1 = glm(z ~ x + y + x:y)
m2 = glm(z ~ y + x:y)
Модели не ведут себя должным образом в отношении условий взаимодействия:
m1:
x:y -0.1565 0.1151 -1.360 0.1770
m2:
y:x -0.2776 0.1416 -1.961 0.0528 .
Я понимаю, что может быть способ использовать функцию взаимодействия () с аргументом lex.order, но я не могу понять, как и действительно ли это лучший путь. Совет?