Python ARIMA 代码 arima模型 python

1.项目背景

      当今世界正处于一个数据信息时代，随着后续互联网的发展各行各业都会产生越来越多的数据，包括不限于商店、超市、便利店、餐厅等等。那么这里面很多数据都是随着时间产生的，这就形成了时间序列数据，而且很多时间序列数据都是非平稳时间序列数据。目前对非平稳时间序列分析应用最多的模型就是ARIMA模型，本项目也是通过Python程序来进行数据探索性分析、数据预处理、构建ARIMA时序模型以及如何把模型预测的结果应用于日常的实际生活当中。

2.数据获取

      本次建模数据来源于某公司餐厅销售数据，其主要的数据基本统计概况如下：

数据项：2项(日期、销量)，如下图所示

 通过观察数据，发现一些空值和一些重复数据，后续章节进行处理。

3.数据预处理

       真实数据中可能包含了大量的缺失值和噪音数据或人工录入错误导致有异常点存在，非常不利于算法模型的训练。数据清洗的结果是对各种脏数据进行对应方式的处理，得到标准的、干净的、连续的数据，提供给数据统计、数据挖掘等使用。数据预处理通常包含数据清洗、归约、聚合、转换、抽样等方式，数据预处理质量决定了后续数据分析挖掘及建模工作的精度和泛化价值。

3.1数据描述性统计与清洗

（1）导入程序库：将所用到的程序库导入到Python程序中，如图所示。

（2）读取数据：使用Pandas库中read_excel方法读取Excel数据，并转为DataFrame类型。读取数据代码如图所示：

3）查看统计信息和空值，代码如下：

 

       从上图可以看到，数据总数据量40条(不包括空值)，以及未去除重复值的平均值、标准差、最小值、最大值、分位数。另外可以看到销量有2条空值。

（4）查看是否有重复数据以及清理重复数据，代码如下：

从上图中的部分展示，可以看到是有重复数据，所以使用drop_duplicates()函数进行了重复数据清理。

（5）空值清理，本项目涉及2条空值数据 直接进行丢弃操作。空值的处理也可以用均值填充、分位数填充、拉格朗日插值填充等等(根据具体数据处理的需要进行即可)。代码如下：

（6）针对清洗后的数据进行统计分析，代码如下：

 从上图可以看出，清洗后数据量为37条，平均值为3542.49、标准差为543.92。

4.探索性数据分析

4.1数据分析

销量时间序列分析：

 

关键代码如下：

5.构建ARIMA时序模型

5.1 ARIMA模型概念

      在自然界中绝大部分序列都是非平稳的。因而非平稳时间序列的分析更普遍、更重要，创造出来的分析方法也更多。

      非平稳时间序列的分析方法可以分为确定性因素分解的时序分析和随机时序分析两大类。

      确定性因素分解的方法把所有序列的变化都归结为4个因素（长期趋势、季节变动、循环变动和随机波动）的综合影响，其中长期趋势和季节变动的规律性信息通常比较容易提取，而由随机因素导致的波动则非常难以确定和分析，对随机信息浪费严重会导致模型拟合精度不够理想。

     随机时序分析法的发展就是为了弥补确定性因素分解方法的不足。根据时间序列的不同特点，随机时序分析可以建立的模型有ARIMA模型、残差自回归模型、季节模型、异方差模型等。

本项目重点介绍ARIMA模型对非平稳时间序列进行建模。

差分运算：

1. p阶差分

相距一期的两个序列值之间的减法运算称为1阶差分运算。

1. k步差分

相距k期的两个序列值之间的减法运算称为k步差分运算。

       差分运算具有强大的确定性信息提取能力，许多非平稳序列差分后会显示出平稳序列的性质，这时称这个非平稳序列为差分平稳序列。差分平稳序列可以使用ARMA模型进行拟合。ARIMA模型的实质就是差分运算与ARMA模型的组合。

 5.2序列平稳性检验

图5.2-1 原始序列的时序图

此图显示该序列具有明显的单调递增趋势，可以判断为非平稳序列。

 

图5.2-1 原始序列的自相关图

此图的自相关图显示自相关系数长期大于零，说明序列间具有很强的长期相关性。

表5.2-1　原始序列的单位根检验

adf

cValue

p值

1%

5%

10%

1.8138

-3.7112

-2.9812

-2.6301

0.9985

      此表单位根检验统计量对应的p值显著大于0.05，最终将该序列判断为非平稳序列（非平稳序列一定不是白噪声序列）。

5.3对原始序列进行一阶差分，并进行平稳性和白噪声检验

1）对一阶差分后的序列再次做平稳性判断。

 

                     图5.3-1　一阶差分之后序列的时序图 

 

图5.3-2　一阶差分之后序列的自相关图

表5.3-1 一阶差分之后序列的单位根检验

结果显示，一阶差分之后序列的时序图在均值附近比较平稳地波动、自相关图有很强的短期相关性、单位根检验p值小于0.05，所以一阶差分之后的序列是平稳序列。

2）对一阶差分后的序列做白噪声检验

表5.3-2　一阶差分之后序列的白噪声检验

stat

p值

11.3040

0.0007734

 此表输出的p值远小于0.05，所以一阶差分之后的序列是平稳非白噪声序列。

5.4对一阶差分后的序列拟合ARMA模型

下面进行模型定阶。模型定阶就是确定p和q。

第一种方法：人为识别，根据图5-25进行模型定阶

 

图5.4-1　一阶差分后序列的偏自相关图

一阶差分后自相关图显示出1阶截尾，偏自相关图显示出拖尾性，所以可以考虑用MA(1)模型拟合1阶差分后的序列，即对原始序列建立ARIMA(0，1，1)模型。

第二种方法：相对最优模型识别。

计算ARMA(p，q) 当p和q均小于等于3的所有组合的BIC信息量，取其中BIC信息量达到最小的模型阶数。

计算完成BIC矩阵如下：

当p值为0、q值为1时，最小BIC值为422.510082。p、q定阶完成。

        用AR(1)模型拟合一阶差分后的序列，即对原始序列建立ARIMA(0，1，1)模型。虽然两种方法建立的模型是一样的，但模型是非唯一的。ARIMA(1，1，0)和ARIMA(1，1，1)这两个模型也能通过检验。下面对一阶差分后的序列拟合AR(1)模型进行分析:

1.参数检验和参数估计见下表:

表5.4-1 模型参数

Parameter

Coef.

Std.Err.

t

const

49.9561

20.1390

2.4806

ma.L1.D.销量

0.6710

0.1648

4.0712

2.模型检验之残差检验：

 

                                        图5.4-2　残差自相关图 

                                      图5.4-3　残差偏自相关图 

                                            图5.4-4　残差偏自相关图

D-W检验：

DW值显著的接近于0或４时，则存在自相关性，而接近于２时，则不存在（一阶）自相关性。

D-W检验的结果为：1.9735。因此，模型预测的残差不存在自相关性性，这说明拟合的模型预测效果很好。

残差序列的白噪声检验结果为：(array([0.00390439]), array([0.95017654]))，可以看出p值为0.9502， P>0.05.说明原假设成立。

5.5 ARIMA模型预测

应用ARIMA(0，1，1)对表某餐厅的销售数据做为期5天的预测，结果如下表所示：

表5.5-1 预测未来5天的销售额

需要说明的是，利用模型向前预测的时期越长，预测误差将会越大，这是时间预测的典型特点。

6.实际应用

通过上面的预测结果分析得到，此模型能够在误差允许的范围内进行预测未来的销量，可以将此模型应用于实际生产中，帮助企业更好地备货。