经济时间序列的季节调整、分解和平滑方法-2.ppt

第二章第二章 经济时间序列的经济时间序列的 季节调整、分解与平滑季节调整、分解与平滑 本章主要介绍经济时间序列的分解和平滑方 法。时间序列分解方法包括季节调整和趋势分解 ，指数平滑是目前比较常用的时间序列平滑方法 。 1 经济指标的月度或季度时间序列包含4种变动要素 ： 长期趋势要素T 循环要素C 季节变动要素S 不规则要素I 经济时间序列的分解经济时间序列的分解 2 图图1 1 我国工业总产值的时间序列我国工业总产值的时间序列 Y Y 图形图形 图图2 2 工业总产值的趋势工业总产值的趋势 循环要素循环要素 TC TC 图形图形 图图3 3 工业总产值的季节变动要素工业总产值的季节变动要素 S S 图形图形 图图4 4 工业总产值的不规则要素工业总产值的不规则要素 I I 图形图形 3 季节调整的概念季节调整的概念 季节性变动的发生，不仅是由于气候的直接影响，而 且社会制度及风俗习惯也会引起季节变动。经济统计中的 月度和季度数据或大或小都含有季节变动因素，以月份或 季度作为时间观测单位的经济时间序列通常具有一年一度 的周期性变化，这种周期变化是由于季节因素的影响造成 的，在经济分析中称为季节性波动。经济时间序列的季节 性波动是非常显著的，它往往遮盖或混淆经济发展中其他 客观变化规律，以致给经济增长速度和宏观经济形势的分 析造成困难和麻烦。因此，在进行经济增长分析时，必须 去掉季节波动的影响，将季节要素从原序列中剔除，这就 是所谓的“季节调整” (Seasonal Adjustment)。 4 2 2.2 .2.1 .1 X-11X-11季节调整方法季节调整方法 2 2.2 .2 经济时间序列的季节调整方法经济时间序列的季节调整方法 X-11方法是基于移动平均法的季节调整方法。它的特 征在于除了能适应各种经济指标的性质，根据各种季节调 整的目的，选择计算方式外，在不作选择的情况下，也能 根据事先编入的统计基准，按数据的特征自动选择计算方 式。在计算过程中可根据数据中的随机因素大小，采用不 同长度的移动平均，随机因素越大，移动平均长度越大。 X-11方法是通过几次迭代来进行分解的，每一次对组成因 子的估算都进一步精化。 5 2 2.2.2.2.2 X12X12季节调整方法季节调整方法 美国商务部国势普查局的X12季节调整程序是在X11方 法的基础上发展而来的，包括X11季节调整方法的全部功 能，并对X11方法进行了以下3方面的重要改进： (1) 扩展了贸易日和节假日影响的调节功能，增加了季 节、趋势循环和不规则要素分解模型的选择功能； (2) 新的季节调整结果稳定性诊断功能； (3) 增加X12-ARIMA模型的建模和模型选择功能。 6 X12季节调整方法的核心算法是扩展的X11季节调整程序 。共包括4种季节调整的分解形式：乘法、加法、伪加法和对 数加法模型。注意采用乘法、伪加法和对数加法模型进行季节 调整时，时间序列中不允许有零和负数。 加法模型 (2.2.1) 乘法模型： (2.2.2) 对数加法模型： (2.2.3) 伪加法模型： (2.2.4) 7 例例2.1 2.1 利用利用X12X12加法模型进行季节调整加法模型进行季节调整 图图2.12.1a a 社会消费品零售总额原序列社会消费品零售总额原序列 图图2.12.1b b 社会消费品零售总额的社会消费品零售总额的TCI TCI 序列序列 8 图图2.12.1d d 社会消费品零售总额 社会消费品零售总额 I I 序列序列 图图2.12.1c c 社会消费品零售总额的社会消费品零售总额的TCTC序列序列 9 TRAMO(Time Series Regression with ARIMA Noise, Missing Observation, and Outliers)用来估计和预测具有缺失 观测值、非平稳ARIMA误差及外部影响的回归模型。它能够 对原序列进行插值，识别和修正几种不同类型的异常值，并 对工作日变化及复活节等特殊回归因素及假定为ARIMA过程 的误差项的参数进行估计。 SEATS(Signal Extraction in ARIMA Time Series)是基 于ARIMA模型来对时间序列中不可观测成分进行估计。 这两个程序往往联合起来使用，先用TRAMO对数据进 行预处理，然后用SEATS将时间序列分解为趋势要素、循环 要素、季节要素及不规则要素4个部分。 2 2.2.3 .2.3 TRAMO/SEATSTRAMO/SEATS方法方法 10 本节主要介绍利用EViews软件对一个月度或季度时间序 列进行季节调整的操作方法。在EViews工作环境中，打开一 个月度或季度时间序列的工作文件，双击需进行数据处理的 序列名，进入这个序列对象，在序列窗口的工具栏中单击 Proc按钮将显示菜单： 2 2. .2. 2.4 4 季节调整相关操作季节调整相关操作 ( (EViewsEViews软件软件) ) 11 1. 1. X11X11方法方法 X-11法是美国商务部标准的季节调整方法(乘法模型、加法 模型)，乘法模型适用于序列可被分解为季节调整后序列（趋势 循环不规则要素项）与季节项的乘积，加法模型适用于序列 可被分解为季节调整后序列与季节项的和。乘法模型只适用于 序列值都为正的情形。 12 2. 2. Census X12Census X12方法方法 EViews是将美国国势调查局的X12季节调整程序直接 安装到EViews子目录中，建立了一个接口程序。 EViews进 行季节调整时将执行以下步骤： 1给出一个被调整序列的说明文件和数据文件； 2利用给定的信息执行X12程序； 3返回一个输出文件，将调整后的结果存在EViews工 作文件中。 X12的EViews接口菜单只是一个简短的描述，EViews 还提供了一些菜单不能实现的接口功能，更一般的命令接口 程序。 13 调用X12季节调整过程，在序列窗口选择Procs/Seasonal Adjustment / Census X12，打开一个对话框： 14 3. 3. 移动平均方法移动平均方法 X-11法与移动平均法的最大不同是：X-11法中季节 因子年与年有可能不同，而在移动平均法中，季节因子 被假设为是一样的。 15 Tramo(Time Series Regression with ARIMA Noise, Missing Observation, and Outliers)是对具有缺失观测值， ARIMA误差、几种外部影响的回归模型完成估计、预测和插 值的程序。 Seats(Signal Extraction in ARIMA Time Series)是基于 ARIMA模型的将可观测时间序列分解为不可观测分量的程序 。这两个程序是有Victor Gomez 和Agustin Maravall 开发的 。 当选择了Pross/Seasonal Adjustment/Tramo Seats 时， EViews执行外部程序，将数据输给外部程序，然后将结果返 回EViews。 4. 4. tramotramo/Seats/Seats方法方法 16 22.3 .3 趋势分解趋势分解 本章第2节介绍的季节调整方法可以对经济时间序列进 行分解，但在季节调整方法中，趋势和循环要素视为一体不 能分开。本节专门讨论如何将趋势和循环要素进行分解的方 法。测定长期趋势有多种方法，比较常用的方法有回归分析 方法、移动平均法、阶段平均法(phase average，PA方法)、 HP滤波方法和频谱滤波方法（frequency (band-pass) filer， BP滤波）。本节主要介绍HP滤波方法和BP滤波方法。 17 22.3.1 .3.1 HodrickHodrick-Prescott-Prescott（HPHP）滤波滤波 在宏观经济学中，人们非常关心序列组成成分中的长期 趋势，Hodrick-Prescott滤波是被广泛使用的一种方法。该 方法在Hodrick and Prescott(1980) 分析战后美国经济周期的 论文中首次使用。我们简要介绍这种方法的原理。 设Yt是包含趋势成分和波动成分的经济时间序列，YtT是 其中含有的趋势成分， YtC是其中含有的波动成分。则 (2.3.1) 计算HP滤波就是从Yt中将YtT 分离出来 。 18 一般地，时间序列Yt中的不可观测部分趋势YtT常被定 义为下面最小化问题的解： (2.3.2) 其中：c(L)是延迟算子多项式 (2.3.3) 将式(2.3.3)代入式(2.3.2)，则HP滤波的问题就是使下面损 失函数最小，即 (2.3.4) 19 最小化问题用c(L)YtT2 来调整趋势的变化，并随着 的增 大而增大。这里存在一个权衡问题，要在趋势要素对实际序列 的跟踪程度和趋势光滑度之间作一个选择。 = 0 时，满足最小 化问题的趋势等于序列Yt； 增加时，估计趋势中的变化总数 相对于序列中的变化减少，即 越大，估计趋势越光滑； 趋 于无穷大时，估计趋势将接近线性函数。一般经验地， 的取 值如下： 20 使用Hodrick-Prescott滤波来平滑序列，选择Procs/ Hodrick Prescott Filter出现下面的HP滤波对话框： 首先对平滑后的序列给一个名字，EViews将默认一个名字， 也可填入一个新的名字。然后给定平滑参数的值，年度数据取100 ，季度和月度数据分别取1600和14400。不允许填入非整数的数据 。点击OK后，EViews与原序列一起显示处理后的序列。注意只 有包括在当前工作文件样本区间内的数据才被处理，平滑后序列 区间外的数据都为NA。 21 例例2.3 2.3 利用利用HPHP滤波方法求经济时间序列的趋势项滤波方法求经济时间序列的趋势项T T 利用HP滤波方法求中国社会消费品零售总额月度时间序 列和中国GDP季度时间序列的趋势项。 图图2.4 2.4 蓝线表示蓝线表示GDPGDP序列、序列、 红线表示趋势红线表示趋势T T序列序列 图图2.5 2.5 蓝线表示社会消费品零售总额、蓝线表示社会消费品零售总额、 红线表示趋势红线表示趋势T T序列序列 22 例例2.4 2.4 利用利用HPHP滤波方法求潜在产出和产出缺口滤波方法求潜在产出和产出缺口 设Yt为我国的季度GDP指标，利用季节调整方法将GDP 中的季节因素和不规则因素去掉，得到GDP_TC序列。本例的 潜在产出Y*，即趋势利用HP滤波计算出来的YtT来代替， GDP的循环要素Yt序列由式(2.3.6)计算： (2.3.6) 图图2.6 2.6 蓝线表示蓝线表示 GDP_TC GDP_TC 、 红线表示趋势序列红线表示趋势序列GDP_TC_HPGDP_TC_HP 图 图2.7 2.7 GDPGDP的循环要素的循环要素 序列序列 23 图2.7显示的GDP的循环要素YtC序列实际上就是围绕趋 势线上下的波动，称为GDP缺口序列。它是一个绝对量的产 出缺口。也可以用相对量表示产出缺口，本例用Gapt来表示相 对产出缺口，可由下式计算得到： (2.3.7 ) 图图2.8 2.8 通货膨胀率通货膨胀率( (红线红线) ) 产出缺口产出缺口Gap (Gap (蓝线蓝线) ) 24 2.3.2 2.3.2 频谱滤波（频谱滤波（BPBP滤波）方法滤波）方法 20世纪以来，利用统计方法特别是时间序列分析方法研 究经济时间序列和经济周期的变动特征得到越来越广泛的应 用。自时间序列分析产生以来，一直存在两种观察、分析和 解释时间序列的方法。第一种是直接分析数据随时间变化的 结构特征，即所谓时域（time domain）分析法，使用的工 具是自相关（或自协方差）函数和差分方程；另一种方法是 把时间序列看成不同谐波的叠加，研究时间序列在频率域（ frequency domain）里的结构特征，由于这种分析主要是用 功率谱的概念进行讨论，所以通常称为谱分析。 25 谱分析的基本思想是：把时间序列看作是互不相关 的周期（频率）分量的叠加，通过研究和比较各分量的 周期变化，以充分揭示时间序列的频域结构，掌握其主 要波动特征。因此，在研究时间序列的周期波动方面， 它具有时域方法所无法企及的优势。 26 BPBP滤波滤波的操作的操作 在EViews中，可以使用 Band-Pass 滤波对经济时间 序列进行趋势循环分解。在序列对象的菜单中选择 Proc/ Frequency Filter，显示如下所示的对话框。 27 为了使用Band-Pass滤波，首先要选择一种滤波 类型。共有3种类型： （1） BK固定长度对称滤波（Fixed length symmetric (Baxter-King，BK)）； （2）CF固定长度对称滤波（Fixed length symmetric (Christiano-Fitzgerald，CF)）； （3）全样本长度非对称滤波（Full sample asymmetric(Christiano-Fitzgerald)）。 EViews默认的是BK 固定长度对称滤波。如果 使用固定长度对称滤波，还必须指定先行/滞后（ Lead/lag）项数n。 28 用户必须选择循环周期（Cycle periods）的区间以计 算Band-Pass滤波的频率响应函数的权重序列。这个区间 由一对数据（PL，PU）描述，PL、PU 由Band-Pass滤波要 保留的循环波动成分所对应的周期来确定。月度数据填月 数；季度数据填季度的个数。EViews将根据数据类型填入 了默认数值。例如，例2.6认为中国社会消费品零售总额的 增长周期大约在1年半（18个月）到5年（60个月），如果 保留在这个区间内的循环要素，则区间的下界是18，上界 是60。因此，设定PL=18，PU=60（相当于例2.6中的 p和q ）。 29 在Band-Pass滤波的输出结果中，左侧的图描述了原序 列、趋势序列和循环序列。对于BK和CF固定长度对称滤波 而言，Eviews 画出频率响应函数w()，频率 的区间是0， 0.5，右面的图描述了频率响应函数。但是，对于时变的CF 滤波，并没有画出频率响应函数，因为滤波的频率响应函数 随数据和观测值个数变化。 用户需要输入希望保存的结果（循环成分、趋势成分 ）对象的名字。循环序列（Cycle series）是包含循环要素的 序列对象；趋势序列(Non-cyclical series)是实际值和循环序 列的差。用户还能得到在滤波中所用的Band-Pass滤波频率 响应函数的权序列，它将存储在矩阵对象中。 30 例例2.6 2.6 利用利用BPBP滤波对经济时间序列进行趋势、循环分解滤波对经济时间序列进行趋势、循环分解 中国社会消费品零售总额月度时间序列（SL）的取值 范围从1980年1月至2004年8月（附录E表E.5）。由于带通 （BP）滤波的两端各欠n项，为了近期的分解结果没有缺 失值，本例利用ARIMA模型将序列外推到2006年2月。然 后对SL进行季节调整去掉季节和不规则要素，得到只包 含趋势循环要素的序列SL_TC。根据增长率周期波动分析 ，我国社会消费品零售总额的增长率大约存在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滤波的频率响应函数。 31 图图2.12 2.12 红线表示红线表示BPBP n n (p,q)(p,q)滤波频率响应函数滤波频率响应函数 蓝线表示带通滤波的频率响应函数蓝线表示带通滤波的频率响应函数 32 图图2.13 2.13 红线表示红线表示SLSL的原序列的原序列 蓝线表示趋势要素序列蓝线表示趋势要素序列SL_T SL_T 33 图图2.14 2.14 循环要素序列循环要素序列SL_C SL_C 34 2 2.3.2 .3.2 指数平滑指数平滑 指数平滑是可调整预测的简单方法。当只有少数观 测值时这种方法是有效的。与使用固定系数的回归预测 模型不同，指数平滑法的预测用过去的预测误差进行调 整。 1 1. .单指数平滑（一个参数）单指数平滑（一个参数） 2 2. .双指数平滑（一个参数）双指数平滑（一个参数） 3 3. .Holt-Winters Holt-Winters 无季节趋势（两个参数）无季节趋势（两个