时间序列之自相关函数

  • 相关系数
  • 自相关函数(Autocorrelation Function,ACF)[^1]
    • 样本自相关系数
  • 参考

相关系数

两个随机变量XXX和YYY的相关系数定义如下:
ρx,y=Cov(X,Y)Var(X)Var(Y)=E[(X−μx)(Y−μy)]E(X−μx)2E(Y−μy)2\begin{aligned} \rho_{x,y} &= \frac{Cov(X,Y)}{\sqrt{Var(X)Var(Y)}} \\ &= \frac{E[(X-\mu_x)(Y-\mu_y)]}{\sqrt{E(X-\mu_x)^2E(Y-\mu_y)^2}} \end{aligned} ρx,y​​=Var(X)Var(Y)​Cov(X,Y)​=E(X−μx​)2E(Y−μy​)2​E[(X−μx​)(Y−μy​)]​​
其中,μx,μy\mu_x,\mu_yμx​,μy​分别表示X,YX,YX,Y的均值,并假定方差是存在的。相关系数度量的是XXX和YYY线性相关的程度。相关系数具有以下几个性质:

  • −1≤ρx,y≤1-1 \leq \rho_{x,y} \leq 1−1≤ρx,y​≤1,说明相关系数处于−1-1−1到111之间;
  • ρx,y=±1\rho_{x,y} = \pm1ρx,y​=±1的充分必要条件是,XXX和YYY几乎处处有线性关系,即,存在a(≠0)a(\neq0)a(​=0)与bbb,使得P(Y=aX+b)=1P(Y=aX+b)=1P(Y=aX+b)=1。

当我们有样本{(xt,yt)}t=1T\{(x_t,y_t)\}_{t=1}^T{(xt​,yt​)}t=1T​,相关系数可以由样本相关系数估计出来,如下所示:
ρ^x,y=∑t=1T(xt−x‾)(yt−y‾)∑t=1T(xt−x‾)2∑t=1T(yt−y‾)2\hat{\rho}_{x,y} = \frac{\sum_{t=1}^T(x_t-\overline{x})(y_t-\overline{y})}{\sqrt{\sum_{t=1}^T(x_t-\overline{x})^2 \sum_{t=1}^T(y_t-\overline{y})^2}} ρ^​x,y​=∑t=1T​(xt​−x)2∑t=1T​(yt​−y​)2​∑t=1T​(xt​−x)(yt​−y​)​
其中,x‾,y‾\overline{x},\overline{y}x,y​分别是XXX和YYY的样本均值。

自相关函数(Autocorrelation Function,ACF)1

对于时间序列rtr_trt​,当我们考虑rtr_trt​与它的过去值rt−ir_{t-i}rt−i​的线性相依关系时,可以把相关系数的概念推广到自相关系数rtr_trt​与rt−lr_{t-l}rt−l​的相关系数称为rtr_trt​的间隔为lll的自相关系数,通常记为ρl\rho_lρl​,在弱平稳性的假定条件下,它只是关于lll的函数:
ρl=Cov(rt,rt−l)Var(rt)Var(rt−l)=Cov(rt,rt−l)Var(rt)=γlγ0\begin{aligned} \rho_l &= \frac{Cov(r_t,r_{t-l})}{\sqrt{Var(r_t)Var(r_{t-l})}} \\ &= \frac{Cov(r_t,r_{t-l})}{Var(r_t)} \\ &= \frac{\gamma_l}{\gamma_0} \end{aligned} ρl​​=Var(rt​)Var(rt−l​)​Cov(rt​,rt−l​)​=Var(rt​)Cov(rt​,rt−l​)​=γ0​γl​​​
这里用到了弱平稳性的性质Var(rt)=Var(rt−l)Var(r_t)=Var(r_{t-l})Var(rt​)=Var(rt−l​)。由定义2,我们有:

  • ρ0=1\rho_0=1ρ0​=1
  • ρl=ρ−l\rho_l=\rho_{-l}ρl​=ρ−l​
  • −1≤ρl≤1-1 \leq \rho_l \leq 1−1≤ρl​≤1

当且仅当对所有的l>0l>0l>0,都有ρl=0\rho_l=0ρl​=0,则称一个弱平稳序列是序列不相关的。

样本自相关系数

对一个给定的样本序列{rt}t=1T\{r_t\}_{t=1}^T{rt​}t=1T​,设 r‾\overline{r}r 是样本均值,则rtr_trt​的间隔为1的样本自相关系数
ρ^1=∑t=2T(rt−r‾)(rt−1−r‾)∑t=1T(rt−r‾)2\hat{\rho}_1 = \frac{\sum_{t=2}^T (r_t-\overline{r})(r_{t-1}-\overline{r})}{\sum_{t=1}^T (r_t-\overline{r})^2} ρ^​1​=∑t=1T​(rt​−r)2∑t=2T​(rt​−r)(rt−1​−r)​
在某些一般性条件下,ρ^1\hat{\rho}_1ρ^​1​是ρ1\rho_1ρ1​的相合估计。rtr_trt​的间隔为lll的样本自相关系数可定义为:
ρ^l=∑t=l+1T(rt−r‾)(rt−l−r‾)∑t=1T(rt−r‾)2,0≤l<T−1\hat{\rho}_l = \frac{\sum_{t=l+1}^T (r_t-\overline{r})(r_{t-l}-\overline{r})}{\sum_{t=1}^T (r_t-\overline{r})^2},0 \leq l < T-1 ρ^​l​=∑t=1T​(rt​−r)2∑t=l+1T​(rt​−r)(rt−l​−r)​,0≤l<T−1
若{rt}\{r_t\}{rt​}是一个独立同分布序列,满足E(rt2)<∞E(r_t^2)<\inftyE(rt2​)<∞,则对任意固定的正整数lll,ρ^l\hat{\rho}_lρ^​l​渐进的服从均值为0,方差为1T\frac{1}{T}T1​的正态分布。

