我的因变量为y，自变量分别为x1，x2，x3，x4（实现情况大约有二十多个变量，这里简化了一下）。其中x1，x2是0、1变量，其余两个为连续性变量。问题是这样的：1.我想生成每一个变量的平方项（0、1变量的平方项还是他自己就不生成了），有什么便捷的命令可以一次性生成吗？
最终效果：y=x1+x2+x3+x4+x3*x3+x4*x4    式（1）


2.然后生成i和j的交互项，i不等于j，有什么便捷的命令可以一次性生成吗？效果如下：
y=x1+x2+x3+x4+x3*x3+x4*x4+x1*x2+x1*x3+x1*x4+x2*x3+x2*x4+x3*x4     式（2）


3.最后在式（2）回归的基础上计算x1的系数乘以x1的均值+x2的系数乘以x2的均值+x3的系数乘以x3的均值+x4的系数乘以x4的均值（我知道predict可以生成拟合值，但是我是在式2上进行回归，predict生成中包含了平方项的系数与均值的乘积，交互项的系数与均值的乘积，这些不是我想要的）