1.4估计CPI的ARIMA模型(1 / 1)

由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中不必再处理数据。

【实验步骤】