由11.3知CPI序列为1阶单整,于是可以对1阶差分后的CPI序列估计其ARMA(p,q)模型,即得原序列CPI的ARIMA(p,1,q)模型。
ARMA(p,q)的定阶是通过看样本的自相关函数跟偏自相关函数图的,于是我们先要画出1阶差分后CPI序列的样本的自相关函数、偏自相关函数图。点击主界面Quick→Series Statistics→Correlogram...,在弹出的Series对话框中输入CPI,点击OK,就会弹出Correlogram对话框,在其中Correlogram of栏中选择1st(表示1阶差分),Lags栏中就默认为12阶(表示输出阶数为12阶),点击OK,即可得到1阶差分后CPI序列的样本的自相关函数、偏自相关函数图,如图1110所示。
图1110
我们可以看到偏自相关函数图中1阶跟2阶是有明显的尖柱的,而自相关函数图中1阶是有明显的尖柱的,不妨初步定阶为p=1,q=2。即使用ARIMA(2,1,1)对原序列进行识别。
点击主界面菜单Quick→Estimate Equation...,在弹出的对话框中输入:D(CPI) C D(CPI(-1)) D(CPI(-2)) MA(1),点击确定,即可得到估计结果如图1111所示。
得到估计后,还要对序列残差进行1阶LM检验,在图1111中点击View→Residual Tests→Serial Correlation LM Test...,在弹出的对话框中输入1,点击OK,即得到序列残差的1阶LM检验结果如图1112所示。
图1111
图1112
根据图1111得到模型的估计结果为:
ΔCPIt=3.963+0.468ΔCPIt-1-0.018ΔCPIt-2+εt+0.981εt-1
AIC=6.011132SC=6.204685LM(1)=0.361366
图1113
但是发现,在5%的显著性水平下,参数都没有通过检验,特别是ΔCPIt-2项的系数,不妨去掉ΔCΡΙt-2项再对序列进行建模,即采用ARIMA(1,1,1)对CPI序列进行估计。
点击主界面菜单Quick→Estimate Equation...,在弹出的对话框中输入:D(CPI) C D(CPI(-1)) MA(1),点击确定,即可得到估计结果如图1112所示。
得到估计后,还要对序列残差进行1阶LM检验,在图1112中点击View→Residual Tests→Serial Correlation LM Test...,在弹出的对话框中输入1,点击OK,即得到序列残差的1阶LM检验结果,如图1114所示。
图1114
根据图1112数据,得到模型的估计结果为:
ΔCPIt=4.517831+0.430897ΔCPIt-1+εt+0.997487εt-1
AIC=5.889310SC=6.033292LM(1)=0.511508
可以看到两个模型的残差都是不相关的,即可以认为残差为随机序列,即模型均是基本符合的。但是比较发现,ARIMA(2,1,1)的模型中存在参数不显著,但是ARIMA(1,1,1)模型中的参数均显著。还有可以通过赤池准则,比较两个模型的AIC值跟SC值,ARIMA(1,1,1)模型的AIC值跟SC值均比ARIMA(2,1,1)的模型的AIC值跟SC值小,所以得到ARIMA(1,1,1)模型为较佳的模型的结论。
综上所述,最后我们选择ARIMA(1,1,1)模型对CPI序列进行估计,得到CPI序列其模型的口径为:
CPIt=4.517831-0.569103CPIt-1-0.430897CPIt-2+εt+0.997487εt-1
【实验内容二】
对中国进出口贸易进行协整分析,如果存在协整关系,则建立ECM模型。
通过查询中国统计年鉴,可以获取1951—2003年中国进口(im)、出口(ex)和商品零售物价指数(pt)。因为该期间物价变化大,特别是改革开放以后变化更为激烈,所以物价指数也作为一个解释变量加入模型中。为消除物价变动对进出口数据的影响以及消除进出口数据中存在的异方差,定义三个变量如下,数据如表112所示:
lnex=ln(ex/price)(自然对数的出口额,不变价价格1990=1)
lnim=ln(im/price)(自然对数的进口额,不变价价格1990=1)
lnpt=ln(price)(自然对数的商品价格指数)
表112中国进口、出口与商品零售物价指数
YearlnexlnimlnptYearlnexlnimlnpt
19514.1084.485-0.921
19524.2264.551-0.926
19534.4414.772-0.892
19544.5594.670-0.870
19554.7464.973-0.860
19564.8804.831-0.860
19574.8424.756-0.844
19585.0464.964-0.842
19595.1905.098-0.832
19604.9494.977-0.801
19614.5174.413-0.652
19624.4674.135-0.614
19634.5874.251-0.675
19644.7284.453-0.713
19654.8854.753-0.740
19664.9324.855-0.742
19674.8254.729-0.751
19684.8644.681-0.751
19694.8504.614-0.759
19704.8034.791-0.764
19714.9994.731-0.772
19725.1924.933-0.774
19735.5295.408-0.768
19745.6995.791-0.761
19755.7245.755-0.761
19765.6615.619-0.757
19775.6785.627-0.738
19785.8515.963-0.730
19796.0676.204-0.711
19806.2556.352-0.652
19816.5366.537-0.629
19826.6366.490-0.611
19836.6796.641-0.596
19846.9316.998-0.567
19857.1797.620-0.483
19867.4117.737-0.425
19877.6477.740-0.354
19887.6627.813-0.185
19897.6007.717-0.021
19908.0027.8530.000
19918.2228.1030.028
19928.3698.3180.081
19938.3688.4920.205
19948.8518.8050.401
19958.8918.7710.539
19968.8418.7570.598
19979.0208.7700.606
19989.0518.7810.580
19999.1418.9780.549
20009.4009.2990.534
20019.4749.3850.526
20029.6889.5900.513
20039.9879.9270.512
数据来源:《中国统计年鉴》
说明:表中数据已为不变价的对数序列,所以在EViews中不必再处理数据。
【实验步骤】