7021 字
35 分钟
Econometrics-Lesson5:一元线性回归:拟合优度、显著性检验、置信区间与预测

课程说明与考核提醒#

  • 记录日期:2026-10-14;对应课堂:2026-10-13,第11—13节;任课教师:孔伟杰。
  • 涉及教材:李子奈、潘文卿《计量经济学》(第五版),第2章:§2.2 拟合优度(后半部分)→ §2.4 一元线性回归模型的统计检验 → §2.5 一元线性回归分析的应用:预测问题。课件中的假设检验推导还用到 §2.3 基本假设与OLS估计量的统计性质。
  • 作业 / 测试 / 考试说明:现有PPT未出现可确认的具体作业题号、截止时间、测试安排或计分规则。所提供的录音DOCX可提取内容主要是时间戳与课堂画面,缺少可读取的口头转写,因此无法据此核实老师有没有口头布置作业。下面不会把自主练习冒充为老师指定作业。
  • 材料核对:课件第1—29页为知识内容,第30页为教室电脑提示,第31—33页为电脑画面;例题数值以表格原始数据和软件输出为准。部分课件文字与教材中的计算采用不同精度,相关位置已注明。

概述#

这一节课的核心逻辑是:先判断模型对样本解释得好不好,再判断估计出的系数是否显著,最后量化参数和预测结果的不确定性。

  1. 拟合优度(R2R^2):样本点离回归线有多近?YY 的样本波动有多大比例被线性回归解释?
  2. 显著性检验(tt 检验与 pp 值):样本得到的斜率可能来自抽样波动吗?有没有足够证据表明总体斜率不等于零?
  3. 参数置信区间:仅给出一个估计值还不够,总体参数可能处在哪个范围?
  4. 统计预测:给定新的收入 X0X_0,预测平均消费或某个家庭的消费,两种区间为什么宽度不同?

一条主线:

样本回归直线 Y^i=β^0+β^1Xi\hat Y_i=\hat\beta_0+\hat\beta_1 X_i → R2R^2 衡量拟合 → se(β^j)se(\hat\beta_j) 衡量估计精度 → tt / pp 判断显著性 → 参数置信区间 → 均值与个值预测区间。


目录#


一、拟合优度与可决系数(教材 §2.2)#

1. 为什么估计完回归直线后,还要衡量拟合优度?#

已知一元线性回归模型:

Yi=β0+β1Xi+μiY_i=\beta_0+\beta_1X_i+\mu_i

用普通最小二乘法(OLS)得到样本回归函数:

Y^i=β^0+β^1Xi,ei=Yi−Y^i\hat Y_i=\hat\beta_0+\hat\beta_1X_i,\qquad e_i=Y_i-\hat Y_i

其中:

  • YiY_i:真实观测值;Y^i\hat Y_i:回归线给出的拟合值;eie_i:残差。
  • OLS 选择 β^0,β^1\hat\beta_0,\hat\beta_1,使残差平方和 ∑ei2\sum e_i^2 最小。
  • 不同样本、不同变量关系都可能得到一条回归线;还要评价这条线在给定样本中的拟合程度。

2. 总离差分解:TSS=ESS+RSSTSS=ESS+RSS#

令 Yˉ\bar Y 为样本均值,yi=Yi−Yˉy_i=Y_i-\bar Y,y^i=Y^i−Yˉ\hat y_i=\hat Y_i-\bar Y。则:

Yi−Yˉ⏟总离差=(Y^i−Yˉ)⏟回归解释部分+(Yi−Y^i)⏟残差部分\underbrace{Y_i-\bar Y}_{\text{总离差}} =\underbrace{(\hat Y_i-\bar Y)}_{\text{回归解释部分}} +\underbrace{(Y_i-\hat Y_i)}_{\text{残差部分}}

教材图2.2.1:总离差分解为回归解释的离差和残差

直观地看,在同一个 XiX_i 处,YiY_i 到样本均值 Yˉ\bar Y 的竖直距离,可以拆成“回归拟合值到均值的距离”与“观测值到拟合值的距离”。

两边平方再加总:

∑(Yi−Yˉ)2=∑(Y^i−Yˉ)2+∑ei2+2∑(Y^i−Yˉ)ei\sum (Y_i-\bar Y)^2 =\sum(\hat Y_i-\bar Y)^2+\sum e_i^2 +2\sum(\hat Y_i-\bar Y)e_i

对于包含截距项的OLS回归,正规方程保证:

∑ei=0,∑Xiei=0\sum e_i=0,\qquad \sum X_ie_i=0

由于 Y^i−Yˉ=β^1(Xi−Xˉ)\hat Y_i-\bar Y=\hat\beta_1(X_i-\bar X),交叉项为零:

∑(Y^i−Yˉ)ei=β^1[∑Xiei−Xˉ∑ei]=0\sum(\hat Y_i-\bar Y)e_i =\hat\beta_1\left[\sum X_ie_i-\bar X\sum e_i\right]=0

因此:

TSS=ESS+RSS\boxed{TSS=ESS+RSS}
缩写名称公式直观含义
TSSTotal Sum of Squares,总离差平方和∑(Yi−Yˉ)2\sum(Y_i-\bar Y)^2YY 原本有多少波动
ESSExplained Sum of Squares,回归平方和∑(Y^i−Yˉ)2\sum(\hat Y_i-\bar Y)^2被回归线解释的波动
RSSResidual Sum of Squares,残差平方和∑(Yi−Y^i)2\sum(Y_i-\hat Y_i)^2回归后仍未解释的波动
TIP

为什么不能直接比较两个模型的 RSS 大小?

RSS 受 YY 的量纲和原始波动大小影响。将它与同一样本、同一个被解释变量的 TSS 相比,能获得更便于解释的相对拟合指标。改变被解释变量或样本时,比较 R2R^2 也需要谨慎。

3. 可决系数 R2R^2#

R2=ESSTSS=1−RSSTSS\boxed{R^2=\frac{ESS}{TSS}=1-\frac{RSS}{TSS}}

