多变量时间序列模型

第7章 多变量时间序列模型§1 Granger 因果检验判断一个变量的变化是否是另一个变量变化的原因，是经济计量学中的常见问题。

Granger 提出一个判断因果关系的一个检验，这就是Granger 因果检验。

一．Granger 因果检验的思想如果x 影响y ，或者x 是y 的原因，此时x 的变化必然先于y 的变化，此时就须满足两个条件：1）x 可以预测y ，即根据y 的过去值对y 进行回归时，如果加上x 的过去值，能显著增强回归的解释能力。

2）不能根据y 预测x ，因为如果根据x 预测y ，又能根据y 预测x ，很可能x 和y 都是由第三个或其他变量决定。

二．检验步骤1）首先检验零假设“x 不影响y ”即x 不是y 的granger 原因。

首先根据x 和y 的滞后值对y 回归（无限制回归），然后用y 的滞后值对y 进行回归（有限制回归）。

即：无限制回归：0111mmt i t ij t j t i i y yx ααβε--===+++∑∑有限制回归：021mt i t it i y yααε-==++∑用F 检验来判断x 是否显著了了无限制回归的解释能力。

此时统计量()()()*,11RSS RSS rF F r n k RSS n k -=----*RSS 是有限制回归的残差平方和，RSS 使无限制回归的残差平方和。

n 是样本容量，k使无限制回归变量的个数。

r 是限制回归模型中的变量个数。

2）检验y 不影响x ，即x 不是y 的granger 原因。

此时调换模型中的变量x 和y 的位置，利用F 统计量来检验。

三．如果一对时间序列是协整的，则至少在某一方面存在granger 原因。

§2 伪回归一．现在考虑非平稳序列回归出现的问题设{}{} t t y x 是两个无关的随机游走过程，()()2121 0, 0,t t t t t t t t uy y IID x x u u IID εεεσσ--=+=+且()0s t E u ε=，任意的,t s 。

此时t x 既不影响t y ，也不受t y 的影响。

建立模型t t t y c x e β=++ 对于这种情况，0pβ⇒，20R → 。

Banerjee 利用蒙特卡罗模拟得出结论：对于相互独立的单整序列{}{} t t y x ，且001,1,0,0u x y εσσ====，进行回归时，t 统计量显示了比正常检验临界值水平还高。

也就是在相互独立的序列进行的实际回归中，经常伴随着高的2R ，并且β系数显著。

这种现象就称为为回归现象。

二．伪回归的产生原因1．伪回归现象产生的根本原因就是序列的非平稳性。

设{}{} t t y x 是两个无关的随机游走过程，并假设模型为t t t y x e αβ=++，()20,t u v IN σ ，此时建立统计推断：原假设0H ：0β=此时根据OLS 估计，2ˆt ttx yxβ=∑∑。

因为{}{} t t y x 是两个无关的随机游走，则ˆβ的分子，分母均为一种维纳过程泛函。

（参见hendry 和秦朵的《动态经济计量学》第三章）。

2．()W r 和()V r 表示定义在[]0,1r ∈上的两个相互独立的维纳过程，则[]()1Tr tt W r ε=⇒[]()1Tr tt uV r =⇒于是可以证明：()2122201T tt TyW r dr εσ-=⇒⎡⎤⎣⎦∑⎰()2122201Ttu t TxV r dr σ-=⇒⎡⎤⎣⎦∑⎰ ()()121Tt t u t Tx y W r V r dr εσσ-=⇒∑⎰于是()()()()()()11221120ˆ0u u u W r V r drW r V r drV r drV r drεεσσσβσσ⇒⇒→⎡⎤⎡⎤⎣⎦⎣⎦⎰⎰⎰⎰但在零假设0β=条件下，1222222221111T T T T t t t t t t t t t t t T y x T x T y T y x -----====⎛⎫⎤⎡⎤⎡⎤⎡⎤=- ⎪⎥⎢⎥⎢⎥⎢⎥⎪⎦⎣⎦⎣⎦⎣⎦⎝⎭∑∑∑∑ 此时有两个特点：1）发散。

2）分布不是常规的分布。

从而利用t 统计量不能得到正确的检验。

直观地看一下：t t t y x e αβ=++，如果0β=，则t e 也是()1I 过程非平稳。

从而最小二乘估计不可用。

三．伪回归检验我们利用残差的平稳性检验来判断是否存在伪回归。

如果残差非平稳，则是伪回归。

四．伪回归的纠正方法：1）在回归模型中包含自变量和因变量的一阶滞后变量，即11t t t t t y x y x u αβφδ--=++++此时能够消除伪回归。

即当y 和x 不相关，则ˆβ和ˆδ依概率收敛于零。

2）对t y 和t x 先做一阶差分，从而使得{}t y ∆和{}t x ∆变成平稳过程，建立模型t t t y x u αβ∆=+∆+此时能够消除伪回归。

即当y 和x 不相关，则ˆβ依概率收敛于零。

3）Cochrane-Orcutt 方法（自相关问题） 如果t t t y x u αβ=++，这里1t t t u u ρε-=+ 则根据广义差分法，建立模型()()111t t t t t y y x x ραρβρε---=-+-+进行迭代估计，可以证明ˆβ依概率收敛于零。

