poly(*,*) in lm() (PR#8972)
I used the function lm() to calculate the coefficients of a polynome. If I used the function poly(t,2) to denote a polynome of form 1 + x + x^2, the coefficients are wrong.
If you take a look at the manual page for poly, you will find that poly (by default) will calculate _orthogonal_ polynomials, which are not the same as 1, x, x^2, etc.
Bj?rn-Helge Mevik