在含截距项的OLS模型且 TSS>0TSS>0 时:

  • 0≤R2≤10\le R^2\le1;R2=1R^2=1 为样本完全拟合;越接近 1,样本拟合通常越紧密。
  • 例如 R2=0.80R^2=0.80,表示样本中 YY 围绕均值的平方离差约有80%被该回归解释。
  • R2R^2 是样本统计量,换一个样本可能得到不同结果。

需要区分:R2R^2 高不能单独证明因果关系成立,也无法保证模型设定正确、预测一定准确。

课件第1页:可决系数定义及离差形式公式

4. R2R^2 与样本相关系数 rr 的关系#

一元含截距项OLS中,令 xi=Xi−Xˉx_i=X_i-\bar X,yi=Yi−Yˉy_i=Y_i-\bar Y,有:

β^1=∑xiyi∑xi2,y^i=β^1xi\hat\beta_1=\frac{\sum x_iy_i}{\sum x_i^2}, \qquad \hat y_i=\hat\beta_1x_i

于是:

R2=∑y^i2∑yi2=β^12∑xi2∑yi2=(∑xiyi)2(∑xi2)(∑yi2)=rXY2\begin{aligned} R^2 &=\frac{\sum\hat y_i^2}{\sum y_i^2} =\frac{\hat\beta_1^2\sum x_i^2}{\sum y_i^2}\\ &=\frac{(\sum x_iy_i)^2}{(\sum x_i^2)(\sum y_i^2)} =\boxed{r_{XY}^2} \end{aligned}

其中:

rXY=∑xiyi∑xi2∑yi2r_{XY}=\frac{\sum x_iy_i}{\sqrt{\sum x_i^2\sum y_i^2}}

注意 R2R^2 丢失了方向信息:正相关或负相关都可能得到较大的 R2R^2;关系方向要看 β^1\hat\beta_1 的正负。

WARNING

课件公式 R2=r2R^2=r^2 的适用条件是一元、包含截距项、用OLS估计。在多元回归、无截距回归或其他估计方法中,不能不加条件直接套用。

二、统计推断的准备:模型假设与抽样分布(教材 §2.3、§2.4)#

1. 为什么 β^\hat\beta 的数值还不足以判断显著性?#

例如某次样本估计得到 β^1=0.67\hat\beta_1=0.67,这只能说明这组样本中的回归斜率是 0.670.67。换样本后,斜率可能变化。为了判断总体中的 β1\beta_1 是否可能为零,必须知道估计量的抽样分布和标准误。

  • 总体参数 β0,β1\beta_0,\beta_1:未知、固定的真实参数。
  • 估计量 β^0,β^1\hat\beta_0,\hat\beta_1:随抽样变化的随机变量。
  • 估计值:从这一组实际样本计算出的具体数字(如 0.670.67)。
  • 标准误 se(β^j)se(\hat\beta_j):估计量抽样波动大小的估计值;越小通常代表估计越精确。

2. 本堂课用到的基本假设#

课件从OLS估计量服从正态分布出发,其依据来自教材§2.3。常规经典线性回归讨论包括:

  1. 模型函数形式与解释变量的选取恰当;
  2. XX 在样本中有变异(∑(Xi−Xˉ)2>0\sum(X_i-\bar X)^2>0);
  3. 随机干扰项条件均值为零:E(μi∣X)=0E(\mu_i\mid X)=0;
  4. 随机干扰项具有条件同方差 Var⁡(μi∣X)=σ2\operatorname{Var}(\mu_i\mid X)=\sigma^2,不同观测间无条件相关(在给定解释变量后);
  5. 做经典小样本精确 tt 推断时还使用条件正态性 μi∣X∼N(0,σ2)\mu_i\mid X\sim N(0,\sigma^2)。

这里 XX 表示给定的解释变量样本向量。前述假设成立时,β^j\hat\beta_j 对总体参数无偏,且具有相应的抽样分布;有关假设失效会影响传统标准误和显著性检验的可靠性。

3. OLS估计量的方差#

记 Sxx=∑(Xi−Xˉ)2=∑xi2S_{xx}=\sum(X_i-\bar X)^2=\sum x_i^2,则条件于解释变量样本:

β^1∣X∼N(β1,σ2Sxx)\hat\beta_1\mid X\sim N\left(\beta_1,\frac{\sigma^2}{S_{xx}}\right)β^0∣X∼N(β0,σ2[1n+Xˉ2Sxx])\hat\beta_0\mid X\sim N\left(\beta_0,\sigma^2\left[\frac1n+\frac{\bar X^2}{S_{xx}}\right]\right)

所以真实标准差为:

sd⁡(β^1∣X)=σSxx,sd⁡(β^0∣X)=σ1n+Xˉ2Sxx\operatorname{sd}(\hat\beta_1\mid X)=\frac{\sigma}{\sqrt{S_{xx}}},\qquad \operatorname{sd}(\hat\beta_0\mid X)=\sigma\sqrt{\frac1n+\frac{\bar X^2}{S_{xx}}}

总体误差方差 σ2\sigma^2 未知,用残差估计:

σ^2=∑ei2n−2=RSSn−2\boxed{\hat\sigma^2=\frac{\sum e_i^2}{n-2}=\frac{RSS}{n-2}}

自由度 n−2n-2 的来源:样本中估计了截距 β0\beta_0 与斜率 β1\beta_1 两个参数,残差因此满足两项OLS约束。

将 σ\sigma 替换为 σ^\hat\sigma,得到实际可算的标准误:

se(β^1)=σ^2Sxx\boxed{se(\hat\beta_1)=\sqrt{\frac{\hat\sigma^2}{S_{xx}}}}se(β^0)=σ^2(1n+Xˉ2Sxx)\boxed{se(\hat\beta_0)=\sqrt{\hat\sigma^2\left(\frac1n+\frac{\bar X^2}{S_{xx}}\right)}}
TIP

记住斜率标准误的结构:

