案例三时间序列分析

.PAGE.案例三：时间序列分析学习目的通过本案例的学习，旨在使同学们到达以下几个方面的学习目标：1、培养学习利用多种时间序列分析方法解决实际问题的能力2.掌握时间序列平滑方法：移动平均、加权移动平均等方法。了解这些平滑方法在处理时时序列数据时各自的优点和缺点，学会用这些方法来处理不同类型和特点的数据。3、掌握时问序列的构成分析方法。影响时间序列的因素大体上可以分为四种，即长期趋势(T)、季节变动(s)、循环波动(c)和不规那么波动(J)，通过本案例的学习，学会如何将各种影响因素分别从时问序列中别离出来并用数量加以测定。4、掌握利用模型对时间序列进展分析的方法。在假定现象未来的开展趋势能够与过去保持一致的前提下，同学们要学会利用以上建立的模型对未来进展预测。数据文件本案例的数据为我国1978至2006年居民收入年度数据资料。数据文件存为ExceI格式，使用时可直接点击cash3对于本案例，可以用各种软件包括Excel、Eview、SPSS或SAS等进展分析。但我们建议使用Eview，因为这个软件在处理时间序列问题时更加方便也更为专业化。案例分析所需统计知识洁明"统计学原理"〔第四版〕复旦大学第130-168页案例分析过程提示首先用软件做出我国1978-2006居民收入时序图。观察数列按时间顺序变化的特点。由图中可以看出，收入有明显的向上开展趋势。在此根底上，我们可用进展以下方法进展分析。1、对城乡居民的收入差异进展描述性分析.2、时间序列平滑法。利用移动平均法对时间数列进展平滑，观察长期趋势3、建立长期趋势模型，进展预测需要讨论的问题1、以上各种分析方法分别适用于什么特点的时间序列数据？2、各种分析方法在分析过程中提供的信息有什么不同？3、除了以上这几种时间序列分析方法之外，你还可以找到其他种类的时间序列分析方法吗？你认为还有哪种分析方法适用于本案例的数据特点？4、利用上述案例资料以小组为单位写一分析报告，题目自拟。附录：阅读案例全国城镇居民收入差异的数量分析本案例来源：朱建平..全国城镇居民收入差异的数量分析随着改革开放的不断深入，社会主义市场经济体制的进一步确立，我国城镇居民的收入普遍提高，人们生活水平明显改善。但是，在开展中另一种趋势也明显起来，即我国城镇居民收入两极分化的程度在加剧，本文试图借助计量经济学的有关理论，从全国城镇居民人均收入的差异着手讨论，通过模型从静态和动态的角度探讨城镇居民收入两极分化的形成过程、现状及开展趋势，进而讨论如何合理确定并及时调整我国城镇居民收入贫困线的标准。一、城镇居民人均收入差异分析几年来我国城镇居民生活有了可喜的变化，居民生活水平有了明显提高。1990年--1996年我国城镇居民收入状况表时间〔年〕人均年收入〔元〕最高收入〔元〕困难户收入〔元〕人均年收入与困难户收入比例19901522.792675.64782.932.05819911713.102956.81928.862003663.001032.002.03219932583.164905.771239.351.99219943502.316837.811566.331.80819954288.098231.311984.921.86219964844.789250.442242.921.858〔如表〕1996年我国城镇居民家庭人均年收入到达4844.78元，比1990年增长218%，年均递增速度为21.28%，各年增长速度均超过了同期各项反项指标，如物价指数、通货膨胀指数的增长速度。可见，全国城镇居民不仅在名义货币收入上有了大幅度的提高，在实际水平上也同样实现了增长。全国困难户人均收入水平由1990年的782.93元提高到1996年的2242.92元，增长了1459.99元；最高收入户人均收入水平由1990年的2675.64元提高为9250.44元，增长2.46倍。这一切都充分说明，社会主义市场经济的逐步确立，使我国城镇居民的收入水平有了明显的提高。但是，应该看到，在全国城镇居民收入水平整体上得到提高的同时，收入的差距被拉大了。七年中，我国困难户与最高收入户居民人均年收入差异从1990年的1892.71元扩大到1996年7007.52元，扩大了2.70倍。均增长速度为24.38%，超过人均收入水平的增长速度。这一结果清楚地说明：七年来全国城镇居民平均收入水平两极分化的程度加剧了。这并不是我们建立社会主义市场经济体制，全面振兴经济的初衷，我们不希望在国家经济明显趋好的大环境下出现更多的"穷人〞，但这又是一个我们不得不承受的现实。进一步的分析我们可以看到，这种差距的拉大还伴随着收入中非工资性收入所占比重增大、灰色收入和资本收入增加的趋势。二、我国城镇居民收入水平及差异的数量分析在上面讨论的基上，根据1989年─1997年"中国统计年鉴"的有关资料，对困难户与最高收入户居民人均年收入的差异及全国城镇居民人均年收入与困难户人均年收入比例分别进展了时间序列分析，建立模型为其中─最高收入户与困难户人均年收入之差；─时间。模型〔1〕均通过了总体与个体的检验显著性检验。其中─全国城镇居民人均年收入与困难户人均年收入比例；─时间。模型〔2〕虽然判定系数=0.725，但是个体检验相当显著，并且标准差和残差平方和都很小，这说明该比例值受时间变化的影响不大。对模型〔1〕求二阶导数，即令模型〔3〕等于零，便得到模型〔1〕所描述的曲线在=4.51处有拐点，如下图。上面的模型及图形清楚地说明：1、全国城镇居民最高收入户与困难户人均年收入差距越拉越大，但是近年来，困难户人均年收入根本上是每年全国平均水平的一半。由于全国城镇居民人均年收入逐年提高，因此划分困难户的标准随之变化。由此可见，人均年收入低于全国城镇居民人均年收入一半的居民户为困难户。2、全国城镇居民最高收入户与困难户人均年收入差距在这七年之的变化可以分为三个阶段：第一阶段是1990年─1992年，从1989年治理整顿后到1992年，收入差距拉大的速度不快；第二阶段是1992年─1994年，在1992年南巡讲话之后，全国经济出现高速开展，收入差距拉大的速度增加。由于模型〔1〕描述的曲线在=4.51处有拐点，那么说明在1993年中间速度最快，但从此之后，收入差距拉大的速度将趋于缓和；第三阶段是1994年─1996年，随着整个经济开展出现软着陆，全国城镇居民最高收入户与困难户人均年收入差距拉大的速度出现了缓和。三、结论1、随着我国城镇居民收入水平的继续提高，最高收入户与困难户人均年收入的差距进一步扩大的趋势将持续下去，这符合收入增长的"马太效应〞理论，是一种正常的变动趋势。2、城镇居民最高收入户与困难户年均收入差距扩大的速度将趋于缓和。这是全社会收入水平普遍提高，收入将逐步趋于规化，社会再分配功能日益发挥作用的必须结果。随着城镇居民收入水平的不断提高，社会再分配手段的作用将日益增大，特别是对高收入阶层来说，政府将通过征收所得税的手段对其高额收入加以适当调节，使其与低收入水平的差距不致过大。另外，随着我国社会保障制度的日益完善，也能在对不同收入水平进展适度调节的前提下缩小收入差距，并提高全社会成员的生活质量。附录：案例阅读基于SARIMA模型的我国入境旅游人数时间序列分析王丽英，X后平，"基于SARIMA模型的我国入境旅游人数时间序列分析"，统计与咨询，2008〔4〕。[摘要]时间序列是一种按照时间顺序取得的一组数据，分析时间序列的常用方法为Box-Jenkins模型。Box-Jenkins模型不以经济理论为指导，依据时间序列自身构造特点建立模型，并利用外推进展预测。本文搜集了2001年1月至2007年9月的入境旅游人数，在此根底上根据Box-Jenkins建模的方法，建立了入境旅游人数带的SARIMA模型，对模型进展了适应性检验，比拟了预测值与观测值的差异，证明模型是较合理的。[关键词]入境人数时间序列SARIMA模型自相关函数偏自相关函数一、关于本文时间序列模型的说明时间序列是一种是按照时间顺序取得的一组数据，大多数的时间序列存在惯性，通过对这种惯性的分析就可以由现在值和过去值对未来值进展预测。时间序列分析是一种根据动态数据提醒系统动态构造和规律的统计方法，其根本思想是根据随机的时间序列建立能够比拟准确的反映时间序列中所包含的动态依存关系的数学模型，并借以对未来进展预测。分析时间序列的方法很多，本文主要讨论Box-Jenkins模型。Box-Jenkins模型不以经济理论为指导，依据时间序列自身构造特点建立模型，并利用外推进展预测。建立时间序列模型的前提条件时如果时间序列是平稳的，就可以用ARMA模型来刻划它。但通常经济时间序列都存在一定的趋势，是不平稳的时间序列，不能直接建立ARMA(p,q)模型，这时差分运算就是一种较好的处理方式，许多非平稳的时间序列差分后会显示出平稳序列的性质，我们称这个非平稳序列为差分平稳序列。对差分平稳序列可以使用ARIMA模型进展拟合。ARIMA〔p,d,q〕模型称为求和自回归移动平均模型。其根本构造为式中：B为滞后算子为平稳可逆ARMA〔p,q〕模型的自回归系数。为平稳可逆的ARMA〔p,q〕模型和移动平滑系数多项式分别表示自回归阶数、差分阶数、移动平均阶数当时，ARIMA〔〕模型就是ARMA〔p,q〕模型。当时,ARIMA〔〕模型可以简记为IMA(d,q)模型当时，ARIMA〔〕模型可以简记为ARI(p,d)模型ARIMA模型可以对具有季节效应的序列建模。乘积季节模型是随机模型与ARIMA模型的结合，其形式为：式中:D为周期步长，d为提取趋势信息所用的差分阶数为白噪声序列。该模型简记为ARIMA〔p,d,q〕×(P,D,Q)二、我国入境旅游人数SARIMA模型的建立自从改革开放以来，我国的旅游事业蓬勃开展，入境旅游人数逐年递增。本文选取了2001年1月至2007年9月入境旅游人数共81个数据，我们利用2001年1月至2007年6月数据进展建模，为检验模型的效果，将2007年7月至9月的3个观测值留意出，作为评价预测精度的参照对象。数据的分析与处理均采用Eviews3．1软件。表1200年1月至2007年9月入境旅游人数〔单位：万人〕时间1月2月3月4月5月6月7月8月9月10月11月12月2001717.38611.6753.26779.7707.87711.8745.4809.17417587508152002740.7711839.54807.77787.92795.8849.4890.78658708428662003848.43737.6785.13564.92543.83652.6776.9884.48088558288772004808.73753.3855.16954.71877.47893.4959971.69179889359902005938.06855.11027.81024.8995.14989.410761067989105599110182006998.85871.410031097.11002.4100010901115104411381042109320071023933.11089.81151.71072.61067114911571123数据来源：../wyzl/more.asp"newsid=004001。2003年3-5月数据是作者根据资料推算。〔一〕数据的平稳性及正态性检验Box-Jenkins时序建模是基于平稳时间序列，因此首先检验数据的平稳性。1、绘制观察值序列时序图图1：入境旅游人数序列时序图时序图显示该序列随时间的推移具有明显的递增趋势，又含有周期为12个月的季节波动。2、图2为根据中国入境旅游人数所作的自相关及偏自相关分析图。从自相关图中可以发现，自相关系数衰减很慢，没有很快衰减到零，因此，该序列含有一定的趋势性。图2：入境旅游人数的自相关分析图3、正态性检验图3：Q-Q图从Q-Q图中我们可知该序列具有正态性〔二〕数据的预处理和模型的识别为了消除异方差，对原数列作对数处理得数列lnrjrs，为消除数列的趋性，对lnrjrs做一阶差分，得到序列dlnrjrs，其时序图和自相关和偏见自相关图如图4-5。图4：序列dlnrjrs的时序图图5：dlnrjrs的自相关和偏自相关图从图4中可以看出，作对数差分后的序列dlnrjrs，其均值在零点附近，原序列的线性递增趋势已被根本消除，该序列是平稳的。从图5中发现，当滞后期K=12时，该序列的自相关系数和偏见自相关系数与零有显著差异，这说明序列具有周期为12个月的季节波动。对序列进展二阶季节差分后发现季节性并没有得到改善，故只做一阶季节差分。经过对数一阶差分，序列的递增趋势根本消除，故d=1,自相关系数和偏见自相关系数均显示出不截尾的性质，同时存在明显的季节效应，可考虑建立乘积季节效应模型〔p,d,q〕×(P,D,Q)，由于实际建模时常用高阶的AR模型代替相应的MA和ARMA。综合考虑，可供选取择的(p,q)组合有：〔1，1〕；〔2，2〕〔2，0〕，〔3，0〕，由于k=12，样本自相关系数和偏自相关系数显著不为0，故P=Q=1。由于是对对数一阶差分建模，我们选择不带截距项的模型〔三〕模型的选择及参数估计将四个模型的相关检验汇入表2。表2各模型的检验结果模型AICSC参数显著性检验〔0.57000.5485-2.6720-2.5371没有通过显著性检验0.57040.5327-2.6417-2.4376没有通过显著性检验0.14140.1046-2.0603-1.9357，没有通过显著性检验0.56840.5385-2.6686-2.4985没有通过显著性检验0.56840.5385-2.6686-2.4985没有通过显著性检验从表2中可见，几个模型的AIC和SC均差异不大，其调整后的决定系数除模型三外其余均在0.56附近。综合考虑可选择模型二进展优化，由于不显著，剔除AR〔2〕重新进展估计。其估计结果均不显著，经反复方试验，我们认为虽然模型二的没有通过检验，但其拟合效果是最优的，因此选择作为估计模型，其估计参数如下:表3：参数估计表VariableCoefficientStd.Errort-StatisticProb.AR(2)0.0715980.1320910.5420370.5899AR(1)0.5864840.1434294.0890200.0001SAR(12)0.8990820.05676315.839220.0000MA(2)-0.4093990.082450-4.9654170.0000MA(1)-0.5763540.082410-6.9937360.0000SMA(12)-0.7239970.127404-5.6827060.0000R-squared0.570418Meandependentvar0.003801AdjustedR-squared0.532736S.D.dependentvar0.090305S.E.ofregression0.061730Akaikeinfocriterion-2.641709Sumsquaredresid0.217202Schwarzcriterion-2.437601Loglikelihood89.21383Durbin-Watsonstat2.020444模型的调整可决系数为0.5704，AIC为-2.6417，SC为-2.4376，除外，各系数均通过显著性检验。因此可以认为不带截距项的ARIMA模