居民消费价格指数的时间序列分析讲解

1、中南大学数学学院大学生课程设计课程设计名称：时间序列分析专业班级：统计091题目：居民消费价格指数的时间序列分析姓名：杨茹岚学号：1304090109指导教师：唐立时间序列分析课程设计中南大学数学院统计系居民消费价格指数的时间序列分析杨茹岚统计09011304090109摘要：时间序列分析是一种根据动态数据揭示系统动态结构和规律的统计方法。本文以我国2007年1月至2011年4月居民消费价格指数为研究对象，基于居民消费价格指数存在明显的非平稳性和季节性特征，运用自回归移动平均季节模型进行建模分析，并利用SPS&4立了居民消费价格指数时间序列的相关关系模型，并对其进行预测，取得较好的效果

2、。关键词：居民消费价格指数ARIMABox-JenkinsPandit-Wu一、引言（一）问题的基本情况及背景居民消费价格指数的调查范围和内容是居民用于日常生活消费品的全部商品和服务项目价格。包括食品、烟酒及用品、衣着、家庭设备用品及维修服务、医疗保健和个人用品、交通和通讯、娱乐教育文化用品及服务、居住等八大类商品及服务项目价格。既包括居民从商店、工厂、集市所购买商品的价格，也包括从餐饮行业购买商品的价格。该指数以实际调查的综合平均单价和根据住户调查有关资料确定的权数，按加权算术平均公式计算。全国居民消费价格指数是反映居民家庭购买生活消费品和支出服务项目费用价格变动趋势和程度的相对数。其目的在

3、于观察居民生活消费品及服务项目价格的变动对城乡居民生活的影响，为各级党政领导掌握居民消费状况，研究和制定居民消费价格政策、工资政策以及为新国民经济核算体系中有消除价格变动因素的不变价格核算提供科学依据。居民消费价格指数还是反映通货膨胀的重要指标。当居民消费价格指数上升时，表明通货膨胀率上升，消费者的生活成本提高，货币的购买能力减弱；相反，当居民消费价格指数下降时，表明通货膨胀率下降，亦即消费者的生活成本降低，货币的购买能力增强。时间序列分析课程设计中南大学数学院统计系居民消费价格指数的高低直接影响居民的生活水平，因此，准确的分析并及时的对居民消费价格指数做出合理的预测，对国家制定相应的经济政策

4、，实行宏观调控，稳定物价，保证经济的增长平稳发展具有重要意义。(二)问题的提出时间序列是指同一种现象在不同时间上的相继观察值排列而成的一组数字序列。时间序列预测方法的基本思想是：预测一个现象的未来变化时，用该现象的过去行为来预测未来。即通过时间序列的历史数据揭示现象随时间变化的规律，将这种规律延伸到未来，从而对该现象的未来做出预测。对此希望建立相关居民消费价格指数的数学模型并预测居民消费价格指数的走势。(三)问题分析居民消费价格指数是一个滞后性的数据，根据居民消费价格指数的这一特点，我们可以运用时间序列分析的方法对居民消费价格指数进行拟合，从而对未来的居民消费价格指数走势做出合理的预测。二、模

5、型的介绍及说明(一)时间序列模型的介绍时间序列是按时间顺序取得的一系列数据，时间序列分析方法有很多，本文主要讨论ARMA1型即自回归移动平均模型的方法。ARMA1型是一类常用的随机时序模型，由博克斯(Box)、詹金斯(Jenkins)创立，简称B-J方法。在B-J方法中，只有平稳的时间序列才能直接建立ARMA®型，这就要求时间序列满足假设条件：(1)对任意时间t,其均值恒为常数；(2)对任意时间t和s,其自相关系数只与时间间隔t-s有关，而与t和s的起始点无关。这样时间序列的统计特征不随时间推移而变化，称为平稳时间序列。时间序列建模基本步骤是：(1)用观测、调查、统计、抽样等方法取得

