13季节ARIMA模型.docx

1、13季节ARIMA模型 季节时间序列模型在某些时间序列中，存在明显的周期性变化。这种周期是由于季节性变化包括季度、月度、周度等变化或其他一些固有因素引起的。这类序列称为季节性序列。比方一个地区的气温值序列每隔一小时取一个观测值中除了含有以天为周期的变化，还含有以年为周期的变化。在经济领域中，季节性序列更是随处可见。如季度时间序列、月度时间序列、周度时间序列等。处理季节性时间序列只用以上介绍的方法是不够的。描述这类序列的模型之一是季节时间序列模型(seasonal ARIMA model)，用SARIMA表示。较早文献也称其为乘积季节模型multiplicative seasonal model

2、。设季节性序列月度、季度、周度等序列都包括其中的变化周期为s，即时间间隔为s的观测值有相似之处。首先用季节差分的方法消除周期性变化。季节差分算子定义为， s = 1- Ls 假设季节性时间序列用yt表示，那么一次季节差分表示为 s yt = (1- Ls) yt = yt - yt - s 对于非平稳季节性时间序列，有时需要进行D次季节差分之后才能转换为平稳的序列。在此根底上可以建立关于周期为s的P阶自回归Q阶移动平均季节时间序列模型注意P、Q等于2时，滞后算子应为(Ls)2 = L2s。 AP (Ls) sDyt = BQ (Ls) ut (2.60)对于上述模型，相当于假定ut是平稳的、非

3、自相关的。当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) AP(Ls) (dsDyt) = q(L) BQ(Ls) vt (2.63)其中下标P, Q, p, q分别表示季节与非季节自回归、移动平均算子的最大滞后阶数，d, D分别表示非季节和季节性差分次数。上式称作

4、(p, d, q) (P, D, Q)s 阶季节时间序列模型或乘积季节模型。保证dsDyt具有平稳性的条件是p(L)AP(Ls) = 0的根在单位圆外；保证dsDyt具有可逆性的条件是q (L)BQ (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- 1 L) (1- 1 L12) 12 yt = (1+1 L) (1

5、+1 L12) vt 12 yt具有平稳性的条件是 | 1 | 1，| 1 | 1， 12 yt具有可逆性的条件是 | 1 | 1，| 1 | 1。设log(Yt) = yt，变量 12 yt在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模型表达为 12 yt = (1+ 1 L) (1+ 1 L12) vt (2.64)(2.64) 式的EViews估计命令是 DLOG(Y,

6、1,12) MA(1) SMA(12) 由(2.64) 式得 12 yt = (1+1 L) (1+1 L12 ) vt = vt +1 L vt +1 L12vt + 1 1 L13vt = vt +1 vt 1 +1 vt 12 + 1 1 vt 13 上式对应的EViews估计命令是DLOG(Y,1,12) MA(1) MA(12) MA(13)模型表达式是 12 yt = vt +1 vt 1 +12 vt 12 + 13 vt 13这是一个非季节模型表达式。以上两个EViews估计命令是等价的，都是估计MA(13)模型。注意：唯一不同点是上式对vt 13的系数没有约束，而对季节模型来

7、说，相当于增加了一个约束条件，13 =1 1。进一步化简 (yt yt - 12) = vt +1 vt 1 +1 vt 12 + 1 1 vt 13 yt yt - 12 = vt +1 vt 1 +1 vt 12 + 1 1 vt 13用于预测的模型型式是 yt = yt -1 + yt - 12 yt 13 + vt +1 vt 1 +1 vt 12 + 1 1 vt 13 (2.65)从上式可以看出SARIMA模型可以展开为ARIMA模型。 对乘积季节模型的季节阶数，即周期长度s的识别可以通过对实际问题的分析、时间序列图以及时间序列的相关图和偏相关图分析得到。以相关图和偏相关图为例，如

8、果相关图和偏相关图不是呈线性衰减趋势，而是在变化周期的整倍数时点上出现绝对值相当大的峰值并呈振荡式变化，就可以认为该时间序列可以用SARIMA模型描述。 建立SARIMA模型，1首先要确定d, D。通过差分和季节差分把原序列变换为一个平稳的序列。令 xt = dsD yt 2然后用xt 建立 p (L) AP (Ls) xt = q (L) BQ (Ls) vt模型。注意：1用对数的季节时间序列数据建模时通常D不会大于1，P和Q不会大于3。2乘积季节模型参数的估计、检验与前面介绍的估计、检验方法相同。利用乘积季节模型预测也与上面介绍的预测方法类似。2.9 季节时间序列建模案例 案例1：文件名：

