基于向量误差修正模型的两区制门限协整检验

1、基于向量误差修正模型的两区制门限协整检验摘要这篇论文检验了一个具有单协整向量的两区制向量误差修正模型以及误差修正中的门限效应。我们采用了一个相关的单一算法，能够给出二元情况下完全门限协整模型的最大似然估计。我们对门限采用SupLM方法进行检验。我们派生出一个空渐近分布，演示模拟临界值的方式，以及给出一个bootstrap近似。我们研究了使用蒙特卡洛模拟的检验的有效性，发现这种检验方式十分有效。通过使用这样的方法对利率期限结构模型的研究，我们发现，存在着非常强的门限效应。I引言门限协整是由Balke和Fomby（1997）提出的用于研究非线性协整关系的可行的工具。特别值得注意的是，这个模型允许长

2、期均衡的非线性调整。这个模型被广泛应用：BalkeandWohar（1998）、LoandZivot（2001）、Martensetal（1998）、Michaeletal（1997）、O'Connell（1998）、O'ConnellandWei（1997）、ObstfeldandTaylor（1997）以及Taylor（2001）、LoandZivot（2001）对这些方法展开了广泛的回顾。这一系列模型最重要的统计问题在于检验门限效应的存在性（空的线性？）。Balke和Fomby（1997）建议采用Hansen（1996）和Tsay（1989）用来检验误差修正（协整残差）的

3、单因素检验方法。众所周知，在协整向量已知的情况下这种检验方法是有效的，但是Balke-Fomby并没有对将此方法用于估计的协整向量的情况的检验做出理论上的说明。Lo和Zivot（2001）将Balke和Fomby的方法扩展到了具有已知协整向量的多元门限协整模型，采用了Tsay（1998）以及Hansen（1996）对多元扩展模型的检验方法。在本文中，我们将把这些方法扩展到检验未知协整向量的情况中。正如Balk-Fomby所研究的，我们的模型是一个具有单协整向量以及误差修正中存在门限效应的向量误差修正模型（VECM）。然而，与Balke-Fomby只关注于单变量估计以及检验方法不同，我们的估计与

4、检验将关注于完全多变量门限模型。事实上，将误差修正作为门限变量并不是本文分析的必要条件，我们这里所讨论的方法能够非常简单的合并到其他模型中，只要这些模型中的门限变量是前定变量的固定的变形形式。本文有两点贡献。第一，我们提出了一个采用最大似然估计方法估计门限模型的方法。这种算法包含了一个搜索门限以及协整向量的联合网络结构。这个算法在二元情况下非常容易实现，但是在更高纬度的情况下却可能非常困难。此外，正因为如此，我们无法提供对于一致性以及最大似然估计分布理论上的证明。第二，我们发展出一个检验门限效应的方法。原假设不存在门限效应，所以模型简化为一个常规的线性VECM。在此原假设下的估计非常容易，简化

5、为一个常规的降秩回归。这就表明其检验可以基于拉格朗日乘子（LM）准侧，只要求在原假设下进行估计。由于门限参数并非在原假设下决定，我们在SupLM检验的基础上进行扩展推论（参照Davies（1987）、Andrews（1993）以及AndrewsandPloberger（1994）对这种检验思路的动机以及理由）。我们的检验采用了与Seo（1998）为了检验误差修正模型中结构变异所派生出来的方法所具有的相似的代数形式。我们派生出SupLM检验的渐近原分布，发现它与Hansen（1996）用来平稳数据进行门限检验的形式一致。一般来说，渐近分布基于数据的协方差结构，预先排除制表。我们建议采用Hanse

6、n（1996,2000b）的平稳回归元bootstrap方法，或者参数残差bootstrap算法，来近似样本的分布。Section2将介绍门限模型以及从Gaussianquasi-MLE方法中派生出来的用于检验该模型的方法。Section3将采用LM检验方法检验门限协整以及它的渐近分布,并给出两种计算P值的方法。Section4将给出一个关于该检验的大小以及有效性的模拟结果。Section5给出一个在利率期限结构上的应用。渐近分布理论的证明将在附录中给出。本文估计、检验的Gauss程序以及本文中主要工作的副本发布在/bhansen.2估计2.1. 线性协整设叫为

