《传统时间序列分析》由会员分享,可在线阅读,更多相关《传统时间序列分析(24页珍藏版)》请在金锄头文库上搜索。
1、第九章第九章 传统时间序列分析传统时间序列分析时间序列的变动主要是由长期趋势、循环波动、季节变动及不规则变动而形成的,其中前三种变动有一个共同的特点,就是依一定的规则而变化,不规则变动则在综合中可以消除。基于这种认识,本章主要是介绍设法消除不规则变动,拟合确定型趋势,因而形成了一系列确定型时间序列分析方法。实验一实验一 季节模型季节模型实验目的:实验目的:掌握季节调整的方法。实验内容:实验内容:对时间序列进行季节调整。知识准备:知识准备:经济时间序列的变化受许多因素的影响,概括地讲,可以将影响时间序列变化的因素分为四种,即长期趋势(T,随着时间的变化,按照某种规律稳步地增长、下降或保持在某一水
2、平上) 、季节变动因素(S,在一个年度内依一定周期规则性变化) 、周期变动因素(C,以若干年为周期的波动变化)和不规则变动因素(I,许多不可控的偶然因素共同作用的结果) 。传统时间序列分析应是设法消除不规则变动,指拟合确定性趋势,因而形成了长期趋势分析、季节变动分析和循环波动测定等一系列确定型时间序列分析方法。季节变动是一种较为普遍的现象,其按照一定的周期循环进行,而且每个周期变化强度大体一致。研究季节变动的目的在于了解季节变动的规律,并进行季节预测。分析季节变动的方法有很多,其中常用的方法有两类:一是不考虑长期趋势的影响;二是考虑长期趋势的影响,运用时间序列模型分解的方法来计算季节指数。谓季
3、节调整,就是将某一统计指标的时间序列中的季节性因素和偶然性因素剔除,从而使经过季节调整的时间序列能够较为准确地反映出社会经济运行基本态势。本章主要介绍 X11 方法、Census X12 方法和移动平均比率法等季节调整方法。一、X11 方法X11 的全称是“X11”变量的第二类调查统计方法季节调整方案,通常简称为 X11 方案。其基本思想是利用一系列处理技术将不可比因素如季节、节假日、各月(季)的星期数量等分离,大大提高数据的可比性,以便于对系统作出正确的分析和客观的评价;同时,通过分离,获得关于系统动态结构和规律的大量信息。X11 季节调整方法包括乘法模型和加法模型。乘法模型将时间序列分解为
4、上述四个因素变动的乘积;加法模型则将序列分解为上述四个因素的和无约束样本回归模型。乘法模型: tttttISCTY(1)加法模型: tttttISCTY(2)对于一个时间序列,采用哪种模型分析取决于各成分之间的关系。一般来说,若 4 种成分是相互独立的应用加法模型;若相互有关联则用乘法模型。相对而言,乘法模型应用比较广泛。在乘法模型中,时间序列值和长期趋势用绝对数表示,季节变动、周期变动和不规则变动用相对数(百分数)表示。X11 是 X12 方法的核心,其目标旨在将经济时序 Yt 分解为趋势循环成分TCt 、季节成分 S t 和不规则分成 I t ,其中月度数据还需要分解出移动节日、交易日数等
5、日历效应分量 Dt 。X11 可以估计趋势和季节因素,可以用自动过滤选择模式选择季节过滤和趋势过滤(对序列结尾有专门的过滤) ,还能够调整异常值。设为原时间序列,现以乘法模型为例,说明 X11 季节调整方法的核心算tY法,共分为三个阶段:第一阶段,初步估计(1)序列趋势循环的粗略估计,即对原序列作中心化移动平均12/ )21.21(6556ttttttYYYYYTC(3)4)5 . 05 . 0(2112ttttttYYYYYTC(4)(3)式适用于月度数据, (4)式适用于季度数据。该移动平均能够保留线性趋势,消除季节性并减小不规则因素的方差。(2)趋势因素过滤,计算得到仅包含季节和不规则因
6、素变化新序列 SItttTCYSI/)1((5)(3)对序列进行 3*3 项移动平均,得到初步的季节因素估计序列 S)1( tSI)(9/ )232()1( 3*3)1( 24)1( 12)1()1( 12)1( 24)1(ttttttt SIMSISISISISIS(6)对每个月的观测值分别应用 3*3 项移动平均,从而保留线性趋势。(4)先将季节因素估计序列除以 2*12 的移动平均值,得到初步的无偏季节因子,然后从中消除季节因子从而得到不规则因子 )1( 12*2)1()1( 6)1( 5)1( 5)1( 6)1()1(24/ )2.2(tttttttt SMSSSSSSS(7)通过对季
7、节因子进行标准化,使得因子之和在每个连续的 12 个月内都近似为 0。(5)将原序列除以季节因子,从而得到初步季节调整后的序列)1()1(/tttSYTCI(8)这是初次估计的季节调整后序列,所包含的季节性因素已经很少。接下来将基于序列继续进行移动平均。)1( tTCI第二阶段 进一步估计和过滤(1)利用长度为 2H+1Henderson 过滤来计算趋势(循环)成分 HHjjtH jTCIhTCt)1()12()2((9)其中 H 值由信噪比标注 I/C(C 代表趋势)确定,为各平滑数据的权数。jhX-12 选择 H 的标准:当 I/C1.0 时,用 9(即 H=4)次过滤;当 1.0=I/C
8、3.5,用 13(即 H=6)次过滤;当 3.5=I/C 时,用 23(即 H=12)次过滤。(2)计算最后的季节-不规则因素的比率)2()2(/tttTCYSI(10)(3)通过 3*5 项移动平均计算出最后的季节因素,再对其进行收敛得到无偏季节因素(11) )2( 5*3)2( 36)2( 24)2( 12)2()2( 12)2( 24)2( 36)2(15/ )23332(ttIttttttt SIMSISSISISISISIS对每月(季)的观测值分别应用 3*5 项移动平均,从而保留线性趋势。(4)计算最终的季节因子(12)24/ )2.2()2( 6)2( 5)2( 5)2( 6)2
9、()2( ttttttSSSSSS上式对季节因子进行标准化,使得因子之和在每一个连续的 12 月(4 季度)内都近似为零。(5)第二步季节调整后的序列(13)2()2(/tttSYTCI第三阶段 最终估计与过滤(1)利用 Henderson 移动平均公式计算最终的趋势循环因素(14) HHjjtH jTCIhTCt)2()12()3(上式表示 2H+1 项 Henderson 移动平均,其移动平均的项数由数据决定,可能会第二阶段步骤一有所不同。(2)计算最终的不规则因素(15)3()2()3(/tttTCTCTI从而得到分解序列的乘法模型:)3()2()3( ttttISTCY至此,季节调整完
10、毕。季节调整加法模型等其他季节调整模型的步骤与乘法模型类似。由此可以看出,X11 是基于移动平均的季节调整方法,用户可以根据各种季节调理的目的不同自主地选择计算公式,从而适应序列的特征。系统也会自动地根据随机因素的大小确定移动平均的长度。在 X11 中,趋势过滤包括212(或 24)初步的趋势估计和最终趋势估计的 Handerson 过滤两个环节,季节(自动)过滤包括 33 初步的季节估计和 33、35 或 39 的最终季节过滤两个环节。多数情况下可选择 35,季节模式迅速变化时用 33,季节模式不是正在变化或不规则因素影响很大时用 39。趋势估计或季节估计都采用平滑方法。X11 方法是通过几
11、次迭代来进行分解的,每一次对组成因子的估算都进一步准确化。应当注意,X11 只能作用于含季节因素的序列,季节调整的观测值的个数是有限制的,最少需要 4 个整年的数据,最多能调整 20 年的月度数据及30 年的季度数据。二、Census X12 季节调整方法Census X12 是美国人口普查局关于季节调整的最新研究成果,X11 是 X12方法的核心。 其主要的变化在于可以对数据作更加丰富的预处理,检测和修正不同类型的离群值,估计并消除日历因素的影响,对季节调整的效果进行更严格的诊断检验。其是一种基于给定序列进行成分提取的季节调整方法,是一种基于移动平均的非参数方法提取成分。1、X12 季节调整
12、的模型选择X12 保留了基本的加法和乘法模型,而且还可以计算另一种乘法模型,即把被调整序列取对数后进行加法分解再指数化,这被称为对数加法模型。它适合于对模型的趋势估计进行偏差修正,同时还可对对数正态分布识别的极端值进行修正。伪加法模型是由英国国家统计办公室开发。时序季节调整方法有如下四种:加法模型: ttttISTCY(16)乘法模型: ttttISTCY(17)利用对数变换可将乘法模型转换成加法模型,即对数加法模型:ttttISTCYlnlnlnln(18)伪加法模型: ) 1(ttttISTCY(19)其中,表示趋势循因素;表示季节因素;表示不规则因素。乘法tTCtStI模型应用最为广泛。
13、伪加法模型主要用来对在每年相同月份取值很小、接近于零的非负时间序列进行季节调整。X12 季节调整方法的基本过程和步骤与 X-11 方法的原理相同,在此就不再重述。但其在 X-11 方法的基础上,增加了交易日的调整。2、交易日影响月度(季度)时间序列中循环出现的星期构成效应称为交易日影响。Yong(1965)讨论了流量交易日的影响,Cleveland and Grupe(1983)讨论了存量贸易日的影响。与季节效应类似,交易日效应也使得序列的跨月比较或者两个序列之间运动规律的比较显得困难。因此,在执行季节调整时如果估计的交易日效应应具有统计显著性,通常要将它们从序列中调整出去。为了从月度流量序列
14、中得到交易日效应模型,设一周中第 i 天的效应记为,表示月份 t 中星期 i 的个数,则该月贸易日的累积效应可表示)7,.,2 , 1( iiitd为:(20) 71iitid记为平均日度效应,该月的长度,则上式tN 71, 7/iittidN可以表示为: 6177171)()(ititit iiitititiddNdNd累积的月度效应可以分解为月份长度效应和日度对立变量的净效应,因此式(20)的季节成分和水平成分归于日历月均值表示,其中:tN(21) 411241*kktNNt至于这些成分如何从模型中消除,取决于用以得到不规则成分的季节分解的类型。对于乘法分解模型,在方程(20)两端同时除以
15、以消除对交易日tN效应的去季节化和退势。令,得到1)/(aii(22)*7161*7 *) 1( )(tiitiittit i tt NdNddNN 这就是 Yong(1965)没有推导而直接给出的交易日效应模型。设 ID 表示包括贸易日影响的不规则因素的初步估计,则其估计可通过最小二乘法拟合下面的模型得到回归系数的估计值(其中):)6,.,2 , 1( ii 617 ii(23)ttiti idd )(761同样地,对于加法模型对应于式(23) ,是由式(20)两端同时减去,tN从而得到(24)tt iitittddNN )()(771* 0其中:,;然后与乘法模型类似,建立回归方ii,061 i程t iititttdNNID 71* 0)((25)利用普通最小二乘法估计上式中的参数,从而求得的估计结果)7,.,2 , 1(