se(β^1)=σ^∑(Xi−Xˉ)2se(\hat\beta_1)=\frac{\hat\sigma}{\sqrt{\sum(X_i-\bar X)^2}}

在其他条件相近时,残差波动越小、XX 的离散程度越大,估计斜率通常越精确。


三、变量的显著性检验(教材 §2.4)#

1. 假设检验到底在检验什么?#

回归分析希望判断:解释变量 XX 对被解释变量 YY 的总体线性关系是否显著。对于一元回归,关键是检验总体斜率 β1\beta_1。

原假设(Null hypothesis)与备择假设(Alternative hypothesis):

H0:β1=0,H1:β1≠0H_0:\beta_1=0,\qquad H_1:\beta_1\ne0
  • H0H_0:总体线性回归斜率为零;
  • H1H_1:总体线性回归斜率不为零;
  • 零假设下,样本斜率即使不完全为零,也可能来自抽样误差。

课件第6页:假设检验的逻辑及小概率事件思想

2. 检验思想:反证法与“小概率事件”#

按照课件的表述,假设检验的直觉是:先按 H0H_0 成立来计算概率,再观察样本结果是否过于极端。

若 H0H_0 成立时,如此极端的统计量很难出现,就倾向于拒绝 H0H_0。但要注意:拒绝或不拒绝原假设始终存在统计错误的可能。

  • 第一类错误:H0H_0 实际正确,却拒绝了它;显著性水平 α\alpha 控制此错误的概率。
  • 第二类错误:H0H_0 实际错误,却没有拒绝它。
  • 不拒绝 H0H_0 只表示当前证据不足以拒绝它,不能推出 H0H_0 已被证明为真。

3. 构造 tt 统计量#

若检验一般假设 H0:β1=β1∗H_0:\beta_1=\beta_1^*:

t=β^1−β1∗se(β^1)∼t(n−2)(H0下)\boxed{t=\frac{\hat\beta_1-\beta_1^*}{se(\hat\beta_1)}\sim t(n-2)\quad(H_0\text{下})}

检验是否等于零时,β1∗=0\beta_1^*=0:

t=β^1se(β^1)\boxed{t=\frac{\hat\beta_1}{se(\hat\beta_1)}}

这一比值表示:估计值距离零有几个标准误。绝对值越大,越不容易用抽样波动来解释。

同样可以检验截距项:

t0=β^0−β0∗se(β^0)t_0=\frac{\hat\beta_0-\beta_0^*}{se(\hat\beta_0)}

为什么是 t(n−2)t(n-2)?由于 σ\sigma 未知,需要用 σ^\hat\sigma 代替;在经典正态假设下,标准化估计量便服从自由度为 n−2n-2 的 tt 分布。

课件第8页:斜率的抽样分布、误差方差估计与t统计量

4. 临界值法(双侧检验)#

以 H0:β1=0H_0:\beta_1=0、H1:β1≠0H_1:\beta_1\ne0 为例:

  1. 指定显著性水平 α\alpha,常见取值为 0.100.10、0.050.05、0.010.01。
  2. 计算 tobs=β^1/se(β^1)t_{obs}=\hat\beta_1/se(\hat\beta_1)。
  3. 查 tt 分布表,找到自由度 n−2n-2 下的双侧临界值 tα/2,n−2t_{\alpha/2,n-2}。
  4. 当 ∣tobs∣>tα/2,n−2|t_{obs}|>t_{\alpha/2,n-2} 时,拒绝 H0H_0;否则不拒绝。

课件第9页:双侧t检验与两端拒绝域

图示中,两端各留出 α/2\alpha/2 的拒绝域,中间为概率 1−α1-\alpha 的不拒绝域。

WARNING

课件说“接受 H0H_0”时,按较严谨的统计术语理解为未能拒绝 H0H_0。在论文和作业中建议使用“拒绝/不拒绝原假设”的表述。

5. pp 值法:与临界值法等价的另一种判断#

对于双侧检验:

p=P(∣Tn−2∣≥∣tobs∣∣H0)\boxed{p=P\left(|T_{n-2}|\ge |t_{obs}|\mid H_0\right)}

也就是:在 H0H_0 成立的前提下,观察到至少和本次样本同样极端的统计量的概率。

  • 若 p<αp<\alpha,拒绝 H0H_0;
  • 若 p≥αp\ge\alpha,不拒绝 H0H_0。

课件第10页:p值、临界值及双侧检验示意

例:p=0.013p=0.013。

  • 与 α=0.05\alpha=0.05 比较:0.013<0.050.013<0.05,在5%水平下显著;
  • 与 α=0.01\alpha=0.01 比较:0.013>0.010.013>0.01,在1%水平下未达到显著。

pp 值的常见误解:p=0.013p=0.013 表示“原假设为真的概率只有1.3%”是错误的。pp 值的定义始终以原假设成立为前提。


四、课堂完整例题:家庭可支配收入与消费支出#

1. 题目与原始数据#

本节反复使用教材家庭“可支配收入 XX—消费支出 YY”的10组观测值,构造线性模型:

Yi=β0+β1Xi+μiY_i=\beta_0+\beta_1X_i+\mu_i

家庭 ii收入 XiX_i消费 YiY_i
1800638
21100935
314001155
417001254
520001408
623001650
726001925
829002068
932002266
1035002530

课件第4页:教材表2.3.1所列参数计算用原始数据及离差汇总

由原始数据计算(数据、输出以课件和教材表格为准):

n=10,Xˉ=2150,Yˉ=1582.9≈1583,Sxx=∑xi2=7 425 000,Sxy=∑xiyi=4 974 750,Syy=∑yi2=3 354 954.9.\begin{gathered} n=10,\quad\bar X=2150,\quad\bar Y=1582.9\approx1583,\\ S_{xx}=\sum x_i^2=7\,425\,000,\\ S_{xy}=\sum x_iy_i=4\,974\,750,\\ S_{yy}=\sum y_i^2=3\,354\,954.9. \end{gathered}

