I have a model that I\'m working with that uses a 2nd order orthogonal polynomial as a predictor using the R function poly(). I like the fit that it produces and the fact th