7、一个p维一阶单整时间序列，具有一个px1维的协整向量仔。用叫仔=“叫表示0阶单整的误差修正序列。那么一个滞后2+1阶的线性VECM模型可以简洁的写成：x=Z'X_tt1仔+ut(1)其中，1xab=t_1尸W4Bt1Kxt_1Xt_2X,t_i回归元1仔为一个kx1阶列向量，A为一个kx卩阶矩阵，其中k=pl+2。误差项®为一个具有有限协方差矩阵Y=Euu't鞅差序列（MDS）向量。符号w一1仔和存_仔代表对仔的一般估计值。当估计仔真实值时，我们将分别将其表示为w一1和X一1。我们需要对川故一些标准化处理以使其能够呗识别。由于仅仅存在一个协整向量，那么就可以比较简单的

8、将仔中的一个元素设为1,在系统是双变量（p=2）的情况下不会造成任何优度损失，而在多变量（p>2）的情况下，也只是施加了在协整关系中增加了一个相关系数元素叫的一个约束。参数（仔，A，Y的估计采用最大似然估计的方法，并假设误差项ut服从高斯分布（对仔做一些标准化处理）。估计结果记为（仔，A，Y，回归残差向量为ut=xtA，Xt-i（B。2.2. 门限协整对模型（1）进行扩展为一个两区制门限协整模型：叫=X“B+uifw“1 t-itt-i"B+uifw“2 t1tt_1其中，y为门限参数。可以简化为如下形式:基=Alxt-i仔dit伤卩+码x-1仔d2t伤卩+叫（2）其中，吒B，

9、y=iw1Bsyd2tB，Y=11仔>y1为一个条件状态函数。门限模型（2）存在两区制，有误差修正的值决定。系数矩阵4和码决定了区制机制。模型（2）允许所有系数才两区制之间转换（除了协整向量仔）。在许多的例子中，通过仅允许个别系数在不同区制间能够变换，可以大大的简化模型。这些是模型（2）的特殊形式，对（Z,2施加了约束。例如，仅允许常数项以及误差修正参数w一1可以变动，施加约束使滞后的Xt-j在区制转移时保持不变。门限效应只有当0<Pwt-1<y<1是才存在，否则模型就是简单的一个线性协整模型。我们通过如下假设假设这一约束兀0<PW1<y<1一兀0（3

10、）其中，兀。>0为显著性参数。根据经验性应用，我们设定兀0=0.05。我们通过最大似然估计方法估计模型（2），其前提假设是误差项色服从高斯分布。高斯似然函数为：1 nLn坷，码，丫，伤卩其中，一kuA,A,B,y'1.-1uA,A,B,v2 t12尸/t12尸/t=1叫坷,码，伤v=xtA1xt-1Bd1t仔,Y-A2X-1BB，Y最大似然估计量（坷，码，1，仔，Y为使似然函数几坷，码，1，B，Y达到最大值的参数值。为了计算上的方便，我们首先集中计算出£,码，1。将伤Y固定，只计算坷,码，1的最大似然估计量。这就只是一个特殊的OLS回归：n一1n仔,y=X_1仔X_1仔

11、，d1tB，Yx1Bt_1尸xatn仔,y（4）t=1t=1n一1n仔,y=X-仔X-1仔，d2tB，rx1Bt_1尸X'drt2t仔,y（5）t=1t=1仔,y=ut（A1b，y，A2伤y,B，y以及：1n1,Y=色仔,Yu,y'（6）ntt值得注意的是，式与（5）式分别是当叫一仔<Y以及叫一1仔>y时叫对£一1B的OLS回归结果。这就将似然函数简化为：npTn札伙y=£九（坷伙y，码伤r,丫伙,仔,y=一2log（丫伙y然后，最大似然估计量（仔,Y就是使log（Y伤y最小化的参数值。根据前面对仔的规范化处理的讨论以及施加于式（3）的约束条件1

