On Tue, Feb 03, 2004 at 05:27:40PM +0000, Prof Brian Ripley wrote: >> out <- lm(T ~ poly(X, 4), data = d) > You asked for 4 and got 2? Really? No, in fact I have a 4-power model, I just posted example with power of 2, for simplicity. I forgot to replace the poly argument :) Thanks too all who helped me! -- WBR, Timur.