9、b2c3，5b2c3北京市1978:11989:12社会商品零售额月度数据yt，单位：亿元人民币曲线见图2.32，数据见表2.3。yt与时间呈指数关系且存在递增型异方差。对数的社会商品零售额月度数据Ln yt曲线见图2.33。Lnyt与时间近似呈线性关系异方差问题也得到抑制。 图2.32 yt 图2.33 Lnyt通过Lnyt的相关图和偏相关图见可以看到Lnyt是一个非平稳序列相关图衰减很慢且Lnyt与其12倍数的滞后期存在自回归关系。图2.34 Lnyt的相关图下和偏相关图上对Lnyt进行一阶差分，得Lnyt。Lnyt进行2次一阶差分的结果，序列2Ln yt是过度差分序列。从 Lnyt的相关

10、图和偏相关图可以看到，通过差分 Lnyt的平稳性得到很大改良，但与其12倍数的滞后期存在显著的自相关关系。 图2.35 Ln yt 图2.36 2Ln yt图2.37 Lnyt的相关图下和偏相关图上对Lnyt进行一次季节性差分或12阶差分，得 12 Lnyt。从 12 Lnyt的相关图和偏相关图图2.39可以看到 12 Lnyt仍然是非平稳的。 图2.38 12 Lnyt，EViews：DLOG(Y,0,12) 图2.39 12 Lnyt的相关图下和偏相关图上 对Lnyt进行一阶差分和一阶季节性差分，得12 Lnyt见图2.40。从xt 的相关图和偏相关图见图2.41可以看到12 Lnyt近似

11、为一个平稳过程。图2.40 12 Lnyt = xt，EViews：DLOG(Y,1,12)图2.41 12 Lnyt的相关图下和偏相关图上 用1978:11989:11期间数据，估计yt 的 (1, 1, 1) (1, 1, 0)12阶季节时间序列模型，得结果如下： (1+ 0.5924 L) (1 + 093 L12) 12Lnyt = (1+0.4734 L) vt (2.66)(4.5) (5.4) ()R2 = 0.33, s.e. = 0.146, Q36 = 15.5, 20.05(36-2-1) = 44EViews估计命令是DLOG(Y,1,12) AR(1) SAR(12)

12、 MA(1)EViews输出结果见图2.42。注意：1仔细对照2.66式和图2.42输出结果，不要把自回归系数估计值的符号写错。通过自回归特征根倒数-0.59可知，把表达式中的算子写作(1+ 0.5924 L)是正确的。通过移动平均特征根倒数-0.47可知，把表达式中的算子写作(1+0.4734 L) 是正确的。2表达式中，季节和非季节因子特征多项式之间是相乘关系。 3在EViews估计命令中把变量写作DLOG(Y,1,12)的好处是可以直接对yt和12 Lnyt预测。模型残差序列的相关与偏相关图如图2.43。图2.42 EViews估计结果 对于12 Lnyt来，模型参数全部有显著性，Q36

13、 = 15.5 20.05(36-2-1) = 44。两种检验通过。见输出结果2.42，对于12 Lnyt，模型共有14个特征根。 图2.44 D12DLnyt的实际与预测序列 图2.45 yt的实际与预测序列对1989年第12月份yt进行样本外1期预测，结果如图2.46。图2.46 EViews预测结果预测误差是 = 0.076用1978:11989:12期间数据得EViews (0, 1, 1) (0, 1, 1)12 模型估计结果如下 12 Lnyt = (1- L) (1 - L12) vt (2.67)(- 4.4) (- 9.1) R2 = 0.36, DW = 1.86, F =

14、 71.9, . = 0.038, Q36 = 2, 20.05 (36-2) = 44模型参数全部有显著性，Q36 = 21.88 20.05 (36-2) = 44。两种检验通过。上式变换为， Ln yt Lnyt - 12 = vt - vt 1 vt 12 + 0.2135 vt 13 Lnyt = Lnyt -1 +Ln yt - 12 Lnyt 13 + vt - vt 1 vt 12 + 0.2135 vt 13 (2.68)2.67式也是一个可以选用的模型。表2.3 北京市社会商品零售额yt月度数据单位：亿元人民币，1978:11989:12年：月yt年：月yt年：月yt年：月

15、yt年：月yt1978:011980:061982:111985:041987:091978:021980:071982:121985:051987:10505.21978:031980:081983:011985:06346.01987:111978:041980:091983:021985:071987:121978:051980:101983:031985:081988:011978:061980:111983:041985:091988:021978:071980:121983:05221.01985:101988:031978:081981:011983:061985:111988:

16、041978:091981:021983:071985:12451.01988:051978:101981:031983:081986:01412.01988:061978:111981:041983:091986:021988:071978:121981:05184.01983:10216.01986:03390.01988:081979:011981:061983:11235.01986:04387.01988:091979:021981:071983:12282.01986:051988:101979:031981:08172.01984:011986:061988:111979:041

