机器学习——时间序列ARIMA模型(四)：自相关函数ACF和偏自相关函数PACF用于判断ARIMA模型中p、q参数取值

文章目录

1、自相关函数ACF
2、偏自相关函数PACF
3、ARIMA(p,d,q)的阶数判断
4、代码实现
- 1、引入所需依赖
- 2、数据读取与处理
- 3、一阶差分与绘图
- 4、ACF
- 5、PACF

1、自相关函数ACF

自相关函数反映了同一序列在不同时序的取值之间的相关性。

公式：
ACF(k)=ρk=Cov(yt,yt−k)Var(yt)ACF(k)=\rho_{k}=\frac{Cov(y_{t},y_{t-k})}{Var(y_{t})} ACF(k)=ρk=Var(yt)Cov(yt,yt−k)
其中分子用于求协方差矩阵，分母用于计算样本方差。求出的ACF值为[-1,1]。

但对于一个平稳的AR模型，求出其滞后值的自相关系数后，得到的并不是x(t)与x(t-k)之间单纯的相关关系，其中掺杂了其他随机变量的影响

2、偏自相关函数PACF

偏自相关系数PACF是严格这两个变量(现值和滞后值)之间的相关性，剔除了中间随机变量的干扰。

3、ARIMA(p,d,q)的阶数判断

q看ACF，p看PACF,d看差分阶数

4、代码实现

1、引入所需依赖

%matplotlib inline
import pandas as pd
import datetime
import matplotlib.pylab as plt
import warnings
import seaborn as sns
import statsmodels
from matplotlib.pylab import style
from statsmodels.tsa.arima_model import ARIMA
from statsmodels.graphics.tsaplots import plot_acf, plot_pacfstyle.use('ggplot')
plt.rcParams['font.sans-serif'] = ['SimHei']
plt.rcParams['axes.unicode_minus'] = False

2、数据读取与处理

还是读取道琼斯数据集，这里就不再说啦，这次以周为单位。

详情请看机器学习——时间序列ARIMA模型(三)：AR、MA、ARMA、ARIMA模型定义及公式介绍及股价预测案例代码_天海一直在的博客-CSDN博客

stockFile = 'data2000_2021.csv'
stock = pd.read_csv(stockFile, index_col=0, parse_dates=[0])#将索引index设置为时间，parse_dates对日期格式处理为标准格式。
stock.head(10)
stock_week = stock['close'].resample('W').mean()
stock_train = stock_week['2000':'2019']

3、一阶差分与绘图

stock_train.plot(figsize=(9,5))
plt.legend(bbox_to_anchor=(1.25, 0.5))
plt.title("Stock Close")
sns.despine()
#一阶差分
stock_diff = stock_train.diff()
stock_diff = stock_diff.dropna()
plt.figure()
plt.plot(stock_diff)
plt.title('一阶差分')
plt.show()

4、ACF

acf = plot_acf(stock_diff, lags=20)
plt.title("ACF")
acf.show()

置信区间内q取1，后面突出的两点判断为离群奇异值不管

5、PACF

pacf = plot_pacf(stock_diff, lags=20)
plt.title("PACF")
pacf.show()

在p取9时后面衰减趋于0且在置信区间内，所以p取9
后续可以根据BIC贝叶斯信息准则来验证，结果和此判断一致。

机器学习——时间序列ARIMA模型(四)：自相关函数ACF和偏自相关函数PACF用于判断ARIMA模型中p、q参数取值相关推荐

分类模型-评估指标（2）：ROC曲线、 AUC值（ROC曲线下的面积）【只能用于二分类模型的评价】【不受类别数量不平衡的影响；不受阈值取值的影响】【AUC的计算方式：统计所有正负样本对中的正序对】
评价二值分类器的指标很多,比如precision.recall.F1 score.P-R曲线等.但这些指标或多或少只能反映模型在某一方面的性能.相比而言,ROC曲线则有很多优点,经常作为评估二值分类器 ...
自相关(ACF)与偏自相关(PACF)（4）
§5.偏自相关系数PACF 在AR(1)AR(1)AR(1)模型中,即使yt−2y_{t-2}yt−2没有直接出现在模型中,但是yty_tyt和yt−2y_{t-2}yt−2之间也相关,偏相关系 ...
自相关系数 ACF与偏自相关系数PACF，拖尾和截尾
1.ACF y(t,s)=E(Xt-µt)(Xs-µs) 定义ρ(t,s)为时间序列的自相关系数,为ACF ρ(t,s)=y(t,s)/sqrt(DXt * DXs) E为期望,D为方差 2.PACF ...
自相关ACF与偏自相关PACF
在做机器学习预测时,往往需要对数据进行auto-correlation和partial auto-correlation function(ACF和PACF),根据查阅文献和相关资料,整理出比较容易上 ...
此UML形状所在的绘图页不是UML模型图的一部分。该形状合计用于利用UML模型图模板创建的绘图。解决办法
使用visio编辑类图时出现的问题解决方法说明:当出现上述问题时,检查是否是在UML模式下编辑的.如果不是请重新创建为UML模型即可!!! 检查方法如下解决方案 ①步骤 ②步骤 ③效果图
arma模型谱估计matlab_基于机器学习的心律失常分类(四)——心电信号特征提取[MATLAB]...
目前比较常用的特征提取是提取心电信号的各波形间期长度.波峰高度等,本文是使用ARMA模型对心电信号进行处理,使用其系数来作为特征. 一.心拍划分大多伴有异常波形的心律失常信号,通常都会具体表现在单个 ...
自相关函数和偏自相关函数（ACF、PACF）
前言简单介绍一下自相关函数ACF和偏自相关函数 PACF的计算和图像显示,概念性的东西可以查阅其他资料了解. 实例数值计算 import statsmodels.tsa.api as smttim ...
ARIMA-GRNN模型的发病率预测GUI：基于Matlab编写（ARIMA部分）
该软件首先通过经典的ARIMA模型得出初步的预测数值,生成绝对误差序列,然后输入GRNN模型得出预测的绝对误差数值,最后通过反算生成最终预测数值. 第一步.Matlab支持插件的安装 1. 双击MCR ...
c4.5决策树算法python_Python3实现机器学习经典算法（四）C4.5决策树
一.C4.5决策树概述 C4.5决策树是ID3决策树的改进算法,它解决了ID3决策树无法处理连续型数据的问题以及ID3决策树在使用信息增益划分数据集的时候倾向于选择属性分支更多的属性的问题.它的大部分 ...