12、n冗0兀1x；B5Y<1-兀°t=1坷与码的最大似然估计量为坷=州（仔，Y和码=码（B，y。关键的（7）式并非平滑的，所以简便的梯度爬坡算法（FOB？）并不适合计算（7）式的最大化问题。在前面的7=2的例子中，我们建议采用一种网格搜索的方法在二维空间伤y中求解最大化的解。在更高纬度的情况下，网格搜索方法并不合适，其他的所搜方法（例如DorseyandMayer，1995，提出的遗传算法）可能更为有效。但在本文中，由于仔是一个一直的前定变量，网格搜索可以极大的简化计算过程。使用网格搜索，需要选取一个搜索的值域。我们采用线性模型中仔的一致估计（2.1中所讨论的MLE）校准得出的值域

13、。设竹=叫（仔，令yL,乙表示叫的经验性结果，在yL,沧上构建一个平均分布的空间网格。令仔厶,陶表示一个基于线性估计仔（渐近近似正态分布）的置信间隔，再在仔厶,儿上建立平均分布的空间网格。然后，网格搜索算法将在,沧以及处,陶上检验所有符合兀°<1二11x；B<Y<1一兀°的（y，仔。在图1和2中，我们给出了在Section5中的经验性例子的不可微标准函数（基于一年期和十年期债券利率的应用性研究，设定2=1）。图1将式（7）作为y的函数绘制，图2则将其作为仔的函数绘制。综上所述，我们关于p=2情形的算法为：1、基于线性估计量仔在yL,乙与处,兀上构建网格；2

14、、对于网格上的每一组（y，仔的值带入（4）（5）（6）式分别计算坷仔,y、码仔,y、丫伤y；3、找到网格中使log（Y伤y取最小值的（y，仔作为估计值（仔,y；4、计算得出坷=坷（仔,y、A2=A2（仔,y、E=E（仔,y以及气=（仔,y。值得注意的是，在第3步中，并不能保证最小化估计值（仔,y是唯一的，因为似然函数不是一个凹函数。我们介绍了一个用于MLE的算法，但需要强调的是这并不是一个理论上的推断。我们并没有给出估计结果的一致性的证明或者是分布理论。在线性模型中，仔以速率n靠近仔，在一个平稳模型中，y以速率n靠近*那么似乎在门限协整模型中做这样的猜测是合理的：爭,Y以速率n靠近爭，y。在这

15、种情况下，由于仔和y已知，斜率估计值A1和码服从渐近正态分布。那么，这些估计参数的标准差也可以轻松得出。3门限检验3.1.统计学检验令捡表示表示线性VECM模型（1）,如表示两区制的门限模型（2）。这些模型是嵌套模型，H0模型是对如模型施加约束4=码的特殊形式。我们想要比较检验捡（线性协整）与如（门限协整）。我们将致力于正规模型基础上的统计学检验，以能够直接比较这两个模型，并指出模型之间的不同之处。此外，有人可能会建议考虑采用Tsay（1989）的非参数非线性检验和Tsay（1998）的单变量和多变量情况下的检验。但正如BalkeandFomby（1997）以及LoandZivot（2001）

16、的模拟研究，这些非参数检验一般不如基于可比较模型的检验具有说服力。在本文中我们采用LM统计量检验。我们这样选择基于两点：第一，LM统计量的计算非常容易快捷，使得bootstrap可行；第二，一个似然比率或者wald检验需要无约束模型参数估计的分布理论，而本文并不具备相关理论。虽然没有相关证据，我们猜测这些检验渐近等价于LM检验。所以我们推断能够使用LM统计量检验。假设爭，y已知且平稳，则叫假设下的模型为：兔=A'Xt-1仔+叫（8）以及如假设下的模型为：叫=Aixt-1仔d1t伤y+码x-1仔d2t伤y+叫（9）在给定爭，Y时，这些模型是线性的，那么最大似然估计就是有效地。由于（8）是

17、（9）的一个嵌套模型，那么对于异方差情况下依旧具有稳健性质的类似于LM的统计量就可以从模型（9）的线性回归中计算得出。特别的，令£仔,Y和X2伤Y分别表示X一1仔dltp,y和X仔d2tp,y的行叠矩阵，令和伤y和§2伤y分别表示仔d仔,y和色0X仔d仔,y的行叠矩阵,其中色为Section2.1中所定义的线性模型的残差向量，再定义外籍矩阵为：M仔,y=仔,y,X1仔,yM2仔,7=Ip$X2仔,y迟仔,y以及Q伤Y=§1B，Y£伙YQ2B，y=§2B，y飞2伤y然后我们定义叫仔,y和伤y为vecA1仔,Y以及vecA2仔,Y的Eicker-W

