经济时间序列季节调整分解与平滑PPT课件

1、1 经济指标的月度或季度时间序列包含经济指标的月度或季度时间序列包含4种变动要素：种变动要素： 长期趋势要素长期趋势要素T 循环要素循环要素C 季节变动要素季节变动要素S 不规则要素不规则要素I第1页/共40页20.760.860.961.061.161981198319851987198919911993199519970.890.951.001.061.11198119831985198719891991199319951997 第2页/共40页3 季节性变动的发生，不仅是由于气候的直接影响，而且社会制度及风季节性变动的发生，不仅是由于气候的直接影响，而且社会制度及风俗习惯也会引起季节变动

2、。经济统计中的月度和季度数据或大或小都含有季俗习惯也会引起季节变动。经济统计中的月度和季度数据或大或小都含有季节变动因素，节变动因素，以月份或季度作为时间观测单位的经济时间序列通常具有一年以月份或季度作为时间观测单位的经济时间序列通常具有一年一度的周期性变化，这种周期变化是由于季节因素的影响造成的，在经济分一度的周期性变化，这种周期变化是由于季节因素的影响造成的，在经济分析中称为季节性波动。经济时间序列的季节性波动是非常显著的，它往往遮析中称为季节性波动。经济时间序列的季节性波动是非常显著的，它往往遮盖或混淆经济发展中其他客观变化规律，以致给经济增长速度和宏观经济形盖或混淆经济发展中其他客观变

3、化规律，以致给经济增长速度和宏观经济形势的分析造成困难和麻烦。因此，在进行经济增长分析时，必须去掉季节波势的分析造成困难和麻烦。因此，在进行经济增长分析时，必须去掉季节波动的影响，将季节要素从原序列中剔除，这就是所谓的动的影响，将季节要素从原序列中剔除，这就是所谓的“季节调整季节调整” (Seasonal Adjustment)。 第3页/共40页4 X-11方法是基于移动平均法的季节调整方法。它的特征在于除了能适方法是基于移动平均法的季节调整方法。它的特征在于除了能适应各种经济指标的性质，根据各种季节调整的目的，选择计算方式外，在不作应各种经济指标的性质，根据各种季节调整的目的，选择计算方式

4、外，在不作选择的情况下，也能根据事先编入的统计基准，按数据的特征自动选择计算方选择的情况下，也能根据事先编入的统计基准，按数据的特征自动选择计算方式。在计算过程中可根据数据中的随机因素大小，采用不同长度的移动平均，式。在计算过程中可根据数据中的随机因素大小，采用不同长度的移动平均，随机因素越大，移动平均长度越大。随机因素越大，移动平均长度越大。X-11方法是通过几次迭代来进行分解的，方法是通过几次迭代来进行分解的，每一次对组成因子的估算都进一步精化。每一次对组成因子的估算都进一步精化。第4页/共40页5 美国商务部国势普查局的美国商务部国势普查局的X12季节调整程序是在季节调整程序是在X11方

5、方法的基础上发展而来的，包括法的基础上发展而来的，包括X11季节调整方法的全部功季节调整方法的全部功能，并对能，并对X11方法进行了以下方法进行了以下3方面的重要改进：方面的重要改进： (1) 扩展了贸易日和节假日影响的调节功能，增加了季扩展了贸易日和节假日影响的调节功能，增加了季节、趋势循环和不规则要素分解模型的选择功能；节、趋势循环和不规则要素分解模型的选择功能； (2) 新的季节调整结果稳定性诊断功能；新的季节调整结果稳定性诊断功能； (3) 增加增加X12-ARIMA模型的建模和模型选择功能。模型的建模和模型选择功能。 第5页/共40页6 X12季节调整方法的核心算法是扩展的季节调整方

