多元回归:深入专题¶
本页收拾几个让回归"用得更好"的专题:度量单位与 \(\beta\) 系数(让不同量纲的变量可比)、对数 / 二次 / 交互等函数形式、用调整 \(R^2\) 做模型选择,以及预测与残差分析。
核心结论速览
| 专题 | 要点 |
|---|---|
| 度量单位 | 改单位:系数、标准误、置信区间端点成倍变;\(t\)、\(F\)、\(R^2\) 不变 |
| \(\beta\) 系数 | 标准化回归的斜率:\(x_j\) 增 1 个标准差 \(\Rightarrow y\) 变 \(b_j\) 个标准差,量纲无关 |
| 对数模型 | 因变量取 log 读"%变化";log-level 大变化用精确式 \(\%\Delta y=100(e^{\beta\Delta x}-1)\) |
| 二次项 | 转折点 $x^*=\left |
| 调整 \(R^2\) | \(\bar R^2=1-\dfrac{\mathrm{SSR}/(n-k-1)}{\mathrm{SST}/(n-1)}\),对加变量惩罚,可比非嵌套模型 |
度量单位与 \(\beta\) 系数¶
改变变量的度量单位不改变任何检验结论:\(t\)、\(F\)、\(R^2\) 不变;但系数、标准误、回归标准误、置信区间端点都按比例变化。若把某变量乘以 \(c\),对应系数与标准误同乘 \(c\)(或 \(1/c\)),残差乘 \(c\)、\(\mathrm{SSR}\) 乘 \(c^2\)。
为消除量纲、让系数可比,可做标准化回归(把每个变量减均值除以样本标准差),得到 \(\beta\) 系数 / 标准化系数:
\(\beta\) 系数的含义
\(x_j\) 提高一个标准差,\(y\) 平均变化 \(\hat b_j\) 个标准差。标准化回归的截距为 0,且 \(t\) 统计量与原回归相同(标准化不影响显著性)。它常用于比较不同量纲变量"谁的影响更大"。
函数形式¶
对数模型¶
四种对数组合及其斜率解读(详见简单回归):因变量取 log → 读作 \(y\) 的百分比变化;自变量取 log → 原因读作百分比变化。
log-level 的精确百分比变化
log-level 模型里 \(\%\Delta y\approx 100\beta\,\Delta x\) 只是近似,\(\beta\Delta x\) 较大时误差明显。精确百分比变化为
推导
由 \(\Delta\ln y=\beta\Delta x\):
故 \(\%\Delta y=100(e^{\beta\Delta x}-1)\)。\(\beta\Delta x\) 很小时 \(e^{\beta\Delta x}-1\approx\beta\Delta x\),回到近似式。\(\blacksquare\)
二次项¶
为刻画递增 / 递减的边际效应,常加二次项 \(y=\beta_0+\beta_1 x+\beta_2 x^2+u\)。偏效应随 \(x\) 变化:\(\dfrac{\partial y}{\partial x}=\beta_1+2\beta_2 x\)。令其为零得转折点:
\(\beta_2<0\) 时为先增后减的"倒 U"(峰),\(\beta_2>0\) 时为"正 U"(谷)。注意转折点是否落在数据范围内,影响解读。
交互项¶
\(y=\beta_0+\beta_1 x_1+\beta_2 x_2+\beta_3 x_1 x_2+u\) 中,\(x_1\) 的偏效应是 \(\beta_1+\beta_3 x_2\)——依赖 \(x_2\)。直接看 \(\beta_2\) 容易误读(它是 \(x_1=0\) 处的效应,常无意义)。重参数化为
则 \(\delta_2=\beta_2+\beta_3\mu_1\) 就是在 \(x_1\) 取均值处 \(x_2\) 对 \(y\) 的偏效应,更可解释。
拟合优度与回归元选择¶
调整 \(R^2\)¶
\(R^2\) 随变量增加只升不降,不能用于选模型。调整 \(R^2\) 对加变量施加自由度惩罚:
性质:加入一个新变量后,\(\bar R^2\) 上升 \(\iff\) 该变量的 \(|t|>1\);加入一组变量后上升 \(\iff\) 其联合 \(F>1\)。\(\bar R^2\) 可能为负。
用 \(\bar R^2\) 选非嵌套模型
\(F\) 检验只能比嵌套模型(一个是另一个的特例)。对非嵌套模型(如自变量用 \(\log x\) vs 用 \(x+x^2\)),可比 \(\bar R^2\) 择优。
易错点:因变量形式不同则 \(R^2\)/\(\bar R^2\) 不可比
比较的前提是因变量形式相同。一个模型用 \(y\)、另一个用 \(\log y\),两者拟合的是不同的因变量,\(R^2\) 与 \(\bar R^2\) 都不能直接比。
回归元的选择¶
别为了拼 \(R^2\) 而塞变量:过多控制变量会损害"其他条件不变"的解释。理想的追加变量是既影响 \(y\)、又与所关心的自变量无关的——它能在不引入多重共线性的前提下降低误差方差、提高精度。
预测与残差分析¶
预测值由 OLS 估计量算得,故带抽样波动,需给预测区间(比单点预测的置信区间更宽,因为还要叠加个体误差 \(u\) 的波动)。当因变量为 \(\log(y)\) 时,预测 \(y\) 不能简单取指数 \(e^{\widehat{\log y}}\)(会系统性低估),需乘一个由残差估计的修正因子。
速查表¶
| 专题 | 关键式 |
|---|---|
| \(\beta\) 系数 | \(\hat b_j=(\hat\sigma_j/\hat\sigma_y)\hat\beta_j\) |
| log-level 精确变化 | \(\%\Delta y=100(e^{\beta\Delta x}-1)\) |
| 二次项转折点 | $x^*= |
| 交互项偏效应 | \(\partial y/\partial x_1=\beta_1+\beta_3 x_2\) |
| 调整 \(R^2\) | \(\bar R^2=1-(1-R^2)\dfrac{n-1}{n-k-1}\) |
易错点与考点提醒:
高频陷阱
- 改单位 → 系数 / 标准误 / 置信区间端点成倍变,但 \(t\)、\(F\)、\(R^2\) 不变;
- 有交互项时不要单看主效应系数;有二次项时记得报转折点并看它是否在数据范围内;
- \(\bar R^2\) 加变量才惩罚,\(|t|>1\) 是它上升的临界;
- 因变量取了对数就别和水平模型比 \(R^2\);
- \(\log(y)\) 预测 \(y\) 要做指数修正,不能直接 \(e^{(\cdot)}\)。