У меня есть набор данных с 70 переменными, и я хочу попробовать полиномиальную регрессию на нем. Если бы число столбцов было три / четыре, я мог бы просто написать что-то вроде этого -

 model <- lm(y ~ poly(var1,3) + poly(var2,3) + poly(var4,4)

Как бы мы поступили об этом, если бы у нас было 70 переменных? Должны ли мы вводить вручную имена всех переменных или есть более простой метод?