glm(y~x+I(x^2), family = gaussian) # non-orthogonal polynomial
glm(y~poly(x,2), family = gaussian) # orthogonal