18、hite协方差矩阵：伤y=M仔,y-1Q仔,y仔,y-1(10)%伤y=M2仔,y-Q2仔,yM2仔,y-1（11）得出异方差的似LM稳健统计量的标准表达形式：LM仔,y=vecA1仔,y-码伤Y八匕伤Y+%伤X-11仔,y-码伤y（12）如果仔和y已知，式（12）就可以作为检验统计量。若未知，则LM统计量由在捡假设下点估计得出的估计量通过（12）式计算得到。对仔的原估计为仔（Section2.1）,但在捡假设下不存在对y的估计，以致并不能简单的得出LM统计量。根据交并原理，Davies（1987）提出了统计量：supSupLM=yL<y<丫/M邙，Y（13）在本文的检验中，由于设

19、定了搜索域人,Yu，那么5就是wt-1的百分之叫，%就是其百分之1-兀0。为了进行检验,餐宿兀0不能接近于0,Andrews（1993）的研究表明这样做将降低检验的说服力。Andrews（1993）指出将兀0设定在0.05和0.15之间是检验效果非常不错的选择。对于选择（13）式的检验统计量的进一步的理由Andrew（s1993）和AndrewsandPloberger（1994）中做了详细讨论。AndrewsandPloberger（1994）认为，通过对,Y的指数加权平均能够增加检验的可靠性。然而，由于对加权函数的选择问题，加权方法天生具有随意性，所以我们的分析依旧基于式（13）。由于函数

20、LM邙,y对于y不可微，对于（13）式所定义的最大化问题的求解就需要在yL,yu上进行网格估计。当协整向量仔0为一个已知的前定变量时，可以认为仔固定为一个已知值仔0，则检验采取式（13）的形式。我们将这个检验统计量表示为：sup/、SupLM0=yL<y<y/M仔0,Y（14）尤其需要指出的是，使（13）与（14）式达到最大值的y的值可能与Section2中的MLE估计值y不同。事实上，基于如下两点原因确实如此：第一，（13）与（14）式为LM检验，是基于原假设下的参数估计结果而不是备选假设；第二，这些LM统计量是基于异方差-自相关的协方差矩阵估计结果计算得出的，在这种情况下，就算

21、是SupWald统计量的最大化解也与MLE的结果不一致（仅当估计的是同方差协方差矩阵时二者才等价）。这种差异普遍存在与门限检验以及对回归模型的估计当中，并非只存在于门限协整模型中。3.2.渐近分布首先考虑协整向量真实值仔0已知的情况。回归元是平稳的，检验问题是一个对Hansen（1996）的多元概括。检验中所用的渐近分布服从于那篇文章中的形式。我们将遵循标准弱依赖性条件。假设：“兔,莓为4r阶矩有界的，具有严格的平稳性与规律性，具有一个混合斜率耳机=°m-A，其中，川>vv1且r>v>1。此外，误差项为一个MDS（多维标度法？）误差修正“兔具有一个有限密度函数。匕在

22、原假设（8）的情况下，当ut独立同分布且具有有限密度函数以及为一个4r阶矩有界时，这些前提条件就全部满足；在备选假设（9）的情况下，必须对参数进行进一步的约束才能满足这些假设。令F表示叫1的边缘分布函数，令符号“n”表示关于兀。,1-兀0弱收敛性的同一指标（？）。定义0t=Fwtl以及=严J；J比i才以及X=E1比i“叫"徑兀一£；i定理一：在原假设叫下supSupLM0nT=兀冬厂冬1nTr其中，Tv=S*丫ZQ*r1S*rQ*r=QrMrM11QrQrM11Mr+MrM11Q1M11Mr以及S*丫=SyMyM11S1其中，Sr为一个具有协方差内核ESr1S、'=

