时间序列模型的分析课件(PPT 74页).ppt

1、第五章 时间序列模型的分析,本章讨论时间序列的基本概念与相关模型,介绍这些模型的特点、参数估计方法和检验方法.,第一节 时间序列模型简介,时间系列数据(time series data)是以时间顺序排列的数据序列，是经济学研究中的一种常见数据形式。图5-1给出了一个时间序列的基本例子(不同年份的 GDP 与财政收入的时间序列数据)。,时间序列模型(time series model),研究和分析时间序列数据的模型称为时间序列模型，时间序列方法已成为现代计量经济学的重要内容。 时间序列模型按时间序列的特性，分为不同的种类。例如根据时间序列的平稳性，可分为平稳时间序列和非平稳时间序列；而根据时间序

2、列中变量的相关性，又可分为自回归模型、移动平均模型和自回归移动平均模型等。,时间序列的基本模型,(1)模型的函数形式； (2)滞后变量的个数； (3)误差项u的相关性。,一个时间序列模型的设定中，应包括三个要素：,(5-1),第二节 平稳时间序列模型,当时间序列 满足： (1) ， 即均值为常数； (2) ，即方差为常数 ； (3) ，即协方差与两个时间相距的长度 有关，而与时间的具体位置无关。则 称为平稳时间序列，或弱平稳序列或协方差平稳序列。 当时间序列的所有统计性质都不会随时间而发生变化时，称为严格平稳。本章主要采用弱平稳定义。,时序图,时序图通常由横轴表示时间，纵轴表示时间序列值。根据

3、时序图可以直观地了解时间序列的一些基本分布特征。 图5-1是中国GDP和财政收入的时序图。由于GDP和财政收入有明显的递增趋势, 所以都不是平稳时间序列。,平稳时间序列的意义,传统的统计分析通常有如下的数据结构,一般, 时间序列的每个 只取得1个观察值时, 通常这样的数据是无法进行分析的。但由于序列的平稳性则可以有效地解决这一问题。,例 5-1,一个人的健康指标会随饮食、休息等因素而变化。若每隔一段时间 (如一个月) 检查一次身体，则当时间序列是平稳时，就可以根据样本数据估计健康指标的均值、方差和协方差。,白噪声时间序列(White noise time series),则称时间序列 为白噪声

4、序列，或纯随机序列。 根据平稳时间序列的定义，白噪音序列是一个平稳时间序列。,如果时间序列 满足如下性质：,一、移动平均过程(MA ),其中 和 为常数，而随机干扰 为白噪声。(5-3)式 称为一阶移动平均过程，记为MA(1)。,移动平均过程(moving average process，简称MA)是最基本的平稳时间序列模型，最简单的移动平均过程为,(5-3),由于时间序列 为白噪音序列，则 的均值、方差和协方差 分别为,MA(1)的统计性质,(5-4),所以MA(1)是一个平稳时间序列，并可得相关系数,(5-5),q 阶移动平均过程,q 阶移动平均过程,(5-7),即MA (q)的方差和协方

5、差只与滞后期数q有关，从而MA(q)也是一个平稳时间序列。为了方便, 今后记,无穷阶移动平均过程,q 阶移动平均过程可推广到无穷阶移动平均过程, 记为MA(), 当 绝对收敛(绝对可加)或平方收敛(平方可加)，即,则MA()仍为平稳时间序列。,(5-9),二、自回归过程(AR),一阶自回归模型,AR(1)的统计性质,1 均值,根据假定，由于,于是,2 方差,根据假定，有,则,由此即得AR(1)的平稳条件 。,3 协方差,对于AR(1)，协方差有如下递推公式,于是,而对于AR(1)，已知,得,(5-14),AR(1)模型的相关图,对于AR(1)模型, 自相关系数,例如，当 时，则,这表明对于平稳

