基于WIN10的64位系统演示

一、写在前面

从这一期，我们使用SPSS进行SARIMA模型的构建。

同样，使用某省2005年1月至2016年12月AIDS死亡率的时间序列数据。

二、SPSS建立SARIMA实战

（1）录入数据和格式调整

双击打开IBM SPSS Statistics 27软件，把数据复制进SPSS界面：

数据格式调整：打开“变量视图”，按需调整，如“小数位数”、“数据类型”、“测量”等。

（2）时间序列数据转变

我们需要把“时间”这一列，转变成SPSS软件能识别的时间序列信息。

打开“定义日期和时间”面板，定义“年”和“月”还有“周期”：

点击确认后，可以看到。软件自动生成了“YEAR”、“MONTH”和“DATE”三列，说明时间序列转换完毕。

（3）自相关和偏相关图

虽然说只是看看而已，但是文章中还是得放的，所以还是得展示。

打开“自相关”对话框，先看不拆分的：

由图可知，在自相关图和偏相关图中，在12、24对应的柱子均处于高水平，由此可判断该序列具有季节性（周期为12个月）。另外，从自相关图中可以看到，均超过虚线，因此该序列为不平稳序列，应该进行拆分运算。

进行了拆分运算，分为一般拆分和季节拆分，进行一次一般拆分，则d = 1，进行一次季节拆分，则D = 1。看着图，就先这样吧。

（4）ARIMA建模

策略还是一样的：采用逐个参数尝试的方法，根据一般经验来看，疾病发病率ARIMA模型的p、q、P和Q选取的范围为0、1和2即可，因此需要尝试81（3×3×3×3=81）个ARIMA模型，再根据AIC和SBC准则以及参数是否通过假设检验来选出最佳的模型。

展示软件如何操作，这里以构建ARIMA(1,1,0)(0,1,2)12模型为例子：

结果解析：

由于有一个参数没差异，因此不考虑此模型，也可以填入下表，不然容易混乱。

继续下一个，ARIMA(1,1,1)(0,1,1)12模型：

按照这个流程，把81个模型全部撸一遍，找出最终的模型。

（5）最优ARIMA建模拟合和预测

假设ARIMA(1,1,1)(0,1,1)12是最优模型（我盲猜的，不过大概率是），如何得出拟合值和预测值呢？

先看拟合值，只需要在构建模型的时候勾选“保存”的选项，即可：

打开数据面板，可以发现多出了四列数据，复制到作图软件作图即可：

再看预测值，比如说我们要预测2017年1月到2017年12月的数据。先在数据面板填好这三列：

（6）评估模型的拟合和预测性能

计算mean square error (MSE)(均方差)， mean absolute error (MAE)(平均绝对误差) 和mean absolute percentage error (MAPE)(平均相对误差)公式如下：

其中，Xt为实际值，Xt^为预测值，n为总数。

计算以上误差可以在excel中实现，最后的结果如下：

	拟合误差			预测误差
	MAE	MAPE	MSE	MAE	MAPE	MSE
ARIMA	自己算	自己算	自己算	算不出	算不出	算不出

我们没有对已知数据进行拆分，所以只能算出拟合误差。

一般来说，2005年1月至2016年12月的数据，会用2005年1月至2015年12月的数据建模，去预测2016年1月至2016年12月的数据，这样就会得到拟合误差和预测误差了。

一般认为，MAPE小于0.1则预测效果优秀。

三、写在后面

（1）SPSS不能算单位根检验，而且SPSS的算法跟Eviews也有些区别，所以不同模型在同一个数据的建模，结果是不一样的。甚至不同版本都是这个情况，所以在写文章的时候，要注明软件和版本号！！！

（2）关于杨-博克斯Q(18)统计检验：它对模型中残差误差的随机检验，当其显著性小于0.05时说明残差误差不是随机的。在我们的例子中，它是小于0.05的，说明原始数据中有些信息还没提取出来。可以考虑进一步优化模型，比如调整拆分、自然对数处理等，实在不行，说明ARIMA处理不完全这个数据，那得用组合模型来处理了。毕竟ARIMA模型只能提取数据的线性信息。

四、数据

链接：https://pan.baidu.com/s/1WKNmpFA4ySE09tTSN8uvrw?pwd=wfyt

提取码：wfyt