Suponha que eu tenha os seguintes dados
x<-rnorm(100,1,10000)
y<-rnorm(100,1,10000)+2*x+x^2
Se eu usar a função lm da seguinte maneira:
model1<-lm(y~x+x^2)
O R não entende que é para colocar entre as variáveis independentes o termo x elevado ao quadrado. Ele simplesmente ignora o termo e faz a regressão do modelo como o código abaixo:
model2<-lm(y~x)