6、时间序列，相隔越远，则过去值对现期值得影响越小。图5-2给出了相应的自相关系数序列，称为相关图。 平稳时间序列的自相关系数的衰减速度很快, 而非平稳时间序列的衰减速度则较慢。,(5-15),样本估计,对于平稳时间序列, 序列的均值为常数。表明这时的每个随机变量 的均值都相等, 即 。从而 成为 的样本观察值, 则 的估计为,同样, 可得协方差的估计为,而自相关系数的估计为,例,1964 1999年我国纱年产量的数据如下, 试给出纱年产量的自相关图。,样本自相关图,二阶自回归模型,二阶自回归模型记为AR(2),假定AR(2)为平稳序列，则,即,(5-16),中心化模型,由于AR(2)为平稳序列时

7、， 代入(5-16),记 , 则得中心化AR(2)模型,两边同乘 后求期望，得,(5-17),根据AR(p)的基本假定,AR(2)的方差与协方差,特征方程(characteristic equation),例如，对AR(2)的中心化模型,记,(5-17),(5-19),于是(5-17)可改写为,把(5-19)中的 L 作为变量并用 z 表示，则,称为AR(2)的特征方程。,特征根判别,上述用系数 判别模型的平稳性方法称为平稳域判别, 另一种方法则是特征根判别。 例如，AR(2)平稳的条件是特征方程的根都落在单位圆之外，即,其中 。 对平稳的AR模型，参数估计可直接采用OLS 法，估计量是无偏和

8、一致的。,例,检验AR(2)模型 的平稳性。,(1)由于,根据平稳域判别，该模型时平稳的 (2)由于特征方程,的两个根分别为5和10/9=1.11，即两个根都大于1，从而根据特征根判别，该模型平稳。 在实际应用中，模型参数是未知的,从而无法用上述方法判别模型的平稳性，下一节将讨论有关的统计方法。,AR(p)模型阶数 p 的确定,AR(p)模型阶数 p 通常可由偏自相关系数来选择，也可根据赤池信息准则 (Akaike information criterion ,简记AIC )来选择，在假设模型的阶数为k时计算,其中RSS为残差平方和。AIC准则选择AIC(k)的最小值所对应的正整数 k 为AR

9、(p)模型的阶数。,(5-21),AR(2)的自相关系数方程,三、ARMA 过程,具有如下结构的模型称为自回归移动平均模型，简记为ARMA (p, q),(5-23),其中 。显然，当 q = 0 时，ARMA (p, q)模型成为AR (p) 模型，而当 p = 0 时，ARMA (p, q) 模型就成为 MA (q)模型。所以, ARMA 模型同时包括自回归模型和移动平均模型，而 AR (p) 模型和 MA (q)模型实际上是ARMA (p, q) 的特例。,ARMA (p, q) 模型的平稳条件,由于 q 阶移动平均过程是平稳的，从而 ARMA (p, q)的平稳性完全由自回归部分的平稳

10、性所决定，于是 ARMA (p, q)模型的平稳条件是,的根都落在单位圆之外。,ARMA (p, q) 模型的均值,设ARMA (p, q)模型满足平稳条件，则均值,即,于是(5-23)可转换为,(5-24),ARMA (p, q) 模型的协方差,(5-24)两边同乘 得,根据AR(p)模型的基本假定,于是,得,(5-25),过度参数化,对给定的样本数据，可能会有多个可以选择的模型。在实际工作中，应尽可能采用较简单的模型，即避免过度参数化。 例如，对于ARMA(3, 3)模型，若自回归项和移动平均项有相同的特征根，则删去后就可以减少模型的参数，ARMA(3, 3)则简化为ARMA(2, 2)，

