预测课程论文.doc

1、中国进出口总额时间序列的分析及其预测摘要：本文采用时间序列分析方法对我国19502009年的年度进出口数据进行分析，运用SAS软件构建了ARIMA模型对数据进行拟合，结果表明该模型具有良好的预测结果，故可运用该模型对我国20102014年的进口总额、出口总额进行点预测，同时本文给出了95%的概率意义下的短期预测区间，进而充分提高了该模型的预测精度，预测结果的可靠性得到了有效的保证。关键词：时间序列分析；进出口总额；ARIMA模型；预测一、引言贸易是一国经济增长的引擎机，而进出口贸易则是一国贸易中的重要组成部分1。自加入WTO以来，中国经济融入世界经济一体化进程有所增加，进出口贸易在我国国民经济

2、中的地位也越来越重要。进出口贸易为充分利用国内国际两种资源、两个市场、调剂国内市场余缺，保持国民经济平稳较快发展做出了重要贡献，因此政府需要对进出口贸易进行适当的调节，而正确的调节策略是建立在准确信息基础上的。准确的进出口预测能大大增强政府制定符合国际规则的对外贸易政策的预见性。本研究运用时间序列分析方法对19502009年我国进出口贸易总额的变动进行分析，建立适当的模型对数据进行拟合，并运用模型对进出口总额进行预测，为进行相关的经济决策提供有效的依据。二、文献综述进出口的研究是一个非常有意义的课题，中国科学院预测科学研究中心对此设立了专门的研究课题组，并发表了2010年我国进出口形势分析与预

3、测一文。该文通过对国内外经济形势及海关统计数据的分析，对我国进出口进行了测算。然而该文仅重点分析了2010年我国对外贸易中需要关注的宏观经济问题，并未对数据的分析处理给出具体方法。在对进出口进行分析的过程中，逐步回归法也曾被用于模型的构建。王红、童恒庆等人应用逐步回归理论（2006）对进出口贸易进行建模研究。该模型充分考虑了对进出口有着重大影响的财政支出总额、经济建设支出、社会文教支出、国防支出、行政管理支出等因素的影响，采用逐步回归的方法进行建模。根据逐步回归的结果，再分别进行过三次逐步回归，然后通过综合比较每次逐步回归的模型评价指标和拟合效果，选择误差低，拟合效果好的模型作为进出口的预测模

4、型。模型很好的对进出口进行了拟合，但是由于指标的数目较多，对数据的要求也相应提高，增加了预测的难度。同时，在指标的选择上不可避免的带有主观性，不利于客观的对进出口进行预测。进出口市场是一个受多种因素影响的复杂系统，罗文彬、王海燕等人（2006）通过研究发现其时间序列中具有确定性混沌特征，因此运用在证券等经济时间序列分析中得到的广泛应用的混沌时间序列方法分析进出口时间序列，从而预测未来进出口的具体走势。该方法有效的对进出口作出了预测，但是由于其操作过程较为复杂，不易于理解掌握。在众多关于进出口总额的过程中，建立ARIMA模型无疑是一种易于操作并且行之有效的方法。ARIMA模型在经济预测过程中既考

5、虑了经济现象在时间序列上的依存性，又考虑了随机波动的干扰性，对于经济运行短期趋势的预测准确率较高，因此本文通过对数据进行分析，建立ARIMA模型以对进出口作出预测。在建模的过程中，模型的定阶是最关键的问题。文章运用了SAS软件进行数据分析及建模，SAS系统提供的相对最优模型识别可以对模型进行客观的定阶，一定程度上可避免人为定阶带来的主观误差，并且提高模型预测的准确度。 三、数据来源本文选用我国19502009年年度进出口总额的统计数据进行分析，数据来源于中国统计年鉴2010，经整理后见表1。表1 19502009年年度进出口额 （单位：亿美元） 年份进口出口年份进口出口19505.85.519

