R中多个预测变量的多项式回归

Question

如何在所有数据帧变量上拟合多项式回归？这是一个线性回归，我会这样做：

lm.fit = lm(response~., data=dataset, subset=train)

但如果我想要一个2级（或更高）的等级，我不能使用poly()，如下所示：

lm.fit = lm(reads_avg~poly(., i), data=dataset, subset=train)

特别是，poly似乎不接受点（.）作为输入。