2

不应该

AIC(full) = 275.93

当 step() 函数以低于 -9.86 的完整模型运行时，匹配 AIC 的输出

Start:  AIC=-9.86
y ~ x + x2

       Df Sum of Sq    RSS      AIC
- x2    1   0.03672 85.372 -11.8147
- x     1   1.03869 86.374 -10.6479
<none>              85.336  -9.8578

Step:  AIC=-11.81
y ~ x

       Df Sum of Sq    RSS     AIC
- x     1     1.004 86.376 -12.646
<none>              85.372 -11.815

Step:  AIC=-12.65
y ~ 1


Call:
lm(formula = y ~ 1, data = data)

Coefficients:
(Intercept)  
   -0.03719

这是完整的代码：

set.seed(101)
y = rnorm(100)
x = rnorm(100)
x2 = rnorm(100)
data = data.frame(y = y, x = x, x2 = x2)
null = lm(y~1, data = data)
full = lm(y~x+x2, data = data)
#step(null, scope= list(lower = null, upper = full) , direction="backward", trace = TRUE)
step(full, direction="backward", trace = TRUE)
AIC(full)
4

1 回答 1

3

你需要extractAIC而不是AIC

extractAIC(lm(y~x+x2, data = data), scale=0)
#OR
#extractAIC(full, scale=0)

如果您参考文档?AIC并且?extractAIC它清楚地说明了

对数似然以及因此的 AIC/BIC 仅定义为一个附加常数。不同的常数通常用于不同的目的，因此 extractAIC 和 AIC 可能会给出不同的值（并且对于“lm”类的模型也是如此：参见 extractAIC 的帮助）。


希望这可以帮助！

于 2017-09-06T13:07:44.257 回答