6、被观测系统时间序列动态数据。(2)根据动态数据作相关图，进行相关分析，求自相关函数。相关图能显示出变化时间序列分析课程设计中南大学数学院统计系的趋势和周期。(3)辨识合适的随机模型，进行曲线拟合，即用通用随机模型去拟合时间序列的观测数据。对于短的或简单的时间序列，可用趋势模型和季节模型加上误差来进行拟合。对于平稳时间序列，可用通用ARIMA真型(自回归滑动平均模型)及其特殊情况的自回归模型、滑动平均模型或组合-ARIMA模型等来进行拟合。当观测值多于50个时一般都采用ARIMA真型。对于非平稳时间序列则要先将观测到的时间序列进行差分运算，化为平稳时间序列，再用适当模型去拟合这个差分序列。通常情

7、况下，自回归移动平均模型的建模过程分为以下几个步骤：(1)对原序列进行平稳性检验，若非平稳序列则通过差分消除趋势；(2)判断序列是否具有季节性，若有季节波动，则通过季节差分消除季节性；(3)进行模型识别；(4)进行模型定阶；(5)对模型的参数进行估计；(6)对模型的适合性进行检验，即对残差序列进行白噪声检验。P阶自回归序列记作AR(p),形如称为自回归系数，是模型的待估参数。q阶移动平均序列记作MA(q),形如'=2-°何_。qaq,。为移动平均系数，是模型的待估参数建立平稳时间序列的ARMA(p,q射型，其具体形式如下:其中：帖8为模型的待估参数求和自回归移动平均模型(au

8、toregressiveintegratedmovingaveragemodel)简称ARIMA(p,d,q)模型，其中AR(p)为自回归模型，MA(q)为滑动平均模型，p、q为各自对应阶数，I表示两种模型结合，d为对含有长期趋势、季节变动、循环变动的非平稳时间序列进行差分处理的次数。ARIMA真型的通式如下：中南大学数学院统计系时间序列分析课程设计.（B）Wdxt=G（B）备2E1=0,Var寸=1,Ej=0,s=tEXs.=0,Vs<t式中，力=（-B）,（B）=1可B-HlpBP,为平稳可逆ARMAp,q）模型的自回归系数多项式；0（B）=1-郎-"I-%吟为移动平滑系数

9、多项式，力为零均值白噪声序列10。ARIMA真型的实质就是差分运算与ARMAS型的组合，任何非平稳序列只要通过适当阶数差分实现差分后平稳，就可以对差分后序列进行ARMA1型拟合。（二）模型的说明时间序列分析主要用于：系统描述。根据对系统进行观测得到的时间序列数据，用曲线拟合方法对系统进行客观的描述。系统分析。当观测值取自两个以上变量时，可用一个时间序列中的变化去说明另一个时间序列中的变化，从而深入了解给定时间序列产生的机理。预测未来。一般用ARM梗型拟合时间序列，预测该时间序列未来值。决策和控制。根据时间序列模型可调整输入变量使系统发展过程保持在目标值上，即预测到过程要偏离目标时便可进行必要的

10、控制。拟合好的模型对短期预测有比较好的预测效果，但随着时间的延长，它呈现出较差的预测效果。三、我国居民消费价格指数的时间序列模型拟合（一）数据的选取及说明本文选取的数据主要来源于国家统计局网站，数据已经进行中心化处理，并在原数据基础上减100以简化计算。（二）时间序列模型1 .数据的录入我国2007年1月至2011年4月居民消费价格指数月度数据中南大学数学院统计系时间序列分析课程设计表1我国居民消费价格指数月度数据月份消费者物价指数月份消费者物价指数200701-1.27200903-4.64200702-0.77200904-4.97200703-0.17200905-4.84200704-

