sas分析报告非平稳时间序列

《sas分析报告非平稳时间序列》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、实用标准文案运用SAS对谷物产量进行分析一、摘要利用SAS软件（程序见附录）判断谷物产量数据为平稳序列且为非白噪声序列，然后先后通过模型的识别、参数的估计、模型的优化、残差白噪声检验，确定AR(1)模型拟合时间序列显著有效。由于时间序列之间的相关关系，且历史数据对未来数据有一定的影响，对未来5期的谷物生产量进行预测。二、理论准备首先判断序列的随机性和平稳性。通过随机性检验，判断该序列是否为白噪声序列，如果是白噪声序列，就认为该随机事件没有包含任何值得提取的有用信息，我们就应该终止分析。通过平稳性检验，序列可以分为平稳序列和非平稳序列。如果序列平稳，通过相关计算进行模型拟合，并利用过去

2、行为对将来行为进行预测，达到预测效果。如果序列为非平稳，再确定模型为非平稳序列中四大类模型中的哪种种模型或者几种模型对序列的综合影响，通过把序列转化为平稳序列，再进一步分析。三、数据选取本实验采用某地区连续74年的谷物产量（单位：千吨），如下所示：0.970.451.611.261.371.431.321.230.840.891.181.331.210.980.910.611.230.971.100.740.800.810.800.600.590.630.870.360.810.910.770.960.930.950.650.980.700.861.320.880.680.781.2

3、50.791.190.690.920.860.860.850.900.540.321.401.140.690.910.680.570.940.350.390.450.990.840.620.850.730.660.760.630.320.170.46四、数据进行平稳性与纯随机性的检验与判别（一）序列的纯随机性检验图1序列延迟6阶LB检验结果序列纯随机性检验结果显示延迟6阶LB检验统计量的P值小于1%的显著性水平0.0001，说明序列之间蕴含着很强的相关信息，即该序列是非随机性序列，为非白噪声。精彩文档实用标准文案（二）模型的平稳性检验——ADF检验绘制该序列时序图：图2序列时序图该时

4、序图显示该地区谷物产量在0.8千吨左右，波动比较平稳，但存在略微趋势。运用单位根检验序列的平稳性。序列时序图得到如下图：图3ADF检验图ADF检验结果表示，单位根统计量ADF=-2.682929，在5%的显著性水平下，接受γ=0的原假设，即认为该序列非平稳且带有趋势。为此对序列进行一阶差分（=▽=-），画出的时序图如下精彩文档实用标准文案由时序图可以看出经过一阶差分后的序列平稳。进一步对差分后的序列进行单位根检验，如图：图4一阶差分后的ADF检验图从图中可以看出，在5%的显著性水平下，拒绝γ=0的假设，即通过一阶差分后，序列基于平稳。下面便对一阶差分序列进行分析。五、模型的识别令=▽

5、=-画出差分序列的自相关图和偏自相关图如下：精彩文档实用标准文案图5样本自相关图样本自相关图显示除了延迟一阶自相关系数在二倍标准差范围以外，其他阶数的自相关系数都在2倍标准差范围内波动。根据自相关系数的这一特点可以初步判断该序列具有短期相关性，进一步确定序列平稳。同时，可以认为该序列的自相关系数1阶截尾。图6样本偏自相关图精彩文档实用标准文案样本偏自相关图显示除了延迟1、2、3阶的偏自相关系数显著大于2倍标准差之外，其他的偏自相关系数都在2倍标准差范围内作小值随机波动，所以该偏自相关系数可认为3阶截尾。为了模型定阶的准确性，下面采用BIC准则，对模型进行判断。根据BIC准则，我们知道

6、在所有通过检验的模型中使得BIC函数达到最小的模型为相对最优模型。为此运用SAS输出所有自相关延迟阶数小于等于10，移动平均延迟阶数小于等于5的所有ARMA(p，q)模型的BIC信息量，如下图所示：图7ARMA(p，q)模型的BIC信息量图表显示一阶差分序列的BIC信息量最小的是MA(1)模型。BIC统计量值为-2.8997。可得该模型为MA(1)模型。六、模型的参数估计由BIC准则确定，该模型为MA(1)模型,运用SAS的ESTIMATE命令输出未知参数估计结果，如下图所示：图8模型参数估计参数显著性结果显示常数项的参数t统计量的P值均大于0.05，即参数不显著。然后通过更改命令，

7、去掉常数项，得到的未知参数估计结果，如下图所示：精彩文档实用标准文案图9未知参数估计结果此时参数显著性结果显示参数t统计量的P值小于0.05，即参数显著。七、模型残差的白噪声检验ESTIMATE命令输出的残差白噪声检验图如下所示：图10残差白噪声检验残差白噪声检验显示延迟6阶，延迟12阶，延迟18阶，延迟24阶LB检验统计量的P值均显著大于0.05，说明模型的残差为白噪声序列，即该MA（1）模型显著有效。八、模型的确定（一）、ESTIMATE命令输出的序列