1、上海财经大学统计学系 1多元时间序列分析多元时间序列分析 多元平稳时间序列建模多元平稳时间序列建模 虚假回归虚假回归 单位根检验单位根检验 协整协整 误差修正模型误差修正模型 上海财经大学统计学系 210.1 多元平稳时间序列建模多元平稳时间序列建模 1976年，Box和Jenkins采用带输入变量的ARIMA模型为平稳多元序列建模。构造思想：假设输出变量序列（因变量序列）和输入变量序列（自变量序列），均平稳，首先构建输出序列和输入序列的回归模型，如果有必要，使用ARMA模型继续提取残差序列 中的相关信息。模型形为 tytx1tx2ktxttttkkitliitaBBxBBByi)()()()
2、(1上海财经大学统计学系 3例例10.1 在天然气炉中，输入的是天然气，输出的是CO2，CO2的输出浓度与天然气的输入速率有关。现在以中心化后的天然气输入速率为输入序列，建立CO2的输出百分浓度模型。时序图及样本自相关图直观显示输入序列和输出序列均平稳 上海财经大学统计学系 4-3-2-1012350100150200250XT44485256606450100150200250YT上海财经大学统计学系 5 不考虑输入序列和输出序列之间的关系，将它们分别作为一元时间序列进行分析 天然气输入速率序列 模型为：CO2的输出浓度序列 为AR(1，2，4)疏系数模型：txtytt 1t 2t 3tx1
3、.9755x1.3740 x0.3429xtt 1t 2t 4ty 53.6736 2.1066y1.3394y0.2123y上海财经大学统计学系 6 考虑到输出CO2浓度和输入天然气速率之间的密切关系，将输入天然气速率作为自变量考虑进输出序列的模型中，进一步研究二者之间的关系。滞后k期协方差函数定义为 滞后k期协相关系数为 )(ktktttkExxEyyECovtt kktt kCov(y,x)CVar(y)Var(x)上海财经大学统计学系 7输入序列 和输出序列 的协相关图 txty上海财经大学统计学系 8从协相关图可以看出，输出序列和输入序列的滞后项有显著的相关关系，且滞后阶数比较多，考
4、虑采用ARMA模型结构，以减少待估参数的个数。通过反复尝试，得出以下回归模型tt 1t 2t 3t 6ty 3.6 7 2 9 1.4 4 5 6y0.5 1 4 3y0.5 0 9 3x0.3 0 0 2x上海财经大学统计学系 9再考虑回归残差序列 的性质，从残差序列的时序图和相关图可以看出，残差平稳且不存在序列相关性，说明拟合模型有效。t上海财经大学统计学系 10模型拟合效果图 返回 44485256606450100150200250YTYTF上海财经大学统计学系 1110.2 虚假回归虚假回归 当因变量序列 和输入变量序列（即自变量序列），都平稳时，可以依据Box和Jenkins的理论
5、和方法构建以输入变量为自变量的ARIMAX回归模型来拟合相应序列的变化。当平稳性条件不满足时，我们就不能大胆地构造ARIMAX模型，因为这时容易产生虚假回归的问题。返回 tytx1tx2ktx上海财经大学统计学系 1210.3 单位根检验单位根检验 DF检验检验 ADF检验检验 PP检验检验 上海财经大学统计学系 13DF统计量统计量 考虑1阶自回归序列：单位根检验的原假设和备择假设分别为：t统计量 DF（Dickey-Fuller）检验统计量 时，其极限分布为：tt 1txxi.i.d.2t N(0,)：0H|11H：|1t()S()1DF=S()T 2121/20(1/2)(W(1)1)(
6、W(r)dr)上海财经大学统计学系 14维纳过程具有如下性质：（1）（2）（3）)1,0()1(NW),0()(2rNrW)1(/)(22xrrW上海财经大学统计学系 15DF检验的等价表达检验的等价表达 DF检验可以通过对参数 的检验等价进行：相应的DF检验统计量为：其中，为参数 的样本标准差。0010：HHDFS()S()上海财经大学统计学系 16DF检验方法的三种适用类型检验方法的三种适用类型 第一种类型如式 第二种类型如式 第三种类型如式 tttxx1i.i.d.2t N(0,)t0t 1txx i.i.d.2t N(0,)t0t 1txatx i.i.d.2t N(0,)上海财经大学
7、统计学系 17例例10.2对某国1960年到1993年GNP平减指数 的季度时间序列进行DF单位根检验。1.直观判断：GNP平减指数的季度时间序列绘制时序图，时序图显示序列显著非平稳。4080120160200240255075100125PT上海财经大学统计学系 18 2.对该时间序列进行DF检验 上海财经大学统计学系 19ADF检验检验 DF检验只适用于1阶自回归过程的平稳性检验，但是实际上绝大多数时间序列不会是一个简单的AR(1)过程。为了使DF检验能适用于AR(p)过程的平稳性检验，对DF检验进行了一定的修正，得到增广DF检验（Augmented Dickey-Fuller），简记为A
8、DF检验。上海财经大学统计学系 20ADF检验的原理检验的原理 对任意一个AR(p)过程 AR(p)过程单位根检验的假设：构造ADF检验统计量：其中，为参数 的样本标准差。t1t 1ptptxxx0010：HHADFS()S()上海财经大学统计学系 21ADF检验的三种适用类型检验的三种适用类型 第一种类型 第二种类型 第三种类型 tt 11t 1p-1t p+1txxxx i.i.d.2t N(0,)t0t 11t 1p-1t p+1txxxx i.i.d.2t N(0,)t0t 11t 1p-1t p+1txatxxx i.i.d.2t N(0,)上海财经大学统计学系 22例例10.2续续
9、 对某国1960年到1993年GNP平减指数的季度时间序列进行ADF检验，检验结果入下：上海财经大学统计学系 23PP检验检验 针对序列可能存在高阶相关的情况和可能的异方差情形，Phillips和Perron于1988年对ADF检验进行了非参数修正，提出了Phillips-Perron检验统计量。PP检验统计量既可适用于异方差场合的平稳性检验，又服从相应的ADF检验统计量的极限分布。上海财经大学统计学系 24例例10.2续续 对某国1960年到1993年GNP平减指数的季度时间序列进行PP检验，结果如下：返回 上海财经大学统计学系 2510.4 协整协整 单整与协整单整与协整 协整检验协整检验
10、 上海财经大学统计学系 26单整（单整（integration）的概念）的概念 在单位根检验的过程中，如果检验结果显著拒绝原假设，即说明序列 显著平稳，不存在单位根，这时称序列 为零阶单整序列，简记为 。如果原序列1阶差分后平稳，说明原序列存在一个单位根，这时称原序列为1阶单整序列，简记为如果原序列至少需要进行d阶差分才能实现平稳，说明原序列存在d个单位根，这时称原序列为d阶单整序列，简记 。tx txtxI(0)1(Ixt)1(IxttxI(d)上海财经大学统计学系 27单整序列的性质单整序列的性质 1若 ，对于任意非零实数a与b，有 2若 ，对于任意非零实数a与b，有 3若 ，对于任意非零
11、实数 a与b，有 4.若 ，对于任意非零实数 a与b，有 式中，)0(Ixt)0(Ibxat)(dIxt)(dIbxatttx I(0),y I(0)0(Ibyaxzttt)(dIxt)(cIyt)(kIbyaxztttmax,kd c上海财经大学统计学系 28协整（协整（cointegration）的概念）的概念 协整理论是Engle and Granger在1987年首先提出来的。设随机向量 中所含分量均为d阶单整，记为 。如果存在一个非零向量 ，使得随机向量 ，则称随机向量 具有d,b阶协整关系，记为 向量 被称为协整向量。tXtXI(d)ttYXI db0b tX上海财经大学统计学系
12、29协整检验协整检验 Engle-Granger两步协整检验法两步协整检验法 Johansen协整检验法协整检验法 上海财经大学统计学系 30Engle-Granger两步协整检验法两步协整检验法 1用ADF检验各变量的单整阶数。协整回归要求所有的变量都是一阶单整的，因此，高阶单整变量需要进行差分，以获得序列。2用OLS法估计长期动态回归方程，然后用AD检验残差估计值的平稳性。(1)I上海财经大学统计学系 31Johansen协整检验法协整检验法 Johansen和Juselius提出的一种在VAR（向量自回归）系统下用极大似然估计来检验多变量之间协整关系的方法，通常称为Johansen协整检
13、验。设一个VAR模型如下 t1 t 12t 2ptptYB YB YB YU上海财经大学统计学系 32 向量误差修正模型(VECM)协整关系的个数可通过下面两个统计量来计算：1ptit it ptiYYYUi1iIBB1pIBB1log 1mtraceii rT max1log 1rT上海财经大学统计学系 33 迹检验：即至多有r个协整关系 即有m个协整关系（满秩）最大特征根检验：0:Hrm1:Hrm0:,1,2,H r qqm1:1Hrq上海财经大学统计学系 34例例10.2续续 我们以人均生活费支出对数序列 和可支配收入对数序列 为例来说明如何进行EG检验。时序图 lntylntylntx
14、4.44.85.25.66.06.46.81020304050607080LNYTLNXT上海财经大学统计学系 35 第一步：用ADF检验分别对序列 和 进行单整检验 第二步：用变量 对 进行普通最小二乘回归，得回归模型如下：残差序列时序图 返回 lntxlntytttln y0.25300.9355lnxlntylntx-.15-.10-.05.00.05.10.15.201020304050607080RESID上海财经大学统计学系 3610.5 误差修正模型误差修正模型 误差修正模型(Error Correction Model)简称为ECM，它的主要形式是由Davidson、Hendr
15、y、Srba和Yeo于1978年提出的，也称为DHSY模型。它常常作为协整回归模型的补充模型出现。E-G两步法建立误差修正模型 第一步，先检验两个变量的单整阶数，如果都是1阶单整，紧着着进行回归（OLS法），检验变量间的协整关系，估计协整向量（长期均衡关系参数）；第二步，若协整性存在，则以第一步求得的残差作为非均衡误差项加入到误差修正模型中，并用OLS法估计相应参数。上海财经大学统计学系 37例例10.2续续 对1992年1月到2019年12月经居民消费价格指数调整的中国城镇居民月人均生活费支出对数序列和可支配收入对数序列构建ECM模型 0111lnlntttyxECM111ln0.2530 0.9355lntttECMyx
