课程笔记

9.2一元线性回归

模型

给定解释变量的观测值 x1,…,xnx_1,\ldots,x_n,一元线性回归模型设为

Yi=β0+β1xi+εi,i=1,…,n,Y_i=\beta_0+\beta_1x_i+\varepsilon_i,\qquad i=1,\ldots,n,

其中 β0\beta_0 是截距,β1\beta_1 是斜率,εi\varepsilon_i 是不可观测的随机误差。YiY_i 是随机变量,xix_i 当作已知常数。

计算最小二乘估计本身不需要误差分布假设;讨论它的无偏性、方差和后续推断时,通常假定三条基本条件:E(εi)=0E(\varepsilon_i)=0;每个误差的方差相同,即 D(εi)=σ2D(\varepsilon_i)=\sigma^2;不同误差不相关,即 Cov⁡(εi,εj)=0\operatorname{Cov}(\varepsilon_i,\varepsilon_j)=0(i≠ji\ne j)。若要使用精确的 tt、FF 推断,还要进一步假定 ε1,…,εn\varepsilon_1,\ldots,\varepsilon_n 相互独立且服从正态分布。此外要求 n>2n>2,并且 x1,…,xnx_1,\ldots,x_n 不全相同;否则 Lxx=0L_{xx}=0,斜率无法识别。

斜率 β1\beta_1 表示 xx 增加一个单位时条件均值 E(Y∣x)E(Y\mid x) 的变化量。单个 YY 还受随机误差影响,观测值不会恰好按 β1\beta_1 变化。回归关系说明的是条件均值随 xx 如何变化,本身不自动说明 xx 是 YY 的原因。

最小二乘估计与拟合

最小二乘估计

用直线 y^=b0+b1x\hat y=b_0+b_1x 拟合数据,最小二乘法选择使残差平方和

Q(b0,b1)=∑i=1n(yi−b0−b1xi)2Q(b_0,b_1)=\sum_{i=1}^n(y_i-b_0-b_1x_i)^2

达到最小的 b0,b1b_0,b_1。记样本均值 xˉ=1n∑ixi\bar x=\frac1n\sum_i x_i、yˉ=1n∑iyi\bar y=\frac1n\sum_i y_i,以及

Lxx=∑i=1n(xi−xˉ)2,Lyy=∑i=1n(yi−yˉ)2,Lxy=∑i=1n(xi−xˉ)(yi−yˉ).L_{xx}=\sum_{i=1}^n(x_i-\bar x)^2,\qquad L_{yy}=\sum_{i=1}^n(y_i-\bar y)^2,\qquad L_{xy}=\sum_{i=1}^n(x_i-\bar x)(y_i-\bar y).

当 Lxx>0L_{xx}>0 时,最小二乘估计为

β^1=LxyLxx,β^0=yˉ−β^1xˉ.\hat\beta_1=\frac{L_{xy}}{L_{xx}},\qquad \hat\beta_0=\bar y-\hat\beta_1\bar x.
证明 由正规方程解出截距与斜率

对 QQ 分别求偏导并令其为零:

∂Q∂b0=−2∑i=1n(yi−b0−b1xi)=0,\frac{\partial Q}{\partial b_0} =-2\sum_{i=1}^n(y_i-b_0-b_1x_i)=0,∂Q∂b1=−2∑i=1nxi(yi−b0−b1xi)=0.\frac{\partial Q}{\partial b_1} =-2\sum_{i=1}^n x_i(y_i-b_0-b_1x_i)=0.

第一式给出 b0=yˉ−b1xˉb_0=\bar y-b_1\bar x。代入第二式,并利用 ∑ixi=nxˉ\sum_i x_i=n\bar x 与 nxˉyˉ=yˉ∑ixin\bar x\bar y=\bar y\sum_i x_i,得到

∑ixiyi−nxˉyˉ+b1(nxˉ2−∑ixi2)=0.\sum_i x_iy_i-n\bar x\bar y +b_1\left(n\bar x^2-\sum_i x_i^2\right)=0.

左边第一项正是 LxyL_{xy},而

∑ixi2−nxˉ2=∑i(xi−xˉ)2=Lxx,\sum_i x_i^2-n\bar x^2=\sum_i(x_i-\bar x)^2=L_{xx},

所以 Lxy−b1Lxx=0L_{xy}-b_1L_{xx}=0,即 b1=Lxy/Lxxb_1=L_{xy}/L_{xx}。

拟合值为 y^i=β^0+β^1xi\hat y_i=\hat\beta_0+\hat\beta_1x_i,残差为 ei=yi−y^ie_i=y_i-\hat y_i。含截距时,最小二乘解满足

∑i=1nei=0,∑i=1nxiei=0,\sum_{i=1}^n e_i=0,\qquad \sum_{i=1}^n x_ie_i=0,

因此回归直线经过点 (xˉ,yˉ)(\bar x,\bar y)。

残差 eie_i 由同一组数据拟合得到,彼此通常并不独立;它也不等于不可观测的随机误差 εi\varepsilon_i。不能把残差当作误差的独立样本。

平方和分解与拟合程度

含截距时,总平方和 SSTSST、回归平方和 SSRSSR 与残差平方和 SSESSE 满足