11、避免了模型过度参数化。,例,1950-1998年北京市城乡居民定期储蓄所占比例的数据如下,1 时序图,时序图表明北京市城乡居民定期储蓄占储蓄存款余额的比例比较稳定, 大致在80%左右。,2 自相关图,序列的自相关系数递减速度很快, 在延迟8项后自相关系数在零附近摆动, 可以认为该序列是平稳的。,3 模型识别,由于认为该序列是平稳的,从而可拟合 ARMA 模型。又由于自相关系数递减过程相当连续, 可认为不截尾(拖尾), 应选AR 模型。先选 AR(1)、AR(2)和AR(3), 通过比较后,再确定阶数。,AR(1)模型,AR(2)模型,AR(3)模型,4 选择合适模型,由于,从而可选择AR(2)

12、, 但由于AR(2)中 的系数不显著, 最后选择AR(1)模型, 得,第三节 非平稳时间序列模型,当时间序列 不满足弱平稳条件，即 不是弱平稳序列，则称为非平稳时间序列。 平稳时间序列 的均值和方差为常数，但许多常用的经济时间序列不具有平稳性，而是持续增大或减小。图5-3给出了商品零售价格指数序列和GDP序列。商品零售价格指数序列基本保持平稳，但GDP序列持续上升。这表明仅采用平稳时间序列模型来研究和分析经济时间序列数据是不够的，还需要研究非平稳时间序列。,趋势平稳模型与单位根过程,对非平稳序列，一般采用如下两种方法进行处理 1、趋势平稳(trend-stationary)模型 这种方法是引入

13、一个非随机的时间趋势项，例如,于是均值不再是常数， 但去掉趋势项 ，模型则为平稳序列。 2、单位根过程(unit root process) 这种方法是采用单位根过程,由于(5-27)的特征方程根为1(称为单位根)，从而(5-27)是非平稳序列,(5-26),(5-27),随机游动过程(random walk process),在经济学研究中，单位根过程是一种最常见的非平稳过程。而随机游动过程是一种典型的非平稳序列模型,其中 为白噪声序列。(5-28)是随机游动的基本模型，而(5-27)则是一种扩展。,(5-28),趋势平稳模型与单位根过程的主要区别,(1) 在单位根过程中， 都对 存在影响。

14、而对趋势平稳模型， 只影响 。 (2) 对趋势平稳模型，去掉趋势项 后模型为平稳序列。记,即差分后的模型是平稳过程，从而单位根过程也称为差分平稳过程。,称为 的一阶差分。对于单位根过程,(5-29),一、非平稳时间序列简介,对于非平稳序列，例如随机游动序列，经典回归分析方法不再适用。 假定 和 是两个互不相关的随机游动序列，即,为了方便，设 。则,若采用线性回归模型,不妨设 ，于是,(5-30),(5-31),误差项的方差,根据(5-31)，误差项的方差,即误差项方差与 t 有关，会随 t 的增加而无限增大。由于存在异 方差，相应的 t 检验和 F 检验不再适用。 这种具有单位根过程的非平稳序

15、列时的回归分析，也称为伪回归(spurious regression)。,例5-3,设 和 是互不相关的随机游动过程(图5-4)，采用(5 -30)的线性回归模型，有,T 检验显著，即认为 。 而这与 与 不相关矛盾。图5-5表明误差项很可能是非平稳的，而这时的 t 统计量不再服从 t 分布。 例5-3说明 和 虽然不相关，但伪回归却得到了相反的结果。因此，在时间序列分析时，应重视伪回归问题。,时间序列的单积性,二、单位根时间序列的检验,在计量经济学中，最常见的非平稳序列就是单位根过程。因此，最常见的非平稳序列检验就是单位根检验。 有一些很平稳的序列，可以通过时间序列图做出判断，但对复杂的时间

16、序列，则应采用统计检验方法做出判断。,三种模型,考虑如下的三种模型,(5-33)是一个简单的随机游动模型，(5-34)则带有移动量(截距)a, 而(5-35)还带有趋势项 。,迪基富勒(Dickey-Fuller)检验(简称DF检验),检验,对如下的原假设,其中 为限制条件的个数，n 是样本量，而 k 表示非限制模型中的系数个数，于是n - k 为非限制模型的自由度。查DF检验的 分布表可得临界值 (三个模型的临界值是不同的)。,检验统计量分别为 和 。且,(5-36),拓展迪基富勒检验(简记ADF检验),DF检验设定的三种模型都是一阶自回归序列，而对于多阶自回归序列，则需对模型进行整理转化。