11、0.47200906-5.14200705-0.07200907-5.292007060.93200908-4.682007072.13200909-4.262007083.03200910-4.002007092.73200911-2.912007103.03200912-1.572007113.43201001-1.972007123.03201002-0.772008013.63201003-1.072008025.23201004-0.672008034.83201005-0.372008045.03201006-0.572008054.23201007-0.172008063.632

12、010080.032008072.832010090.132008081.432010100.932008091.132010111.632008100.532010121.13200811-1.072011011.43200812-2.272011021.43200901-2.472011031.93200902-5.072011041.836时间序列分析课程设计中南大学数学院统计系2 .时间序列数据图及平稳性检验图1居民消费价格指数序列图用eviews软件做出数据序列图（图1）并对序列的平稳性进行游程检验。在表2中，ADF>-3.4946,该序列非平稳表2时间序列数据是否平稳的游程检

13、验结果AugmentedDickeyJuHerUnitRootTestonXt-StatisticProb*AugmentedDickey-FuNerteststatistic-231024001731Testcriticalvalues-3677723-2925169-2.6006581%level5%level1。%levelNullHypothesisXhasaunitrootExogenousConstantLagLength3(AutomaticbasedonSICMAXLAG=10)3 .时间序列的预处理为消除序列的趋势同时减少序列的波动，可以对原有时间序列做二阶逐期差分，并绘制差

14、分后的时序图（见图2）0可以看出经过差分处理后的序列趋势基本上消除。为了更好地描述月度数据时间序列并进行模拟，需对该序列再进行季节差分，进一步消除季节性（见图3）。中南大学数学院统计系时间序列分析课程设计图2居民消费价格指数一阶差分后时序图NullHypothesis:YhasaunitrootExogenous:ConstantLagLength:3(AutomaticbasedonSIC,MAXLAG=9)t-StatisticProb.*AugmentedDickey-Fullerteststatistic-2.7463110.0763Testcriticalvalues:1%level

15、-3.6267845%level-2.94584210%level-2.611531*MacKinnon(1996)one-sidedp-values.AugmentedDickey-FullerTestEquationDependentVariable:D(Y)Method:LeastSquaresDate:06/13/12Time:18:50Sample(adjusted):2008M052011M04Includedobservations:36afteradjustmentsVariableCoefficientStd.Errort-StatisticProb.Y(-1)-0.1116

