> 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