luranguo
请帮忙看一下,我用Step回归,得到的回归方程是:q4y ~ q4x1 + q4x2 + q4x4
但用leaps,却发现只引入q4x1和q4x2时的Cp值最小,第二小的Cp值才是q4x1,q4x2和q4x4。
那我应该用哪个回归方程最好呢?为什么?
还有怎么利用condition number选择引入变量的个数呢?(数据和R的output如下)
数据:
Case q4x1 q4x2 q4x3 q4x4 q4y
1 7 26 6 60 78.5
2 1 29 15 52 74.3
3 11 56 8 20 104.3
4 11 31 8 47 87.6
5 7 52 6 33 95.9
6 11 55 9 22 109.2
7 3 71 17 6 102.7
8 1 31 22 44 72.5
9 2 54 18 22 93.1
10 21 47 4 26 115.9
11 1 40 23 34 83.8
12 11 66 9 12 113.3
13 10 68 8 12 109.4
> q4step = step(q4ls, direction=\"both\")
Step: AIC=24.97
q4y ~ q4x1 + q4x2 + q4x4
> X = matrix(c(q4x1, q4x2, q4x3, q4x4), n, 4);
> leaps (X,q4y)
$which
1 2 3 4
1 FALSE FALSE FALSE TRUE
1 FALSE TRUE FALSE FALSE
1 TRUE FALSE FALSE FALSE
1 FALSE FALSE TRUE FALSE
2 TRUE TRUE FALSE FALSE
2 TRUE FALSE FALSE TRUE
2 FALSE FALSE TRUE TRUE
2 FALSE TRUE TRUE FALSE
2 FALSE TRUE FALSE TRUE
2 TRUE FALSE TRUE FALSE
3 TRUE TRUE FALSE TRUE
3 TRUE TRUE TRUE FALSE
3 TRUE FALSE TRUE TRUE
3 FALSE TRUE TRUE TRUE
4 TRUE TRUE TRUE TRUE
$Cp
[1] 138.730833 142.486407 202.548769 315.154284 2.678242 5.495851
[7] 22.373112 62.437716 138.225920 198.094653 3.018233 3.041280
[13] 3.496824 7.337474 5.000000