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

居民消费价格指数的时间序列分析摘要时间序列分析是一种根据动态数据揭示系统动态结构和规律的统计方法本文以我国2007年1月至2011年4月居民消费价格指数为研究对象，基于居民消费价格指数存在明显的非平稳性和季节性特征，运用自回归移动平均季节模型进行建模分析，并利用SPSS建立了居民消费价格指数时间序列的相关关系模型，并对其进行预测，取得较好的效果关键词居民消费价格指数SPSS软件时间序列分析预测、引言大于显著性水平，因此接受零假设，即认为序列随机表3一阶差分和一阶季节差分后数据自相关与偏自相关函数的数据统计Runs TestSDIFF居民物DIFF（居民物价指数_1,1,价指数,1）12Test Value

1.

20.00CasesTest Value2519Cases=Test Value2620Total Cases5139Number ofRuns2620Z-

139.000Asymp.Sig.2-tailed.

8901.000a.Median

4.模型的建立经过一阶差分和一阶季节差分后数据已经平稳化，下面对平稳后的数据进行平稳时间序列的ARMA（p,q）模型的拟合

（1）模型的识别画自相关系数（图4）和偏自相关系数（图5）图居民物价指数□Coefficient------UpperConfidenceLimit------LowerConfidenceLimit图4居民物价指数自相关系数图居民物价指数OCoefficient

1.0-------UpperConfidenceLimit-----LowerConfidenceLimit-

1.0-i i iiiii iiiiiiiiii357911131517192123252729313335Lag Number图5居民物价指数偏自相关系数A A由图4和图5可以看出序列与原序列皆不截尾，但都被负指数pk函数控制收敛到零，此时时间序列有可能为ARMA序列2模型定阶及模型的参数估计通过SPSS软件中的结果对季节差分改进后的时间序列模型ARIMAp,d,qP,D,Q12进行拟合效果的比较，从而最终确定模型的阶数见表4表4各模型参数估计及检验结果3,22,22,11,21,11,0P,q

0.

3270.

7150.

7760.

8610.

8780.37P1dLL.Nrawpo%

0.163-

0.

5850.08———________g.__________“

310.23—————%-

0.503-

0.494-

0.503-

0.496-

0.515-

0.

54010.

1170.

6240.

560.

6690.62一e

0.052-

0.997—-

0.092—一2r

20.

9590.

9560.

9580.

9580.

9580.954BIC-

0.458-

0.498-

0.633-

0.634-

0.729-

0.745RMSE

0.

6310.

6430.

6250.

6240.

6190.638MAPE

71.

8888.

92276.

28275.

70276.

6480.298根据表4中调整后的样本决定素数，以及BIC准则，考察模型的整体拟合效果，力求简洁、有效表6时间序列模型的参数估计ARIMA ModelParametersEstimate SEt Sig.居民物价居民物No ARLag.

370.

1372.

699.010指数-模价指数Transformation1型_1Difference1AR,Lag

540.122-

4.

439.000Seasonal1模型ARMA1,0的BIC值较小，且系数均通过检验见图6,所以最终确定改进后的ARIMAL h01,0,012模型为时间序列凡的最佳预测模型:1—

0.37Bl+

0.54B12l—BX=at t3模型的诊断和检验对模型进行适应性检验，SPSS输出的模型适应性检验的Ljung-Box结果如下见表7表7时间序列模型的检验Model StatisticsModelFit statisticsLjung-Box Ql8ModeNumber ofINumber ofPredictorRMSNormalize StatisticD SigOutlier ssR-square dEMAPE dBIC sF•

0.

954.

6380.

2974537.228I.000居民物价8862指数-模型」P值表明ARIMA1,1,01,0,012模型是合适的残差自相关如图6所示，残差自相关检验也表明ARIMA1,1,01,0,012模型是适合的其图形输出在下一页:Residual ACFResidual PACF口I□□□□D□□□居民物数□，=模口泮ZD□一0□，0-；0-

0.5o1o I-

1.0-

0.

50.

00.

50.

51.0Residual图6时间序列模型的残差自相关图4模型的预测首先，将数据往期的拟合值与实际值对比见图7,可以看到拟合效果比较好Z222211111111114-3-2-L0-9-8-7-6-5-4-3-2-L0-9-8-7-6-5-4-3-2-L.-----Observed——FitForecast居民物」价指数模型-

2.50--