2. 第一步:用OLS计算回归系数#

β^1=SxySxx=4 974 7507 425 000=0.67\hat\beta_1=\frac{S_{xy}}{S_{xx}} =\frac{4\,974\,750}{7\,425\,000}=0.67β^0=Yˉ−β^1Xˉ=1582.9−0.67×2150=142.4\hat\beta_0=\bar Y-\hat\beta_1\bar X =1582.9-0.67\times2150=142.4

因此:

Y^i=142.40+0.67Xi\boxed{\hat Y_i=142.40+0.67X_i}
  • 斜率 0.670.67:样本中,可支配收入每增加1单位,平均拟合消费增加 0.670.67 单位。
  • 截距 142.4142.4:回归线在 X=0X=0 时的拟合值;若样本不覆盖 X=0X=0,其经济解释需要谨慎。

3. 第二步:计算 R2R^2#

ESS=β^12Sxx=0.672×7 425 000=3 333 082.5ESS=\hat\beta_1^2S_{xx}=0.67^2\times7\,425\,000=3\,333\,082.5RSS=TSS−ESS=3 354 954.9−3 333 082.5=21 872.4RSS=TSS-ESS=3\,354\,954.9-3\,333\,082.5=21\,872.4R2=3 333 082.53 354 954.9≈0.9935\boxed{R^2=\frac{3\,333\,082.5}{3\,354\,954.9}\approx0.9935}

解释:约 99.35% 的样本消费支出总平方离差可以由该样本线性回归解释,拟合优度很高。

课件第2页:家庭消费与另一组y1数据的Stata回归输出

WARNING

课件数值的精度差异

第3页把 R2R^2 写成 0.9934,Stata 输出和教材采用 0.9935。根据所列精确平方和计算,四位小数应是 0.9935;复习时以原始平方和及Stata输出为准。

课件第2页同时展示 YY 与另一组 Y1Y_1 各自对 XX 的回归:R2R^2 约分别为 0.9935 与 0.9762。它们用于说明不同数据的拟合紧密程度。由于两次回归的被解释变量不同,不宜仅凭这两个 R2R^2 排名来决定哪个模型更有经济解释力。

4. 第三步:估计随机干扰项的方差#

σ^2=RSSn−2=21 872.48=2734.05\hat\sigma^2=\frac{RSS}{n-2} =\frac{21\,872.4}{8}=2734.05

由此:

se(β^1)=2734.057 425 000≈0.0191891se(\hat\beta_1)=\sqrt{\frac{2734.05}{7\,425\,000}} \approx0.0191891se(β^0)=2734.05(110+215027 425 000)≈44.4467se(\hat\beta_0)=\sqrt{2734.05\left(\frac1{10}+\frac{2150^2}{7\,425\,000}\right)} \approx44.4467

课件第12页:误差方差、斜率和截距标准误及t值的完整计算

5. 第四步:显著性检验#

分别检验:

H0:β1=0与H0:β0=0H_0:\beta_1=0\quad\text{与}\quad H_0:\beta_0=0

用未过早舍入的标准误:

t1=0.670.0191891≈34.92,t0=142.444.4467≈3.20t_1=\frac{0.67}{0.0191891}\approx34.92, \qquad t_0=\frac{142.4}{44.4467}\approx3.20

设 α=0.05\alpha=0.05,df=n−2=8df=n-2=8,双侧临界值为:

t0.025,8=2.306t_{0.025,8}=2.306

  • 斜率检验:∣34.92∣>2.306|34.92|>2.306,拒绝 β1=0\beta_1=0;可支配收入对消费的总体线性影响在5%水平下显著。
  • 截距检验:∣3.20∣>2.306|3.20|>2.306,拒绝 β0=0\beta_0=0;截距在5%水平下显著。

课件第13页:例题t值及5%显著性水平判断

WARNING

两个需要辨认的课件排版 / 舍入问题

  1. 课件第13页有一处将截距 t0t_0 计算式写成了 负号,但 β^0=142.4>0\hat\beta_0=142.4>0,标准误也为正,所以正确结果是 t0≈+3.20t_0\approx+3.20。即便做双侧检验时比较 ∣t0∣|t_0|,也应确保统计量符号正确。
  2. 教材的部分演算先把斜率标准误简写为 0.0190.019,于是出现 0.67/0.019≈35.260.67/0.019\approx35.26。Stata 使用更高精度标准误,得到 34.92。实际报告统计量时,建议保留充分精度再计算。

6. 用p值和回归输出再次核验#

课件Stata回归表给出:

参数系数 Coef.标准误 Std. Err.ttP>∣t∣P>\lvert t\rvert5%水平是否显著
XX0.670.670.01918910.019189134.9234.920.0000.000(显示精度)是
截距142.4142.444.446744.44673.203.200.0130.013是

软件显示 0.000 或 0.0000 表示 p值极小、被四舍五入为零,并不代表数学上恰好等于零。截距的 p=0.013p=0.013 表明其在5%水平下显著,在1%水平下不显著。


五、参数的置信区间(教材 §2.4)#

1. 为什么还要估计区间?#

显著性检验关心能否拒绝某个假设值。参数的**区间估计(interval estimation)**则试图说明:在既定置信水平下,按照一个覆盖率受到控制的构造规则,可能得到多宽的参数区间。

课件先引入形式:

P(β^−δ≤β≤β^+δ)=1−αP(\hat\beta-\delta\le\beta\le\hat\beta+\delta)=1-\alpha

这里应按重复抽样理解概率:总体参数 β\beta 是固定值,随机性来自区间端点随样本变化。比如“95%置信区间”,严格地说,指采用这一构造方法反复抽样,长期来看约有95%的区间覆盖真实参数。

  • 1−α1-\alpha:置信水平(confidence level),常见为95%、99%。
  • α\alpha:显著性水平。
  • tα/2,n−2t_{\alpha/2,n-2}:在双侧区间中使用的临界值。

2. 推导参数置信区间#

根据:

