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

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

《居民消费价格指数的时间序列分析讲解.docx》由会员分享，可在线阅读，更多相关《居民消费价格指数的时间序列分析讲解.docx（27页珍藏版）》请在冰豆网上搜索。

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

中南大学数学学院

大学生课程设计

课程设计名称：

时间序列分析

专业班级：

统计09—1

题目：

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

姓名：

***

学号：

**********

*******

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

杨茹岚统计0901********

摘要：

时间序列分析是一种根据动态数据揭示系统动态结构和规律的统计方法。

本文以我国2007年1月至2011年4月居民消费价格指数为研究对象，基于居民消费价格指数存在明显的非平稳性和季节性特征，运用自回归移动平均季节模型进行建模分析，并利用SPSS建立了居民消费价格指数时间序列的相关关系模型，并对其进行预测，取得较好的效果。

关键词：

居民消费价格指数ARIMABox-JenkinsPandit-Wu

一、引言

（一）问题的基本情况及背景

居民消费价格指数的调查范围和内容是居民用于日常生活消费品的全部商品和服务项目价格。

包括食品、烟酒及用品、衣着、家庭设备用品及维修服务、医疗保健和个人用品、交通和通讯、娱乐教育文化用品及服务、居住等八大类商品及服务项目价格。

既包括居民从商店、工厂、集市所购买商品的价格，也包括从餐饮行业购买商品的价格。

该指数以实际调查的综合平均单价和根据住户调查有关资料确定的权数，按加权算术平均公式计算。

全国居民消费价格指数是反映居民家庭购买生活消费品和支出服务项目费用价格变动趋势和程度的相对数。

其目的在于观察居民生活消费品及服务项目价格的变动对城乡居民生活的影响，为各级党政领导掌握居民消费状况，研究和制定居民消费价格政策、工资政策以及为新国民经济核算体系中有消除价格变动因素的不变价格核算提供科学依据。

居民消费价格指数还是反映通货膨胀的重要指标。

当居民消费价格指数上升时，表明通货膨胀率上升，消费者的生活成本提高，货币的购买能力减弱；相反，当居民消费价格指数下降时，表明通货膨胀率下降，亦即消费者的生活成本降低，货币的购买能力增强。

居民消费价格指数的高低直接影响居民的生活水平，因此，准确的分析并及时的对居民消费价格指数做出合理的预测，对国家制定相应的经济政策，实行宏观调控，稳定物价，保证经济的增长平稳发展具有重要意义。

（二）问题的提出

时间序列是指同一种现象在不同时间上的相继观察值排列而成的一组数字序列。

时间序列预测方法的基本思想是：

预测一个现象的未来变化时，用该现象的过去行为来预测未来。

即通过时间序列的历史数据揭示现象随时间变化的规律，将这种规律延伸到未来，从而对该现象的未来做出预测。

对此希望建立相关居民消费价格指数的数学模型并预测居民消费价格指数的走势。

（三）问题分析

居民消费价格指数是一个滞后性的数据，根据居民消费价格指数的这一特点，我们可以运用时间序列分析的方法对居民消费价格指数进行拟合，从而对未来的居民消费价格指数走势做出合理的预测。

二、模型的介绍及说明

（一）时间序列模型的介绍

时间序列是按时间顺序取得的一系列数据，时间序列分析方法有很多，本文主要讨论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阶自回归序列记作AR（p），形如

称为自回归系数,是模型的待估参数。

q阶移动平均序列记作MA（q），形如

，

为移动平均系数，是模型的待估参数。

建立平稳时间序列的ARMA（p,q）模型，其具体形式如下：

其中：

为模型的待估参数。

求和自回归移动平均模型（autoregressiveintegratedmovingaveragemodel）简称ARIMA（p,d,q）模型，其中AR（p）为自回归模型，MA（q）为滑动平均模型，p、q为各自对应阶数，I表示两种模型结合，d为对含有长期趋势、季节变动、循环变动的非平稳时间序列进行差分处理的次数。

ARIMA模型的通式如下：

式中，

，

，为平稳可逆ARMA（p,q）模型的自回归系数多项式；

，为移动平滑系数多项式，{t}为零均值白噪声序列[10]。

ARIMA模型的实质就是差分运算与ARMA模型的组合，任何非平稳序列只要通过适当阶数差分实现差分后平稳，就可以对差分后序列进行ARMA模型拟合。

（二）模型的说明

时间序列分析主要用于：

①系统描述。

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

②系统分析。

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

③预测未来。

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

④决策和控制。

根据时间序列模型可调整输入变量使系统发展过程保持在目标值上，即预测到过程要偏离目标时便可进行必要的控制。

