知识点

一元线性回归

定义

给定解释变量的观测值 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 只在含截距的一元线性回归中成立;换成多元回归或去掉截距,这个等式不再成立。

例子

由三组数据算出回归直线与误差方差

数据为 (x,y)=(1,2),(2,3),(3,5)(x,y)=(1,2),(2,3),(3,5)。此时 xˉ=2\bar x=2,yˉ=10/3\bar y=10/3,并且

Lxx=2,Lxy=3.L_{xx}=2,\qquad L_{xy}=3.

于是

β^1=32,β^0=13,y^=13+32x.\hat\beta_1=\frac32,\qquad \hat\beta_0=\frac13,\qquad \hat y=\frac13+\frac32x.

三个拟合值为 11/6,20/6,29/611/6,20/6,29/6,残差依次为 1/6,−1/3,1/61/6,-1/3,1/6,所以

SSE=16,s2=SSEn−2=16.SSE=\frac16,\qquad s^2=\frac{SSE}{n-2}=\frac16.

本例中 R2=27/28≈0.964R^2=27/28\approx0.964,样本内拟合看起来很好,但误差方差估计的自由度只有 11。能算出回归直线,并不表示这组数据已经提供了足够证据支持可靠的推断。

探索

一个观测值怎样改变拟合直线?

调整第三个点,观察直线与残差平方和。

正在准备演示…

题集

练习

  1. 对数据 (x,y)=(1,2),(2,3),(3,5)(x,y)=(1,2),(2,3),(3,5),斜率的最小二乘估计为 。

    参考解答

    由 Lxy=3L_{xy}=3、Lxx=2L_{xx}=2,得 β^1=3/2=1.5\hat\beta_1=3/2=1.5。

  2. 对含截距的一元线性回归,决定系数 R2R^2 与样本相关系数 rr 的关系是

    选择答案
    参考解答

    在含截距的一元线性回归中

    R2=SSRSST=Lxy2LxxLyy=r2.R^2=\frac{SSR}{SST}=\frac{L_{xy}^2}{L_{xx}L_{yy}}=r^2.
  3. 含截距的一元线性回归的残差 ei=yi−y^ie_i=y_i-\hat y_i 一定满足

    题 residual-property · choice-residual
    参考解答

    前两条是最小二乘正规方程的直接结果。残差由同一组数据共同拟合得到,通常并不独立,也不等于不可观测的随机误差。