SST=Lyy=SSR+SSE,SST=L_{yy}=SSR+SSE,

其中

SSR=∑i=1n(y^i−yˉ)2=Lxy2Lxx,SSE=∑i=1nei2=Lyy−Lxy2Lxx.SSR=\sum_{i=1}^n(\hat y_i-\bar y)^2=\frac{L_{xy}^2}{L_{xx}},\qquad SSE=\sum_{i=1}^n e_i^2=L_{yy}-\frac{L_{xy}^2}{L_{xx}}.

误差方差 σ2\sigma^2 的无偏估计为

s2=SSEn−2.s^2=\frac{SSE}{n-2}.

自由度是 n−2n-2,因为模型中估计了截距和斜率两个参数。

当 Lyy>0L_{yy}>0 时,决定系数为

R2=1−SSESST=SSRSST=r2,R^2=1-\frac{SSE}{SST}=\frac{SSR}{SST}=r^2,

其中

r=LxyLxxLyyr=\frac{L_{xy}}{\sqrt{L_{xx}L_{yy}}}

是样本相关系数。

R2R^2 度量的是样本内的线性拟合程度。它随样本变化,不能说明模型条件成立,也不保证对新的观测预测准确。R2=r2R^2=r^2 只在含截距的一元线性回归中成立;换成多元回归或去掉截距,这个等式不再成立。

斜率的区间与检验

在一元线性回归中,设 n>2n>2、Lxx=∑i(xi−xˉ)2>0L_{xx}=\sum_i(x_i-\bar x)^2>0,s2=SSE/(n−2)s^2=SSE/(n-2),且误差独立同分布为正态,s>0s>0。

以下结论都在正态误差模型下成立;tt 分布与上侧分位点的定义见上侧分位点。估计量的方差为

D(β^1)=σ2Lxx,D(β^0)=σ2(1n+xˉ2Lxx).D(\hat\beta_1)=\frac{\sigma^2}{L_{xx}},\qquad D(\hat\beta_0)=\sigma^2\left(\frac1n+\frac{\bar x^2}{L_{xx}}\right).

用 s2s^2 代替 σ2\sigma^2 后,斜率的 1−α1-\alpha 置信区间为

β^1±tα/2(n−2)sLxx,\hat\beta_1\pm t_{\alpha/2}(n-2)\frac{s}{\sqrt{L_{xx}}},

其中 s=s2s=\sqrt{s^2},tα/2(n−2)t_{\alpha/2}(n-2) 是自由度为 n−2n-2 的 tt 分布上侧 α/2\alpha/2 分位点。

检验 H0:β1=β1,0H_0:\beta_1=\beta_{1,0} 时,使用统计量

T=β^1−β1,0s/Lxx∼t(n−2),T=\frac{\hat\beta_1-\beta_{1,0}}{s/\sqrt{L_{xx}}}\sim t(n-2),

双侧拒绝域为 ∣T∣>tα/2(n−2)|T|>t_{\alpha/2}(n-2)。

对 H0:β1=0H_0:\beta_1=0,也可以用

F=SSRSSE/(n−2)∼F(1,n−2),F=\frac{SSR}{SSE/(n-2)}\sim F(1,n-2),

并且 F=T2F=T^2。对同一个原假设,单侧 FF 检验与双侧 tt 检验给出相同的结论。

均值区间与个体预测区间

在一元线性回归中,设 n>2n>2、Lxx=∑i(xi−xˉ)2>0L_{xx}=\sum_i(x_i-\bar x)^2>0,s2=SSE/(n−2)s^2=SSE/(n-2),且误差独立同分布为正态,s>0s>0。

在给定的 x0x_0 处,点预测为

y^0=β^0+β^1x0.\hat y_0=\hat\beta_0+\hat\beta_1x_0.

条件均值 E(Y∣x0)E(Y\mid x_0) 的 1−α1-\alpha 置信区间为

y^0±tα/2(n−2)s1n+(x0−xˉ)2Lxx.\hat y_0\pm t_{\alpha/2}(n-2)s\sqrt{\frac1n+\frac{(x_0-\bar x)^2}{L_{xx}}}.

设新个体 Y0Y_0 与已有观测独立,在 x0x_0 处服从同一线性模型且误差方差仍为 σ2\sigma^2,则它的预测区间为

y^0±tα/2(n−2)s1+1n+(x0−xˉ)2Lxx.\hat y_0\pm t_{\alpha/2}(n-2)s\sqrt{1+\frac1n+\frac{(x_0-\bar x)^2}{L_{xx}}}.

预测区间根号内多出的 11 对应新个体自身的随机波动,所以同一位置上的预测区间总比均值响应区间宽。两种区间都在 x0x_0 远离 xˉ\bar x 时变宽。

若 x0x_0 超出观测到的 xx 范围,公式仍按同一条直线外推,但线性关系在范围外是否成立没有任何数据支持;外推结论要比内插结论弱得多。

均值响应区间描述条件均值的不确定性;新个体预测区间还包括个体波动。同一位置上的后者更宽。外推时,即使公式可计算,观测范围外的线性关系仍需另行判断。