6、80199.41180.991951127.61981220.14220.07195211.28.21982192.85223.21195313.510.21983213.9222.26195412.911.51984274.1261.39195517.314.11985422.52273.5195615.616.51986429.04309.42195720.3122.141987432.16394.37195825.0627.261988552.68475.16195928.9231.721989591.42525.38196026.4825.711990533.45620.9119611

7、7.4719.421991637.91719.1196213.7319.131992805.85849.4196314.520.3119931039.59917.44196417.122.519941156.151210.06196522.4625.6319951320.841487.8196624.8226.8119961388.331510.48196721.6923.8819971423.71827.92196820.6823.419981402.371837.12196919.1724.2919991656.991949.31197022.7923.0720002250.942492.

8、03197121.2927.8320012435.532660.98197228.5136.9320022951.73255.96197352.0858.7620034127.64382.28197477.9171.0820045612.295933.26197579.2676.8920056599.537619.53197666.669.4320067914.619689.78197771.4875.220079561.1512200.61978111.3199.55200811325.6714306.931979156.21136.14200910056.8812015.34四、ARIMA

9、模型简介ARIMA模型是由统计学家Box和Jenkins提出的，又称B-J模型，其所依赖的原理是：某些时间序列是依赖于时间t的一组随机变量，构成该时序的单个序列值虽然具有不确定性，但整个序列的变化却有一定的规律性，可以用相应的数学模型(即ARIMA)近似描述。通过对该数学模型的分析研究，能够从本质上认识时间序列的结构与特征，达到最小方差意义下的最优预测。(一)ARIMA(p，d，q)模型的形式2。ARIMA模型有四种基本类型：自回归(AR)模型、移动平均(MA)模型、自回归移动平均(ARMA)模型以及差分自回归移动平均模型(ARIMA)。1AR(p)模型：P阶自回归模型，满足下面的方程： (1

10、)其中：c参数为常数；是自回归系数；p为自回归模型阶数；是均值为0、方差为的白噪声序列。2MA(q)模型：q阶的移动平均模型，满足下面的方程：， (2)其中：参数为常数；参数，是q阶移动平均模型的系数；是均值为0、方差为的白噪声序列。3ARMA(p，q)模型将纯AR(p)与纯MA(q)结合，得到一个一般的自回归移动平均方程ARMA(p，q)：，t=l，2，3，T (3)其中为回归系数，为移动平均系数，是模型的待估参数。当p=0时，ARMA(p，q)=MA(q)；当q=0时，ARMA(p，0)=AR(p)。4ARIMA(p，d，q)模型：通过对不平稳的时间序列进行d阶差分，将其转化为平稳时间序列

11、，然后建立ARMA(p，q)模型。设是d阶单整时间序列，即 I(d)，则： (4)为平稳时间序列，即 I(0)，于是可以对建立ARMA(p，q)模型： (5)公式(5)表明如果一个序列是单整序列，那么该序列可以由其自身的滞后值以及随机扰动项来解释。即如果该序列平稳(它的行为并不会随着时间的推移而变化)，那么就可以通过该序列过去的行为来预测未来。（二）ARIMA(p，d，q)模型建模思路3Box和Jenkins提出了具有广泛影响的建模思想，能够对实际建模起到指导作用，其建模思想可以分为以下五个步骤：1对原时间序列进行平稳性检验，如果序列不满足平稳性条件，可以通过差分变换(如果差分的阶数是d，则进

12、行d阶差分)或者其他变换，如对数差分变换使序列满足平稳性的条件。将非平稳时间序列转化为平稳时间序列，实现短期的均衡，是对非平稳时间序列进行ARMA分析的必要前提。2对原时间序列进行白噪声检验。若序列为白噪声，则分析结束；若序列为非白噪声序列，则可以拟合ARIMA模型。3对模型进行定阶，通过计算能够描述时间序列特征的一些统计量，如自相关(AC)系数和偏自相关(PAC)系数，以确定ARMA模型的阶数p和q。另外运用AIC准则或者SBC准则对模型进行识别，从而达到在初始估计中选择尽可能少的参数，简化建模过程的目的。4对模型的未知参数进行估计，并检验参数的显著性，以及模型本身的合理性。5进行诊断分析，

