普通OLS在时间序列回归中失效，因残差自相关违反误差独立同分布假设，导致标准误低估、t统计量虚高、p值失真；需用Newey-West校正或Prais-Winsten等方法处理，并检验平稳性与内生性。
为什么普通OLS会失效
时间序列回归里，如果残差存在高度自相关（比如 AR(1) 系数接近 0.9），statsmodels.OLS 估计出的系数标准误会严重低估，t 统计量虚高，p 值失真——不是模型拟合得好，而是统计推断完全不可信。
根本原因在于 OLS 假设误差项独立同分布，而自相关直接违反这一假设。常见现象包括：模型 R² 很高但预测不稳定、残差图呈现明显趋势或周期性、durbin_watson 检验值远偏离 2（比如 1.2 或 0.8）。
用 Newey-West 校正标准误（最简实用方案）
不改模型结构，只修正推断——适合快速验证变量是否真有解释力。核心是用异方差自相关稳健（HAC）协方差估计替代默认协方差矩阵。
- 用
statsmodels的get_robustcov_results方法，指定cov_type='HAC'和滞后阶数maxlags -
maxlags建议设为int(4*(nobs/100)**(1/4))（Newey-West 默认公式），也可试1到5观察 t 值变化 - 注意：该方法不提升预测能力，只让 p 值可信；若原始 OLS 系数本身已偏移（如遗漏动态项），校正后仍可能有偏
import statsmodels.api as sm model = sm.OLS(y, X) result = model.fit() robust_result = result.get_robustcov_results(cov_type='HAC', maxlags=3) print(robust_result.summary())
改用 Prais-Winsten 或 Cochrane-Orcutt（显式建模 AR(1) 误差）
当自相关主要由一阶自回归主导时，这两者比 Newey-West 更进一步：先估计误差的 rho，再对数据做广义差分，使变换后误差近似白噪声。
立即学习“Python免费学习笔记（深入）”；
-
statsmodels.tsa.arima_model.ARIMA已弃用；推荐用statsmodels.tsa.statespace.SARIMAX直接建模y ~ X + AR(1) - 若坚持两步法，可用
statsmodels.stats.stattools.cochrane_orcutt（需手动迭代，易收敛失败） - Prais-Winsten 保留首期观测，Cochrane-Orcutt 丢弃它；小样本下前者效率略高，但两者在
n > 50时差异很小 - 务必检查变换后残差的
durbin_watson是否落在1.5–2.5区间，否则说明 AR(1) 假设不足，需尝试 AR(2) 或 ARCH 类模型
警惕“伪回归”和内生性陷阱
高度自相关常伴随单位根或共同趋势，此时即使两个无关序列也可能产生极高 R² 和显著 t 值——这不是模型好，是 spurious regression。
- 先对每个变量做 ADF 检验（
statsmodels.tsa.stattools.adfuller），确认平稳性；非平稳序列必须差分或协整处理 - 若
y和X同阶单整且存在协整关系，应使用statsmodels.tsa.vector_ar.vecm.Vecm或 Engle-Granger 两步法 - 时间序列中
X若含滞后因变量（如y_{t-1}），OLS 会产生小样本偏差，此时必须用工具变量（如statsmodels.sandbox.regression.gmm.GMM）或系统 GMM
实际建模时，自相关程度、平稳性、外生性这三项得逐个排查，漏掉任何一项，结果都可能完全误导决策。
