Skip to content

偏自相关函数

Partial autocorrelation function · PACF

在剔除中间滞后最佳线性解释后衡量两个相隔时点剩余相关的函数。

领域
统计学
条目类型
定义

形式陈述

设零均值弱平稳过程具有自相关函数 ρ,且各有限滞后协方差矩阵正定。对 h2,分别把 XtXth 在线性空间 span(Xt1,,Xth+1) 上投影,记两端投影残差为 et(h)eth(h);偏自相关定义为

α(h)=Corr(et(h),eth(h)),

并约定 α(0)=1α(1)=ρ(1)α(h)=α(h)。它剔除的是中间变量的最佳线性影响,不是一般条件独立关系。

等价地,令 ϕh1,,ϕhh 解 Toeplitz 方程

(1ρ(1)ρ(h1)ρ(1)1ρ(h2)ρ(h1)ρ(h2)1)(ϕh1ϕhh)=(ρ(1)ρ(h)).

α(h)=ϕhh,即用最近 h 个值预测 Xt 时最远滞后的系数。正定性保证解和残差相关有定义;退化过程可能使矩阵奇异,此时普通 PACF 不唯一。

直觉

普通 ρ(2) 同时包含 Xt2Xt 的直接线性联系,以及经 Xt1 逐步传递的联系。PACF 先让中间滞后从两端各自解释能解释的部分,再比较剩余。因而对非退化、因果且为最小阶的 AR(p),超过 p 的旧值对下一期没有新增线性预测信息,总体 PACF 在 p 阶后截尾;普通 ACF 仍会因为反馈链条继续拖尾。

“偏”并不等于因果调整。投影只使用列出的滞后、只消除线性关系,也没有处理潜在共同原因。对联合高斯向量,零偏相关与相应精度矩阵元素为零可对应条件独立;离开高斯族,零偏相关只是一条二阶正交陈述。把 PACF 柱状图解释为某个滞后具有直接因果效应,超出了定义。

例子与边界

考虑 AR(2)

Xt=12Xt1+15Xt2+εt.

其 Yule–Walker 关系给 ρ(1)=5/8ρ(2)=41/80。因此一阶偏自相关是 5/8;二阶用两变量偏相关公式可复算

α(2)=ρ(2)ρ(1)21ρ(1)2=41/8025/64125/64=39/32039/64=15.

它恰好恢复第二个 AR 系数。对 h>2,加入 Xth 不再改善给定最近两期后的最佳线性预测,故总体 α(h)=0。这说明截尾位置是 2,不是说原序列在相隔三期后不相关。

有限样本的 α^(h) 几乎不会精确为零。对独立同分布白噪声常画 ±1.96/n 渐近参考带,但本文允许的仅不相关白噪声不自动满足该校准,逐步估计、同时检查许多滞后和先行选模也都会改变解释。若 Xt=Xt1 几乎处处,则相关矩阵全为 1、投影残差方差为零,PACF 在高阶没有定义;软件通过微小正则化给出的数字是算法约定,不是原过程的总体偏自相关。

推论与应用

Durbin–Levinson 递推从 ρ(1),,ρ(h) 逐阶得到 α(h) 与预测误差方差,利用 Toeplitz 结构把重复线性方程高效更新。正定自相关序列对应 |α(h)|<1;反过来,一列位于 (1,1) 的反射系数可逐阶参数化正定协方差,这在约束平稳 AR 参数和数值优化中很有用。

识别模型时,AR(p) 的 PACF 截尾与 MA(q) 的 ACF 截尾应作为总体基准,再结合信息准则、留出预测和残差诊断。近单位根会使样本 PACF、尤其第一阶接近 1,但不能仅凭该图判定单位根;结构突变也可造成类似持久性。PACF 适合提出阶数候选,不替代根检验和完整模型比较。

参考资料
  • Peter J. Brockwell and Richard A. Davis, Introduction to Time Series and Forecasting, 3rd ed., Springer, 2016,§3.2,PACF and the Durbin–Levinson algorithm。
  • James D. Hamilton, Time Series Analysis, Princeton University Press, 1994,§3.4,partial autocorrelations of autoregressions。
  • George E. P. Box, Gwilym M. Jenkins, Gregory C. Reinsel, and Greta M. Ljung, Time Series Analysis: Forecasting and Control, 5th ed., Wiley, 2015,Chs. 2 and 6,ACF/PACF identification tools。
关系图谱3 个相邻概念 · 1 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组