13、以证实所得模型确实与所观察到的数据特征相符。在ARIMA(p，d，q)建模过程中，需要一些统计量和检验来分析模型形式的选择是否适宜。所需要的统计量和检验包括：检验模型参数显著性水平的t统计量；模型的残差序列应为白噪声序列，可采用检验序列相关的方法或残差单位根检验方法进行检验。五、ARIMA模型的建立（一）进出口序列平稳性检验 本文用x代表中国出口额年度数据序列，用y代表中国进口额年度数据序列，则绘制x与y序列的时序图如图1、图2所示。图1 序列x的时序图图2 序列y的时序图从时序图可以看出x与y序列均呈现指数变化趋势，即两序列为非平稳序列，故需对其取对数再做平稳性检验。本文采用ADF(Augm

14、ented DickeyFuller)方法进行序列单位根检验，对序列lnx与lny的检验结果如下图3、图4所示。图3 序列lnx的ADF检验结果图4 序列lny的ADF检验结果检验结果显示，无论考虑何种类型的模型，检验统计量的P值（第6列）均显著大于（=0.05），所以可以认为两序列均为非平稳序列，需要对序列进行差分运算再做平稳性检验，对差分后的序列 Dlnx与 Dlny的检验结果如下图5、图6所示。图5 序列Dlnx的ADF检验结果图6 序列 Dlny的ADF检验结果检验结果表示所有变量在5的显著水平上均通过ADF检验。（二）进出口序列白噪声检验我们知道如果一个序列是纯随机序列，那它的序列值

15、之间应该没有任何相关关系，也就是说很难根据历史信息预测未来年份的进出口额2，故在建立ARIMA模型之前，必须对 Dlnx与 Dlny进行白噪声检验，结果如下图7、图8所示。图7 序列Dlnx的白噪声检验结果图8 序列 Dlny的白噪声检验结果检验结果显示，在各阶延迟下LB检验统计量的P值(第4列)都非常小，故可判断两序列均为非白噪声。（三）对模型进行最优定阶经过一阶逐期差分后， Dlnx与 Dlny变为平稳时间序列，故d = l。所以选取ARIMA(p，1，q)模型来对两序列进行建模。接下来我们确定是使用AR(p)模型还是MA(q)模型，或者是ARMA(p，q)模型对平稳的时间序列 Dlnx与

16、 Dlny进行估计。本文采用SAS进行数据分析及建模，为了避免因个人经验不足导致的模型识别问题，SAS系统还提供了相对最优模型识别，为了提高分析速度，需要先观察序列自相关图与偏自相关图，以决定在何范围内进行定阶，序列 Dlnx与 Dlny的AC与PAC图如下图9图12所示：图9 序列Dlnx的AC图图10 序列Dlnx的PAC图图11 序列 Dlny的AC图图12 序列 Dlny的PAC图综合各图所示，对序列 Dlnx进行自相关延迟阶数小于等于5，移动平均延迟阶数也小于等于5的最优模型识别，结果如下图13所示；对序列 Dlny进行自相关延迟阶数小于等于7，移动平均延迟阶数也小于等于7的最优模型

17、识别，结果如下图14所示。图13 序列 Dlny的模型识别图图14 序列 Dlny的模型识别图在模型选择的过程中，BIC信息量达到最小的模型的即为最优模型2，从图13最后一条信息显示，在自相关延迟阶数阶数小于等于5，移动平均阶数也小于等于5的所有ARMA（p,q）模型中，BIC信息量相对最小的是ARMA(1,0)模型，即AR(1)模型。故序列 Dlnx的最优拟合模型为AR（1）。同理，序列 Dlny的最优拟合模型则为AR(2)。(四)对模型的未知参数进行估计，并检验参数的显著性，以及模型本身的合理性。确定了拟合模型的阶数之后，可以对模型中的未知参数进行估计，以确定模型的口径，并对拟合好的模型进

