生成lm的所有交互变量

Question

我有一个非常基本的问题： 如何生成所有线性回归的交互作用？

我认为使用y~.可以得到所有主要效果，使用y~. + .^2可以得到所有主要效果以及二阶相互作用。

但是我如何也包括“正常”交互作用？ （例如y~x,y,z,x*y,y*z,x*z）

还是这些已经由y~.生成了？如果是这样，我如何仅获得许多可变项（n> 500）的主要效果（例如y~x,y,z）。

在此先感谢您，对基本问题表示抱歉。