β^j−βjse(β^j)∼t(n−2)\frac{\hat\beta_j-\beta_j}{se(\hat\beta_j)}\sim t(n-2)

两侧各取 α/2\alpha/2,有:

P(−tα/2,n−2≤β^j−βjse(β^j)≤tα/2,n−2)=1−αP\left(-t_{\alpha/2,n-2}\le\frac{\hat\beta_j-\beta_j}{se(\hat\beta_j)}\le t_{\alpha/2,n-2}\right)=1-\alpha

把不等式移项,得到:

βj 的 (1−α)置信区间:β^j±tα/2,n−2 se(β^j)\boxed{\beta_j\text{ 的 }(1-\alpha)\text{置信区间}: \quad\hat\beta_j\pm t_{\alpha/2,n-2}\,se(\hat\beta_j)}

课件第16页:由t统计量推导置信区间

3. 课堂案例:99%参数置信区间#

给定 α=0.01\alpha=0.01,自由度 df=8df=8,查表得:

t0.005,8≈3.355387(表中常简写为 3.355)t_{0.005,8}\approx3.355387\quad(\text{表中常简写为 }3.355)

对斜率:

0.67±3.355387×0.0191891≈(0.6056, 0.7344)0.67\pm3.355387\times0.0191891 \approx\boxed{(0.6056,\ 0.7344)}

对截距:

142.4±3.355387×44.4467≈(−6.736, 291.536)142.4\pm3.355387\times44.4467 \approx\boxed{(-6.736,\ 291.536)}

课件与Stata使用更精确的临界值,故截距区间端点约为 (−6.736,291.536)(-6.736,291.536);若先把临界值舍入成 3.3553.355,区间端点会产生约 0.0170.017 的差异。

课件第17页:99%置信区间的具体数值及软件输出

结论:

  • 斜率的99%置信区间不包含零,相应双侧1%显著性检验拒绝斜率为零。
  • 截距的99%置信区间包含零,相应双侧1%显著性检验不拒绝截距为零。
  • 95%置信区间会比99%区间窄,因为95%对应的临界值更小。
TIP

置信区间与显著性检验的联系:对于检验 H0:βj=0H_0:\beta_j=0 的双侧 tt 检验,若同一置信水平 1−α1-\alpha 的置信区间不包含0,就在显著性水平 α\alpha 下拒绝 H0H_0。

4. 什么因素会影响置信区间宽度?#

从半宽度 tα/2,n−2se(β^j)t_{\alpha/2,n-2}se(\hat\beta_j) 出发:

  • 增大样本容量 nn:在常规可比条件下通常减小标准误,并使tt临界值趋向正态临界值。
  • 降低误差方差 σ2\sigma^2:样本观测更紧贴总体关系,参数标准误通常更小。
  • 增加 XX 的样本离散程度 SxxS_{xx}:斜率标准误随之下降。
  • 降低置信水平(如99%改为95%):临界值下降,区间变窄,但覆盖率要求同时降低。

课件第18页把“拟合程度提高”列为收窄区间的方向。理解时需要补上其他条件相同这一前提:单纯追求高 R2R^2、扩大样本中的变量范围,并不能自动保证估计符合因果解释或经典假设。

课件第18页:置信区间的大小、自由度和估计精度


六、Stata回归结果怎么看?#

课件第19页逐项标出了软件回归表,它将前面的所有统计量汇总在一起。

课件第19页:Stata回归表的各项统计量标注

1. 回归表上半部分:平方和与整体拟合#

Stata字段本例数值(约)意义
Number of obs10样本容量 nn
Model SS3,333,082.5回归平方和 ESSESS
Residual SS21,872.4残差平方和 RSSRSS
Total SS3,354,954.9总离差平方和 TSSTSS
Residual df8残差自由度 n−2n-2
Residual MS2734.05RSS/(n−2)=σ^2RSS/(n-2)=\hat\sigma^2
Root MSE52.288σ^2=σ^\sqrt{\hat\sigma^2}=\hat\sigma
R-squared0.9935R2=1−RSS/TSSR^2=1-RSS/TSS
Adj R-squared0.9927调整后的可决系数
F(1,8)1219.10模型整体显著性检验的F统计量
Prob > F0.0000F检验的p值(显示精度)

关于课件出现但本节没有重点推导的两个字段:

  • 调整 R2R^2:校正解释变量数目对普通 R2R^2 的影响。对于含截距的一元回归,Rˉ2=1−RSS/(n−2)TSS/(n−1)\bar R^2=1-\dfrac{RSS/(n-2)}{TSS/(n-1)}。
  • F检验:在一元含截距回归中,检验 β1=0\beta_1=0 时有 F=t12F=t_1^2,例如 34.922≈1219.434.92^2\approx1219.4;因精度不同会略有差异。更系统的F检验会在多元回归中展开。

2. 回归表下半部分:系数与统计推断#

  • Coef.:β^j\hat\beta_j,参数估计值;
  • Std. Err.:se(β^j)se(\hat\beta_j),标准误;
  • t:系数除以标准误(当 H0:βj=0H_0:\beta_j=0);
  • P>|t|:对应的双侧p值;
  • [95% Conf. Interval]:系数的95%置信区间;
  • _cons:常数项、截距 β^0\hat\beta_0。

本例中,如果只看到 x 的一行 0.67 / 0.0191891 / 34.92 / 0.000,就应立即能读出:估计斜率0.67,估计精度较高,拒绝斜率为0的原假设。


七、一元线性回归的预测问题(教材 §2.5)#

1. 两种预测对象:总体条件均值 vs. 总体个值#

这一节最容易混淆的地方是:给定 X=X0X=X_0,我们可以预测两种对象。

A. 条件均值(conditional mean)

m(X0)=E(Y∣X=X0)=β0+β1X0m(X_0)=E(Y\mid X=X_0)=\beta_0+\beta_1X_0

问的是:所有收入为 X0X_0 的家庭,其平均消费会是多少?

B. 个别值(individual value)