16、350.040649-2.7463110.0099D(Y(-1)-0.0416870.152752-0.2729080.7867D(Y(-2)0.4322300.1386793.1167620.0039D(Y(-3)0.4447120.1589452.7978940.0088C-0.0977790.203069-0.4815070.6335中南大学数学院统计系时间序列分析课程设计R-squared0.442918Meandependentvar-0.083333AdjustedR-squared0.371036S.D.dependentvar1.520432S.E.ofregression1.

17、205813Akaikeinfocriterion3.340431Sumsquaredresid45.07354Schwarzcriterion3.560365Loglikelihood-55.12777F-statistic6.161764Durbin-Watsonstat2.126619Prob(F-statistic)0.000907表3居民消费价格指数一阶差分检验在表3中，ADF=-3.626784<-3.4946因而变换后是平稳序列5.模型的建立1) Pandit-wuDependentVariable:YMethod:LeastSquaresDate:06/13/12Time

18、:19:46Sample(adjusted):2008M032011M03Includedobservations:37afteradjustmentsConvergenceachievedafter86iterationsBackcast:OFF(RootsofMAprocesstoolarge)VariableCoefficientStd.Errort-StatisticProb.AR(1)1.9102010.04124246.316690.0000AR(2)-0.9318600.041397-22.510450.0000MA(1)-1.2908890.152106-8.4867940.0

19、000R-squared0.967035Meandependentvar-0.597297AdjustedR-squared0.965096S.D.dependentvar5.255345S.E.ofregression0.981840Akaikeinfocriterion2.878829Sumsquaredresid32.77635Schwarzcriterion3.009444Loglikelihood-50.25833Durbin-Watsonstat2.658866InvertedARRoots.96+.14i.96-.14iInvertedMARoots1.29EstimatedMA

20、processisnoninvertibleSSE=32.77AIC=2.879SC=3.009再拟合ARMA(4,3澈型：中南大学数学院统计系时间序列分析课程设计DependentVariable:YMethod:LeastSquaresDate:06/13/12Time:19:51Sample(adjusted):2008M052011M03Includedobservations:35afteradjustmentsConvergenceachievedafter39iterationsBackcast:OFF(RootsofMAprocesstoolarge)VariableCoeff

21、icientStd.Errort-StatisticProb.AR(1)2.0580000.18769710.964490.0000AR(2)-1.4384110.380206-3.7832440.0007AR(3)0.5976110.5062351.1805000.2477AR(4)-0.2519980.287915-0.8752540.3889MA(1)-1.6958600.294125-5.7657770.0000MA(2)1.1535670.3645093.1647170.0037MA(3)-0.8683220.338996-2.5614520.0161R-squared0.97503

22、9Meandependentvar-0.911429AdjustedR-squared0.969691S.D.dependentvar5.229039S.E.ofregression0.910356Akaikeinfocriterion2.826894Sumsquaredresid23.20494Schwarzcriterion3.137964Loglikelihood-42.47065Durbin-Watsonstat2.402784InvertedARRootsInvertedMARoots.98+.17i.98-.17i.05-.50i1.32.19-.79i.19+.79iEstima

23、tedMAprocessisnoninvertible.05+.50iSSE=23.205AIC=2.827SC=3.138F检验:=5.049F-.(4,49)32.77-23.205/23.205452-3DependentVariable:YMethod:LeastSquaresDate:06/19/12Time:10:40Sample(adjusted):2008M072011M04Includedobservations:34afteradjustmentsConvergenceachievedafter36iterationsBackcast:2006M082006M12Varia

24、bleCoefficientStd.Errort-StatisticProb.AR(1)0.6200670.2251752.7537150.0113AR(2)0.6682970.2532062.6393400.014710中南大学数学院统计系时间序列分析课程设计AR(3)0.1159680.2898580.4000860.6928AR(4)-0.3089890.289929-1.0657410.2976AR(5)-0.0692330.240974-0.2873040.7765AR(6)-0.1328890.208902-0.6361290.5310MA(1)0.1176960.0774961.

25、5187240.1425MA(2)-0.1177470.073804-1.5953960.1243MA(3)0.0885790.0694141.2761020.2147MA(4)-0.0644420.076148-0.8462710.4061MA(5)-0.9399000.058310-16.118920.0000R-squared0.974955Meandependentvar-1.017647AdjustedR-squared0.964067S.D.dependentvar5.269210S.E.ofregression0.998838Akaikeinfocriterion3.091744

26、Sumsquaredresid22.94657Schwarzcriterion3.585566Loglikelihood-41.55964Durbin-Watsonstat2.074405InvertedARRoots.96-.17i.96+.17i-.04+.50i-.04-.50i-.62-.41i-.62+.41iInvertedMARoots.98.29+.94i.29-.94i-.84-.53i-.84+.53iSSE=22.947AIC;3.092SC=：3.58622.947-23.20522.947F检验：F=/=-0.126<F仃(4,45)452-7'无显著故

27、ARMA(2,1)WARMA(4,3濮型有显著性差异，ARMA6,5)与ARMA4,差异，则选择ARMA(4,3澈型,即yt2.06yt一1.44yt工0.60yt=at1.70at1.15at20.87at上2) Box-Jenkins识别模型Series:cpiAutocorrelationsLagAutocorrelation一.一aStd.ErrorBox-LjungStatisticValuedfSig.b1.253.1363.4541.0632.374.13511.1812.0043.385.13319.5303.0004.206.13221.9804.00011中南大学数学院统计

28、系时间序列分析课程设计5.195.13024.2135.0006.179.12926.1376.0007.088.12826.6147.0008.130.12627.6718.0019-.114.12528.5099.00110-.332.12335.79410.00011.097.12236.43411.00012-.481.12052.46712.00013-.183.11954.84613.00014-.160.11756.72114.00015-.276.11562.44415.00016-.257.11467.54516.00017-.246.11272.35117.00018-.2

29、61.11077.94118.00019-.127.10979.30019.00020-.208.10783.08020.00021-.173.10585.76521.00022.116.10487.01222.00023-.230.10292.12023.00024-.019.10092.15724.00025-.005.09892.16025.00026-.030.09692.25826.00027.005.09492.26127.00028.033.09292.39128.00029.036.09092.54829.00030.095.08893.71630.00031.008.0869

30、3.72631.00032.044.08494.00732.00033.123.08296.26733.00034-.039.07996.50534.00035.036.07796.71935.00036.032.07496.90336.00037.022.07296.99637.00038.034.06997.22938.00039.041.06797.61439.00040.036.06497.94140.00012中南大学数学院统计系时间序列分析课程设计41.023.06198.08341.00042-.015.05898.15542.00043.019.05498.27743.0004

31、4.042.05198.94544.00045.003.04798.95045.00046-.003.04398.95546.00047-.005.03898.97047.00048.009.03399.04048.00049.004.02799.05749.000a. Theunderlyingprocessassumedisindependence(whitenoise).b. Basedontheasymptoticchi-squareapproximation.PartialAutocorrelationsSeries:cpiLagPartialAutocorrelationStd.E

32、rror1.253.1402.332.1403.286.1404.000.1405-.048.1406-.001.1407-.037.1408.042.1409-.247.14010-.508.14011.308.14012-.365.14013.114.14014.057.14015-.004.14016-.031.14017-.042.14018.031.14019-.083.14013中南大学数学院统计系时间序列分析课程设计20-.039.14021.035.14022-.157.14023.121.14024-.243.14025.012.14026.017.14027-.100.14

33、028.002.14029-.063.14030.037.14031.044.14032-.020.14033-.073.14034-.077.14035-.085.14036-.106.14037-.001.14038.076.14039-.017.14040.018.14041-.041.14042.001.14043-.027.14044-.002.14045-.117.14046-.082.14047-.106.14048.042.14049.081.1403）模型的识别画自相关系数（图4）和偏自相关系数（图5）图14中南大学数学院统计系时间序列分析课程设计cpiLagNumber图4

34、居民物价指数自相关系数图15中南大学数学院统计系时间序列分析课程设计5G0,0.LO4豆tEdcpiCoefficientUpperConfidenceLimitLowerConfidenceLmirttrir1357iiiiiiiiiiiir91113151719212325272931333537394143454749LagNumber图5居民物价指数偏自相关系数AA由图4和图5可以看出Pk序列与理卜序列皆是截尾，都被负指数函数控制收敛到零,此时时间序列有可能为ARM前列。16时间序列分析课程设计中南大学数学院统计系CorrelogramofYAutocorrelationPartial

35、CorrelationACPACQ-StatProb11I108380038210660000111120,670-0.11136.3610.00011I1130.505-009644.9460.000i11140366-0023496450000Z111150.239-006651.7330.000"1iI160136-0022524410.0001E170028-011352.4730000匚11匚1a-D.110-0.21152.9810.000匚i11(3-0.217-0026550630.000匚1I110-0.28300Q3587910000匚11111-0.312-00

36、0663.6060.000匚111112-0,2800.12367.7350000匚1I口13-0.146032968.920000011140.00001396B.9200.0001Li150000-045568.9200000|1111160000-0042609200000111170.00000386B.9200.00011-180000-0080689200000|11C1190000-0.10468920000011111200.000-0.08468.9200.0001132100000.19368.9200.000|11n220000031060.9200000111i230.

37、000002568.9200.000i1240000000368.9200000图4相关系数(A。和偏自相关系数(PAC人C则当K>2时，则Qk<L-0.33,即kk呈现1步截尾现象，而J36序列被负指数函数控制收敛于零，呈拖尾现象，故可初步判定序歹1Y适合AR(1)模型。yt=1.37ytat4)模型定阶及模型的参数估计通过SPS球件中的结果对季节差分改进后的时间序列模型ARIMA(pd,q)(P,D,Q)12进行拟合效果的比较，从而最终确定模型的阶数(见表4)。表4各模型参数估计及检验结果(p,q)(3,2)(2,2)(2,1)(1,2)(1,1)(1,0)17叼0.3270.

38、7150.7760.8610.8780.37%0.163-0.5850.08-0.23-*1-0.503-0.494-0.503-0.496-0.515-0.54610.1170.6240.560.6690.62-820.052-0.997-0.092-20.9590.9560.9580.9580.9580.954BIC-0.458-0.498-0.633-0.634-0.729-0.745RMSE0.6310.6430.6250.6240.6190.638MAPE71.8888.92276.28275.70276.6480.298中南大学数学院统计系时间序列分析课程设计根据表4中调整后的样本

39、决定系数r2,以及BIC准则，考察模型的整体拟合效果,力求简洁、有效表6时间序列模型的参数估计EstimateSEtSig.370,1372.6990101-.540.122-4.439.000ARIMAModelParameters居民物价居民物NoARLag指数-模价指数Transformation1型DifferenceAR,LagSeasonal16),所以最终确定改进后模型ARMA(10)的BIC值较小,且系数均通过检验(见图的ARIMA(1,1,0)(1,0,0)12模型为时间序列Xt的最佳预测模型：(1-0.37B)(10.54B12)(1-B)Xt-atLjung-Box结果如

40、下(见表5)模型的诊断和检验对模型进行适应性检验，SPSS俞出的模型适应性检验的18时间序列分析课程设计中南大学数学院统计系7)：表7时间序列模型的检验ModelStatisticsModelNumberofPredictorsModelFitstatisticsLjung-BoxQ(18)NumberofOutliersR-squaredRMSEMAPENormalizedBICStatisticsDFSig.居民物价指数-模型_10.954.63880.298-.74537.22816.0020P值表明ARIMA(1,1,0)(1,0,0)12模型是合适的。残差自相关如图6所示，残差自相关

41、检验也表明ARIMA(1,1,0)(1,0,0)12模型是适合的其图形输出在下一页：ResidualPACFn口JJODnunIDnolmDaID=口-1.045D.OO.B1.&-1.0-OQ.D0.51.DResidual图6时间序列模型的残差自相关图7),可以看到拟合效果比较好6)模型的预测首先，将数据往期的拟合值与实际值对比(见图19时间序列分析课程设计中南大学数学院统计系居民模型115cp冬一JuwoliIf20T1If81Jan片11lz¥xa5宅loaJulasoIf202If2。占It20Ims詈XQWIJ2009I5200ggn289If8Slsp苍20。，

42、£zoom1手20081?200BIfsQrIm=2QCT1LJM007IF2007IE287Dat图7居民消费价格指数实际值与拟合值序列图在建模时特将我国2011年5月至2011年10月的居民消费价格指数的实际观测值留出，作为预测精度的参照对象。利用建立的ARIMA(1,1,0)(1,0,0)12模型对这6个月的CPI指数进行预测，通过SPSSa件可以直接得到数据的短期预测值，预测出我国2011年5月至2011年10月的居民消费价格指数与实际值基本吻合(见表8)。同时给出2011年11月和12月的CPI指数预测值(见表8)。表8模型的预测结果月份消费者物价指数实际值观测值2011055.55.202011066.45.362011076.55.192011086.25.122011096.15.102011105.54.71201111-4.37201112-4.68从表8可以看出