参考


  1. Ruey S. Tsay.金融时间序列分析[D].王远林,王辉,潘家柱 译. ↩︎

  2. 时间序列之相关概念 ↩︎

时间序列之自相关函数相关推荐

  1. R语言绘制时间序列的自相关函数图:使用acf函数可视化时间序列数据的自相关系数图

    R语言绘制时间序列的自相关函数图:使用acf函数可视化时间序列数据的自相关系数图 目录

  2. python 自相关函数_用Python的statsmodels库计算时间序列的自相关函数和画图

    在时间序列分析课程中会需要用到自相关函数的计算,也就是当前期的值和滞后期的值之间的关系,这个指标的计算在计量软件中会比较容易实现,但是如果想要用python做怎么实现呢.代码如下:#导入库 impor ...

  3. x12arima季节调整方法_《时间序列X-12-ARIMA季节调整:原理与方法》

    对时间序列进行季节调整是经济分析的基础性工作.人民银行组织力量对季节调整基本方法进行了研究,结合调整中国特有的移动假日--春节的需要,提出了不同的处理模型,对各国较为通用的季节调整软件X-12-ARI ...

  4. 时间序列分析(2)| ARMA模型的(偏)自相关函数

    自相关函数和偏自相关函数是判断ARMA模型平稳性和阶数的有效工具. 1 自相关函数 根据平稳性的定义,与的方差相同,协方差仅与时间间隔有关,那么二者之间的相关系数也就仅与有关了.因为与对应的是同一个变 ...

  5. R语言ARIMA,SARIMA预测道路交通流量时间序列:季节性、周期性

    最近我们被客户要求撰写关于时间序列的研究报告,包括一些图形和统计输出. 相关视频:在Python和R语言中建立EWMA,ARIMA模型预测时间序列 本文从实践角度讨论了季节性单位根.我们考虑一些时间序 ...

  6. 算法模型---时间序列模型

    文章来源 时间序列 时间序列是时间间隔不变的情况下收集的不同时间点数据集合,这些集合被分析用来了解长期发展趋势及为了预测未来. 时间序列与常见的回归问题的不同点在于: 1.时间序列是跟时间有关的;而线 ...

  7. ARMA 时间序列模型

    更好的理解协方差以及相关系数 ###X因素和Y因素协方差公式: 自相关系数ACF 直观上来说,ACF 描述了一个观测值和另一个观测值之间的自相关,包括直接和间接的相关性信息. 其实自相关系数可以这么理 ...

  8. 【时间序列】怎么理解ACF 和PACF

    目录 1 概述 2 截尾与拖尾 3 Auto regressive (AR) process 4 Moving average(MA) Process 5 总结 1 概述 ACF 是一个完整的自相关函 ...

  9. 时间序列:理解好ACF自相关 与 PACF偏自相关

    ACF和PACF图的直观认识 先不说啥别的概念了,了解世界观不如了解方法论 自回归直观认识(intuition) 由自回归(AR)过程产生的滞后时间为k的时间序列. ACF描述了一个观测值与另一个观测 ...

最新文章

  1. 神经网络:你的大脑是一台计算机吗?
  2. python中变量名有哪些_Python变量范围有哪些?如何在Python中定义变量?
  3. CV:利用cv2自带两步法haarcascade_frontalcatface.xml实现对猫脸检测
  4. 李开复 6000 字长文,揭开 AI 创业的 10 个真相
  5. 【渝粤教育】电大中专就业指导作业 题库
  6. c++语言常用语句,c++语言词汇.doc
  7. 直播预告 | 腾讯千帆区块链云市场发布会震撼来袭
  8. 洛谷 P3375 【模板】KMP字符串匹配
  9. java对象说明 PO,BO,VO,DTO,POJO,DAO,Entity
  10. 职称计算机和英语有效性,2015职称计算机考试突破小诀窍
  11. 全局空间自相关算法:Join Count
  12. 《玩透嵌入式C的角角落落》当你需要循环体至少执行一次时,选择do
  13. 宏碁公司再次遭遇入侵,160GB敏感数据泄露
  14. 计算机共享网络热点,(传输)将win7计算机无线网络变成WiFi热点,使​​手机和笔记本电脑可以共享Internet...
  15. 在Word中引用参考文献
  16. Android集成PayPal
  17. 20 kubeadm高可用部署1.23以上k8s版本
  18. postman需要激活吗_postman接口测试实战例子
  19. 存储高性能、传输如闪电,焱融科技携手瑞云科技给你带来不一样的体验
  20. Openjudge NOI题库1.11编程基础之二分查找04 网线管理

热门文章

  1. 【个人建站心得】wordpress页面真正静态化插件really static
  2. Redhat6.5——解决yum功能不能正常使用
  3. 还在用图片验证码来刁难用户么?教你彻底去除图片验证码!
  4. 几十款免费好用的AI集合免费使用[⭐建议收藏⭐]
  5. 女生美容养颜喝什么茶好
  6. 《哈迪斯:杀出地狱》:超爽的地狱逃生之旅
  7. 【ORB_SLAM3源码解读】Ubuntu 16.04系统下编译运行 ORB_SLAM3
  8. 微信惊现任意代码执行漏洞 360手机卫士提供自检方案
  9. sql去重复操作详解SQL中distinct的用法
  10. cocos2dx经典魔塔游戏改编完整源码