Y0=β0+β1X0+μ0Y_0=\beta_0+\beta_1X_0+\mu_0

问的是:某一个收入为 X0X_0 的家庭,实际消费会是多少?

两者都使用同一个点预测:

Y^0=β^0+β^1X0\boxed{\hat Y_0=\hat\beta_0+\hat\beta_1X_0}

区别在于:预测个别值时,除回归参数的估计误差外,还包含该个体自身的随机扰动 μ0\mu_0。

课件第20页:均值预测和个值预测的学习框架

2. 为什么 Y^0\hat Y_0 是无偏的点预测?#

由OLS参数的无偏性:

E(Y^0∣X)=E(β^0∣X)+X0E(β^1∣X)=β0+β1X0=E(Y∣X0)\begin{aligned} E(\hat Y_0\mid X)&=E(\hat\beta_0\mid X)+X_0E(\hat\beta_1\mid X)\\ &=\beta_0+\beta_1X_0=E(Y\mid X_0) \end{aligned}

在预测新个体时,若该个体随机扰动也满足零均值并与训练样本相应独立,那么:

E(Y^0−Y0∣X,X0)=0E(\hat Y_0-Y_0\mid X,X_0)=0

因此 Y^0\hat Y_0 对条件均值和新个体的真实值都具有无偏预测的性质;个体预测的实际误差仍然可能很大。

课件第21页:条件均值的无偏点预测

3. 条件均值预测的方差:为什么出现 1/n1/n 与 (X0−Xˉ)2(X_0-\bar X)^2?#

先从:

Y^0=β^0+X0β^1\hat Y_0=\hat\beta_0+X_0\hat\beta_1

求方差(X0X_0 给定):

Var⁡(Y^0∣X)=Var⁡(β^0∣X)+X02Var⁡(β^1∣X)+2X0Cov⁡(β^0,β^1∣X)\operatorname{Var}(\hat Y_0\mid X) =\operatorname{Var}(\hat\beta_0\mid X) +X_0^2\operatorname{Var}(\hat\beta_1\mid X) +2X_0\operatorname{Cov}(\hat\beta_0,\hat\beta_1\mid X)

由OLS:

Var⁡(β^0∣X)=σ2(1n+Xˉ2Sxx),Var⁡(β^1∣X)=σ2Sxx,Cov⁡(β^0,β^1∣X)=−Xˉσ2Sxx.\begin{aligned} \operatorname{Var}(\hat\beta_0\mid X)&=\sigma^2\left(\frac1n+\frac{\bar X^2}{S_{xx}}\right),\\ \operatorname{Var}(\hat\beta_1\mid X)&=\frac{\sigma^2}{S_{xx}},\\ \operatorname{Cov}(\hat\beta_0,\hat\beta_1\mid X)&=-\frac{\bar X\sigma^2}{S_{xx}}. \end{aligned}

代入并合并平方:

Var⁡(Y^0∣X)=σ2[1n+Xˉ2+X02−2X0XˉSxx]=σ2(1n+(X0−Xˉ)2Sxx)\begin{aligned} \operatorname{Var}(\hat Y_0\mid X) &=\sigma^2\left[\frac1n+\frac{\bar X^2+X_0^2-2X_0\bar X}{S_{xx}}\right]\\ &=\boxed{\sigma^2\left(\frac1n+\frac{(X_0-\bar X)^2}{S_{xx}}\right)} \end{aligned}

课件第24页:均值预测方差的化简步骤

直观解释:

  • 1/n1/n:样本容量有限,平均水平仍有估计误差;
  • (X0−Xˉ)2/Sxx(X_0-\bar X)^2/S_{xx}:预测点 X0X_0 越远离样本中心 Xˉ\bar X,估计直线的不确定性通常越大;
  • SxxS_{xx} 越大(在其他条件相近时),参数估计误差通常越小。

4. 条件均值的置信区间#

用 σ^2\hat\sigma^2 代替 σ2\sigma^2:

semean(X0)=σ^1n+(X0−Xˉ)2Sxxse_{\mathrm{mean}}(X_0) =\hat\sigma\sqrt{\frac1n+\frac{(X_0-\bar X)^2}{S_{xx}}}

故在 1−α1-\alpha 置信水平下:

Y^0±tα/2,n−2 σ^1n+(X0−Xˉ)2Sxx\boxed{\hat Y_0\pm t_{\alpha/2,n-2}\,\hat\sigma \sqrt{\frac1n+\frac{(X_0-\bar X)^2}{S_{xx}}}}

课件第25页:条件均值的t统计量和置信区间

这条区间描述给定收入水平 X0X_0 的平均消费可能落在哪个区间。

5. 个别值的预测区间:为什么比均值区间多了一个1?#

个别家庭的真实消费是:

Y0=β0+β1X0+μ0Y_0=\beta_0+\beta_1X_0+\mu_0

预测误差:

Y0−Y^0=[(β0+β1X0)−Y^0]⏟回归线估计误差+μ0⏟该家庭自身的随机波动Y_0-\hat Y_0 =\underbrace{\left[(\beta_0+\beta_1X_0)-\hat Y_0\right]}_{\text{回归线估计误差}} +\underbrace{\mu_0}_{\text{该家庭自身的随机波动}}

当新个体误差与估计样本相互独立、同方差时,两项方差可以相加:

Var⁡(Y0−Y^0∣X)=σ2(1+1n+(X0−Xˉ)2Sxx)\operatorname{Var}(Y_0-\hat Y_0\mid X) =\sigma^2\left(1+\frac1n+\frac{(X_0-\bar X)^2}{S_{xx}}\right)

多出来的 11 对应新个体自身的随机误差方差 σ2\sigma^2。

于是其标准误为:

seind(X0)=σ^1+1n+(X0−Xˉ)2Sxxse_{\mathrm{ind}}(X_0)=\hat\sigma\sqrt{1+\frac1n+\frac{(X_0-\bar X)^2}{S_{xx}}}