23、Qr1卜丫2的零均值高斯过程矩阵。定理一中的渐近分布与Hansen（1996）所介绍的一样。一般的，渐近分布不能再进一步的简化。然而，在本节最后我们将讨论一种特殊的简化形式。接下来我们考虑仔需要估计的情况。由于"仔-仔°=0卩1,那么，在仔0的比率为n-1的置信邻域内研究厶M仔,y的变化特征就足够了。定理二：在原假设捡下，LMaa,y=LMB+n,y具有与厶M仔y相同的有限维度渐近分布（fidi's）。另外，如果我们能够说明LMaa,y过程是一个紧凑的集合（？），那么SupLM和SupLM0就具有相同的渐近分布，也就是T。这就表明使用估计值仔而不是真实值仔°

24、;,并没有改变LM检验的原渐近分布。不幸的是，我们无法证明这个想法的正确性。证明时存在双重困难，由于虚拟变量的存在，导致LMa,y在仔上是不连续的，而且LMaa,y也是非平稳变量莓_的函数。仅有很少一部分文献将经验性过程描述成时间序列过程，几乎没有任何一篇文献讨论非平稳数据。再者，非平稳变量基一1出现在虚拟变量函数中，所以Taylor序列方法并不能简化这个问题。我们认为，虽然缺乏完备的证明，定理二的fidi结果足够说明使用渐近分布T计算SupLM统计量的合理性。定理一给出了渐近分布T的一个表达形式。它具有一个随进过程Fr的上确界形式。Tr为一个卡方过程，对每一个r而言，Tr的边缘分布为卡方分布

25、。由于T为这个随机过程的上确界，它的分布形式就由卡方过程的联合分布决定,那么就取决于未知函数Mr和Qr。由于这两个函数可以选取很多种形式，T的临界值无法制成图表。s=sM1。不是一般性，我们标准化M1=1以及Y=/。那么,在一种特殊的情况下，我们可以给出一种重要的简化形式。令模型（2）不具有截距项，且叫没有滞后项，同时满足E珀曾一=2,那么唯一的回归元就是误差修正项叫一1。由于Mr为纯量的单调递增函数，那么就存在函数少s使得M®=Ws就是一个标准布朗运动，S申s=WssW1为布朗桥,而且有T=sup兀0<r<1一兀01supn0<r<1一n0T邰s(Ws-sW

26、12sup=S<s<s2s1s其中，S1=1兀0以及s2=9-11一兀0。这是Andrews（1993）中用来检验未知时间点上结构变异是所用的分布，也是满足s1sS=250s1s的一个函数。（？）汩丄3.3.渐近的p值：固定回归元的bootstrap除去Section3.2最后讨论的特殊情况，定理一和二所描述的渐近分布基于矩函数Mr和Qr，那么临界值就无法制成图表。在本节中，我们将讨论Hansen（1996，2000b）中所使用的固定回归元bootstrap入耳用来计算渐近临界值和p值，得出一阶渐近修正推论。定理二表明，第一阶段对于协整向量仔的估计并不会影响SupLM检验的渐近分布

27、。那么在讨论门限推论时就不需要考虑仔估计方法的影响。然而，由于当仔为需要估计的时候，定理二并不是一个关于SupLM渐近分布的完备证明，我们必须强调这个思路的一些部分只是猜测。接着，我们来描述一个固定回归元的bootstrap。我们令一1（仔以及Xt一1=乂上一1（仔，令叫为Section2中所介绍的降秩回归的残差。对于接下来的讨论，ut、wt1、Xt1以及仔都固定为其样本值。设ebt独立同分布且服从标准正态分布W0,1，令九=utebt。用Xt一1Mybt回归，得到回归残差ubt。用X一1d昇,Y以及X一d2昇,Y对回归，得到回归系数矩阵41yb和码/b，以及残差"眈。定义yb和yb为设定仔=仔，并由屮，Y和§2仔，Y所定义的%代替色之后由式（10）和（11）所计算的值。接着再令SUp/4SUPLM*=y<y<yvecAy厂码y巧Yb+岭-叱41YbA2Hansen（1996）的分析表明，在对叫条件做出局部替代的情况下，SupLM匕T，SupLM*的分布具有对SupLM渐进原分布有效地一阶近似。其中，“巳”表示有Gine和Zinn（1990）定义的依概率弱收敛。SupLM*的分布位置，但是可以是使用模拟的方法计算得出。根据上面所介绍的方法能够得到一个分布形式。根据对于误差项ebt相互独立的模拟结果，能够得出一个新的分布形