交互术语中的变量顺序

Question

我正在尝试安装一些线性模型，如下所示。所有交互术语按字典顺序排序非常重要。请注意，第二个模型缺少x的主要效果。

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参数的interaction（）函数，但我无法弄清楚这是否是最好的方法。建议？