线性回归模型，最大数量为n次幂

Question

我知道，当我使用lm()或glm()函数来拟合R中的回归模型时，可以像这样编写第n级的交互：

fit <- glm(formula=outVar ~ (inVar1 + inVar2 + inVar3)^n,
           data=d)

但是有可能用变量的力量做类似的事情，所以我不必指定I(inVar1^2)，I(inVar1^3)并排除同一变量的不同权力之间的相互作用吗？

修改

我想做这样的事情：

formula=outVar ~ (poly(inVar1 + inVar2 + inVar3, 2))^2

所以我得到了公式

outVar ~ inVar1 + inVar2 + inVar3 + I(inVar1^2) + I(inVar2^2) + I(inVar3^2) + inVar1:inVar2 + inVar1:inVar3 + inVar2:inVar3 + I(inVar1^2):I(inVar2^2) + I(inVar1^2):I(inVar3^2) + I(inVar1^2):I(inVar3^2) + inVar1:I(inVar2^2) + inVar1:I(inVar3^2)...