个别值预测区间(prediction interval):

Y^0±tα/2,n−2 σ^1+1n+(X0−Xˉ)2Sxx\boxed{\hat Y_0\pm t_{\alpha/2,n-2}\,\hat\sigma \sqrt{1+\frac1n+\frac{(X_0-\bar X)^2}{S_{xx}}}}

课件第26页:个别值预测区间及其比均值区间多出的随机误差

TIP

最快记忆方法:

均值: σ^1n+(X0−Xˉ)2Sxx个值: σ^1n+(X0−Xˉ)2Sxx+1\begin{array}{l} \text{均值:}\ \hat\sigma\sqrt{\color{blue}{\frac1n+\frac{(X_0-\bar X)^2}{S_{xx}}}}\\[5pt] \text{个值:}\ \hat\sigma\sqrt{\color{blue}{\frac1n+\frac{(X_0-\bar X)^2}{S_{xx}}}+\color{red}{1}} \end{array}

个值需要覆盖个体自己的不确定性,所以个值预测区间必然比对应均值置信区间更宽。

6. 教材案例的完整预测:收入 X0=1000X_0=1000#

沿用回归结果:

Y^=142.4+0.67X,n=10,Xˉ=2150,Sxx=7 425 000,σ^2=2734.05\hat Y=142.4+0.67X,\quad n=10,\quad \bar X=2150, \quad S_{xx}=7\,425\,000,\quad \hat\sigma^2=2734.05

假设需要预测收入为1000的家庭消费:

Y^0=142.4+0.67×1000=812.4\boxed{\hat Y_0=142.4+0.67\times1000=812.4}

这既可用作该收入组的平均消费预测,也可作某个具体家庭消费的点预测。

(1)95%条件均值置信区间

semean=2734.05[110+(1000−2150)27 425 000]≈27.6\begin{aligned} se_{\mathrm{mean}} &=\sqrt{2734.05\left[\frac1{10}+\frac{(1000-2150)^2}{7\,425\,000}\right]}\\ &\approx27.6 \end{aligned}

因为 t0.025,8=2.306t_{0.025,8}=2.306:

812.4±2.306×27.6≈(748.8, 875.9)\boxed{812.4\pm2.306\times27.6\approx(748.8,\ 875.9)}

(2)95%个别值预测区间

seind=2734.05[1+110+(1000−2150)27 425 000]≈59.1\begin{aligned} se_{\mathrm{ind}} &=\sqrt{2734.05\left[1+\frac1{10}+\frac{(1000-2150)^2}{7\,425\,000}\right]}\\ &\approx59.1 \end{aligned}812.4±2.306×59.1≈(676.1, 948.7)\boxed{812.4\pm2.306\times59.1\approx(676.1,\ 948.7)}

比较:点预测相同,均值置信区间约宽127,个别值预测区间约宽273。单个家庭消费还有自身随机因素,所以可能偏离同收入组平均值更远。

两组数值与教材§2.5的 X0=1000X_0=1000 例子对应。保留的中间位数不同,端点可能存在约0.1的舍入差异。

7. 预测区间的图形与宽度规律#

课件第27页:均值置信带与个别值预测带的对比示意图

图中两条带状区域都围绕样本回归线分布:

  • 中间较窄的带:条件均值的置信区间;
  • 外侧较宽的带:个别值预测区间;
  • 距离样本均值 Xˉ\bar X 越远,带越宽:更远处的预测更依赖估计斜率,估计不确定性增加;
  • 在固定 nn、SxxS_{xx}、σ^\hat\sigma 和置信水平下,X0=XˉX_0=\bar X 时带宽最窄。

课件第29页还列出缩小区间的条件:增加样本容量、降低残差波动、增加解释变量样本分散程度、优先在 Xˉ\bar X 附近预测。对于位于样本观察范围以外的 X0X_0,需要额外小心**外推预测(extrapolation)**可能遇到的函数形式失效问题。

课件第29页:预测区间宽度的影响因素及均值与个值比较


八、公式辨析、易错点与复习总结#

1. 核心公式速查#

问题应用的公式含义
斜率OLS估计β^1=Sxy/Sxx\hat\beta_1=S_{xy}/S_{xx}X变化1单位,Y的拟合变化
截距OLS估计β^0=Yˉ−β^1Xˉ\hat\beta_0=\bar Y-\hat\beta_1\bar X回归线过 (Xˉ,Yˉ)(\bar X,\bar Y)
总平方和分解TSS=ESS+RSSTSS=ESS+RSS需含截距项的OLS
可决系数R2=1−RSS/TSSR^2=1-RSS/TSS样本拟合优度
误差方差估计σ^2=RSS/(n−2)\hat\sigma^2=RSS/(n-2)两个参数占用两个自由度
斜率标准误se(β^1)=σ^/Sxxse(\hat\beta_1)=\hat\sigma/\sqrt{S_{xx}}斜率精度
斜率t检验t=(β^1−β1∗)/se(β^1)t=(\hat\beta_1-\beta_1^*)/se(\hat\beta_1)与假设值相差几个标准误
参数置信区间β^j±tα/2,n−2se(β^j)\hat\beta_j\pm t_{\alpha/2,n-2}se(\hat\beta_j)参数的不确定范围
点预测Y^0=β^0+β^1X0\hat Y_0=\hat\beta_0+\hat\beta_1X_0均值与个值同一个中心
均值区间SEσ^1/n+(X0−Xˉ)2/Sxx\hat\sigma\sqrt{1/n+(X_0-\bar X)^2/S_{xx}}平均响应估计精度
个值区间SEσ^1+1/n+(X0−Xˉ)2/Sxx\hat\sigma\sqrt{1+1/n+(X_0-\bar X)^2/S_{xx}}另加新个体扰动

2. 最容易混淆的五组概念#

(1)拟合优度与显著性

  • R2R^2:当前样本中,回归线解释了多少 YY 的平方离差;
  • tt / pp:关于总体斜率的假设是否被当前样本证据拒绝。

