R 语言经典实例（原书第 2 版）

340

｜

第

11

章

11.13.4 另请参阅

参见 11.12 节和 11.16 节。

11.14 回归系数的置信区间

11.14.1 问题

正在执行线性回归，需要回归系数的置信区间。

11.14.2 解决方案

将回归模型保存在一个对象中；然后使用 confint 函数提取置信区间：

load(file = './data/conf.rdata')

m <- lm(y ~ x1 + x2)

confint(m)

#> 2.5 % 97.5 %

#> (Intercept) -3.90 6.47

#> x1 -2.58 6.24

#> x2 4.67 5.17

11.14.3 讨论

该解决方案使用模型

y

=

β

0

+

β

1

(

x

1

)

i

+

β

2

(

x

2

)

i

+

ε

i

。函数 confint 返回截距（

β

0

）、

x

1

的

系数（

β

1

）和

x

2

的系数（

β

2

）的置信区间：

confint(m)

#> 2.5 % 97.5 %

#> (Intercept) -3.90 6.47

#> x1 -2.58 6.24

#> x2 4.67 5.17

默认情况下，confint 使用 95% 的置信度。可以使用 level 参数选择不同的置信

水平：

confint(m, level = 0.99)

#> 0.5 % 99.5 %

#> (Intercept) -5.72 8.28

#> x1 -4.12 7.79

#> x2 4.58 5.26

11.14.4 另请参阅

arm 包的 coefplot 函数可以绘制回归系数的置信区间。