湖北省人均GDP的ARIMA模型及预测
- 格式:pdf
- 大小:308.22 KB
- 文档页数:6
湖北省人均GDP的ARIMA模型及预测
张鸽
摘要:本文介绍求和自回归移动平均模型ARIMA(p,d,q)的模型拟合方法,利用SAS统计软件实现模型
的拟合。采用时间序列分析方法,对湖北省1978-2013年人均GDP的数据进行分析。通过对数
据平稳性检验、模型参数检验、白噪声检验等分析,建立了ARIMA(1,1,0)时间序列模型,并对
未来十年的湖北省人均GDP数据进行预报。
关键词:人均GDP;ARIMA模型;预报
人均GDP是衡量一个国家或者地区经济发展水平和综合经济实力的重要指标,其计算公式是某地
区一定时期国内生产总值除以同期平均人口。在国际上它是划分一个区域经济发展阶段的依据之一,
人均GDP不仅考虑了经济总量的大小,而且结合了人口的因素,用它来反映区域经济增长和发展情况
会更加准确和富有现实意义。湖北省是我国中部崛起战略的重要支点,其经济体系在中部六省崛起战
略中具有典型性和引导性。因此对湖北省人均GDP的历史数据,进行研究分析,掌握其统计规律性和
变化趋势,为地方政策的制定提供经济依据。
近年来,已经有一些学者研究了人均GDP的发展规律及预测,SpencerLJames等[1]研究了210个
国家,从1950-2015年的人均GDP数据,建立了一个综合的时间序列模型;张毅[2]应用了ARIMA(1,1,2)
模型,对我国1952-2001年的人均GDP数据进行建模,并用过去人均GDP的值和过去误差来预测未来
人均GDP的走势;对湖北省人均DGP的分析也有一些研究,如:张波[3]利用Eviews软件,对1978-2004
年湖北省人均GDP历史数据,采用Box-Jenkins法建立ARMA(3,2)模型,并对此模型进行短期预报;
候大强[4]以2000-2010年的资料为依据,采用灰色系统理论方法[5-6],建立了GM(1,1)模型,揭示了湖
北省人均GDP增长变化的规律性。大多数学者选择使用SPSS或者Eviews软件进行分析,使用SAS来
分析的文章很少,本文采用1978-2013年的最新数据,应用权威的SAS软件,应用时间序列的方法,建
立了ARIMA(1,1,0)模型,对湖北省人均GDP数值进行分析,通过建立相关模型来探讨湖北省经济发展
状态和趋势,并对2014-2023年的人均GDP进行预测。
1、),,(qdpARIMA模型
ARIMA模型[7](IntegratedAutoregressiveMovingAverageModel)是研究时间序列的重要方法,
由自回归模型(AR模型)和滑动平均模型(MA模型)构成。简记为ARIMA(p,d,q),其中p,d,q分别表
示自回归阶数,差分阶数和移动平均阶数。具有如下结构
,,0)(,,0)(,)(,0)(,)(
2
tsXEtsEVarEBXB
tsstttttd
)(
其中,----1,----1221221qqppBBBBBBBB)()(t为零均值的白噪声序
列。
ARIMA模型的基本思想是:某些时间序列是依赖时间T的一组随机变量,构成时间序列的单个序列
虽然有不确定性,但整个序列的变化却有一定的规律性。该方法利用外推机制描述时间序列的变化,
达到最小方差下的最优预测,是一种精度较高的时间序列短期预测方法。
2、ARIMA模型的建模过程
任何非平稳时间序列只要通过适当阶的差分实现平稳,就可以对平稳后的序列进行ARMA(p,q)拟
合。
2.1对数据进行平稳性处理和检验
首先对给定的数据画出散点图(折线图)、自相关系数(ACF)图,判断序列的平稳性。如果序列
是非平稳的,可能有递增、递减或者季节性趋势,需要对数据进行差分或者对数差分,重复以上步骤,
使序列平稳,差分的次数就是ARIMA中d的值。对平稳的序列还需进行白噪声检验,白噪声序列则没
有分析的意义,对于平稳的非白噪声序列则进行ARMA(p,q)模型的拟合。
2.2模型的定阶
对于处理后的平稳序列进行分析,应用SAS程序的PROCARIMA过程,分析其自相关系数(ACF)
和偏自相关系数(PACF)的变化趋势并估计自相关阶数p和移动平均阶数q的值。选择模型的原则如
表:模型AR(p)MA(q)ARMA(p,d)
ACF拖尾截尾拖尾
PACF截尾拖尾拖尾
2.3模型的参数估计
估计方法有三种:矩估计、极大似然估计和最小二乘估计,在SAS中默认的是条件最小二乘估计。
参数检验就是要检验每个参数是否显著不为零,如有不显著的参数,则剔除不显著的参数对应的自变
量重新拟合模型。
假设条件:mjHHjj1,0:010:若拒绝原假设认为参数显著,否则认为参数
不显著。
2.4模型的检验
模型的检验是检验拟合的模型是否提取充分的信息,即残差的白噪声检验。如果不通过白噪声检
验,则重新拟合模型;如果是白噪声,则说明模型拟合有效。模型的有效性检验是LB统计量对残差序
列进行卡方检验。
2.5模型的预测
通过ARIMA过程拟合模型,利用SAS程序的FORECAST过程,可以预报后期的数据并分析发展趋势。
3、实例分析
下面给出湖北省1978-2013年人均GDP的数据(来自国家统计局湖北省人均GDP统计公报)年份人均GDP/元年份人均GDP/元年份人均GDP/元
19783221990155620028319
19794091991168920039011
198042819921986200410190
198146619932565200511554
198250619943341200613360
198354919954162200716386
198467619965122200819858
198580819975899200922677
198689119986289201027906
1987103119996514201134197
1988122820007188201238572
1989138320017813201342613
在SAS中,使用GPLOT过程作出人均GDP序列的散点图,
图1:人均GDP序列图图2:人均GDP的对数差分序列图由图1可以看出人均GDP随时间有明显的增长趋势,说明该序列不平稳。对序列取对数并作一次
差分,再作变换后的数据随时间图(图2),可以初步判断变换后的序列为平稳序列。
使用ARIMA过程对对数差分后的人均GDP作纯随机性及平稳性检验,见下表3.
表3.对数差分后人均GDP序列的纯随机和平稳性检验
Autocorrelations
LagCovarianceCorrelation-198765432101234567891
00.00373461.00000||********************|0
10.00177770.47601|.|**********|0.169031
20.000497850.13331|.|***.|0.203762
30.000074760.02002|.|.|0.206239
4-0.0007412-.19846|.****|.|0.206294
5-0.0009891-.26484|.*****|.|0.211679
6-0.0009629-.25782|.*****|.|0.220943
7-0.0011399-.30523|.******|.|0.229378
8-0.0007310-.19572|.****|.|0.240703
9-0.0003456-.09255|.**|.|0.245208
10-0.0002713-.07265|.*|.|0.246204
11-0.0001902-.05092|.*|.|0.246816
12-0.0004363-.11682|.**|.|0.247116
AutocorrelationCheckforWhiteNoise
ToChi-Pr>
LagSquareDFChiSq--------------------Autocorrelations--------------------
616.9860.00930.4760.1330.020-0.198-0.265-0.258
1224.74120.0161-0.305-0.196-0.093-0.073-0.051-0.117
分析可知延迟6阶,延迟12阶后的0.05p,故序列为非白噪声。该序列的自相关系数迅速
衰减为零,故该序列是平稳的。
下面拟合模型的阶数,应用ARIMA过程的IDENTIFY命令的)5:0()5:0(minicqp进行最优模
型定阶。
表4.相对最优模型定阶
MinimumInformationCriterion
LagsMA0MA1MA2MA3MA4MA5
AR0-5.89756-6.23114-6.16712-6.09183-6.06788-6.31539
AR1-6.40211-6.30316-6.41791-6.3177-6.36466-6.45062
AR2-6.47271-6.43031-6.35197-6.27357-6.27784-6.36434
AR3-6.38285-6.35229-6.25276-6.1758-6.1785-6.26672
AR4-6.45781-6.39069-6.31545-6.22816-6.13903-6.16744
AR5-6.38408-6.29215-6.22969-6.16055-6.05905-6.12469
Errorseriesmodel:AR(7)
MinimumTableValue:BIC(2,0)=-6.47271
上表可以看出在2p时,BIC函数值最小。应用ARIMA过程的2pestimate做参数检验,
结果如下
表5.参数的显著性检验
ConditionalLeastSquaresEstimation
StandardApprox
ParameterEstimateErrortValuePr>|t|Lag
MU0.143870.015949.03<.00010
AR1,10.547960.176523.100.00401
AR1,2-0.125780.17673-0.710.48182
分析上表,存在参数不显著,说明未通过参数检验。下面检验1p和1q情形,
表6.)1(),1(MAAR模型参数检验与模型检验结果
ConditionalLeastSquaresEstimation
StandardApprox
ParameterEstimateErrortValuePr>|t|Lag
MU0.143450.017518.19<.00010
AR1,10.486850.153663.170.00331
ConstantEstimate0.07361
AIC-101.498
ConditionalLeastSquaresEstimation
StandardApprox
ParameterEstimateErrortValuePr>|t|Lag
MU0.143520.0140010.25<.00010
MA1,1-0.530880.14804-3.590.00111
AIC-101.941从上面2个模型的检验可以看到,它们均为有效模型,本文为方便,根据AIC准则,AR(1)模型的
AIC比MA(1)的值更小,所以选择前者作为预测模型。