请帮忙解读结果。

daiwenyan · 2009年1月5日

> seedrates.lm1<-lm(grain~rate,data=seedrates)

> seedrates.lm2<-lm(grain~rate+I(rate^2),data=seedrates)

> anova(seedrates.lm2,seedrates.lm1)

Analysis of Variance Table

Model 1: grain ~ rate + I(rate^2)

Model 2: grain ~ rate

Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)

1 2 0.026286

2 3 0.187000 -1 -0.160714 12.228 0.07294

什么意思没明白，是将两次拟合结果做对比，看有无显著性差异吗？还是什么别的？

yihui · 2009年1月6日

I(rate^2)这一项在alpha=0.1的情况下不能舍弃

jhfu · 2009年1月7日

[quote]引用第1楼谢益辉于2009-01-06 20:54发表的“”:

I(rate^2)这一项在alpha=0.1的情况下不能舍弃[/quote]

那么I(rate^2)语句是什么意思呢，这个语句的作用是什么呢

jhfu · 2009年1月7日

I(rate^2)与rate^2有什么区别啊

yihui · 2009年1月7日

有区别，如果要在回归方程中加入自变量的某种变换（如平方项），则必须用I()“保护”起来。

> x = iris[, 1]<br />
> y = iris[, 2]<br />
> lm(y ~ x)<br />
<br />
Call:<br />
lm(formula = y ~ x)<br />
<br />
Coefficients:<br />
(Intercept)            x  <br />
    3.41895     -0.06188  <br />
<br />
> lm(y ~ x + x^2)<br />
<br />
Call:<br />
lm(formula = y ~ x + x^2)<br />
<br />
Coefficients:<br />
(Intercept)            x  <br />
    3.41895     -0.06188  <br />
<br />
> lm(y ~ x + I(x^2))<br />
<br />
Call:<br />
lm(formula = y ~ x + I(x^2))<br />
<br />
Coefficients:<br />
(Intercept)            x       I(x^2)  <br />
     6.4158      -1.0856       0.0857  <br />

jhfu · 2009年1月8日

[quote]引用第4楼谢益辉于2009-01-07 22:57发表的“”:

有区别，如果要在回归方程中加入自变量的某种变换（如平方项），则必须用I()“保护”起来。

> x = iris[, 1]<br />
> y = iris[, 2]<br />
> lm(y ~ x)<br />
<br />
.......[/quote]<br />
<br />
<br />
明白了，谢谢