17、考虑 p 阶自回归模型,(5-37),可转化为如下拓展模型,对(5-38)可采用DF检验 。也可对模型增加截距项或趋势项。ADF检验统计量 ，并采用DF检验的临界值。,(5-38),例,利用1970至1991年期间美国的GDP数据说明DF检验。,1 带有移动量的随机游动模型,带有移动量随机游动模型的DF检验,查DF统计量分布表得,由于,从而接受 , 认为GDP序列存在单位根。,2 带有移动量和线性时间趋势量的随机游动模型,带有移动量和线性时间趋势量随机游动模型的DF检验,查DF统计量分布表得,由于,从而接受 , 仍认为GDP序列存在单位根。,第四节 向量自回归模型,把自回归模型拓展到向量，称为

18、向量自回归模型(vector autoregression model, 简称VAR模型)。向量自回归模型非常适合于研究各种变量之间的关系, 基本形式为,(5-42),其中 为 向量，即包含n个变量， 。(5-42)常称为VAR模型的标准形式，可认为是结构模型的简约型。 本节假定(5-42)中 都是平稳序列，分别讨论VAR模型的识别、估计和检验以及一些重要的应用。,例 5-4,设双变量自回归模型,(5-43),(5-44),其中 和 为白噪声且不相关。(5-43)和(5-44)往往建立在经济学理论的基础上，常称为结构型。通过变换，可将结构型转换为简约型。记,则(5-43)和(5-44)可联合表

19、示为,从而简约型为,(5-45),(5-46),其中 。,一、向量自回归模型的识别与估计,根据模型识别方法，例5-4的模型不可识别，从而不能由简约型参数估计值给出结构型参数的估计值。 若 ，则,(5-47) 和 (5-48) 成为递归模型。即 (5-48) 只包含一个内生变量 ，其余都是前定变量，而(5-47)包含两个内生变量，其中 为前一个方程的内生变量。这种递归模型的内生变量之间并不相互影响，而只是单向影响。例如(5-47)和(5-48)表明 影响 ，而 不影响 。,(5-47),(5-48),递归模型的参数估计,根据递归模型(5-47)、(5-48), 得简约型模型,(5-49),(5-

20、50),其中,即(5-50)中系数与(5-48)完全相同。 对于递归模型, 各方程可分别采用OLS估计。,似乎不相关回归模型(简称SUR模型),对简约型模型(5-49)和(5-50), 记,由于 和 不相关, 得,(5-51),即协方差矩阵,因而, 每个方程的误差项都不相关, 但不同方程之间的误差项相关, 这种模型称为似乎不相关回归模型( seemingly unrelated regression model), 即SUR模型。,SUR模型的参数估计,对SUR模型, 有两种参数估计方法。 (1) 广义最小二乘估计(GLS估计): 当协方差矩阵末知时, 则需要先估计协方差, 即采用可行广义最小二乘估计(FGLS估计)。 (2)最大似然估计(ML估计): 由于VAR模型存在滞后变量, 在采用ML估计时, 似然函数是在以前的 p 个观察值(相当于解释变量观察值)条件下, 其余观察值(相当于被解释变量观察值)的条件似然函数。 GLS估计和ML估计都是一致和有效的估计。,二、向量自回归模型的检验(1) 滞后期长度的确定,在VAR模型中, 一个重要的问题是滞后期长度的确定, 即确定应包含多少个滞后变量。通常可采用对数似然比检验。,例,例如对季度数据 取8个滞后变量 取12个滞后变量 LR检验统计量,其中 和 分别是取8个滞后变