13季节ARIMA模型

12.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)-dq(L)vt(2.62)把(2.62)式代入(2.60)式，于是得到季节时间序列模型的一般表达式。p(L)P(Ls)(dsDyt)=q(L)Q(Ls)vt(2.63)其中下标P,Q,p,q分别表示季节与非季节自回归、移动平均算子的最大滞后阶数，d,D分别表示非季节和季节性差分次数。上式称作(p,d,q)(P,D,Q)s阶季节时间序列模型或乘积季节模型。保证（dsDyt）具有平稳性的条件是p(L)P(Ls)=0的根在单位圆外；保证（dsDyt）具有可逆性的条件是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估计命令是2DLOG(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+11L13vt=vt+1vt–1+1vt–12+11vt–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=11。进一步化简(yt–yt-12)=vt+1vt–1+1vt–12+11vt–13yt–yt-12=vt+1vt–1+1vt–12+11vt–13用于预测的模型型式是yt=yt-1+yt-12–yt–13+vt+1vt–1+1vt–12+11vt–13(2.65)从上式可以看出SARIMA模型可以展开为ARIMA模型。对乘积季节模型的季节阶数，即周期长度s的识别可以通过对实际问题的分析、时间序列图以及时间序列的相关图和偏相关图分析得到。以相关图和偏相关图为例，如果相关图和偏相关图不是呈线性衰减趋势，而是在变化周期的整倍数时点上出现绝对值相当大的峰值并呈振荡式变化，就可以认为该时间序列可以用SARIMA模型描述。建立SARIMA模型，（1）首先要确定d,D。通过差分和季节差分把原序列变换为一个平稳的序列。令xt=dsDyt（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与时间近似呈线性关系（异方差问题也得到抑制）。302004006008001000787980818283848586878889Y4.55.05.56.06.57.0787980818283848586878889LNY图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倍数的滞后期存在显著的自相关关系。-0.4-0.20.00.20.4787980818283848586878889DLNY-0.4-0.20.00.20.4787980818283848586878889D2LNY图2.35Lnyt图2.362Lnyt4图2.37Lnyt的相关图（下）和偏相关图（上）对Lnyt进行一次季节性差分（或12阶差分），得12Lnyt（图2.38）。从12Lnyt的相关图和偏相关图（图2.39）可以看到12Lnyt仍然是非平稳的。-0.4-0.20.00.20.4787980818283848586878889SDLNY图2.3812Lnyt，（EViews：DLOG(Y,0,12)）图2.3912Lnyt的相关图（下）和偏相关图（上）对Lnyt进行一阶差分和一阶季节性差分，得12Lnyt（见图2.40）。从xt的相关图和偏相关图（见图2.41）可以看到12Lnyt近似为一个平稳过程。5-0.4-0.20.00.20.4787980818283848586878889DSDLNY图2.4012Lnyt=xt，（EViews：DLOG(Y,1,12)）图2.4112Lnyt的相关图（下）和偏相关图（上）用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)=44EViews估计命令是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。6图2.42EViews估计结果图2.43模型残差序列的相关与偏相关图对于12Lnyt来，模型参数全部有显著性，Q36=15.520.05(36-2-1)=44。两种检验通过。见输出结果（2.42），对于12Lnyt，模型共有14个特征根。7图2.44D12DLnyt的实际与预测序列图2.45yt的实际与预测序列对1989年第12月份yt进行样本外1期预测，结果如图2.46。图2.46EViews预测结果预测误差是=1.7341.7349.789=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–13Lnyt=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年：月yt1978:01134.31980:06168.21982:11205.81985:04343.41987:09499.51978:02119.41980:07163.51982:12248.21985:05341.21987:10505.21978:03128.31980:08161.61983:01243.21985:06346.01987:11518.71978:04126.41980:09172.91983:02217.51985:07329.91987:12617.91978:05128.81980:10166.51983:03226.21985:08328.11988:01570.71978:06127.81980:11175.21983:04223.51985:09358.21988:02561.31978:07121.11980:12197.71983:05221.01985:10358.41988:03570.41978:08118.41981:01212.11983:06220.51985:11376.61988:04567.91978:09125.71981:02177.91983:07205.81985:12451.01988:05570.91978:10123.61981:03182.91983:08206.91986:01412.01988:06603.91978:11128.51981:04184.21983:09218.81986:02374.51988:07591.