Tengo un conjunto de datos con 70 variables, y quiero probar la regresión polinómica en él. Si el número de columnas fuera tres / cuatro, podría simplemente codificar algo como esto:

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

¿Cómo haríamos esto, si tenemos 70 variables? ¿Deberíamos escribir manualmente los nombres de todas las variables o hay un método más fácil?