17、981:091984:021986:071988:121979:051981:101984:031986:081989:011979:061981:111984:04247.01986:091989:021979:071981:121984:051986:101989:031979:081982:011984:061986:111989:041979:091982:02182.01984:071986:121989:051979:101982:031984:081987:011989:061979:11149.01982:041984:091987:021989:071979:121982:0

18、51984:101987:031989:081980:011982:061984:111987:041989:091980:021982:071984:121987:051989:101980:031982:081985:011987:061989:111980:041982:091985:021987:071989:121980:051982:101985:031987:08 案例2 香港季节GDPt数据的拟合季节时间序列模型，file：HongKong1980:12002:4年香港季度GDPt序列曲线见图2.27数据见表2.4，单位：港元。19801997年GDPt随时间呈指数增长。199

19、7年由于遭受东南亚金融危机的影响，经济开展处于停滞状态，19982002年底GDPt总量几乎没有增长。另一个特征是GDPt随时间呈递增型异方差。所以，用对数的季度GDPt数据LnGDPt，曲线见图2.48建立季节时间序列模型。 图2.47 GDPt8 LnGDPt通过LnGDPt的相关图和偏相关图图2.49可以看到LnGDPt是一个非平稳序列相关图衰减得很慢。图2.49 LnGDPt的相关图和偏相关图对LnGDPt进行一阶差分，得 DLnGDPt见图2.50。DLnGDPt的平稳性得到很大改良，但其季节因素影响还很大。从 DLnGDPt的相关图和偏相关图图2.51也可以明显地看到这个特征。假设

20、对LnGDPt直接进行一次季节差分四阶差分，得D4LnGDPt见图2.52。其波动性也很大。相关图和偏相关图见图2.53。D2LnGDPt显然是过度差分序列图2.54。 图2.50 DLnGDPt ,s.d. = 0.062 图2.51 DLnGDPt的相关图和偏相关图 2 D4LnGDPt,s.d. = 0.076 图2.53 D4LnGDPt的相关图和偏相关图 4 D2LnGDPt ,s.d. = 0.062 在DLnGDPt的根底上进行一阶季节差分，或在D4LnGDPt根底上进行一阶非季节差分，得 D4DLnGDPt图2.55。其相关图和偏相关图见图2.56。D4DLnGDPt中已经根本

21、消除了季节变化因素。在D4DLnGDPt的根底上建立时间序列模型。 图2.55 D4DLnGDPt,s.d. = 0.029 图2.56 D4DLnGDPt的相关和偏相关图通过对D4DLnGDPt的相关和偏相关图分析，应该建立(2, 1, 2) (1, 1, 1)4 模型。EViews估计命令是DLOG(GDP,1,4) C AR(1) AR(2) SAR(4) MA(1) MA(2) SMA(4)用1980:12002:3的数据得估计结果如下EViews输出结果如图：D4DLnGDPt = - 0.0023 + ut 1980:12002:3 ()(1-1.20 L+0.66 L2) (1

22、- 3 L4) ut = (1 - 1.16 L+ 0.97 L2) (1 - 0.95 L4) vt () (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.57 EViews估计结果()误差项的相关和偏相关图注意：1不要把自回归系数估计值的符号写错。不要把均值- 0.0023项表达错。EViews仍然是对(D4DLnGDPt + 0.0023)建立(2, 1, 2) (1, 1, 1)4 阶季节时间序列模型，而不是对D4

23、DLnGDPt建立季节时间序列模型。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.58 D4DLnGDPt的实际与预测序列 图2.59 GDPt的实际与预测序列对2002年第4季度GDPt进行样本外1期预测，结果如下：预测误差是 = 0.006表2.4 香港季度GDPt数据1980:120

24、02:4，单位：港元年：月GDPt /1011年：月GDPt /1011年：月GDPt /1011年：月GDPt /10111980:11985:41991:301997:21980:21986:11991:41997:31980:31986:21992:11997:41980:41986:31992:21998:11981:101986:41992:31998:21981:21987:11992:41998:31981:31987:21993:11998:41981:41987:31993:21999:11982:11987:41993:31999:21982:21988:11993:419

25、99:31982:31988:21994:101999:41982:41988:31994:22000:11983:11988:41994:32000:201983:21989:11994:42000:3001983:31989:21995:12000:41983:41989:31995:22001:11984:11989:41995:32001:201984:201990:11995:402001:31984:301990:21996:12001:41984:4001990:301996:22002:11985:11990:41996:32002:21985:21991:11996:42002:31985:31991:21997:12002:40