拟合好的模型对短期预测有比较好的预测效果，但随着时间的延长，它呈现出较差的预测效果。

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

（一）数据的选取及说明

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

（二）时间序列模型

1.数据的录入

我国2007年1月至2011年4月居民消费价格指数月度数据

表1我国居民消费价格指数月度数据

月份

消费者物价指数

月份

消费者物价指数

200701

-1.27

200903

-4.64

200702

-0.77

200904

-4.97

200703

-0.17

200905

-4.84

200704

-0.47

200906

-5.14

200705

-0.07

200907

-5.29

200706

0.93

200908

-4.68

200707

2.13

200909

-4.26

200708

3.03

200910

-4.00

200709

2.73

200911

-2.91

200710

3.03

200912

-1.57

200711

3.43

201001

-1.97

200712

3.03

201002

-0.77

200801

3.63

201003

-1.07

200802

5.23

201004

-0.67

200803

4.83

201005

-0.37

200804

5.03

201006

-0.57

200805

4.23

201007

-0.17

200806

3.63

201008

0.03

200807

2.83

201009

0.13

200808

1.43

201010

0.93

200809

1.13

201011

1.63

200810

0.53

201012

1.13

200811

-1.07

201101

1.43

200812

-2.27

201102

1.43

200901

-2.47

201103

1.93

200902

-5.07

201104

1.83

2.时间序列数据图及平稳性检验

图1居民消费价格指数序列图

用eviews软件做出数据序列图（图1）并对序列的平稳性进行游程检验。

在表2中，

，该序列非平稳。

表2时间序列数据是否平稳的游程检验结果

3.时间序列的预处理

为消除序列的趋势同时减少序列的波动，可以对原有时间序列做二阶逐期差分，并绘制差分后的时序图（见图2）。

可以看出经过差分处理后的序列趋势基本上消除。

为了更好地描述月度数据时间序列并进行模拟，需对该序列再进行季节差分，进一步消除季节性（见图3）。

图2居民消费价格指数一阶差分后时序图

NullHypothesis:

Yhasaunitroot

Exogenous:

Constant

LagLength:

3（AutomaticbasedonSIC,MAXLAG=9）

t-Statistic

  Prob.*

AugmentedDickey-Fullerteststatistic

-2.746311

 0.0763

Testcriticalvalues:

1%level

-3.626784

5%level

-2.945842

10%level

-2.611531

*MacKinnon（1996）one-sidedp-values.

AugmentedDickey-FullerTestEquation

DependentVariable:

D（Y）

Method:

LeastSquares

Date:

06/13/12Time:

18:

50

Sample（adjusted）:

2008M052011M04

Includedobservations:

36afteradjustments

Variable

Coefficient

Std.Error

t-Statistic

Prob.  

Y（-1）

-0.111635

0.040649

-2.746311

0.0099

D（Y（-1））

-0.041687

0.152752

-0.272908

0.7867

D（Y（-2））

0.432230

0.138679

3.116762

0.0039

D（Y（-3））

0.444712

0.158945

2.797894

0.0088

C

-0.097779

0.203069

-0.481507

0.6335

R-squared

0.442918

    Meandependentvar

-0.083333

AdjustedR-squared

0.371036

    S.D.dependentvar

1.520432

S.E.ofregression

1.205813

    Akaikeinfocriterion

3.340431

Sumsquaredresid

45.07354

    Schwarzcriterion

3.560365

Loglikelihood

-55.12777

    F-statistic

6.161764

Durbin-Watsonstat

2.126619

    Prob（F-statistic）

0.000907

表3居民消费价格指数一阶差分检验

在表3中，

因而变换后是平稳序列。

5.模型的建立

1）Pandit-wu

DependentVariable:

Y

Method:

LeastSquares

Date:

06/13/12Time:

19:

46

Sample（adjusted）:

2008M032011M03

Includedobservations:

37afteradjustments

Convergenceachievedafter86iterations

Backcast:

OFF（RootsofMAprocesstoolarge）

Variable

Coefficient

Std.Error

t-Statistic

Prob.  

AR

（1）

1.910201

0.041242

46.31669

0.0000

AR

（2）

-0.931860

0.041397

-22.51045

0.0000

MA

（1）

-1.290889

0.152106

-8.486794

0.0000

R-squared

0.967035

    Meandependentvar

-0.597297

AdjustedR-squared

0.965096

    S.D.dependentvar

5.255345

S.E.ofregression

0.981840

    Akaikeinfocriterion

2.878829

Sumsquaredresid

32.77635

    Schwarzcriterion

3.009444

Loglikelihood

