Fitting Polynomial Regression Models In R