r - 交互项中的变量顺序

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() 函数，但我不知道如何，或者实际上，这是否是最好的方法。建议？