My questions is, that I see people using R´s lm() (linear regression model) with Y ~ X^2
e.g. here: Simple non-linear regression problem
But I dont see how and why it works, hence lm() is clearly stated as a linear estimator (using OLS), why does it work for quadratic, clearly non-linear estimation, as well?
Thank you.