6、法的核心算法是扩展的X11季节调整程序。共包括季节调整程序。共包括4种季节调种季节调整的分解形式：乘法、加法、伪加法和对数加法模型。注意采用乘法、伪加法和整的分解形式：乘法、加法、伪加法和对数加法模型。注意采用乘法、伪加法和对数加法模型进行季节调整时，时间序列中不允许有零和负数。对数加法模型进行季节调整时，时间序列中不允许有零和负数。 加法模型加法模型 (2.2.1) 乘法模型：乘法模型： (2.2.2) 对数加法模型：对数加法模型： (2.2.3) 伪加法模型：伪加法模型： (2.2.4)ttttISTCYttttISTCYttttISTCYlnlnlnln) 1(ttttISTCY第6页/

7、共40页7 第7页/共40页8 第8页/共40页9 TRAMO(Time Series Regression with ARIMA Noise, Missing Observation, and Outliers)用来估计和预测用来估计和预测具有缺失观测值具有缺失观测值、非平稳、非平稳ARIMA误误差及外部影响的回归模型。它能够对原序列进行插值，识别和修正几种不同类型差及外部影响的回归模型。它能够对原序列进行插值，识别和修正几种不同类型的异常值，并对工作日变化及复活节等特殊回归因素及假定为的异常值，并对工作日变化及复活节等特殊回归因素及假定为ARIMA过程的误过程的误差项的参数进行估计。差项的

8、参数进行估计。 SEATS(Signal Extraction in ARIMA Time Series)是基于是基于ARIMA模型来对模型来对时间序列中不可观测成分进行估计。时间序列中不可观测成分进行估计。 这两个程序往往联合起来使用，先用这两个程序往往联合起来使用，先用TRAMO对数据进行预处理，然后用对数据进行预处理，然后用SEATS将时间序列分解为趋势要素、循环要素、季节要素及不规则要素将时间序列分解为趋势要素、循环要素、季节要素及不规则要素4个部分。个部分。第9页/共40页10 本节主要介绍利用本节主要介绍利用EViews软件对一个月度或季度时间序列进行季节调整的软件对一个月度或季度

9、时间序列进行季节调整的操作方法。在操作方法。在EViews工作环境中，打开一个月度或季度时间序列的工作文件，工作环境中，打开一个月度或季度时间序列的工作文件，双击需进行数据处理的序列名，进入这个序列对象，在序列窗口的工具栏中单双击需进行数据处理的序列名，进入这个序列对象，在序列窗口的工具栏中单击击Proc按钮将显示菜单：按钮将显示菜单：第10页/共40页11 X-11法是美国商务部标准的季节调整方法法是美国商务部标准的季节调整方法(乘法模型、加法模型乘法模型、加法模型)，乘法模型适，乘法模型适用于序列可被分解为季节调整后序列（趋势用于序列可被分解为季节调整后序列（趋势循环循环不规则要素项）与季

10、节项的乘积，项）与季节项的乘积，加法模型适用于序列可被分解为季节调整后序列与季节项的和。乘法模型只适用于加法模型适用于序列可被分解为季节调整后序列与季节项的和。乘法模型只适用于序列值都为正的情形。序列值都为正的情形。 第11页/共40页12 EViews是将美国国势调查局的是将美国国势调查局的X12季节调整程序直接安装到季节调整程序直接安装到EViews子目子目录中，建立了一个接口程序。录中，建立了一个接口程序。 EViews进行季节调整时将执行以下步骤：进行季节调整时将执行以下步骤： 1给出一个被调整序列的说明文件和数据文件；给出一个被调整序列的说明文件和数据文件； 2利用给定的信息执行利用

11、给定的信息执行X12程序；程序； 3返回一个输出文件，将调整后的结果存在返回一个输出文件，将调整后的结果存在EViews工作文件中。工作文件中。 X12的的EViews接口菜单只是一个简短的描述，接口菜单只是一个简短的描述，EViews还提供了一些菜单还提供了一些菜单不能实现的接口功能，更一般的命令接口程序。不能实现的接口功能，更一般的命令接口程序。 第12页/共40页13 调用调用X12季节调整过程，在序列窗口选择季节调整过程，在序列窗口选择Procs/Seasonal Adjustment / Census X12，打开一个对话框：，打开一个对话框： 第13页/共40页14 X-11法与移

12、动平均法的最大不同是：法与移动平均法的最大不同是：X-11法中季节因子年与年有可能法中季节因子年与年有可能不同，而在移动平均法中，季节因子被假设为是一样的。不同，而在移动平均法中，季节因子被假设为是一样的。 第14页/共40页15 Tramo(Time Series Regression with ARIMA Noise, Missing Observation, and Outliers)是对具有缺失观测值，是对具有缺失观测值，ARIMA误差、几种外部影响误差、几种外部影响的回归模型完成估计、预测和插值的程序。的回归模型完成估计、预测和插值的程序。 Seats(Signal Extracti

13、on in ARIMA Time Series)是基于是基于ARIMA模型的将可模型的将可观测时间序列分解为不可观测分量的程序。这两个程序是有观测时间序列分解为不可观测分量的程序。这两个程序是有Victor Gomez 和和Agustin Maravall 开发的。开发的。 当选择了当选择了Pross/Seasonal Adjustment/Tramo Seats 时，时，EViews执行外部执行外部程序，将数据输给外部程序，然后将结果返回程序，将数据输给外部程序，然后将结果返回EViews。 第15页/共40页16 本章第本章第2节介绍的季节调整方法可以对经济时间序列进行分解，但在季节节介绍

14、的季节调整方法可以对经济时间序列进行分解，但在季节调整方法中，趋势和循环要素视为一体不能分开。本节专门讨论如何将趋势和调整方法中，趋势和循环要素视为一体不能分开。本节专门讨论如何将趋势和循环要素进行分解的方法。测定长期趋势有多种方法，比较常用的方法有回归循环要素进行分解的方法。测定长期趋势有多种方法，比较常用的方法有回归分析方法、移动平均法、阶段平均法分析方法、移动平均法、阶段平均法(phase average，PA方法方法)、HP滤波方法滤波方法和频谱滤波方法（和频谱滤波方法（frequency (band-pass) filer， BP滤波）。本节主要介绍滤波）。本节主要介绍HP滤波方法和

15、滤波方法和BP滤波方法。滤波方法。 第16页/共40页17 在宏观经济学中，人们非常关心序列组成成分中的长期趋势，在宏观经济学中，人们非常关心序列组成成分中的长期趋势，Hodrick-Prescott滤波是被广泛使用的一种方法。该方法在滤波是被广泛使用的一种方法。该方法在Hodrick and Prescott(1980) 分析战后美国经济周期的论文中首次使用。我们简要介绍这种方法的原理。分析战后美国经济周期的论文中首次使用。我们简要介绍这种方法的原理。设设Yt是包含趋势成分和波动成分的经济时间序列，是包含趋势成分和波动成分的经济时间序列，YtT是其中含有的趋势成是其中含有的趋势成分，分， Y

16、tC是其中含有的波动成分。则是其中含有的波动成分。则 (2.3.1) 计算计算HP滤波就是从滤波就是从Yt中将中将YtT 分离出来分离出来 。ctTttYYYTt,2, 1第17页/共40页18 一般地，时间序列一般地，时间序列Yt中的不可观测部分趋势中的不可观测部分趋势YtT常被定义为下面最小化问题常被定义为下面最小化问题的解：的解： (2.3.2)其中：其中：c(L)是延迟算子多项式是延迟算子多项式 (2.3.3) 将式将式(2.3.3)代入式代入式(2.3.2)，则，则HP滤波的问题就是使下面损失函数最小，即滤波的问题就是使下面损失函数最小，即 (2.3.4) TtTtTttYLcYY1

17、22min LLLc111 TtTtTtTtTtTtTttYYYYYY121112min第18页/共40页19 最小化问题用最小化问题用c(L)YtT2 来调整趋势的变化，并随着来调整趋势的变化，并随着 的增大而增大。这里存在的增大而增大。这里存在一个权衡问题，要在趋势要素对实际序列的跟踪程度和趋势光滑度之间作一个选择。一个权衡问题，要在趋势要素对实际序列的跟踪程度和趋势光滑度之间作一个选择。 = 0 时，满足最小化问题的趋势等于序列时，满足最小化问题的趋势等于序列Yt； 增加时，估计趋势中的变化总数增加时，估计趋势中的变化总数相对于序列中的变化减少，即相对于序列中的变化减少，即 越大，估计趋

18、势越光滑；越大，估计趋势越光滑； 趋于无穷大时，估计趋趋于无穷大时，估计趋势将接近线性函数。一般经验地，势将接近线性函数。一般经验地， 的取值如下：的取值如下： 月度数据，季度数据，年度数据144001600100第19页/共40页20 使用使用Hodrick-Prescott滤波来平滑序列，选择滤波来平滑序列，选择Procs/ Hodrick Prescott Filter出现出现下面的下面的HP滤波对话框：滤波对话框： 首先对平滑后的序列给一个名字，首先对平滑后的序列给一个名字，EViews将默认一个名字，也可填入一个新的将默认一个名字，也可填入一个新的名字。然后给定平滑参数的值，年度数据

19、取名字。然后给定平滑参数的值，年度数据取100，季度和月度数据分别取，季度和月度数据分别取1600和和14400。不允许填入非整数的数据。点击不允许填入非整数的数据。点击OK后，后，EViews与原序列一起显示处理后的序列。注意与原序列一起显示处理后的序列。注意只有包括在当前工作文件样本区间内的数据才被处理，平滑后序列区间外的数据都为只有包括在当前工作文件样本区间内的数据才被处理，平滑后序列区间外的数据都为NA。 第20页/共40页21 利用利用HP滤波方法求中国社会消费品零售总额月度时间序列和中国滤波方法求中国社会消费品零售总额月度时间序列和中国GDP季度时季度时间序列的趋势项。间序列的趋势

20、项。第21页/共40页22 设设Yt为我国的季度为我国的季度GDP指标，利用季节调整方法将指标，利用季节调整方法将GDP中的季节因素和不规中的季节因素和不规则因素去掉，得到则因素去掉，得到GDP_TC序列。本例的潜在产出序列。本例的潜在产出Y*，即趋势利用，即趋势利用HP滤波计算出来滤波计算出来的的YtT来代替，来代替，GDP的循环要素的循环要素Yt序列由式序列由式(2.3.6)计算：计算：(2.3.6)TttctYYYTt,2, 1 第22页/共40页23 图图2.7显示的显示的GDP的循环要素的循环要素YtC序列实际上就是围绕趋势线上下的波动，序列实际上就是围绕趋势线上下的波动，称为称为G

21、DP缺口序列。它是一个绝对量的产出缺口。也可以用相对量表示产出缺口，缺口序列。它是一个绝对量的产出缺口。也可以用相对量表示产出缺口，本例用本例用Gapt来表示相对产出缺口，可由下式计算得到：来表示相对产出缺口，可由下式计算得到： (2.3.7) TtTtttYYYGap100第23页/共40页24 20世纪以来，利用统计方法特别是时间序列分析方法研究经济时间序列和世纪以来，利用统计方法特别是时间序列分析方法研究经济时间序列和经济周期的变动特征得到越来越广泛的应用。自时间序列分析产生以来，一直经济周期的变动特征得到越来越广泛的应用。自时间序列分析产生以来，一直存在两种观察、分析和解释时间序列的方

22、法。第一种是直接分析数据随时间变存在两种观察、分析和解释时间序列的方法。第一种是直接分析数据随时间变化的结构特征，即所谓时域（化的结构特征，即所谓时域（time domain）分析法，使用的工具是自相关）分析法，使用的工具是自相关（或自协方差）函数和差分方程；另一种方法是把时间序列看成不同谐波的叠（或自协方差）函数和差分方程；另一种方法是把时间序列看成不同谐波的叠加，研究时间序列在频率域（加，研究时间序列在频率域（frequency domain）里的结构特征，由于这种分）里的结构特征，由于这种分析主要是用功率谱的概念进行讨论，所以通常称为谱分析。析主要是用功率谱的概念进行讨论，所以通常称为谱

23、分析。第24页/共40页25 谱分析的基本思想是：把时间序列看作是互不相关的周期（频率）分谱分析的基本思想是：把时间序列看作是互不相关的周期（频率）分量的叠加，通过研究和比较各分量的周期变化，以充分揭示时间序列的频量的叠加，通过研究和比较各分量的周期变化，以充分揭示时间序列的频域结构，掌握其主要波动特征。因此，在研究时间序列的周期波动方面，域结构，掌握其主要波动特征。因此，在研究时间序列的周期波动方面，它具有时域方法所无法企及的优势。它具有时域方法所无法企及的优势。 第25页/共40页26 在在EViews中，可以使用中，可以使用 Band-Pass 滤波对经济时间序列进行趋势循环滤波对经济时

24、间序列进行趋势循环分解。在序列对象的菜单中选择分解。在序列对象的菜单中选择 Proc/ Frequency Filter，显示如下所示的，显示如下所示的对话框。对话框。第26页/共40页27 为了使用为了使用Band-Pass滤波，首先要选择一种滤波类型。共有滤波，首先要选择一种滤波类型。共有3种种类型：类型： （1） BK固定长度对称滤波（固定长度对称滤波（Fixed length symmetric (Baxter-King，BK)）；）； （2）CF固定长度对称滤波（固定长度对称滤波（Fixed length symmetric (Christiano-Fitzgerald，CF)）；）

25、； （ 3 ） 全 样 本 长 度 非 对 称 滤 波 （） 全 样 本 长 度 非 对 称 滤 波 （ F u l l s a m p l e asymmetric(Christiano-Fitzgerald)）。）。 EViews默认的是默认的是BK 固定长度对称滤波。如果使用固定长度对固定长度对称滤波。如果使用固定长度对称滤波，还必须指定先行称滤波，还必须指定先行/滞后（滞后（Lead/lag）项数）项数n。第27页/共40页28 用户必须选择循环周期（用户必须选择循环周期（Cycle periods）的区间以计算）的区间以计算Band-Pass滤滤波的频率响应函数的权重序列。这个区间由

26、一对数据波的频率响应函数的权重序列。这个区间由一对数据（PL，PU）描述，描述，PL、PU 由由Band-Pass滤波要保留的循环波动成分所对应的周期来确定。月度数据滤波要保留的循环波动成分所对应的周期来确定。月度数据填月数；季度数据填季度的个数。填月数；季度数据填季度的个数。EViews将根据数据类型填入了默认数值。将根据数据类型填入了默认数值。例如，例例如，例2.6认为中国社会消费品零售总额的增长周期大约在认为中国社会消费品零售总额的增长周期大约在1年半（年半（18个月）个月）到到5年（年（60个月），如果保留在这个区间内的循环要素，则区间的下界是个月），如果保留在这个区间内的循环要素，则

27、区间的下界是18，上界是上界是60。因此，设定。因此，设定PL=18，PU=60（相当于例相当于例2.6中的中的 p和和q）。）。第28页/共40页29 在在Band-Pass滤波的输出结果中，左侧的图描述了原序列、趋势序列和循滤波的输出结果中，左侧的图描述了原序列、趋势序列和循环序列。对于环序列。对于BK和和CF固定长度对称滤波而言，固定长度对称滤波而言，Eviews 画出频率响应函数画出频率响应函数w()，频率频率 的区间是的区间是0，0.5，右面的图描述了频率响应函数。但是，对于时变的，右面的图描述了频率响应函数。但是，对于时变的CF滤波，并没有画出频率响应函数，因为滤波的频率响应函数随

28、数据和观测值滤波，并没有画出频率响应函数，因为滤波的频率响应函数随数据和观测值个数变化。个数变化。 用户需要输入希望保存的结果（循环成分、趋势成分）对象的名字。循用户需要输入希望保存的结果（循环成分、趋势成分）对象的名字。循环序列（环序列（Cycle series）是包含循环要素的序列对象；趋势序列）是包含循环要素的序列对象；趋势序列(Non-cyclical series)是实际值和循环序列的差。用户还能得到在滤波中所用的是实际值和循环序列的差。用户还能得到在滤波中所用的Band-Pass滤滤波频率响应函数的权序列，它将存储在矩阵对象中。波频率响应函数的权序列，它将存储在矩阵对象中。 第29

29、页/共40页30 中国社会消费品零售总额月度时间序列（中国社会消费品零售总额月度时间序列（SL）的取值范围从）的取值范围从1980年年1月月至至2004年年8月（附录月（附录E表表E.5）。由于带通（）。由于带通（BP）滤波的两端各欠）滤波的两端各欠n项，为了项，为了近期的分解结果没有缺失值，本例利用近期的分解结果没有缺失值，本例利用ARIMA模型将序列外推到模型将序列外推到2006年年2月。月。然后对然后对SL进行季节调整去掉季节和不规则要素，得到只包含趋势循环要素的进行季节调整去掉季节和不规则要素，得到只包含趋势循环要素的序列序列SL_TC。根据增长率周期波动分析，我国社会消费品零售总额的

30、增长率。根据增长率周期波动分析，我国社会消费品零售总额的增长率大约存在大约存在1.5年年5年之间的波动。年之间的波动。 取取 p = 18 ( p = 1/18)，q = 60 ( q = 1/60)，利用带通滤波方法希望得到利用带通滤波方法希望得到只保留只保留1.5年年5年周期成分的滤波序列。而取年周期成分的滤波序列。而取n =18的的BPn(p，q) 滤波中滤波中2年年3.5年周期成分的权重最大，可以近似地作为中国社会消费品零售总额年周期成分的权重最大，可以近似地作为中国社会消费品零售总额的循环要素序列的循环要素序列SL_C，同时可以从同时可以从SL_TC中去掉中去掉SL_C，得到趋势要素序列得到趋势要素序列SL_T。图图2.12是是BP滤波的频率响应函数。滤波的频率响应函数。 第30页/共40页31第31页/共40页32第32页/共40页33第33页/共40页34 指数平滑是可调整预测的简单方法。当只有少数观测值时这种方法指数平滑是可调整预测的简单方法。当只有少数观测值时这种方法是有效的。与使用固定系数的回归预测模型不同，指数平滑法的预测用是有效的。与使用固定系数的回归预测模型不同，指数平滑法的预测用过去的预测误差进行调整。过去的预