第七章季节性时间序列分析方法由于季节性时间序列在经济生活中大量存在,故将季节时间序列从非平稳序列中抽出来,单独作为一章加以研究,具有较强的现实意义。
本章共分四节:简单随机时间序列模型、乘积季节模型、季节型时间序列模型的建立、季节调整方法X-11程序。
本章的学习重点是季节模型的一般形式和建模。
§1 简单随机时序模型在许多实际问题中,经济时间序列的变化包含很多明显的周期性规律。
比如:建筑施工在冬季的月份当中将减少,旅游人数将在夏季达到高峰,等等,这种规律是由于季节性(seasonality)变化或周期性变化所引起的。
对于这各时间数列我们可以说,变量同它上一年同一月(季度,周等)的值的关系可能比它同前一月的值的相关更密切。
一、季节性时间序列1.含义:在一个序列中,若经过S个时间间隔后呈现出相似性,我们说该序列具有以S为周期的周期性特性。
具有周期特性的序列就称为季节性时间序列,这里S为周期长度。
注:①在经济领域中,季节性的数据几乎无处不在,在许多场合,我们往往可以从直观的背景及物理变化规律得知季节性的周期,如季度数据(周期为4)、月度数据(周期为12)、周数据(周期为7);②有的时间序列也可能包含长度不同的若干种周期,如客运量数据(S=12,S=7)2.处理办法:(1)建立组合模型;(1)将原序列分解成S个子序列(Buys-Ballot 1847)对于这样每一个子序列都可以给它拟合ARIMA模型,同时认为各个序列之间是相互独立的。
但是这种做法不可取,原因有二:(1)S 个子序列事实上并不相互独立,硬性划分这样的子序列不能反映序列{}t x 的总体特征;(2)子序列的划分要求原序列的样本足够大。
启发意义:如果把每一时刻的观察值与上年同期相应的观察值相减,是否能将原序列的周期性变化消除?(或实现平稳化),在经济上,就是考查与前期相比的净增值,用数学语言来描述就是定义季节差分算子。
定义:季节差分可以表示为S t t t S t S t X X X B X W --=-=∇=)1(。
二、 随机季节模型1.含义:随机季节模型,是对季节性随机序列中不同周期的同一周期点之间的相关关系的一种拟合。
AR (1):t t S t S t t e W B e W W =-⇔+=-)1(11ϕϕ,可以还原为:t t S S e X B =∇-)1(1ϕ。
MA (1):t S t S t t t e B W e e W )1(11θθ-=⇔-=-,可以还原为:t S t S e B X )1(1θ-=∇。
2.形式:广而言之,季节型模型的ARMA 表达形式为t S t S e B V W B U )()(= (1)这里,⎪⎩⎪⎨⎧----=----=∇=qSq S S S pS P S S S td S t B V B V B V B V B U B U B U B U X W ΛΛ2212211)(1)()(平稳。
注:(1)残差t e 的内容;(2)残差t e 的性质。
§2 乘积季节模型一、 乘积季节模型的一般形式由于t e 不独立,不妨设),,(~m d n ARIMA e t ,则有t t d a B e B )()(Θ=∇φ (2)式中,t a 为白噪声;n n B B B B ϕϕϕφ----=Λ22111)(;m m B B B B θθθ----=ΘΛ22111)(。
在(1)式两端同乘d B ∇)(φ,可得:t S t d S t DS d S t d S a B B V e B B V X B U B W B U B )()()()()()()()(Θ=∇=∇∇=∇φφφ (3)注:(1)这里t DSS X B U ∇)(表示不同周期的同一周期点上的相关关系;t d X B ∇)(φ则表示同一周期内不同周期点上的相关关系。
二者的结合就能同时刻划两个因素的作用,仿佛是显像管中的电子扫描。
(2)从结构上看,它是季节模型与ARIMA 模型的结合形式,称之为乘积季节模型,阶数用S q D p m d n ),,(),,(⨯来表示。
(3)将乘积季节模型展开便会得到一般的ARIMA 模型。
例如:t S t a B V B X B )1)(1()1(11--=-θ,可以展开为t S S t a B V B V B X B )1()1(11111++--=-θθ,此时也有)1,1,0(~+S ARIMA X t ,并且其中有许多系数为0。
但其参数并不独立。
所以尽管模型的阶数可能很高,然而真正独立的参数不多,我们称这类模型为疏系数模型(带有一定约束条件的疏系数模型)。
二、 常用的两个模型1.t t a B B X B B )1)(1()1)(1(1212112θθ--=-- 类型为:S )1,1,0()1,1,0(⨯ (4) 2.t t a B B X B )1)(1()1(1212112θθ--=- 类型为:S )1,1,0()1,0,0(⨯ (5)三、 乘积季节模型与ARIMA 模型的关系我们可以将乘积季节模型t S t d S t DS d S t d S a B B V e B B V X B U B W B U B )()()()()()()()(Θ=∇=∇∇=∇φφφ (3)展成ARIMA 模型形式。
例如,t S t a B V B y B )1)(1()1(11--=-θ是)1,0,0()1,1,0(⨯季节模型,将式子的右边展成:t S j jj t S S t a B a BV B V B y B )1()1()1(11*11111∑+=+-=+--=-θθθ (6)这是一个)1,1,0(+S 阶ARIMA 模型,但是其参数不是独立的,有下面的约束关系11*11**1*21*1,,0,V V S S S θθθθθθθ-======+-Λ (7)尽管模型的阶数很高,然而真正独立的参数并不多,有许多参数取值为零§3 季节性时间序列模型的建立季节性时间序列模型的建立也包含这样几个过程:模型的识别、模型的定阶、参数估计、诊断检验等。
基本上采用的是BOX-JENKINS 方法,也就是立足于考察数据序列的样本自相关、偏自相关函数。
如果样本自相关、偏自相关函数既不截也不拖尾,而且也不呈线性衰减趋势,相反地,在相应于周期S 的整数倍点上,自相关(或偏自相关)函数出现绝对值相当大的峰值并呈现振荡变化,我们就可以判明原数据序列适合于乘积季节模型。
一、 季节性MA 模型的自相关函数{}t X 是一个季节性时间序列,如果S t MA X )1(~,则t S S t e B X )1(θ-= (6)t e 不平稳,设)1(~MA e t ,则t t a B e )1(1θ-= (7)我们就能得到一个乘积季节模型t S S t a B X )1)(1(1θθ--= (8)1111----+--=S t S S t S t t t a a a a X θθθθ (9)当S=12时,有)13(~131********MA a a a a X t t t t t ---+--=θθθθ (10)可以计算出: 因此有:注:(1)1ρ为t t a B e )1(1θ-=的一阶自相关系数,12ρ为t S S t e B X )1(θ-=的一阶自相关系数; (2)1θ与12θ比较容易求解; (3)可以推广到更一般的形式。
二、 季节性AR 模型的偏自相关函数 {}t X 是一个季节性时间序列,如果S t AR X )1(~,则t t S S e X B =-)1(ϕ (11)t e 不平稳,设)1(~AR e t ,则t t a e B =-)1(1ϕ (12)我们就能得到一个乘积季节模型t t S S a X B B =--)1)(1(1ϕϕ (13) t t S S S S a X B B B =+--+)1(111ϕϕϕϕ (14)当S=12时,有)13(~131********AR a X X X X t t t t t =+-----ϕϕϕϕ (15)可以根据YULE-WORK 方程求出偏自相关函数。
注:(1)根据它在周期点上的偏自相关函数的截尾性和拖尾性识别模型的类型和定阶; (2)可以推广到更一般的形式。
三、 季节性时间序列模型的建模方法利用B-J 建模方法:判别周期性,即S 的取值;根据SACF 和SPACF 提供的信息识别模型类型和阶数,最后进行估计和诊断检验。
具体做法:第一步:对时间序列{}t X 进行普通差分∆和季节差分S ∆,以得到平稳的序列{}t W ,t D S d t X W ∆∆=; 第二步:计算差分后序列的SACF 和SPACF ,选择一个暂定的模型;第三步:由SACF 和SPACF 函数的值,利用矩估计法得到的值作为初始值,对模型参数作最小二乘估计;第四步:模型的诊断与检验。
注:(1)关于差分阶数d 和季节差分阶数D 的选取可采用试探的方法1;也可使用差分后序列均方差的大小挑选;(2)季节差分算子的阶数不宜过高。
四、 应用实例【例6-1】试用1987年到1996年甲地某商品各月销售量资料为例建立季节性时间序列模型2。
建模型过程: 1.时间序列图明显存在着季节性变化,并且以12为周期。
2.SACF 和SPACF 函数图1 详见备课笔记。
2 资料来源王振龙:《时间序列分析》,中国统计出版社,P189。
再次证明,时间序列存在着以S=12为周期的季节性变动。
3.进行差分变换需要进行一阶普通差和以12为周期的季节差分,得到t t X B Y )1(-= (17) t t t t W X B B Y B X =--=-=)1)(1()1(1212 (16)计算其自相关系数。
一阶普通差分图一阶普通差分和一阶季节差分序列图 4.模型的识别与定阶 5.参数估计 6.诊断检验7.模型应用预测结果【例6-2】表显示了我国1990年1月至1997年12月工业总产值的月度资料(1990年不变价格),记作IP t,共有96个观测值,对序列IP t建立ARMA模型3,在建模过程中将1997年12个月的观测值留出作为评价预测精度的参照对象。
1990年1月至1997年12月我国工业总产值单位:亿元1.时间序列图表明数据或者序列是非平稳的。
2.进行相应的差分变换为消除趋势同时减小序列的波动,对原序列做一阶自然对数并逐期差分,即是差分运算与对数运算的结合。
由时间序列图可以看到,序列的趋势已经基本消除,但可能存在着季节性变化,这一点可以从序列的自相关图看出。
由图形可以看出,在12的整数倍上,样本的偏自相关系数显着不为零,因此需要做季节差分处理。
此时差分后序列的自相关图为可以对序列进行零均值的检验,详见易丹辉:《数据分析与EVIEWS应用》,P128。