五．总结1）伪回归现象：对于任何两个（或两个以上）不相关的单位根过程，只要样本量足够大，检验他们相关性的统计量一定呈显著性，这就是伪回归现象。

2）回归分析将平稳过程当作非平稳过程来处理是十分危险的。

因此回归中必须分清平稳过程和非平稳过程。

3）伪回归的本质问题是变量的非平稳性。

§3 协整一．定义1．零阶整()0I ：对于n 维线性过程0t i t ii X C ε∞-==∑，如果0ii C∞=≠∑，i C 为系数，()20,t IID εσ ，则称t X 为零阶整的。

2．1阶整()1I ：如果t x ∆为()0I 过程，则n 维随机过程t X 称为1阶整的。

()1I 过程都是非平稳的，但非平稳性可以通过差分消除。

这里任何形式为1t t i i t ii i X C C εε∞-===+∑∑，当0C ≠时，都是()1I 过程，3．协整：如果t X 为1阶整的，但对于某些线性组合，只要恰当的选择0X β'，0β≠，t X β'变为平稳的，则t X 称为协整，β称为协整向量，线性不相关的协整向量数目称为协整的秩。

协整向量张成的空间称为协整空间。

即表示为()0t X I β' 。

二．说明1）协整向量β至少有两个非零元素。

2）对于常数C ，C β都是协整向量。

3）对于k 个线性无关的t X ，()1,...,k βββ=满秩，则t X β'平稳。

三．协整的意义1．协整可以用来描述某些经济变量水平值间存在稳定的长期关系，也就是某些经济变量之间不能相互分离太远。

也就是如果变量不存在协整性，他们可以任意分离。

协整向量可以认为是经济系统中对变量长期行为的限制条件。

2．协整响亮多好还是少好协整向量多，经济系统就越稳定，即在尽可能多的方向上存在稳定的关系。

协整向量少，经济模型存在唯一的稳定状态均衡，这样能够获得协整向量的精确估计。

四．协整于长期均衡的联系如果t t t y x u β-=，()0t u I ，(),1t t y x I ，则此时t y 和t x 协整。

这里t u 可以测量“系统”的失衡程度，可以称为均衡误差。

如果是协整的，()0t u I ，均衡误差具有均值回返特性。

如果不是协整的，则()1t u I ，则均衡误差将会出现剧烈波动，不具有均值回返特性。

五．协整和伪回归的关系对于()()2121 0, 0,t t t t t t t t uy y IID x x u u IID εεεσσ--=+=+ ，t t t y x v β=+，Philips 证明伪回归的长期协方差矩阵22u u S u εεεσσσσ⎡⎤=⎢⎥⎣⎦∑非奇异。

当0v ερ=时，此时系统残差非平稳，为伪回归。

而对于协整方程，长期协方差矩阵22u u S uεεεσσσσ⎡⎤=⎢⎥⎣⎦∑奇异。

此时220u u u εεεσσσσ=，从而1uu u εεεσρσσ==。

残差矩阵奇异是系统为协整过程的必要条件。

六．协整检验1．对于回归模型t t t y x u β=+我们检查t u 的生成过程：1t t t u u ρε-=+是否是单位根过程。

此种方法就是E-G 两步法。

2．Machinnon 检验方法：小样本方法 是一种从一般到特殊的检验方法：对于k 个()1I 序列12,,....,t t kt x x x ，1k ≥，1,2,....,t T = 可以建立三种协整回归方程 （1） 11kt j jtt j x xβμ==+∑（2） 102kt jjtt j x xαβμ==++∑（3） 1012kt jjtt j x t xααβμ==+++∑t μ为扰动项。

协整检验的基本方法就是检验其扰动项的平稳性，利用ADF 检验：ˆt u∆()1ρ=-1ˆt u -1ˆli t i t i u γε-=+∆+∑ ˆt u∆()01δρ=+-1ˆt u -1ˆli t i t i u γε-=+∆+∑ ˆt u∆()011t δδρ=++-1ˆt u -1ˆli t i t i u γε-=+∆+∑ 原假设：1ρ=，备则假设1ρ<。

具体检验步骤参见单位根检验步骤。

§3 误差修正模型ECM （E-G 两步法）作用：模型取决于解释变量与因变量长期关系得偏离以及对这些因变量的调整。

一．设两个经济序列t y 和t x 是()1I 过程，并且是协整的，协整向量为()1,β-，则设()*t t t y f x x β==，于是()11t t t t t y x y x u αθβ--∆=∆+-+，则这种模型称为误差修正模型ECM 。

这个模型需要t y ∆对t x ∆，11,t t y x --进行回归。

因此误差修正项11t t y x β---的系数可作为前期的y 对x β之间偏离程度的测定及偏离的短期调整。

从而ECM 不再是使用变量的水平值或变量的差分变量来建模，而是把两者有机的结合起来。

短期看，y ∆是由较稳定的长期趋势和短期波动所决定，短期内系统对于均衡状态的偏离程度的大小直接导致波动振幅的大小。

长期看，协整关系式起到吸引力的作用，将非均衡状态拉回到均衡状态。

ECM 能清楚地显示出于长期均衡偏离的程度，所以立刻显示出关于这种偏离的调整信息。

二．E-G 两步法的估计步骤方法1．建立静态模型：t t t y x u αβ=++，进行最小二乘估计并进行协整检验。