> init.h = c(600,700,800,950,1100,1300,1500) > h.d = c(253, 337, 395, 451, 495, 534, 573) > res.lm3 = lm(h.d ~ init.h + I(init.h^2) + I(init.h^3)) > res.lm4 = update(res.lm3, . ~ . + I(init.h^4)) > anova(res.lm3, res.lm4) Analysis of Variance Table Model 1: h.d ~ init.h + I(init.h^2) + I(init.h^3) Model 2: h.d ~ init.h + I(init.h^2) + I(init.h^3) + I(init.h^4) Res.Df RSS Df Sum of Sq F Pr(>F) 1 3 48.3 2 2 12.7 1 35.5 5.58 0.14