На этот вопрос уже есть ответ:

У меня есть датафрейм с 16 переменными. Когда я делаю множественную линейную регрессию, я делаю следующее:

fit <- lm(y ~ .,data=data)

Теперь я знаю, как добавить член второй степени к одной из переменных:

fit2 <- lm(y ~ poly(x1,2) + .,data=data)

Но сейчас я не хочу писать это для всех моих 16 переменных. Как я могу сделать это простым способом для всех моих переменных?