居民消费价格指数的时间序列分析讲解.docx
- 文档编号:5010242
- 上传时间:2022-12-12
- 格式:DOCX
- 页数:27
- 大小:290.90KB
居民消费价格指数的时间序列分析讲解.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
.
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 居民消费 价格指数 时间 序列 分析 讲解
![提示](https://static.bdocx.com/images/bang_tan.gif)