5.00-———————————————————一一一Date图7居民消费价格指数实际值与拟合值序列图在建模时特将我国2011年5月至2011年10月的居民消费价格指数的实际观测值留出,作为预测精度的参照对象利用建立的ARIMA1,1,01,0,012模型对这6个月的CPI指数进行预测，通过SPSS软件可以直接得到数据的短期预测值，预测出我国2011年5月至2011年10月的居民消费价格指数与实际值基本吻合见表8同时给出2011年11月和12月的CPI指数预测值见表8模型的预测结果消费者物价指数月份实际值观测值

2011055.

55.

202011066.

45.

362011076.

55.

192011086.

25.

122011096.

15.

102011105.

54.71201111—

4.37201112—

4.68DqEnN52O-O-O-050Jan2007-Mar2007-May27Ju-2007Sep2007Nov2007Jan2008Mar28May28Ju28Sep2008Nov28Jan2009-Mar29-May29Ju2009Sep2009lzo29Jan2010-Mar2010-May201Ju-2010Sep21Nov21Jan2011Mar221May2011Ju-21Sep20111从表8可以看出，滞后一期的预测效果较好，之后的2步、3步等等预测得到的预测值效果不是太好而当到了t+1的时候，滞后一期已经成为已知,我们习惯上利用这一最新的信息，对预测值进行修正，那么能否原来的预测得到新的预测呢？下面就预测值的适时修正进行讨论

五、预测值的适时修正对于预测模型1-

0.37Bl+

0.54B12l-BX=a,即t tX=

1.37X-

0.17X_+

0.74X_+

0.2X_+aJIt t12t13t14t首先由Pi、pi213和Pi4可求出格林函数G0=1Gi—PiG=1,370G=PiG]=

1.8772P12GOG12=S1G11+G13=W1G12+P12Gl+P13Go■■对于一个ARMA系统，我们有00XtQ+1=W Ga_j1+1+]tj=0ooXt+iD=W Gjja_j+t+1j=0因而有Xt+iQ=XtQ+1+G]at+1其中at+i=Xt+i-X lot当我们已知观测值Xo5=Yo5-丫=

2.03,那么a=X-X l=

2.03-

1.73=

0.3050504X l=X2+Gia5=

1.89+

0.3x

1.37=

2.30504X2=X3+G2a05=

1.72+

0.3x

1.877=

2.28O504将新预测值加上均值可得2011年6月和7月的新预测值为

5.77和

5.75,比之前效果好了很多由此对模型评价如下预测时短期预测有比较好的预测效果，但随着时间的延长，它呈现出较差的预测效果随着时间的推进，可以根据观测的实际值进一步修正模型对10月份以后的预测,得到精确度更高的预测值预测值在不断修正之后变得比较准确应用时间序列分析的方法对居民消费价格指数（2007年-2011年）的变化建立模型，可以较好的模拟和预测价格指数的变化规律，对数据的预报有一定的参考价值

六、基于预测值看我国对经济发展的决策与控制居民消费价格指数表示对普通家庭的支出来说，购买具有代表性的一组商品，在今天要比过去某一时间多花费多少在日常中我们更关心的是通货膨胀率，它被定义为从一个时期到另一个时期价格水平变动的百分比，公式为T=（Pl—P0）/P0,式子中T为1时期的通货膨胀率，P1和P0分别表示1时期（代表报告期）和0时期（代表基期）的价格水平如果用上面介绍的消费价格指数来衡量价格水平，则通货膨胀率就是不同时期的消费价格指数变动的百分比如此可用预测的居民消费价格指数来估计2011年11月、12月的通货膨胀率分别为

4.37%和