(2)估计值与标准误

  • β^1=0.67\hat\beta_1=0.67:估计的影响大小;
  • se(β^1)=0.0192se(\hat\beta_1)=0.0192:斜率估计的不确定性。

(3)显著性水平与p值

  • α\alpha:分析者预先给出的判断阈值;
  • pp:按原假设和实际样本统计量算出的尾部概率。

(4)参数置信区间与预测区间

  • 参数区间:针对未知的 β0,β1\beta_0,\beta_1;
  • 均值区间:针对 E(Y∣X0)E(Y\mid X_0);
  • 个值区间:针对未来某一个 Y0Y_0。

(5)统计显著与经济显著

  • 统计显著说明样本证据有助于排除某个假设值;
  • 经济显著还需结合系数大小、单位、理论机制和识别条件判断。

3. 做题时的标准步骤#

针对一元线性回归的综合题,推荐按以下顺序:

  1. 写出模型及 n,Xˉ,Yˉ,Sxx,Sxy,Syyn,\bar X,\bar Y,S_{xx},S_{xy},S_{yy};
  2. 求 β^0,β^1\hat\beta_0,\hat\beta_1,写出样本回归方程;
  3. 求 ESS,RSS,R2ESS,RSS,R^2,解释拟合优度;
  4. 求 σ^2\hat\sigma^2 与 se(β^0),se(β^1)se(\hat\beta_0),se(\hat\beta_1);
  5. 列出 H0,H1,αH_0,H_1,\alpha,计算 tt 或 pp,用“拒绝/不拒绝”给结论;
  6. 根据置信水平计算参数区间;
  7. 如需预测,先判断目标是均值还是个体,再选择正确的标准误公式。
WARNING

课件及公式使用边界

  1. R2∈[0,1]R^2\in[0,1] 和 R2=r2R^2=r^2 均有特定的含截距OLS前提。
  2. 小样本 tt 检验与预测区间的精确结论依赖经典模型假设;大样本推断还需相应的渐近条件。
  3. p<αp<\alpha 意味着在给定检验下拒绝原假设,不直接证明因果关系。
  4. R2R^2 高、斜率显著,都不能排除遗漏变量、内生性、函数形式错误。
  5. 预测个别值时,公式根号内必须有额外的 11。

九、自主巩固(非教师布置作业)#

目前提供的材料无法核实老师是否口头指定了教材习题,因此本节仅给出直接对应课堂内容的自测与解答,供复习;不将其列为课程正式作业。

自测1|教材本章练习“证明一元回归的 R2=r2R^2=r^2”#

题目(依据教材第2章练习有关要求整理):证明一元含截距OLS回归中的可决系数 R2R^2 等于 X,YX,Y 样本相关系数 rr 的平方。

解答:

β^1=SxySxx,ESS=β^12Sxx=Sxy2Sxx\hat\beta_1=\frac{S_{xy}}{S_{xx}},\quad ESS=\hat\beta_1^2S_{xx} =\frac{S_{xy}^2}{S_{xx}}

从而:

R2=ESSTSS=Sxy2SxxSyy=(SxySxxSyy)2=r2R^2=\frac{ESS}{TSS} =\frac{S_{xy}^2}{S_{xx}S_{yy}} =\left(\frac{S_{xy}}{\sqrt{S_{xx}S_{yy}}}\right)^2=r^2

要点:用到离差形式、样本平方和分解和包含截距的OLS性质。

自测2|同时判断显著性与区间结论#

题目:样本量 n=10n=10,得到某回归截距 β^0=142.4\hat\beta_0=142.4,其标准误为 44.446744.4467。H0:β0=0H_0:\beta_0=0,给定 α=0.05\alpha=0.05,t0.025,8=2.306t_{0.025,8}=2.306;若改为99%置信水平,t0.005,8=3.355t_{0.005,8}=3.355。问:是否在5%水平下显著?99%参数区间是否包含零?

解答:

t=142.4/44.4467≈3.20>2.306t=142.4/44.4467\approx3.20>2.306

因此在5%水平下拒绝 H0H_0。

99%区间为(使用未提前舍入的临界值):

142.4±3.355387×44.4467≈(−6.736,291.536)142.4\pm3.355387\times44.4467\approx(-6.736,291.536)

包含零,因此在1%水平下不拒绝 H0H_0。两个结论完全相容,因为检验标准不同。

自测3|区分两类预测#

题目:在同一个 X0X_0 处,为什么均值置信区间与个值预测区间具有相同中心,但个值区间更宽?

解答:两者均用 Y^0=β^0+β^1X0\hat Y_0=\hat\beta_0+\hat\beta_1X_0 做点预测,区别在于个值还包含新个体的 μ0\mu_0。若新个体误差和训练样本相应独立、具有方差 σ2\sigma^2,预测误差方差就比均值估计误差方差多 σ2\sigma^2,所以标准误公式根号中多一个 11,区间更宽。


资料定位#

  • 课堂资料:孔伟杰,计量经济学第5次课,2026-10-13;课件PPT第1—29页(第30—33页是课堂电脑画面)。
  • 教材:李子奈、潘文卿,《计量经济学》(第五版),高等教育出版社,2020,第2章:§2.2、§2.3、§2.4、§2.5。尤其参考教材图2.2.1、家庭可支配收入—消费支出例题、§2.5 X0=1000X_0=1000 预测例题。
  • 使用说明:正文配图均采用 Astro 静态绝对路径 /images/econometrics/lesson5/...。应将随附 public/images/econometrics/lesson5/ 放入 Astro 项目根目录的对应位置。本文公式采用 Markdown 中的 KaTeX/MathJax 常用写法,需网站已启用相应数学公式渲染支持。
Econometrics-Lesson5:一元线性回归:拟合优度、显著性检验、置信区间与预测
https://www.sleepyfish2031.top/posts/课程笔记/计量经济学/lesson5/
作者
Sleepyfish
发布于
2026-10-14
许可协议
CC BY-NC-SA 4.0