-50.25833

    Durbin-Watsonstat

2.658866

InvertedARRoots

 .96+.14i

     .96-.14i

InvertedMARoots

      1.29

EstimatedMAprocessisnoninvertible

再拟合ARMA（4,3）模型：

DependentVariable:

Y

Method:

LeastSquares

Date:

06/13/12Time:

19:

51

Sample（adjusted）:

2008M052011M03

Includedobservations:

35afteradjustments

Convergenceachievedafter39iterations

Backcast:

OFF（RootsofMAprocesstoolarge）

Variable

Coefficient

Std.Error

t-Statistic

Prob.  

AR

（1）

2.058000

0.187697

10.96449

0.0000

AR

（2）

-1.438411

0.380206

-3.783244

0.0007

AR（3）

0.597611

0.506235

1.180500

0.2477

AR（4）

-0.251998

0.287915

-0.875254

0.3889

MA

（1）

-1.695860

0.294125

-5.765777

0.0000

MA

（2）

1.153567

0.364509

3.164717

0.0037

MA（3）

-0.868322

0.338996

-2.561452

0.0161

R-squared

0.975039

    Meandependentvar

-0.911429

AdjustedR-squared

0.969691

    S.D.dependentvar

5.229039

S.E.ofregression

0.910356

    Akaikeinfocriterion

2.826894

Sumsquaredresid

23.20494

    Schwarzcriterion

3.137964

Loglikelihood

-42.47065

    Durbin-Watsonstat

2.402784

InvertedARRoots

 .98+.17i

     .98-.17i

   .05-.50i

 .05+.50i

InvertedMARoots

      1.32

     .19-.79i

   .19+.79i

EstimatedMAprocessisnoninvertible

F检验：

DependentVariable:

Y

Method:

LeastSquares

Date:

06/19/12Time:

10:

40

Sample（adjusted）:

2008M072011M04

Includedobservations:

34afteradjustments

Convergenceachievedafter36iterations

Backcast:

2006M082006M12

Variable

Coefficient

Std.Error

t-Statistic

Prob.  

AR

（1）

0.620067

0.225175

2.753715

0.0113

AR

（2）

0.668297

0.253206

2.639340

0.0147

AR（3）

0.115968

0.289858

0.400086

0.6928

AR（4）

-0.308989

0.289929

-1.065741

0.2976

AR（5）

-0.069233

0.240974

-0.287304

0.7765

AR（6）

-0.132889

0.208902

-0.636129

0.5310

MA

（1）

0.117696

0.077496

1.518724

0.1425

MA

（2）

-0.117747

0.073804

-1.595396

0.1243

MA（3）

0.088579

0.069414

1.276102

0.2147

MA（4）

-0.064442

0.076148

-0.846271

0.4061

MA（5）

-0.939900

0.058310

-16.11892

0.0000

R-squared

0.974955

    Meandependentvar

-1.017647

AdjustedR-squared

0.964067

    S.D.dependentvar

5.269210

S.E.ofregression

0.998838

    Akaikeinfocriterion

3.091744

Sumsquaredresid

22.94657

    Schwarzcriterion

3.585566

Loglikelihood

-41.55964

    Durbin-Watsonstat

2.074405

InvertedARRoots

 .96-.17i

     .96+.17i

  -.04+.50i

-.04-.50i

-.62-.41i

    -.62+.41i

InvertedMARoots

      .98

     .29+.94i

   .29-.94i

-.84-.53i

-.84+.53i

F检验：

故ARMA（2,1）与ARMA（4,3）模型有显著性差异，ARMA（6,5）与ARMA（4，3）无显著差异，则选择ARMA（4,3）模型，即

2）Box-Jenkins识别模型

Autocorrelations

Series:

cpi

Lag

Autocorrelation

Std.Errora

Box-LjungStatistic

Value

df

Sig.b

1

.253

.136

3.454

1

.063

2

.374

.135

11.181

2

.004

3

.385

.133

19.530

3

.000

4

.206

.132

21.980

4

.000

5

.195

.130

24.213

5

.000

6

.179

.129

26.137

6

.000

7

.088

.128

26.614

7

.000

8

.130

.126

27.671

8

.001

9

-.114

.125

28.509

9

.001

10

-.332

.123

35.794

10

.000

11

.097

.122

36.434

11

.000

12

-.481

.120

52.467

12

.000

13

-.183

.119

54.846

13

.000

14

-.160

.117

56.721

14

.000

15

-.276

.115

62.444

15

.000

16

-.257

.114

67.545

16

.000

17

-.246

.112

72.351

17

.000

18

-.261

.110

77.941

18

.000

19

-.127

.