4.68%此数据距公认的5%恶性通胀是咫尺之遥O高通胀意味着价格的高涨新世纪以来，中央政府曾先后两次在全国范围内启动价格临时干预措施，对稳定市场价格发挥了重要作用国务院近日下发《关于稳定消费价格总水平保障群众基本生活的通知》，称必要时对重要的生活必需品和生产资料可以实行价格临时干预措施社会的稳定是我国的头等大事，是社会和谐、经济发展的前提，而生活必需品和生产资料价格的稳定对社会的稳定有着重要作用参考文献[11王振龙.应用时间序列分析[M].北京中国统计出版社，2010

[2]张丽，牛惠芳.时间序列分析方法在居民消费价格指数预测中的应用.洛阳师范学院学报，2008年第二期

[3]严健标，李强，基于ARIMA模型的我国农村居民消费水平预测分析【4】中国统计年鉴2010

（一）问题的基本情况及背景居民消费价格指数的调查范围和内容是居民用于日常生活消费品的全部商品和服务项目价格包括食品、烟酒及用品、衣着、家庭设备用品及维修服务、医疗保健和个人用品、交通和通讯、娱乐教育文化用品及服务、居住等八大类商品及服务项目价格既包括居民从商店、工厂、集市所购买商品的价格，也包括从餐饮行业购买商品的价格该指数以实际调查的综合平均单价和根据住户调查有关资料确定的权数，按加权算术平均公式计算全国居民消费价格指数是反映居民家庭购买生活消费品和支出服务项目费用价格变动趋势和程度的相对数其目的在于观察居民生活消费品及服务项目价格的变动对城乡居民生活的影响，为各级党政领导掌握居民消费状况，研究和制定居民消费价格政策、工资政策以及为新国民经济核算体系中有消除价格变动因素的不变价格核算提供科学依据居民消费价格指数还是反映通货膨胀的重要指标当居民消费价格指数上升时，表明通货膨胀率上升，消费者的生活成本提高，货币的购买能力减弱；相反，当居民消费价格指数下降时，表明通货膨胀率下降，亦即消费者的生活成本降低，货币的购买能力增强居民消费价格指数的高低直接影响居民的生活水平，因此，准确的分析并及时的对居民消费价格指数做出合理的预测，对国家制定相应的经济政策，实行宏观调控，稳定物价，保证经济的增长平稳发展具有重要意义

（二）问题的提出时间序列是指同一种现象在不同时间上的相继观察值排列而成的一组数字序列时间序列预测方法的基本思想是预测一个现象的未来变化时，用该现象的过去行为来预测未来即通过时间序列的历史数据揭示现象随时间变化的规律，将这种规律延伸到未来，从而对该现象的未来做出预测对此希望建立相关居民消费价格指数的数学模型并预测居民消费价格指数的走势三问题分析居民消费价格指数是一个滞后性的数据，根据居民消费价格指数的这一特点，我们可以运用时间序列分析的方法对居民消费价格指数进行拟合,从而对未来的居民消费价格指数走势做出合理的预测

二、模型的介绍及说明一时间序列模型的介绍时间序列是按时间顺序取得的一系列数据，时间序列分析方法有很多,本文主要讨论ARMA模型即自回归移动平均模型的方法ARMA模型是一类常用的随机时序模型，由博克斯Box、詹金斯Jenkins创立，简称B—J方法在B—J方法中，只有平稳的时间序列才能直接建立ARMA模型，这就要求时间序列满足假设条件1对任意时间t,其均值恒为常数；2对任意时间t和s,其自相关系数只与时间间隔t-s有关，而与t和s的起始点无关这样时间序列的统计特征不随时间推移而变化，称为平稳时间序列时间序列建模基本步骤是1用观测、调查、统计、抽样等方法取得被观测系统时间序列动态数据2根据动态数据作相关图，进行相关分析，求自相关函数相关图能显示出变化的趋势和周期3辨识合适的随机模型，进行曲线拟合，即用通用随机模型去拟合时间序列的观测数据对于短的或简单的时间序列，可用趋势模型和季节模型加上误差来进行拟合对于平稳时间序列，可用通用ARIMA模型自回归滑动平均模型及其特殊情况的自回归模型、滑动平均模型或组合-ARIMA模型等来进行拟合当观测值多于50个时一般都采用ARIMA模型对于非平稳时间序列则要先将观测到的时间序列进行差分运算，化为平稳时间序列,再用适当模型去拟合这个差分序列通常情况下，自回归移动平均模型的建模过程分为以下几个步骤1对原序列进行平稳性检验，若非平稳序列则通过差分消除趋势；2判断序列是否具有季节性，若有季节波动，则通过季节差分消除季节性；3进行模型识别；4进行模型定阶；5对模型的参数进行估计；6对模型的适合性进行检验，即对残差序列进行白噪声检验P阶自回归序列记作ARp,形如凡=+…+@pXjp,6称为自回归系数,是模型的待估参数q阶移动平均序列记作MAq,形如凡=网-91a.i9a.,为移动平均系数，q tq是模型的待估参数建立平稳时间序列的ARMAp,q模型，其具体形式如下:Xt_^iX-i@pXt_p=a-0a_0q^t-qt t1t1其中6与为模型的待估参数求和自回归移动平均模型autoregressive integratedmoving averagemodel简称ARIMA p,d,q模型，其中AR p为自回归模型，MA q为滑动平均模型，p、q为各自对应阶数，I表示两种模型结合，d为对含有长期趋势、季节变动、循环变动的非平稳时间序列进行差分处理的次数ARIMA模型的通式如下O BVdx,=05^,£/==0,s wtExs=0,Vst、t式中，DK—,中姻-为平稳可逆ARMA p,q模型的自回归系数I B$=1--AB,多项式；田-_纥日,为移动平滑系数多项式，E｝为零均值白噪声序列皿3=1-4ARIMA模型的实质就是差分运算与ARMA模型的组合，任何非平稳序列只要通过适当阶数差分实现差分后平稳，就可以对差分后序列进行ARMA模型拟合二模型的说明时间序列分析主要用于

①系统描述根据对系统进行观测得到的时间序列数据，用曲线拟合方法对系统进行客观的描述

②系统分析当观测值取自两个以上变量时，可用一个时间序列中的变化去说明另一个时间序列中的变化，从而深入了解给定时间序列产生的机理

③预测未来一般用ARMA模型拟合时间序列，预测该时间序列未来值

④决策和控制根据时间序列模型可调整输入变量使系统发展过程保持在目标值上，即预测到过程要偏离目标时便可进行必要的控制拟合好的模型对短期预测有比较好的预测效果，但随着时间的延长，它呈现出较差的预测效果

三、我国居民消费价格指数的时间序列模型拟合

（一）数据的选取及说明本文选取的数据主要来源于国家统计局网站，数据已经进行中心化处理，并在原数据基础上减100以简化计算

（二）时间序列模型

1.数据的录入我国2007年1月至2011年4月居民消费价格指数月度数据表1我国居民消费价格指数月度数据月份消费者物价指数月份消费者物价指数200701-

1.27200903-

4.64200702-

0.77200904-

4.97200703-

0.17200905-

4.84200704-

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.

632010080.

032008072.

832010090.

132008081.

432010100.

932008091.

132010111.

632008100.

532010121.13200811-

1.

072011011.43200812-

2.

272011021.43200901-

2.

472011031.93200902-

5.

072011041.

835.00-1■——电

2.时间序列数据图及平稳性检验图1居民消费价格指数序列图

2.5OH用SPSS软件做出数据序列图（图1）并对序列的平稳性进行游程检验在表2中，概率的P值为

0.000,如果显著性水平为

0.05,由于概率P值小于显著性水平，因此拒绝零假设，即认为序列非随机其检验的SPSS输出结果如下表2时间序列数据是否平稳的游程检验结果Runs Test居民物价指数Test Value1-.02CasesTest Value

26.JAN2007•MAR20072007MAYJUL2007-SEP2007-NOV2007JAi8-MAR2008-MAY20082008JULJEP2008-NOV2008JAi9Dt-MAR2009a-MAY2009JUL289-S2009-NOV2009—JAN200-MAR20-0¥200JUL220-SEP200-NOV20-0JAN2011-MAR20MCases=Test Value26Total Cases52Number ofRuns4Z-

6.443Asymp.Sig.2-tailed.

0003.时间序列的预处理为消除序列的趋势同时减少序列的波动，可以对原有时间序列做二阶逐期差分，并绘制差分后的时序图（见图2）可以看出经过差分处理后的序列趋势基本上消除为了更好地描述月度数据时间序列并进行模拟，需对该序列再进行季节差分，进一步消除季节性（见图3）

4.00-居民物价指数-

2.00-—————袈————Transforms:difference

（2）图2居民消费价格指数二阶差分后时序图-

2.50--

5.00-MARMAYJULSEPNOVJANMARMAYJULSEPNOVJANMARMAYJULSEPNOVJANMAR2008200820082008200820092009200920092009200920102010201020102010201020112011DateaQ-MAR28720072875^2007JUL287SEP2008-NOV2008JAN288-MAR2008-MAY2008JUL2008-SEP2009-NOV-MAR289JAN2009dUL2009SEP2009-Nov2009JAN2010MAR200MAY2010200dUL200200SEPJAN2011IIR2011Transforms:difference2,seasonal difference1,period12图3居民消费价格指数一阶差分和一阶季节差分后时序图在表3中，概率的P值为L000,如果显著性水平为

0.05,由于概率P值。