基于ARIMA模型对某种重要化工原料价格变动的分析.docVIP

基于ARIMA模型对某种重要化工原料价格变动的分析 　　摘要：在分析过程中，首先利用SAS软件绘制时序图，从整体上把握数据信息。再进行数据的预处理，发现该化工原料价格序列具有非平稳非白噪声的属性，说明该序列有值得提取的相关信息，同时考虑差分方法使序列平稳。 　　然后分别采取1阶差分和2阶差分并将得到的结果进行对比，从中选出相对最优模型为ARIMA（2，1，0）。接着对所选取的模型进行参数估计以及模型检验，发现ARIMA（2，1，0）模型对该化工原料价格序列建模成功。 　　最后根据建立的模型对该化工原料的价格进行未来五个月的预测以期得到有价值的信息。 　　关键词：ARIMA模型；SAS软件；差分平稳；序列预测 　　一、描绘时序图 　　针对某种重要化工原料1990年至2004年的月度价格作出时序图，从图中可以发现，该序列波动范围较大，并不在某一个常数值附近波动，同时显示出一定的趋势性。 　　二、时序特征分析 　　（一）描述性统计 　　由SAS编程可得到化工原料月度价格序列描述性统计的值，其中均值为215.3616，标准差为60.29508，该序列共有177个观察值。 　　（二）自相关图检验 　　作出化工原料月度价格序列的自相关图，其中横轴表示自相关系数，纵轴表示延迟时期数，用水平方向的垂线表示自相关系数的大小。可以发现，序列的自相关系数递的速度相当缓慢，直至延迟24期自相关系数也未衰减到零。同时，自相关图显示序列自相关系数长期位于零轴的一边，说明序列具有单调趋势。由以上的分析可知该序列为非平稳序列。 　　（三）纯随机性检验 　　由SAS可得到序列的纯随机性检验结果： 　　延迟阶数LB统计量的值 P值 　　6554.08.0001 　　12731.53.0001 　　18783.1.0001 　　24801.66.0001 　　检验结果显示，在各阶延迟下LB检验统计量的P值都非常小（99.999%）认为该序列属于非白噪声序列。结合前面的平稳性检验结果，说明该序列为非平稳序列，同时还蕴含着值得提取的相关信息。 　　三、差分平稳及模型确定 　　由于该化工原料月度数据表现出非平稳性，为了进一步的分析，下面对序列进行差分处理。 　　（一）1阶差分平稳 　　1、1阶差分平稳序列分析 　　作图可以看出，1阶差分已经达到了很好的效果。从时序图可以发现，序列不存在明显的趋势性或周期性，没有显著的非平稳特征。 　　白噪声检验结果显示，在各阶延迟下LB检验统计量的P值都很小（98%）认为该序列属于非白噪声序列。再做出样本自相关图和样本偏自相关图，进一步确定模型平稳性并给拟合模型定阶。 　　从自相关图也可以看出自相关系数很快地衰减到零，延迟1阶后，自相关系数基本都在2倍标准差范围之内，故进一步确定该序列具有短期相关性，1阶差分序列平稳。再进一步考察自相关系数衰减到零的过程，可以看到自相关系数衰减到零并不是一个突然的过程，而是有一个连续渐变的过程，故认为自相关系数拖尾。 　　最后考察偏自相关系数衰减到零的过程，除了1～2阶偏自相关系数在2倍标准差范围之外，其他阶数的偏自相关系数都在2倍标准差范围内，这是偏自相关系数2阶截尾的典型特征。根据自相关系数拖尾，偏自相关系数2阶截尾的属性，初步确定拟合模型为AR（2）模型。 　　2、1阶差分序列定阶 　　（1）相对最优定阶。考虑人为定阶的主观性过强，利用SAS进行相对最优定阶，得到结果。最后一条信息显示，在自相关延迟阶数小于等于5，移动平均延迟阶数也小于等于5的所有ARIMA（p，1，q）模型中，BIC信息量最小的是ARIMA（3，1，0）。 　　（2）参数检验。由检验结果可知，常数项和φ3项无法通过显著性检验，去除常数项得到新的检验结果。 　　结果发现，去除常数项后φ3项仍无法通过显著性检验，而去除φ3后参数均通过显著性检验，故认为此模型为ARIMA（2，1，0）。 　　（二）2阶差分平稳 　　未避免差分方法选择不当对模型建立造成的影响，下面对序列进行2阶差分，并将得到的结果与1阶差分的结果进行对比，从而筛选出较优模型。 　　1、2阶差分平稳序列分析 　　对比1阶差分时序图，可以看到2阶差分与1阶差分得到的结果并没有显著的区别。从时序图可以发现，序列不存在明显的趋势性或周期性，没有显著的非平稳特征。 　　白噪声检验显示，在各阶延迟下LB检验统计量的P值都很小（99.999%）认为该序列属于非白噪声序列。再做出样本自相关图和样本偏自相关图，进一步确定模型平稳性并给拟合模型定阶。从自相关图也可以看到，延迟3阶后，自相关系数基本都在2倍标准差范围之内，故进一步确定该序列具有短期相关性，2阶差分序列平稳。 　　再进一步考察自相关系数衰减到零