13季节ARIMA模型.docx
- 文档编号:23322297
- 上传时间:2023-05-16
- 格式:DOCX
- 页数:19
- 大小:233.88KB
13季节ARIMA模型.docx
《13季节ARIMA模型.docx》由会员分享,可在线阅读,更多相关《13季节ARIMA模型.docx(19页珍藏版)》请在冰豆网上搜索。
13季节ARIMA模型
2.8季节时间序列模型
在某些时间序列中,存在明显的周期性变化。
这种周期是由于季节性变化(包括季度、月度、周度等变化)或其他一些固有因素引起的。
这类序列称为季节性序列。
比如一个地区的气温值序列(每隔一小时取一个观测值)中除了含有以天为周期的变化,还含有以年为周期的变化。
在经济领域中,季节性序列更是随处可见。
如季度时间序列、月度时间序列、周度时间序列等。
处理季节性时间序列只用以上介绍的方法是不够的。
描述这类序列的模型之一是季节时间序列模型(seasonalARIMAmodel),用SARIMA表示。
较早文献也称其为乘积季节模型(multiplicativeseasonalmodel)。
设季节性序列(月度、季度、周度等序列都包括其中)的变化周期为s,即时间间隔为s的观测值有相似之处。
首先用季节差分的方法消除周期性变化。
季节差分算子定义为,
s=1-Ls
若季节性时间序列用yt表示,则一次季节差分表示为
syt=(1-Ls)yt=yt-yt-s
对于非平稳季节性时间序列,有时需要进行D次季节差分之后才能转换为平稳的序列。
在此基础上可以建立关于周期为s的P阶自回归Q阶移动平均季节时间序列模型(注意P、Q等于2时,滞后算子应为(Ls)2=L2s。
P(Ls)sDyt=Q(Ls)ut(2.60)
对于上述模型,相当于假定ut是平稳的、非自相关的。
当ut非平稳且存在ARMA成分时,则可以把ut描述为
p(L)dut=q(L)vt(2.61)
其中vt为白噪声过程,p,q分别表示非季节自回归、移动平均算子的最大阶数,d表示ut的一阶(非季节)差分次数。
由上式得
ut=p-1(L)-dq(L)vt(2.62)
把(2.62)式代入(2.60)式,于是得到季节时间序列模型的一般表达式。
p(L)P(Ls)(dsDyt)=q(L)Q(Ls)vt(2.63)
其中下标P,Q,p,q分别表示季节与非季节自回归、移动平均算子的最大滞后阶数,d,D分别表示非季节和季节性差分次数。
上式称作(p,d,q)(P,D,Q)s阶季节时间序列模型或乘积季节模型。
保证(dsDyt)具有平稳性的条件是p(L)P(Ls)=0的根在单位圆外;保证(dsDyt)具有可逆性的条件是q(L)Q(Ls)=0的根在单位圆外。
当P=D=Q=0时,SARIMA模型退化为ARIMA模型;从这个意义上说,ARIMA模型是SARIMA模型的特例。
当P=D=Q=p=q=d=0时,SARIMA模型退化为白噪声模型。
(1,1,1)(1,1,1)12阶月度SARIMA模型表达为
(1-1L)(1-1L12)12yt=(1+1L)(1+1L12)vt
12yt具有平稳性的条件是1<1,1<1,12yt具有可逆性的条件是1<1,1<1。
设log(Yt)=yt,变量12yt在EViews中用DLOG(Y,1,12)表示(这样表示的好处是EViews可以直接预测到Y),上式的EViews估计命令是
DLOG(Y,1,12)AR
(1)SAR(12)MA
(1)SMA(12)
(0,1,1)(0,1,1)12阶月度SARIMA模型表达为
12yt=(1+1L)(1+1L12)vt(2.64)
(2.64)式的EViews估计命令是
DLOG(Y,1,12)MA
(1)SMA(12)
由(2.64)式得
12yt=(1+1L)(1+1L12)vt=vt+1Lvt+1L12vt+11L13vt
=vt+1vt–1+1vt–12+11vt–13
上式对应的EViews估计命令是
DLOG(Y,1,12)MA
(1)MA(12)MA(13)
模型表达式是
12yt=vt+1vt–1+12vt–12+13vt–13
这是一个非季节模型表达式。
以上两个EViews估计命令是等价的,都是估计MA(13)模型。
注意:
唯一不同点是上式对vt–13的系数没有约束,而对季节模型来说,相当于增加了一个约束条件,13=11。
进一步化简
(yt–yt-12)=vt+1vt–1+1vt–12+11vt–13
yt–yt-12=vt+1vt–1+1vt–12+11vt–13
用于预测的模型型式是
yt=yt-1+yt-12–yt–13+vt+1vt–1+1vt–12+11vt–13(2.65)
从上式可以看出SARIMA模型可以展开为ARIMA模型。
对乘积季节模型的季节阶数,即周期长度s的识别可以通过对实际问题的分析、时间序列图以及时间序列的相关图和偏相关图分析得到。
以相关图和偏相关图为例,如果相关图和偏相关图不是呈线性衰减趋势,而是在变化周期的整倍数时点上出现绝对值相当大的峰值并呈振荡式变化,就可以认为该时间序列可以用SARIMA模型描述。
建立SARIMA模型,
(1)首先要确定d,D。
通过差分和季节差分把原序列变换为一个平稳的序列。
令
xt=dsDyt
(2)然后用xt建立p(L)P(Ls)xt=q(L)Q(Ls)vt模型。
注意:
(1)用对数的季节时间序列数据建模时通常D不会大于1,P和Q不会大于3。
(2)乘积季节模型参数的估计、检验与前面介绍的估计、检验方法相同。
利用乘积季节模型预测也与上面介绍的预测方法类似。
2.9季节时间序列建模案例
案例1:
(文件名:
b2c3,5b2c3)北京市1978:
1~1989:
12社会商品零售额月度数据(yt,单位:
亿元人民币)曲线见图2.32,数据见表2.3。
yt与时间呈指数关系且存在递增型异方差。
对数的社会商品零售额月度数据(Lnyt)曲线见图2.33。
Lnyt与时间近似呈线性关系(异方差问题也得到抑制)。
图2.32yt图2.33Lnyt
通过Lnyt的相关图和偏相关图(见图2.34)可以看到Lnyt是一个非平稳序列(相关图衰减很慢)且Lnyt与其12倍数的滞后期存在自回归关系。
图2.34Lnyt的相关图(下)和偏相关图(上)
对Lnyt进行一阶差分,得Lnyt(图2.35)。
图2.36是对Lnyt进行2次一阶差分的结果,序列2Lnyt是过度差分序列。
从Lnyt的相关图和偏相关图(图2.37)可以看到,通过差分Lnyt的平稳性得到很大改进,但与其12倍数的滞后期存在显著的自相关关系。
图2.35Lnyt图2.362Lnyt
图2.37Lnyt的相关图(下)和偏相关图(上)
对Lnyt进行一次季节性差分(或12阶差分),得12Lnyt(图2.38)。
从12Lnyt的相关图和偏相关图(图2.39)可以看到12Lnyt仍然是非平稳的。
图2.3812Lnyt,(EViews:
DLOG(Y,0,12))
图2.3912Lnyt的相关图(下)和偏相关图(上)
对Lnyt进行一阶差分和一阶季节性差分,得12Lnyt(见图2.40)。
从xt的相关图和偏相关图(见图2.41)可以看到12Lnyt近似为一个平稳过程。
图2.4012Lnyt=xt,(EViews:
DLOG(Y,1,12))
图2.4112Lnyt的相关图(下)和偏相关图(上)
用1978:
1~1989:
11期间数据,估计yt的(1,1,1)(1,1,0)12阶季节时间序列模型,得结果如下:
(1+0.5924L)(1+0.4093L12)12Lnyt=(1+0.4734L)vt(2.66)
(4.5)(5.4)(1.9)
R2=0.33,s.e.=0.146,Q36=15.5,20.05(36-2-1)=44
EViews估计命令是
DLOG(Y,1,12)AR
(1)SAR(12)MA
(1)
EViews输出结果见图2.42。
注意:
(1)仔细对照(2.66)式和图2.42输出结果,不要把自回归系数估计值的符号写错。
通过自回归特征根倒数-0.59可知,把表达式中的算子写作(1+0.5924L)是正确的。
通过移动平均特征根倒数-0.47可知,把表达式中的算子写作(1+0.4734L)是正确的。
(2)表达式中,季节和非季节因子(特征多项式)之间是相乘关系。
(3)在EViews估计命令中把变量写作DLOG(Y,1,12)的好处是可以直接对yt和12Lnyt预测。
模型残差序列的相关与偏相关图如图2.43。
图2.42EViews估计结果
图2.43模型残差序列的相关与偏相关图
对于12Lnyt来,模型参数全部有显著性,Q36=15.5<20.05(36-2-1)=44。
两种检验通过。
见输出结果(2.42),对于12Lnyt,模型共有14个特征根。
图2.44D12DLnyt的实际与预测序列图2.45yt的实际与预测序列
对1989年第12月份yt进行样本外1期预测,结果如图2.46。
图2.46EViews预测结果
预测误差是
=
=0.076
用1978:
1~1989:
12期间数据得EViews(0,1,1)(0,1,1)12模型估计结果如下
12Lnyt=(1-0.35L)(1-0.61L12)vt(2.67)
(-4.4)(-9.1)
R2=0.36,DW=1.86,F=71.9,s.e.=0.038,Q36=21.88,20.05(36-2)=44
模型参数全部有显著性,Q36=21.88<20.05(36-2)=44。
两种检验通过。
上式变换为,
Lnyt–Lnyt-12=vt-0.35vt–1-0.61vt–12+0.2135vt–13
Lnyt=Lnyt-1+Lnyt-12–Lnyt–13+vt-0.35vt–1-0.61vt–12+0.2135vt–13(2.68)
(2.67)式也是一个可以选用的模型。
表2.3北京市社会商品零售额(yt)月度数据(单位:
亿元人民币,1978:
1~1989:
12)
年:
月
yt
年:
月
yt
年:
月
yt
年:
月
yt
年:
月
yt
1978:
01
134.3
1980:
06
168.2
1982:
11
205.8
1985:
04
343.4
1987:
09
499.5
1978:
02
119.4
1980:
07
163.5
1982:
12
248.2
1985:
05
341.2
1987:
10
505.2
1978:
03
128.3
1980:
08
161.6
1983:
01
243.2
1985:
06
346.0
1987:
11
518.7
1978:
04
126.4
1980:
09
172.9
1983:
02
217.5
1985:
07
329.9
1987:
12
617.9
1978:
05
128.8
1980:
10
166.5
1983:
03
226.2
1985:
08
328.1
1988:
01
570.7
1978:
06
127.8
1980:
11
175.2
1983:
04
223.5
1985:
09
358.2
1988:
02
561.3
1978:
07
121.1
1980:
12
197.7
1983:
05
221.0
1985:
10
358.4
1988:
03
570.4
1978:
08
118.4
1981:
01
212.1
1983:
06
220.5
1985:
11
376.6
1988:
04
567.9
1978:
09
125.7
1981:
02
177.9
1983:
07
205.8
1985:
12
451.0
1988:
05
570.9
1978:
10
123.6
1981:
03
182.9
1983:
08
206.9
1986:
01
412.0
1988:
06
603.9
1978:
11
128.5
1981:
04
184.2
1983:
09
218.8
1986:
02
374.5
1988:
07
591.8
1978:
12
145.2
1981:
05
184.0
1983:
10
216.0
1986:
03
390.0
1988:
08
636.2
1979:
01
164.7
1981:
06
182.4
1983:
11
235.0
1986:
04
387.0
1988:
09
674.5
1979:
02
126.2
1981:
07
175.6
1983:
12
282.0
1986:
05
389.8
1988:
10
647.7
1979:
03
143.7
1981:
08
172.0
1984:
01
268.4
1986:
06
397.7
1988:
11
640.5
1979:
04
143.7
1981:
09
184.9
1984:
02
227.6
1986:
07
381.4
1988:
12
804.2
1979:
05
145.5
1981:
10
184.7
1984:
03
248.6
1986:
08
386.9
1989:
01
694.3
1979:
06
143.7
1981:
11
195.1
1984:
04
247.0
1986:
09
429.8
1989:
02
673.8
1979:
07
138.4
1981:
12
224.8
1984:
05
249.9
1986:
10
428.8
1989:
03
718.7
1979:
08
136.7
1982:
01
233.6
1984:
06
253.1
1986:
11
444.4
1989:
04
690.3
1979:
09
145.5
1982:
02
182.0
1984:
07
245.5
1986:
12
527.7
1989:
05
676.6
1979:
10
150.7
1982:
03
206.6
1984:
08
249.6
1987:
01
478.3
1989:
06
665.8
1979:
11
149.0
1982:
04
202.2
1984:
09
272.3
1987:
02
442.4
1989:
07
642.2
1979:
12
164.7
1982:
05
201.7
1984:
10
278.7
1987:
03
461.4
1989:
08
638.9
1980:
01
190.3
1982:
06
202.6
1984:
11
299.4
1987:
04
458.2
1989:
09
674.1
1980:
02
174.9
1982:
07
192.8
1984:
12
366.3
1987:
05
458.2
1989:
10
652.7
1980:
03
163.2
1982:
08
186.2
1985:
01
364.8
1987:
06
468.5
1989:
11
641.9
1980:
04
168.4
1982:
09
199.3
1985:
02
349.1
1987:
07
454.5
1989:
12
734.1
1980:
05
168.6
1982:
10
198.2
1985:
03
359.1
1987:
08
458.9
案例2香港季节GDPt数据的拟合(季节时间序列模型,file:
HongKong)
1980:
1~2002:
4年香港季度GDPt序列曲线见图2.27(数据见表2.4,单位:
港元)。
1980~1997年GDPt随时间呈指数增长。
1997年由于遭受东南亚金融危机的影响,经济发展处于停滞状态,1998~2002年底GDPt总量几乎没有增长。
另一个特征是GDPt随时间呈递增型异方差。
所以,用对数的季度GDPt数据(LnGDPt,曲线见图2.48)建立季节时间序列模型。
图2.47GDPt图2.48LnGDPt
通过LnGDPt的相关图和偏相关图(图2.49)可以看到LnGDPt是一个非平稳序列(相关图衰减得很慢)。
图2.49LnGDPt的相关图和偏相关图
对LnGDPt进行一阶差分,得DLnGDPt(见图2.50)。
DLnGDPt的平稳性得到很大改进,但其季节因素影响还很大。
从DLnGDPt的相关图和偏相关图(图2.51)也可以明显地看到这个特征。
若对LnGDPt直接进行一次季节差分(四阶差分),得D4LnGDPt见图2.52。
其波动性也很大。
相关图和偏相关图见图2.53。
D2LnGDPt显然是过度差分序列(图2.54)。
图2.50DLnGDPt,(s.d.=0.062)图2.51DLnGDPt的相关图和偏相关图
图2.52D4LnGDPt,(s.d.=0.076)图2.53D4LnGDPt的相关图和偏相关图
图2.54D2LnGDPt,(s.d.=0.062)
在DLnGDPt的基础上进行一阶季节差分,或在D4LnGDPt基础上进行一阶非季节差分,得D4DLnGDPt(图2.55)。
其相关图和偏相关图见图2.56。
D4DLnGDPt中已经基本消除了季节变化因素。
在D4DLnGDPt的基础上建立时间序列模型。
图2.55D4DLnGDPt,(s.d.=0.029)图2.56D4DLnGDPt的相关和偏相关图
通过对D4DLnGDPt的相关和偏相关图分析,应该建立(2,1,2)(1,1,1)4模型。
EViews估计命令是
DLOG(GDP,1,4)CAR
(1)AR
(2)SAR(4)MA
(1)MA
(2)SMA(4)
用1980:
1~2002:
3的数据得估计结果如下(EViews输出结果如图2.57):
D4DLnGDPt=-0.0023+ut(1980:
1~2002:
3)
(-2.4)
(1-1.20L+0.66L2)(1-0.33L4)ut=(1-1.16L+0.97L2)(1-0.95L4)vt(2.69)
(14.4)(-8.8)(2.8)(55.9)(86.1)(-32.8)
R2=0.57,DW=2.0,F=16.1,Q36=19.3,20.05(36-3-3-1)=42.6
图2.57EViews估计结果
图2.57模型(2.99)误差项的相关和偏相关图
注意:
(1)不要把自回归系数估计值的符号写错。
不要把均值(-0.0023)项表达错。
EViews仍然是对(D4DLnGDPt+0.0023)建立(2,1,2)(1,1,1)4阶季节时间序列模型,而不是对D4DLnGDPt建立季节时间序列模型。
(2)季节和非季节因子之间是相乘关系。
(3)在EViews估计命令中把变量写作DLOG(GDP,1,4),好处是预测时可直接预测GDPt,也可以预测D4DLnGDPt。
模型参数全部有显著性,Q36=19.6<20.05(36-3-3-1)=42.6。
两种检验通过。
依据输出结果,对于D4DLnGDPt,模型共有12个特征根。
4个实根,8个复根。
图2.58D4DLnGDPt的实际与预测序列图2.59GDPt的实际与预测序列
对2002年第4季度GDPt进行样本外1期预测,结果如下:
预测误差是
=
=0.006
表2.4香港季度GDPt数据(1980:
1~2002:
4,单位:
港元)
年:
月
GDPt/1011
年:
月
GDPt/1011
年:
月
GDPt/1011
年:
月
GDPt/1011
1980:
1
0.31489
1985:
4
0.70951
1991:
3
1.80580
1997:
2
3.32335
1980:
2
0.34505
1986:
1
0.69108
1991:
4
1.86722
1997:
3
3.51501
1980:
3
0.37641
1986:
2
0.73102
1992:
1
1.76661
1997:
4
3.53431
1980:
4
0.38566
1986:
3
0.83372
1992:
2
1.89129
1998:
1
3.09548
1981:
1
0.38830
1986:
4
0.88419
1992:
3
2.10027
1998:
2
3.18978
1981:
2
0.40632
1987:
1
0.84318
1992:
4
2.15501
1998:
3
3.26421
1981:
3
0.44307
1987:
2
0.89861
1993:
1
2.05568
1998:
4
3.24905
1981:
4
0.47473
1987:
3
1.05262
1993:
2
2.18604
1999:
1
2.90525
1982:
1
0.45042
1987:
4
1.06886
1993:
3
2.41488
1999:
2
3.05798
1982:
2
0.45973
1988:
1
1.00261
1993:
4
2.47149
1999:
3
3.19225
1982:
3
0.50703
1988:
2
1.06898
1994:
1
2.36370
1999:
4
3.30586
1982:
4
0.51368
1988:
3
1.22224
1994:
2
2.49299
2000:
1
3.08457
1983:
1
0.47199
1988:
4
1.27865
1994:
3
2.68834
2000:
2
3.12780
1983:
2
0.50214
1989:
1
1.18206
1994:
4
2.75271
2000:
3
3.30500
1983:
3
0.55978
1989:
2
1.26058
1995:
1
2.53532
2000:
4
3.36602
1983:
4
0.60088
1989:
3
1.39551
1995:
2
2.66221
2001:
1
3.0
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 13 季节 ARIMA 模型