18、行显著性诊断，拟合结果如下图15、图16所示。图15 序列 Dlny的模型拟合结果图图16 序列 Dlny的模型拟合结果图图15的结果表明，序列 Dlnx的拟合模型各参数均通过显著性检验，但是序列 Dlny的二阶自相关系数却未能通过显著性检验，可考虑将其模型调整为AR(1)，重新进行拟合，结果如下图17所示。图17 序列 Dlny的模型再拟合结果图从检验结果可知，系数显著性比AR(2)模型要有所改善，故最终决定对序列 Dlny采用AR(1)模型进行拟合。(五)进行诊断分析，以证实所得模型确实与所观察到的数据特征相符。在建模的最后，为了确定模型的拟合效果，还需对残差进行白噪声检验，结果如下图18

19、、图19所示。图18 序列 Dlnx的残差白噪声检验结果图图19 序列 Dlny的残差白噪声检验结果图检验结果显示，在各阶延迟下LB检验统计量的P值均大于0.05，故可判断两序列均为白噪声，基本没有未提取的信息了，模型已经提取了有规律的信息，模型拟合和效果较好，ARIMA(1,1,0)即为两序列最佳预测模型。五、模型的预测分析根据图15、图17所示，序列 Dlnx与 Dlny的ARIMA(1,1,0)模型分别为： （6） （7）由 Dlnx、 Dlny模型可推导序列lnx、lny的预测公式为： （8） （9） 于是，得到原时间序列x、y的预测公式为： (10) (11)根据预测公式(10)和公

20、式(11)对20102014年中国进出额进行事后预测，预测结果见表2。表2 中国出口、进口额预测值（20102014）年份点估计值95%置信区间下限95%置信区间上限201011869.6410790.49252.4117330.59715227.2115883.12201112652.8712115.157490.8645529.27721374.1926518.88201213977.9913740.636251.2674212.50331261.444775.31201315703.5115623.25298.293217.34746548.0175796.87201417780.181

21、7774.364524.8492458.7669873.34128388六、结论从以上分析可知，中国进出口序列是受多种制约因子影响的非平稳时间序列，是一组依赖于时间变化的随机变量，可用ARIMA模型予以近似描述。经过模型参数估计与诊断检验以及实证检验发现，可利用其对中国未来几年进出口进行预测分析。从上表可以看出：继20082009年的短暂下降之后2010一2014年中国进出口额继续增加。2010年的进、出口总额分别为10790.4、11869.64亿美元，在95%置信水平下，进口额和出口额的区间估计分别为（9252.411，15227.21）亿美元、（7330.597，15883.12）亿美元

22、。同理，20112014年的预测值也可从上表查得，预测区间的存在也更有利于准确地把握进出口额的变动。同时，通过对预测结果的观察，2010一2014年中国对外贸易顺差仍然存在。但顺差额呈递减趋势，顺差额分别为1079.24亿美元、537.72亿美元、237.36亿美元、80.31亿美元、5.82亿美元。根据该模型预测出来的进出口额较为准确，可为国家科学制定对外贸易政策和合理规划外贸进出口提供重要的理论依据。参考文献：【1】 N.格里高利.曼昆.宏观经济学M.北京：中国人民大学出版社，2005.【2】 汪艳涛，王记志.中国农产品贸易ARIMA模型的建立及预测：20092012年J.经济与管理，2009,23（7）.【3】 王燕.应用时间序列分析M.北京：中国人民大学出版社，2005.【4】 程兰芳，陆敏.中国服务贸易进出口规模的ARIMA模型和预测J.黑龙江对外贸易，2010（12）.【5】 徐国祥.统计预测和决策M.上海：上海财经大学出版社，2008.11.【6】 沈汉溪，林坚.基于ARIMA模型的中国外贸进出口预测：20062010J.国际贸易问题，2007（6）.【7】 栾惠德，张晓峒.协整还是协变：来自中国进出口时间序列的经验证据（19502004）J.南开经济研究，2007（2）【8】 美达摩达尔.N.古扎拉蒂.计量经济学基础M.北京：中国人民大学出版社，2004.10