跳到内容
公开可用 发布人:德古意特出版社 2015年5月21日

一般确定性成分的递归调整和改进的协整秩检验

  • 本杰明·博恩 马泰·德米特里斯库 电子邮件徽标

摘要

本文讨论了当序列递归调整为确定性成分时,积分向量自回归的协整秩的检验。为此,在两种不同的互补情况下,分析了去除一般可加确定性成分的递归或自适应过程的渐近性质。当被检验时间序列的随机分量是弱平稳的(平衡误差也是如此)时,递归调整的影响随着样本量的增加而消失。当适当规范化的随机分量弱收敛到具有可积路径的极限连续时间过程时(就像常见的随机趋势一样),递归调整具有永久性的效果,甚至是渐近的:归一化递归调整过程弱收敛到极限过程的递归调整版本。协整秩检验的零极限分布可以用递归调整的布朗运动表示。此外,在经验相关性的情况下,检验了带有递归调整的协整秩检验的有限样本性质:所考虑的确定性成分分别是常数、常数和线性趋势。与似然比检验或广义最小二乘平差检验相比,有限样本在零点下的经验抑制频率有所改善;在备选方案下也发现了改进,随着初始条件的大小增加,似然比检验的性能越来越好。关于等级选择,三种测试程序的简单组合和不同的调整效果最好。

1引言

如何解释数据中确定性成分的问题可能与时间序列计量经济学本身一样古老。将随机成分从确定性成分中分离出来的更流行的方法似乎仍然是使用普通最小二乘法(OLS)进行调整,即使在广义最小二乘法(GLS)可能是有效解决方案的情况下也是如此。

递归删除确定性趋势分量,或递归调整,[1]被引入作为替代方案So and Shin(1999)用于在协方差和自回归参数估计的背景下贬低。反复贬低相当于从每次观察中减去t吨平均值ˉt吨=t吨负极1j个=1t吨j个; 它具有自适应性,例如在自回归模型中,与通常的OLS降阶不同,它具有保持回归因子和扰动乘积的鞅差性质的优点。这样,用递归降级计算的统计数据的偏差(如果不是方差的话)可以减少;看见So and Shin(1999)了解详细信息。偏差减少对于处理可能的集成时间序列非常有用。在这种情况下,递归贬低被用于单位根的Dickey–Fuller测试Shin and So(2001);泰勒(2002)在季节单位根检验的更一般设置中,以递归方式解释线性趋势;另请参见罗德里格斯(2006)用于带有递归趋势删除和库津(2005)针对递归反理性问题。由于减少了自回归参数的估计偏差(在自回归根接近单位的情况下,估计偏差很大),因此可以获得更强大的单位根检验。这由以下人员记录Leybourne、Kim和Newbold(2005)他们表明,使用递归调整的测试与使用局部到单位GLS调整的单位根测试具有可比性(Elliott、Rothenberg和Stock,1996年). 同时,递归过程保持了OLS方法固有的简单性,记住GLS调整需要指定微调参数。但单位根测试并不是递归调整的唯一应用领域:Sul、Phillips和Choi(2005)使用递归降级来减少预白化长期方差估计的自回归估计的偏差。

在多元非平稳模型的背景下,当测试(无)协整或协整秩时,很自然会问递归调整是否比其他调整方法带来优势。

因此,本文的贡献是双重的。首先,它研究了一般确定性分量的递归调整。这是分析带有递归调整的协整秩检验的必要步骤;但是这个主题本身就很有趣,我们在比向量自回归(VAR)过程分析所需的更一般的条件下处理它。具体地,区分了随机分量的两种可能情况。在第一种情况下,假设弱平稳性——例如,通过协整VAR的平衡误差满足的条件。在第二种情况下(非平稳),当使用适当的归一化时,应存在具有可积路径的极限连续时间过程;这适用于VAR的随机趋势。但是非平稳分数积分过程也可以通过这个假设来适应,弱极限是分数布朗运动。在这两种情况下,对递归调整序列的影响是不同的,并由此导出了精确的结果。在平稳情况下,随着时间指数增长到无穷大,效应消失,但在样本开始时仍然不可忽略。在非平稳情况下,极限行为是一个连续的递归调整过程。

第二个贡献是在处理协整VAR的确定性成分时使用递归调整,以测试关于协整秩的假设。由提出的估计集成VAR的协整秩的程序约翰森(1995)事实证明,在指定协整VAR模型以进行后续推断方面,它是一个极有价值的工具。[2]但是,当数据生成过程显示确定性组件时,对于如何从数据中删除它们,没有黄金标准。约翰森(1995)建议保持在似然比(LR)框架内;Saikkonen和Lütkepohl(2000a、2000b)吕特凯波尔和赛科宁(2000)另一方面,建议在应用LR或拉格朗日乘数测试之前,通过GLS调整去除常量或线性趋势。这两种程序都不统一地支配另一种程序。由于LR秩检验可以被视为Dickey–Fuller检验的推广,因此预计递归调整在多元设置中也表现良好。本文研究应用于递归调整数据的LR型秩检验的渐近性质和小样本性质。渐近分布依赖于递归调整的Wiener过程,这一点不足为奇。在小样本中,递归调整被发现在大多数情况下在零和备选方案下的经验拒绝频率方面表现更好;初始条件也起作用。当关注于选择协整秩而不是测试特定的零(例如无协整)时,将测试与递归、OLS和GLS类型调整简单结合似乎是更好的选择。

在进入论文的主要部分之前,让我们介绍一些符号。让粗体符号表示列向量,““表示定义于[0,1]空间被赋予适当的范数,“I(0)”和“I(1)”分别代表0阶和1阶的积分。此外,表示欧几里德向量范数和相应的诱导矩阵范数,以及朗道符号()第页()有他们通常的意思。d日=d日1,,表示×具有元素的对角矩阵在主对角线上。最后,C类是一个常量,其值在不同的情况下可能不同。

2模型和试验

2.1假设

让观测序列根据

[1]t吨=x个t吨+d日t吨,t吨=1,,T型,

哪里x个t吨是一个K(K)-维度集成VAR(第页+1)过程(可能是协整的)和向量d日t吨=d日1t吨,,d日K(K)t吨是的确定函数t吨.

让的纠错表示x个t吨

[2]Δx个t吨=Πx个t吨负极1+j个=1第页A类j个Δx个t吨负极j个+εt吨,

哪里x个0固定和创新εt吨是不相关的。什么时候?Π=0,x个t吨集成,并且,当Π拥有全军衔,x个t吨是静止的。什么时候?第页=第页n个k个Π非零,但小于K(K),矩阵可以分解为Π=αβ,αβ两者兼而有之K(K)×第页矩阵。表示方式A类正交补码w.r.t。K(K)矩阵的A类具有K(K)行和0<第页<K(K)列。以下假设确保了x个t吨I(1)具有协整秩吗第页.

假设1

让特征多项式

det(探测)(1负极z(z))K(K)负极Πz(z)负极j个=1第页A类j个(1负极z(z))z(z)j个=0
在单位圆内没有根K(K)负极第页单位圆上的根,如果第页>0,让矩阵αK(K)负极j个=1第页A类j个β拥有完整的等级.

创新εt吨被用来满足协整文献中的条件标准:

假设2

εt吨要么是具有有限方差的iid序列,要么是具有一致有界四阶矩的弱平稳鞅差分序列。表示方式Σε的协方差矩阵εt吨,并让Σε是正定的.

Granger表示定理意味着x个t吨,

[3]x个t吨=Ξj个=1t吨εj个+第页(1),

其中矩阵¦Β表示为Ξ=βαK(K)负极j个=1第页A类j个β负极1α假设2,紧接着

1T型x个[T型]B类(),

其中K(K)-维布朗运动B类()ΞΣε¦Β并满足第页n个k个ΞΣεΞ=第页n个k个Ξ=K(K)负极第页.

关于模型的随机部分,假设事实上相当标准;看见约翰森(1995)和Lütkepohl(2005). 本文更多关注的是确定性组件。

让每个K(K)确定性组件d日t吨是简单趋势函数的线性组合,所有趋势函数都相同K(K)的元素t吨,

[4]d日t吨==1L(左)c(c)(f)(t吨).

此外,让趋势发挥作用(f)()为人所知。当然,系数可能在K(K)的元素t吨; 特别是,有些可能是0,从而适应某些趋势函数特定于的某些元素的情况t吨.

对于渐近分析,关于趋势函数行为的假设(f)(t吨)作为T型需要增长到无穷大。

假设3

D类T型=d日T型α1,,T型αL(左)对于一些实际常数α,=1,,L(左); 趋势函数(f)(t吨),=1,,L(左),线性独立,满足

(f)(t吨)=D类T型τt吨T型,t吨=1,,T型,
哪里τ()=τ1(),,τL(左)()是上的确定向量函数[0,1],具有有界导数的分段光滑0τ(第页)τ(第页)d日第页具有完整等级0<1.

极限中满秩的条件不等于趋势函数的线性独立性,因此是必需的。该假设允许例如常数(其中α1=0τ1=1),对于多项式趋势(其中α1=k个τ1=k个),或用于平滑过渡(其中τ包含所谓的sigmoid函数,如双曲正切)。分段线性趋势或多项式样条也包含在假设中。假设在这一点上,一个常数总是包含在确定性函数集合中:例如,如果没有非零截距,线性趋势就根本不可能成立。[3]实际上,如果没有线性趋势,二次趋势也不可能成立,依此类推。在前面提到的包含常数之前,我们留下(not)选项对从业者来说,考虑高阶多项式趋势函数中的所有幂:该假设包括限制和不限制的情况,并允许人们在认为适合手头数据集的情况下进行递归调整。

罗德里格斯岛(2013)通过递归调整截距、趋势和中断来分析单位根测试。通过涉及常见虚拟变量的交互项,很容易处理样本中固定位置的结构断裂。根据假设3,break的行为几乎与常量相似α1=0; τ1=1(>λ)τ1=1(<λ),使用1()指示功能,用于某些(相对)中断时间0<λ<1.

在存在中断的情况下,递归调整过程实际上有一个有趣的计算方面:一些代数表明递归调整(如以下小节所述)对于趋势组件及其与虚拟变量的交互,相当于只要所有趋势组件同时出现中断,就可以分别调整中断前和中断后的子样本。这为假设3中的问题提供了一个直接的解决方案,而步骤虚拟模型定义为τ1=1(>λ)违反了关于0τ(第页)τ(第页)d日第页对于<λ。处理该问题的另一种方法是使用定义为τ1=1(<λ).

该假设还意味着断裂时间已知;但这是一个似是而非的猜想,即中断时间可以直接内生化。例如,在单位根测试的上下文中,可以使用inf公司-键入统计信息齐沃特和安德鲁斯(1992); 看见Lütkepohl、Saikkonen和Trenkler(2004)对于协整秩检验的相应讨论。[4]

请注意,假设3并不直接涵盖季节决定论;但是泰勒(2002)库津(2005)指出如何将季节性成分减少到相应的非季节性情况。最后,有些情况下0τ(第页)τ(第页)d日第页不是一个合适的积分,不同的方法适用;参见命题1后的讨论。

2.2递归调整

到目前为止,文献中只考虑了递归贬低和贬低。但对一般加性趋势成分的扩展是直接的。该程序也是基于OLS的,但未知系数不是从整个样本中估计的;更确切地说,每个人都有一个估计t吨基于样本的未知趋势参数1,,t吨然后使用它们来调整观察结果t吨多元时间序列的递归调整被理解为逐项调整。

表示方式C类系数矩阵c(c)=1,,K(K),=1,,L(左)和依据(f)(t吨)向量(f)1(t吨),,(f)L(左)(t吨),一个有矩阵表示法

d日t吨=C类(f)(t吨).

现在让我们将注意力转向递归调整序列的相应过程。正式来说,有T型负极L(左)+1最小二乘回归,各一个t吨=L(左),,T型:

j个=C类ˆ(t吨)(f)(j个)+x个ˆj个(t吨),j个=1,,t吨,

C类ˆ(t吨)=j个=1t吨(f)(j个)(f)(j个)负极1j个=1t吨(f)(j个)j个.

表示方式˜t吨递归调整序列,上面的定义为其提供了

[5]˜t吨=t吨负极j个=1t吨(f)(j个)(f)(j个)负极1j个=1t吨(f)(j个)j个(f)(t吨)

对于t吨=L(左),,T型.在样品开始时,获得˜L(左)=0; 设置˜t吨=0也适用于t吨=1,,L(左)负极1.

为每个样本拟合OLS回归t吨之间L(左)T型可能比简单的OLS调整更昂贵;但简单的封闭形式表达式˜t吨例如,在递归贬低或递归减损的情况下(参见示例3),可以使用Sherman–Morrisson公式简化C类ˆ(t吨)通过更新C类ˆ(t吨负极1).

等式[5]暗示用于t吨L(左)那个

˜t吨=x个t吨负极j个=1t吨(f)(j个)(f)(j个)负极1j个=1t吨(f)(j个)x个j个(f)(t吨),

也,d日˜t吨=0为所有人t吨更一般地说,递归调整与通常的OLS程序一样是线性的:

[6]A类~t吨=A类˜t吨

对于共形矩阵A类,因此递归调整t吨不影响潜在的协整结构。具体来说,如果x个t吨是I(0)线性组合,x个t吨~行为渐近稳定,参见命题1,并且,作为等式的结果[6],也会x个˜t吨另一方面,如果x个t吨是I(1)线性组合,x个˜t吨也是随机趋势。

2.3递归调整测试

协整性质对递归调整的不变性表明,这种调整可以与Johansen的(1995)协整秩检验。即使创新是高斯的,结果也不是LR检验,因为趋势参数没有集中在似然函数之外。

如果没有确定性,LR检验可以解释为差异的典型相关性问题Δt吨和滞后水平t吨负极1(如果考虑短期动态,则进行预美白)。因此,本小节讨论了如何t吨负极1Δt吨可以针对其各自的确定性组件进行调整。

e(电子)c(c)τ表示具有递归调整的LR型检验统计量。要计算它,滞后水平t吨负极1按等式处理[5],提供递归调整的˜t吨负极1.

确定性分量是否只包括截距,例如差分序列Δt吨根本不需要调整。但一般情况并非如此;此外,差分序列的确定性成分,d日t吨负极d日t吨负极1通常不同于级别的确定性组件;请参见下面的引理。

引理1

表示方式S公司点集,其中τ或其导数具有跳跃不连续性。τ˙()=(f)1,,(f)L(左)对于[0,1]S公司,为了方便起见(f)(负极1)=(f)(0)在假设3的条件下

啜饮[0,1]S公司T型D类T型负极1(f)([T型])负极(f)([T型]负极1)负极τ˙()0.
证明

明显且省略.■

不连续点S公司对应于差异中的附加异常值。它们对这里研究的测试统计数据的影响取决于它们的大小:平均值的差异偏移是渐近可忽略的;但是,如果线性趋势的系数出现中断,而平均值没有同时移动,则不会消除不连续性,并且在通过脉冲虚拟调整差异时必须考虑到这一点。然而,这种程度的不连续性是不合理的。所以,考虑到一个常数差的情况,让我们把向量τ˙()不包含零元素或“离群值”。[5]引理的条件排除了诸如t吨第页对于0<第页<1,因为它的导数有极点;排除的原因是技术性的,因为τ可能导致在最高权力下出现不确定的形式,结果不会立即出现。

引理1清楚地表明了Δt吨需要进行适当的调整,但(分段)常量确定性组件的情况除外。在递归去趋势上下文中,罗德里格斯(2006)列出了三种可能的方法。两种是调整水平t吨和滞后水平t吨负极1分别,以及差异的样本平均值Δt吨从调整后的水平中减去。这相当于OLS调整差异Δt吨直接。第三个,由提出泰勒(2002),最终相当于使用不同的调整序列。人们还可以考虑递归地调整不同决定论的差异。无论是哪种方式,差异和滞后水平的调整程序的相互作用留下了不可忽视的条件,如下例所示。

示例1

为简单起见,考虑单变量情况t吨=c(c)0+c(c)1t吨+x个t吨哪里x个t吨是一个综合过程。递归去趋势水平由下式给出

˜t吨=x个t吨+2t吨j个=1t吨x个j个负极6t吨(t吨+1)j个=1t吨j个x个j个,
它们的差异由下式给出
Δ˜t吨=Δx个t吨负极4t吨负极2t吨(t吨+1)x个t吨负极2t吨(t吨负极1)j个=1t吨负极1x个j个+12t吨(t吨2负极1)j个=1t吨负极1j个x个j个.
为了进行比较,递归贬低的差异如下所示
Δ˜t吨=Δx个t吨负极1t吨x个t吨.
使用x个t吨一个综合的过程,紧接着.1t吨x个t吨=第页(t吨负极0.5)以及其他术语,如2t吨(t吨负极1)j个=1t吨负极1x个j个可以显示出具有相同的顺序。现在检查样本的叉积矩˜t吨负极1Δ˜t吨(或Δ˜t吨)如Dickey–Fuller型统计中所示:而由调整程序相互作用产生的术语,例如.t吨=2T型1t吨x个t吨x个˜t吨负极1,可以证明是不可忽略的,确定它们的极限行为并非没有挑战.

我们选择对不同序列进行OLS调整,以确定其自身的确定成分。与其他方案相比,这有两个优点。一方面,相互作用项的复杂性得到了控制;关于限制中的后果,请参见命题3。另一方面,从技术角度来看,它使滞后增加更容易处理(t吨由于不清楚简单的滞后(例如递归调整的差异)是否会导致关键的测试统计,或者对于每个滞后差异是否需要不同的调整方案。

通过以下方式表示Δt吨系列Δt吨带OLS调整,即。

[7]Δt吨=Δt吨负极t吨=2T型Δt吨Δ(f)(t吨)t吨=2T型Δ(f)(t吨)Δ(f)(t吨)负极1Δ(f)(t吨),

哪里Δ(f)(t吨)仅包含与差异相关的确定性组件。尽管T型-一致的是,这一步确实会影响极限分布;见提案3。

在调整其特定确定性成分的水平和差异后,回归Δt吨来自等式[7]和˜t吨负极1来自等式[5]关于滞后z(z)t吨负极1=((Δt吨负极1),,(Δt吨负极第页))获得预白残差第页0t吨第页1t吨使用它们,可以建立常用的矩矩阵,

S公司j个=1T型t吨=第页+1T型第页t吨第页j个t吨,
,j个=0,1。表示方式λ^1,,λ^K(K)实变量方程的降阶解λ
λS公司11负极S公司10S公司00负极1S公司01=0,

即的特征值S公司11负极1S公司10S公司00负极1S公司01。(跟踪)e(电子)c(c)τ统计数据是

[8]e(电子)c(c)τ(第页)=T型=第页+1K(K)k个λ^.

或者,可以使用最大特征值测试,

[9]e(电子)c(c)最大值τ(第页)=λˆ第页+1.

下一节讨论跟踪和最大特征值统计的极限行为。在此过程中,我们分析了递归调整对平稳序列(非零协整秩下的均衡误差)和非平稳序列(随机趋势)的影响。

3主要成果

我们首先检查递归删除确定性趋势组件的效果d日t吨根据观测序列。它取决于随机分量的性质。因此,我们区分了两种不同的情况。首先,我们处理弱平稳过程。其次,适当的归一化过程被认为具有a.s.可积路径的极限连续时间过程。基于这些结果,我们给出了带有递归调整的协整秩检验的渐近性。

3.1平稳性下的递归调整

让我们首先处理平稳性的情况。均衡误差序列属于这一类,在我们的模型中实际上是I(0)。可以建立以下结果来扩展So and Shin(1999); 虽然它的公式化程度略高于要求,但在证明复杂性方面的额外成本可以忽略不计。

提议1

假设3成立并假设w个t吨具有自方差矩阵的零均值弱平稳过程Γ小时这样,就像T型,ϕ(T型)=1T型小时=1T型Γ小时0.然后C类+不依赖t,因此

E类w个˜t吨负极w个t吨2C类ϕ(t吨).
证明:

见附录. ■

命题的假设允许在短记忆下通常理解的内容,因此(即使是间接的)适用于w个t吨=βx个t吨负极1即平衡误差。但是,只要样本均值的方差消失,就可以允许偏离弱平稳性T型,参见。菲利普斯(1987); 然而,这里没有讨论这个主题,就像本小节没有讨论具有无限方差的过程一样。[6]注意,还涵盖了具有长记忆的弱平稳过程;这将允许使用递归调整对顺序进行部分集成的过程d日,只要d日<0.5.

备注1

对于样本开始处的固定t,递推平差方案不一致地过滤随机分量,而只是无偏地过滤;这并不奇怪,因为它只使用了时间t时可用的有限信息.t吨=[T型]对于一些固定的0<<1),有一个w个˜t吨负极w个t吨第页0切比雪夫不等式用一种可能不太严格的符号表示w个˜t吨负极w个t吨=第页ϕ(t吨)均匀(单位:t).

备注2

结果要求在任何情况下w个t吨:即使对于iid序列,如果样本平均值不一致,也不成立。利率w个˜t吨收敛到w个t吨与样本平均值的收敛速度有关w个ˉ,更一般地说,与级数的级数依赖性有关w个t吨.

如果α对某些人来说是负面的1L(左),趋势函数(f)(t吨)消失为T型此类趋势函数通常不会影响基于数据的统计的渐近行为没有解释了这一趋势成分(另请参见示例2)。然而,在小样本中,即使这样的趋势函数也会产生差异,因此建议将其考虑在内。在这种情况下,可能无法一致地估计其系数,但仍会消除干扰趋势参数。假设3中的可逆性要求也可能对α; 例如,如果(f)(t吨)=t吨α,要求α>负极0.5.

对于以下情况(f)(t吨)=t吨α具有α<负极0.5必须采用不同的方法。这是因为0τ(第页)τ(第页)d日第页不是一个真积分。分析此类情况时应使用哪种方法取决于各自数据生成过程的特殊性。例如,假设向量(f)(t吨)具有绝对可和的元素,j个1|(f)(j个)|<=1,,L(左)、和限制t吨1(f)(t吨)(f)(t吨)为了简单起见,是非奇异的。然后,j个=1t吨w个j个(f)(j个)方差一致有界。因此,C类ˆ(t吨)负极C类概率一致有界,并且C类ˆ(t吨)负极C类(f)(t吨)将按要求消失t吨自从(f)(t吨)0.

递归确定性成分在样本开始时并没有“完全”删除,这一事实并不妨碍递归调整数据的使用:例如,协方差矩阵是一致估计的,如下例所示。

示例2

假设流程w个t吨成为协变-遍历,T型负极1小时+1T型w个t吨w个t吨负极小时第页Γ小时。表示方式Γˆ小时具有递归调整的样本自方差,Γˆ小时=1T型t吨=小时+1T型w个˜t吨w个˜t吨负极小时和依据1城市街区规范。然后,

E类(血管内皮细胞(Γ^小时负极1T型t吨=小时+1T型w个t吨w个t吨负极小时)1)
1T型t吨=小时+1T型E类(血管内皮细胞(w个t吨(w个˜t吨负极小时负极w个t吨负极小时))1)
+1T型t吨=小时+1T型E类(血管内皮细胞((w个˜t吨负极w个t吨)w个t吨负极小时)1)
+1T型t吨=小时+1T型E类(血管内皮细胞((w个t吨˜负极w个t吨)(w个˜t吨负极小时负极w个t吨负极小时))1)
C类K(K)2T型t吨=小时+1T型ϕ(t吨负极小时)+t吨=小时+1T型ϕ(t吨)+t吨=小时+1T型ϕ(t吨)ϕ(t吨负极小时).
(为了建立第二个不等式,Cauchy–Schwarz不等式被应用于元素,事实是w个t吨利用了一致有界二阶矩。)ϕ(t吨)0,收敛如下。注意Γˆ小时
Γˆ小时负极Γ小时=第页最大值1T型t吨=1T型ϕ(t吨),E类v(v)e(电子)c(c)Γ小时负极1T型t吨=小时+1T型w个t吨w个t吨负极小时1.

递归调整对(共)积分VAR的共同趋势的影响是截然不同的。我们将在下一小节中使用更通用的设置对其进行研究。

3.2非平稳条件下的递推平差

在非平稳情况下,比如在共同趋势空间中工作时,调整后序列的行为会发生变化w个t吨=βx个t吨也就是说w个˜t吨w个t吨不会渐近消失,即使在样本末尾也不会。为了使术语非平稳性更精确,我们假设适当的归一化和插值过程存在一个极限连续时间过程w个t吨有疑问。所以让我们

[10]d日¦ΑT型负极1w个[T型]V(V)(),[0,1]

对一些人来说K(K)-量纲向量¦ΑT型规范化函数,以便¦ΑT型作为T型,其中K(K)-尺寸加工V(V)()是路径可积的,并且在原点几乎肯定是连续的,0V(V)()=V(V)(0)a.s.为了形成结果,定义递归调整的限制过程,

V(V)˜τ()=V(V)()负极V(V)(0)负极0V(V)(第页)负极V(V)(0)τ(第页)d日第页0τ(第页)τ(第页)d日第页负极1τ()

对于(0,1].减法V(V)(0)无关紧要,因为w个~t吨对于一个恒定的变化是不变的第2.1节常数总是要计算的。同时,这确保了V(V)˜τ()原产地;见命题2的证明。我们现在可以证明适当规范化过程的弱收敛性˜t吨递归调整的过程V(V)˜τ().

提议2

假设3成立,同时方程收敛.[10]. 然后

d日¦ΑT型负极1w个˜[T型]V(V)˜τ()
V(V)˜τ()几乎可以肯定是连续的=0
证明:

见附录. ■

例如,命题的假设w个t吨是一个I(1)过程,弱收敛为多元布朗运动(可能具有奇异协方差矩阵)¦ΑT型=T型0.5,,T型0.5也适用于收敛到分数布朗运动(I型或II型)且具有归一化因子的I(d)过程¦ΑT型取决于分数积分参数d日.流程聚合到α-也允许稳定的Lévy运动,因此不需要w个t吨具有有限方差;自然,¦ΑT型取决于索引α在这种情况下。也包括非零起始值;关键要求是路径的连续性=0.

等式[6]当然也适用于非平稳情况,导致

d日¦ΑT型负极1A类w个˜[T型]A类V(V)˜τ()

和至

d日¦ΑT型负极1A类w个˜[T型]A类V(V)˜τ()

对于某些共形矩阵A类.

示例3

K(K)=1.如果T型负极0.5w个[T型]弱收敛到布朗运动B类()并且递归地移除一个常数(L(左)=1具有(f)1(t吨)=1α1=0)或者是一个恒定的线性趋势(L(左)=2,使用(f)1(t吨)=1α1=0,以及(f)2(t吨)=t吨α2=1),恢复了文献中先前讨论的两个特殊案例。在通常的情况下,它认为

[11]w个˜t吨=w个t吨负极1t吨j个=1t吨w个j个;
在恒常和线性趋势案例中,它认为
[12]w个˜t吨=w个t吨+2t吨j个=1t吨w个j个负极6t吨(t吨+1)j个=1t吨j个w个j个;
w个˜[T型]归一化到递归降阶布朗运动后弱收敛B类˜1(),
B类˜1()={B类()负极10B类(第页)d日第页,>00..,=0,
以及递归的有偏布朗运动B类˜1;(),
B类˜1;()={B类()+20B类(第页)d日第页负极620第页B类(第页)d日第页,>00..,=0,
分别;看见Shin and So(2001),泰勒(2002),或Chang(2002).

3.3限制分布

我们现在可以给出迹和最大特征值检验统计量的渐近分布;请参阅以下建议。正如所料,它们具有多元Dickey–Fuller分布的典型结构,并依赖于假定的确定性趋势项。

提案3

W公司()成为(K(K)负极第页)-独立标准Wiener过程和W公司~τ()相应的递归调整版本,以及d日W公司τ˙()表示调整后的差速器d日W公司()负极01d日W公司τ˙01τ˙τ˙d日第页负极1τ˙()d日哪里τ˙()=d日d日τ()假设1-3和T型,

1.适用于e(电子)c(c)τ来自等式的统计[8]用计算t吨来自等式[1]那个

e(电子)c(c)τ(第页)d日tr公司[01d日W公司τ˙(W公司˜τ)(01W公司˜τ(W公司˜τ)d日)负极101W公司˜τ(d日W公司τ˙)],

2.e(电子)c(c)最大值τ第页来自等式[9]分布收敛于括号内矩阵的最大特征值.

证明:

见附录. ■

测试的一致性第页=第页0反对第页>第页0基于e(电子)c(c)τ第页0e(电子)c(c)最大值τ第页0给出:证明中认为第页非零特征值具有相同的用不可观测值计算它们时的极限x个t吨不适应趋势。基于命题3的证明,可以分析局部备选方案对序列的影响。虽然这里没有讨论这个话题,但从证据中可以明显看出,权力在1/T型null的邻域。

备注3

如果假设的协整排名,第页0,大于实际值e(电子)c(c)τ第页0不是如上所述,而是由K(K)负极第页特征值问题的最小解

|ρ01W公司˜τ(W公司˜τ)d日负极01W公司˜τ(d日W公司τ˙)01d日W公司τ˙(W公司~τ)|=0;
看见约翰森(1995)第158页顶部.

如果差异中的某些确定性成分与协整关系正交,则命题3仍然适用,但确定性成分将被过度指定——与正确指定的模型相比,这种情况下可能会出现局部替代序列的功率损失;看见Doornik、Hendry和Nielsen(1998)Saikkonen和Lütkepohl(1999)对于线性趋势的情况。这可以通过构建只考虑协整关系中确定性的测试来抵消。[7]缺点是这样的测试只是渐近相似;看见尼尔森和拉赫贝克(2000)讨论其含义。此外,正如在Demetrescu、Lütkepohl和Saikkonen(2009年)本文采用了简单的方法,只讨论了相似的测试统计数据。

4与LR和GLS调整试验的比较

本小节介绍了基于递归调整数据的协整秩检验的小样本性质的检验结果。考虑了实证工作中最常见的两种情况:常数、常数和线性趋势。所有模拟均在Matlab公司.[8]我们比较了e(电子)c(c)τ使用原始LR测试进行测试约翰森(1995)以及使用GLS调整的测试(此处用SL表示);看见Saikkonen和Luukkonen(1997)Saikkonen和Lütkepohl(1999)对于平均调整版本,以及Saikkonen和Lütkepohl(2000年b)Lütkepohl、Saikkonen和Trenkler(2001)用于趋势调整版本。

带有递归调整的测试统计数据的计算如所示第2.3节特别是,对测定值进行了如下调整。滞后水平t吨负极1按等式递归调整(元素方向)[11]以及[12]. 差异Δt吨在常数的情况下不需要调整。在线性趋势的情况下,通过减去样本平均值来调整非零平均值。

但在继续检查协整秩检验的小样本性质之前,需要使用临界值。这些是通过使用100000次重复和400个积分点对相关随机泛函进行离散近似得到的。

表1

递归贬低和递归贬低测试的临界值

痕迹测试最大特征值测试
常量恒定趋势(&T)常量恒定趋势(&T)
K(K)负极第页10%5%1%10%5%1%10%5%1%10%5%1%
12.783.956.772.293.305.892.763.896.632.283.285.79
29.4611.2815.327.669.3613.188.5410.2814.136.958.5712.07
19.3721.9527.2615.5617.8722.7414.3616.5120.9912.0213.9718.28
432.7335.8542.3426.2729.0735.2620.2722.6327.4717.3519.6724.28
550.1453.9961.3640.2043.7450.9326.1928.703422.7825.2930.26
671.2775.7184.5957.6161.8070.3032.0634.7840.3128.2330.8736.45
796.14101.2111.278.1782.9792.3538.0740.9546.8533.8336.6742.39
8124.9130.7142.4102.4107.9118.943.9846.9753.3239.4442.4448.39
9158164.6177.1130.3136.3148.650.0253.1759.7945.2048.2954.86
10194.9202215.8161.5168.3181.755.9659.2666.1250.9554.2260.86

跟踪和最大特征值统计的临界值如下所示表1对于K(K)负极第页范围为1至10;但在下面的小样本比较中,只给出了跟踪测试的结果,因为它似乎是应用工作中比较流行的变体。LR测试的临界值取自约翰森(1995)和SL测试的测试特伦克勒(2003).

获得所需的临界值后,现在可以检查e(电子)c(c)τ测试。由于渐近分析侧重于类似的测试统计,因此生成的蒙特卡罗样本没有确定性成分,t吨=x个t吨随机成分生成如下:

[13]x个t吨=ψ第页00K(K)负极第页x个t吨负极1+u个t吨.

对于ψ=1(或第页=0),过程x个t吨是综合的,但不是协整的;对于|ψ|<10<第页<K(K),的满秩线性变换x个t吨与协整秩协整第页。我们生成了具有初始值的序列

[14]x个0=γ/1负极ψ2第页000K(K)负极第页,

哪里γ为了研究初始条件对测试行为的影响,进行了更改。已知初始条件会影响单位根测试的性能;看见缪勒和埃利奥特(2003)Elliott和Müller(2006)此外,初始值预计也会对协整秩检验产生影响。例如,Ahlgren和Juselius(2012)量化LR和SL测试对初始值的不同反应。

这个K(K)流程的要素u个t吨从等式中驱动(共同)集成VAR[13]被认为是稳定的单变量AR(1)过程,

[15]u个t吨=ϕu个t吨负极1+εt吨,

具有|ϕ|<1iid公司创新εt吨; 更准确地说,

[16]εt吨N个0,Σε=第页ΘΘK(K)负极第页.

在方程式中使用VAR(2)数据生成过程[13]–[15]允许直接区分长期和短期属性。通过设置ϕ=0,x个t吨成为一个特定的VAR(1)过程,数据生成过程也由托达(1994)随后在其他一些模拟研究中,探索了协整秩检验的性质;看看其他人Demetrescu、Lütkepohl和Saikkonen(2009年)Toda认为,这种数据生成过程对于研究协整秩检验的性质特别有用,因为其他VAR(1)过程可以通过测试不变的全秩线性变换从中获得。同时,这是一种节省的参数化,可以降低蒙特卡罗实验的复杂性。仅搜索结果K(K)=报告;模拟也使用K(K)=5,但结果相似,我们不在此处报告它们以节省空间。

矩阵Θ是一个第页×(K(K)负极第页)矩阵。这个矩阵要么为零,

Θ={(0.8,0.4)(f)o个第页K(K)=,第页=1,(0.8,0.4)(f)o个第页K(K)=,第页=2

Θ={(负极0.8,0.4)(f)o个第页K(K)=,第页=1,(负极0.8,0.4)(f)o个第页K(K)=,第页=2

参数ψ设置为0.5。我们通过设置γ=0γ=5.生成的样本大小相同T型=100对每种情况进行10000次重复进行模拟。

对于每种情况,基于e(电子)c(c)τ、LR和SL统计数据,以及根据约翰森(1995)显示。

此外,我们报告了每种情况下通过组合程序选择的等级的相对频率,并选择使用三种测试中的每种测试估计的等级的四舍五入平均值作为协整等级的估计值。除了提供协整秩估计(这是向量误差修正模型估计的关键输入)之外,使用组合程序的原因如下。本节中的蒙特卡罗实验表明,对比试验的行为存在一些差异,尤其是在备选方案下的抑制频率方面。但差异取决于数据生成过程,所以事先不知道这三种测试中哪一种性能更好。因此,将这三项测试结合起来似乎是一种更好的策略,例如。哈维、莱伯恩和泰勒(2009)拜耳和汉克(2013)然而,直接组合测试(例如,通过拒绝组合)有一个严重的缺点,因为必须为每种组合方法专门获得临界值。与此相反,如上所述,将所分析VAR的协整秩的估计值组合起来是很简单的。

小样本检查集中在两个病例上。首先,ϕ设置为0,VAR(1)进程的顺序x个t吨被假定为已知。这有助于分离协整秩检验问题和选择适当的滞后次序。第二,ϕ设置为-0.5,无限制VAR的滞后顺序是通过Akaike信息标准(AIC)选择的,因为它更接近于实际应用中可能会执行的操作。

这里的讨论重点是短期动态规范,即。ϕ=负极0.5.[9]案例的结果γ=0在上半部分报告表24.在以下情况下Θ=0递归调整检验对零下的经验抑制频率具有良好的控制;在偏离标称5%的情况下,除一种情况外,经验抑制频率在4.4%和5.9%之间。相反,在大多数情况下,LR测试在null下会更频繁地拒绝。[10]与LR测试相比,SL测试更接近标称5%,尤其是在减损时,但不如e(电子)c(c)τ贬低时进行测试。此外,当降低等级时,LR和SL测试在备选方案下倾向于具有较低的经验拒绝频率,尽管在零下具有较高的经验拒绝次数;特别是SL测试总是拒绝较少,而LR测试这次排名居中。当偏离时,LR测试在某些情况下比e(电子)c(c)τ测试,但这可以归因于零位下较高的经验抑制频率。查看所选列组的相对频率时Θ=0,只有一种情况是e(电子)c(c)τ测试优于其他两个单独测试;组合程序略好于所有三个单一测试。在以下情况下Θ是如上所示的非零矩阵,整体图像会发生变化,但不是本质上的。现在e(电子)c(c)τ测试不再那么清晰了。当关注等级选择时e(电子)c(c)τ对于常量的情况,测试比其他两个单一测试更好,但LR优于e(电子)c(c)τ在常量和趋势的情况下进行测试,至少针对第页=1第页=2.组合程序在选择正确等级方面总是更好。非零的选择Θ似乎并不重要,因为这两种规范的结果非常相似。

现在让我们检查初始值非零的情况,即。γ=5在我们的模拟中。对于零位下的经验抑制频率,图片与以下情况类似γ=0关于备选方案下的经验拒绝频率,SL测试是三个单一测试中受影响最大的测试。在备选方案下,其经验拒绝频率大幅下降,其选择正确协整秩的能力也大幅下降第页=1和第页= 2. 通过构造,当三个测试中的一个测试显示其选择正确协整秩的能力下降时,组合程序会受到影响,因此并不总是优于本规范中的三个单一测试。但如果底层数据生成过程未知,这仍然是一个有效的选择。Meng等人(2013)对于单变量情况,当初始值较大时,与LR测试相比,递归调整过程的能力下降。[11]这尤其适用于Θ=0.何时Θ是一个非零矩阵,差异要小得多,这可以用不同的总功率不同来解释Θ.

总结我们的发现e(电子)c(c)τ在大多数研究案例中,test在null和alternative两种情况下的经验拒绝频率都优于其竞争对手。类似地,估计协整秩的组合过程表现良好。在大多数情况下,LR测试的性能比e(电子)c(c)τ测试,但它似乎对初始值的增加更稳健,至少在功率方面。考虑到初始值的增加,SL测试是受影响最大的测试:在某些规范中,它显示在以下替代方案下,经验拒绝频率急剧下降γ=5简而言之e(电子)c(c)τ在许多情况下,试验在经验拒绝频率方面优于LR和SL试验,包括零拒绝频率和替代拒绝频率,但并非全部。没有一个测试能够统一地支配其他两个测试。组合过程包括选择由所选等级的平均值表示的等级作为协整等级[12]在估计协整秩时,似乎是更好的选择。

表2

相对拒绝频率和所选等级的频率τ=1τ=(1,)

第页0拒绝选择
常量持续和趋势常量持续和趋势
e(电子)c(c)τL(左)S公司L(左)e(电子)c(c)τL(左)S公司L(左)e(电子)c(c)τL(左)S公司L(左)梳子e(电子)c(c)τL(左)S公司L(左)梳子
07.18.87.35.89.3792.991.392.793.494.290.79394.5
第页= 016.57.76.56.15.58.66.55.4
20.50.90.70.50.30.70.40.2
094.390.582.182.984.976.45.79.517.97.217.115.123.615.1
γ=0第页= 115.86.66.14.46.35.588.583.976.287.678.578.671.181.4
25.26.15.24.945.953.3
099.999.999.699.399.898.50.10.10.40.10.70.21.50.4
第页= 219895.584.288.68978.81.94.415.52.910.710.820.110.1
25.95.564.65.45.492.19078.792.28483.673.386.1
07.18.87.35.89.3792.991.392.793.494.290.79394.5
第页= 016.57.76.56.15.58.66.55.4
20.50.90.70.50.30.70.40.2
079.994.426.272.490.629.420.15.673.820.927.69.470.626.8
γ=5第页= 115.16.86.446.75.274.887.522.87468.583.925.969.7
24.56.23.24.73.66.33.33.4
099.710095.998.299.985.70.304.10.31.80.114.41.3
第页= 2183.796.842.679.292.643.6163.254.914.318.97.345.817.4
25.25.75.84.15.3578.59138.181.275.287.336.678.2
表3

相对拒绝频率和所选等级的频率τ=1τ=(1,),Θ=(0.8,0.4)

第页0拒绝选择
常量持续和趋势常量持续和趋势
e(电子)c(c)τL(左)S公司L(左)e(电子)c(c)τL(左)S公司L(左)e(电子)c(c)τL(左)S公司L(左)梳子e(电子)c(c)τL(左)S公司L(左)梳子
07.18.87.35.89.3792.991.392.793.494.290.79394.5
第页= 016.57.76.56.15.58.66.55.4
20.50.90.70.50.30.70.40.2
099.799.997.796.699.897.30.30.12.30.23.40.22.70.8
γ=0第页= 116.37.66.75.37.55.493.492.39194.291.392.391.994.9
25.76.96.15.357.15.24.2
010010099.999.210099.7000.100.800.30.1
第页= 2198.998.287.493.495.885.311.812.51.35.94.214.44.6
25.65.75.64.45.54.793.392.582.394.58990.380.992.4
07.18.87.35.89.3792.991.392.793.494.290.79394.5
第页= 016.57.76.56.15.58.66.55.4
20.50.90.70.50.30.70.40.2
096.199.969.29299.874.33.90.130.92.88.10.125.73.8
γ=5第页= 115.77.86.74.37.85.490.392.163.491.987.79269.492.4
25.27.25.24.94.17.44.73.8
099.510099.297.710096.50.500.80.22.303.50.7
第页= 219298.274.386.296.173.87.51.825.44.511.53.923.88
25.15.75.63.65.84.286.992.569.691.982.690.369.188.7
表4

相对拒绝频率和所选等级的频率τ=1τ=(1,),Θ=(负极0.8,0.4)

第页0拒绝选择
常量持续和趋势常量持续和趋势
e(电子)c(c)τL(左)S公司L(左)e(电子)c(c)τL(左)S公司L(左)e(电子)c(c)τL(左)S公司L(左)梳子e(电子)c(c)τL(左)S公司L(左)梳子
第页= 007.18.87.35.89.3792.991.392.793.494.290.79394.5
16.57.76.56.15.58.66.55.4
20.50.90.70.50.30.70.40.2
γ=0第页= 1099.799.997.596.699.897.50.30.12.50.23.40.22.50.8
16.28.16.85.37.95.693.591.890.994.191.491.891.894.9
25.67.36.15.34.97.35.44.3
第页= 2010010099.999.310099.8000.100.800.20.1
199.198.387.793.396.185.30.91.712.21.25.93.914.54.5
25.25.95.64.85.74.993.992.482.594.588.590.480.692
第页= 007.18.87.35.89.3792.991.392.793.494.290.79394.5
16.57.76.56.15.58.66.55.4
20.50.90.70.50.30.70.40.2
γ=5第页= 1096.399.968.792.299.8743.70.131.32.87.80.2264
15.387.14.38.15.89191.962.79287.891.768.691.9
24.87.25.54.947.55.14
第页= 2010010099.998.710098.70.100.101.301.30.3
186.899.318.183.798.320.913.20.781.812.9151.77814.7
24.86.15.54.46.24.38293.216.983.579.492.119.182.2

5结论

递归贬低和去趋势已成功用于时间序列推断,例如,减少自回归参数估计的偏差或增加单位根检验的威力。接下来是多元的、潜在的非平稳时间序列,我们讨论了在推断协整秩时,这些改进是否可以实现的问题。

因此,本文分析了所谓的e(电子)c(c)τ用于检验协整VAR的协整秩的统计数据。它们是通过在应用测试之前递归调整确定性组件的滞后水平而产生的。

为此,本文还分析了一般趋势分量的递归删除(递归调整)对时间序列的影响。讨论了两种情况:第一,假设序列是弱平稳的(对应于平衡误差),第二,假设序列为非平稳的,并在适当归一化时弱收敛到连续过程(对应于共同趋势)。在第一种情况下,递归调整的效果被证明是渐近消失的,而在第二种情况中,限制过程被证明是以可预测的方式受到影响的。因此,协整秩检验的零极限分布可以用递归调整的布朗运动来给出。

针对常数和趋势的经验相关情况,通过蒙特卡罗模拟发现e(电子)c(c)τ测试程序在有限样本中表现良好,并且在零和备选方案下的经验拒绝频率方面通常优于其竞争对手。然而,没有一项测试比这更好。当关注协整秩的估计时,一个非常简单的程序将e(电子)c(c)τLR和SL测试的表现非常有希望。

附录

命题1的证明

推导w个˜t吨w个t吨,请注意

w个˜t吨=w个t吨负极1t吨j个=1t吨w个j个D类t吨负极1(f)(j个)1t吨j个=1t吨D类t吨负极1(f)(j个)D类t吨负极1(f)(j个)负极1D类t吨负极1(f)(t吨)

具有D类t吨负极1假设3中定义。那就让我们

j个,t吨=D类t吨负极1(f)(j个)1t吨j个=1t吨D类t吨负极1(f)(j个)D类t吨负极1(f)(j个)负极1D类t吨负极1(f)(t吨)

并注意,对于[0,1],

[t吨],t吨τ()01τ()τ()d日负极1τ(1)

作为t吨这样的话

啜饮[0,1][t吨],t吨负极τ()01τ()τ()d日负极1τ(1)0;

上确界序列是收敛的,它必须是有界的。此外,矩阵01τ()τ()d日根据假设3,是可逆的。因此[t吨],t吨收敛为t吨连续函数,定义在紧区间上[0,1],也是有界的。这意味着j个,t吨自身有界(一致在):

啜饮[0,1][t吨],t吨啜饮[0,1][t吨],t吨负极τ()01τ()τ()d日负极1τ(1)
+啜饮[0,1]τ()01τ()τ()d日负极1τ(1)
C类

现在,j个,t吨是标量,所以它可以容纳

C类o个v(v)w个˜t吨负极w个t吨=C类o个v(v)1t吨j个=1t吨w个j个j个,t吨=1t吨2=1t吨j个=1t吨,t吨j个,t吨E类w个w个j个

因此

C类o个v(v)w个˜t吨负极w个˜t吨1t吨2=1t吨j个=1t吨,t吨j个,t吨Γ负极j个啜饮t吨T型,j个t吨j个,t吨21t吨2=1t吨j个=1t吨Γ负极j个.

在命题的假设中,利用r.h.s.上的二重和受和的平方限制的事实,我们得出

C类o个v(v)w个˜t吨负极w个t吨(ϕ(t吨))2C类ϕ(t吨)

对一些人来说C类不依赖于t吨(我们使用的事实是ϕ(t吨)0意味着存在一个有限的t吨0对于其中|ϕ(t吨)|<1为所有人t吨>t吨0,因此(ϕ(t吨))2ϕ(t吨)为所有人t吨>t吨0). 进一步注意,由于柯西-施瓦兹不等式,它认为

E类w个~t吨负极w个t吨2=t吨第页C类o个v(v)w个˜t吨负极w个t吨K(K)C类o个v(v)w个˜t吨负极w个t吨K(K),

从而确定结果。

命题2的证明

我们从建立连续性开始=0属于V(V)˜τ().让D类τ=d日(τ())并表示为ι一个向量。立即写入

0V(V)(第页)负极V(V)(0)τ(第页)d日第页0τ(第页)τ(第页)d日第页负极1τ()
=10V(V)(第页)负极V(V)(0)D类τ负极1τ(第页)d日第页10D类τ负极1τ(第页)D类τ负极1τ(第页)d日第页负极1D类τ负极1τ().

τ持续打开[0,1],D类τ负极1τ(第页)ι作为00<第页<因此,D类τ负极1τ(第页)1=(ι)此外,0D类τ负极1τ(第页)D类τ负极1τ(第页)d日第页对于任何情况都是可逆的>0使用第一个积分元素的中值定理可以看出10D类τ负极1τ(第页)D类τ负极1τ(第页)d日第页远离0。总结一下,

0V(V)(第页)负极V(V)(0)τ(第页)d日第页0τ(第页)τ(第页)d日第页负极1τ()=10V(V)(第页)负极V(V)(0)d日第页

沿着每个路径V(V)().再次应用第一个积分中值定理,它适用于每个元素1k个K(K)那个0<ξk个<这样的话

10V(V)k个(第页)负极V(V)k个(0)ιd日第页=V(V)k个(ξk个)负极V(V)k个(0)ι.

随着V(V)()在原点,我们有0(因此ξk个0)那个

V(V)˜τ()0

几乎可以按要求确定。

为了建立期望的弱收敛性,假设w.l.o.gτ0并考虑功能

F类V(V)={V(V)()负极V(V)(0)负极0V(V)(第页)负极V(V)(0)τ第页d日第页0τ第页τ第页d日第页负极1τ>0,0=0

结果遵循连续映射定理(CMT),如果F类是一个连续函数,即。

啜饮0,1V(V)1负极V(V)20第页e(电子)啜饮0,1F类V(V)1负极F类V(V)20

使用写入V(V)(0)取消

F类V(V)1负极F类V(V)2=V(V)1负极V(V)2
负极0V(V)1负极V(V)2τ第页d日第页0τ第页τ第页d日第页负极1τ.

然后,啜饮0,1F类V(V)1负极F类V(V)2以为界

啜饮0,1V(V)1负极V(V)2
+0啜饮0,1V(V)1负极V(V)2τ第页d日第页0τ第页τ第页d日第页负极1τ,

哪里啜饮0,1V(V)1负极V(V)2是矢量啜饮0,1V(V)1,负极V(V)2,1k个K(K)。对于任何q个,1,第二次对皇家海军的指控在以下情况下消失

啜饮q个,1V(V)1,负极V(V)2,啜饮q个,1V(V)1负极V(V)2啜饮0,1V(V)1负极V(V)20

因此,我们只需表明该属性扩展到0,q个.如果

0ιτ第页d日第页0τ第页τ第页d日第页负极1τ=1

作为0反过来,有界性是从连续性的证明中隐含的F类V(V)为0,结果如下。

我们现在陈述并证明命题3证明所需的两个初步引理。

引理2

1.在假设1-3下,它保持为T型那个啜饮t吨T型Δt吨负极Δx个t吨=啜饮t吨T型Δx个t吨负极Δx个t吨=第页(T型负极0.5).

2.出租v(v)t吨是一些I(1)向量过程,可能是协积分的,这样T型负极0.5v(v)[T型]B类(),[0,1],对于一些布朗运动B类使用合适的协方差矩阵。然后,t吨=第页+1T型v(v)˜t吨负极1Δv(v)t吨负极j个=第页(T型).

引理2的证明

1.在假设1-3下,

T型t吨=2T型Δt吨Δ(f)(t吨)t吨=2T型Δ(f)(t吨)Δ(f)(t吨)负极1Δ(f)(t吨)
01d日V(V)第页τ˙(第页)01τ˙(第页)τ˙(第页)d日第页负极1τ˙()

(请参阅导致方程式的参数[19]在引理3.5的证明中)。τ˙()限定于[0,1],r.h.s.为第页(1)统一,我们有

啜饮t吨T型|Δt吨负极Δx个t吨|=t吨T型啜饮|(t吨=2T型Δt吨Δ(f)(t吨))(t吨=2T型Δ(f)(t吨)Δ(f)(t吨))负极1Δ(f)(t吨)|=第页(T型负极0.5)

根据需要。

2.在无偏向的情况下,对滞后水平和滞后差异的样本交叉积矩的行为有很好的理解;参见示例。菲利普斯和杜劳夫(1986)菲利普斯(1988).具体而言,t吨=第页+1T型v(v)t吨负极1Δv(v)t吨负极j个=第页(T型),1j个第页在涵盖我们假设的相当一般的条件下。对于递归调整的级别,它认为

t吨=第页+1T型v(v)˜t吨负极1Δv(v)t吨负极j个=t吨=第页+1T型v(v)t吨Δv(v)t吨负极j个负极t吨=第页+1T型(j个=1t吨v(v)j个(f)(j个))(j个=1t吨(f)(j个)(f)(j个))负极1(f)(t吨)Δv(v)t吨负极j个.

关于r.h.s.的第一个摘要是第页(1); 利用事实v(v)t吨负极1=第页(T型0.5)再加上与命题1和命题2的证明中使用的类似的论点,我们可以得出结论,r.h.s.上的第二个求和的行为如下T型0.5t吨=第页+1T型τ([T型])Δv(v)t吨负极j个.但是T型负极0.5t吨=第页+1T型τ([T型])Δv(v)t吨负极j个01τ()d日B类,所以

t吨=第页+1T型v(v)t吨负极1Δv(v)t吨负极j个=第页(T型)

根据需要。

引理3

(这是的引理10.3的类似物约翰森1995.)让β是…的正交补β,γˉ=βββ负极1,以及B类T型=γˉ,以保持约翰森的原始符号。还定义

C类o个v(v)|Δx个t吨βx个t吨负极1Δx个t吨负极1,,Δx个t吨负极第页=Σ00Σ0βΣβ0Σββ.

然后保持为T型那个

1S公司00第页Σ00,

2βS公司11β第页Σββ,

三。βS公司10第页Σβ0,

4T型负极1B类T型S公司11B类T型d日01B类˜τB类˜τd日,

5B类T型S公司10αd日01B类˜τd日V(V)τ˙α,

6B类T型S公司11β=第页(1),

哪里V(V)是具有协方差矩阵的K维布朗运动C类o个v(v)V(V)=Σε,d日V(V)τ˙是命题3中定义的微分,以及B类˜τ是递归调整的(K(K)负极第页)-维布朗运动B类对于其中C类o个v(v)B类=γˉΞΣεΞγˉ.

引理3的证明

1.让z(z)t吨负极1=Δx个t吨负极1,,Δx个t吨负极第页。这得益于标准OLS代数

第页0t吨=Δt吨负极1T型t吨=第页+1T型Δt吨z(z)t吨负极11T型t吨=第页+1T型z(z)t吨负极1z(z)t吨负极1负极1z(z)t吨负极1.

召回,啜饮t吨T型Δt吨负极Δx个t吨=第页T型负极0.5因此啜饮t吨T型z(z)t吨负极1负极z(z)t吨负极1=第页T型负极0.5.

然后,

1T型t吨=第页+1T型z(z)t吨负极1z(z)t吨负极1=1T型t吨=第页+1T型z(z)t吨负极1+z(z)t吨负极1负极z(z)t吨负极1z(z)t吨负极1+z(z)t吨负极1负极z(z)t吨负极1
=1T型t吨=第页+1T型z(z)t吨负极1z(z)t吨负极1+1T型t吨=第页+1T型z(z)t吨负极1z(z)t吨负极1负极z(z)t吨负极1
+1T型t吨=第页+1T型z(z)t吨负极1负极z(z)t吨负极1z(z)t吨负极1+1T型t吨=第页+1T型z(z)t吨负极1负极z(z)t吨负极1z(z)t吨负极1负极z(z)t吨负极1

所以,对于一个合适的常数C类,

1T型t吨=第页+1T型z(z)t吨负极1z(z)t吨负极1负极1T型t吨=第页+1T型z(z)t吨负极1z(z)t吨负极1
C类啜饮t吨z(z)t吨负极1负极z(z)t吨负极11T型t吨=第页+1T型z(z)t吨负极1+啜饮t吨z(z)t吨负极1负极z(z)t吨负极12.

给定力矩条件z(z)t吨负极1,1T型t吨=第页+1T型z(z)t吨负极1=第页1所以

1T型t吨=第页+1T型z(z)t吨负极1z(z)t吨负极1=1T型t吨=第页+1T型z(z)t吨负极1z(z)t吨负极1+第页T型负极0.5

同样,

1T型t吨=第页+1T型Δt吨负极1z(z)t吨负极1=1T型t吨=第页+1T型Δx个t吨负极1z(z)t吨负极1+第页T型负极0.5.

因此,使用啜饮t吨T型z(z)t吨负极1负极z(z)t吨负极1=第页T型负极0.5再一次,一个人马上就知道了

啜饮t吨T型第页0t吨负极Δx个t吨负极1T型t吨=第页+1T型Δx个t吨z(z)t吨负极11T型t吨=第页+1T型z(z)t吨负极1z(z)t吨负极1负极1z(z)t吨负极1=第页(T型负极0.5),

第页0t吨=Δx个t吨负极1T型t吨=第页+1T型Δx个t吨z(z)t吨负极11T型t吨=第页+1T型z(z)t吨负极1z(z)t吨负极1负极1z(z)t吨负极1+第页(T型负极0.5)

均匀地。然后我们使用了第页(T型负极0.5)术语

S公司00=1T型Δx个t吨负极1Δx个t吨负极1
负极1T型Δx个t吨负极1z(z)t吨负极11T型z(z)t吨负极1z(z)t吨负极1负极11T型z(z)t吨负极1Δx个t吨负极1+o个第页(1)

这样的话

S公司00第页C类o个v(v)Δx个t吨负极1负极C类o个v(v)Δx个t吨负极1,z(z)t吨负极1C类o个v(v)z(z)t吨负极1负极1C类o个v(v)z(z)t吨负极1,Δx个t吨负极1
=C类o个v(v)Δx个t吨负极1|z(z)t吨负极1=Σ00.

2.我们有

第页1t吨=˜t吨负极1负极1T型t吨=第页+1T型˜t吨负极1z(z)t吨负极1(1T型t吨=第页+1T型z(z)t吨负极1z(z)t吨负极1)负极1z(z)t吨负极1

具有˜t吨负极1=x个˜t吨负极1.在β,x个t吨负极1为I(0),因此假设力矩条件表明1T型t吨=第页+1T型βx个˜t吨负极1=第页(1)根据引理2.1的证明

β第页1t吨=βx个˜t吨负极1负极1T型t吨=第页+1T型βx个˜t吨负极1z(z)t吨负极11T型t吨=第页+1T型z(z)t吨负极1z(z)t吨负极1负极1z(z)t吨负极1+第页(T型负极0.5)

均匀地。由于第页(T型负极0.5)项的样本协方差β第页1t吨等于的样本协方差,直至消失项βx个˜t吨负极1负极1T型t吨=第页+1T型βx个˜t吨负极1z(z)t吨负极11T型t吨=第页+1T型z(z)t吨负极1z(z)t吨负极1负极1z(z)t吨负极1.所以βS公司11β是递归调整的I(0)过程的样本协方差,直到消失项为止,它沿着示例2的行保持不变

βS公司11β第页C类o个v(v)βx个t吨负极1
负极C类o个v(v)βx个t吨负极1,z(z)t吨负极1C类o个v(v)z(z)t吨负极1负极1C类o个v(v)z(z)t吨负极1,x个t吨负极1β
=C类o个v(v)βx个t吨负极1|z(z)t吨负极1
=Σββ.

类似于引理3.1和3.2并省略。

4.沿γ¯,x个t吨负极1是I(1);利用CMT确定

[17]1T型1.5t吨=第页+1T型γ¯'x个˜t吨负极1d日01B类˜τd日,

并再次利用这个事实啜饮t吨T型z(z)t吨负极1负极z(z)t吨负极1=第页(T型负极0.5)得出结论

γ¯'第页1t吨=γ¯'x个~t吨负极1负极1T型t吨=第页+1T型γ¯'x个˜t吨负极1z(z)t吨负极1'1T型t吨=第页+1T型z(z)t吨负极1z(z)t吨负极1'负极1z(z)t吨负极1+第页(1)

均匀地。现在使用引理2.2得出以下结论

1T型t吨=第页+1T型γ¯'x个˜t吨负极1z(z)t吨负极1'=第页(1),

导致

啜饮t吨T型γ¯'第页1t吨负极γ¯'x个˜t吨负极1=第页(1),

以及

[18]啜饮t吨T型1T型γ¯'第页1t吨1T型γ¯'x个˜t吨负极1=第页(T型负极0.5).

限制T型负极1B类T型'S公司11B类T型以下是等式[18]、2号提案和CMT。

5.表示方式x个t吨系列x个t吨通过以下方式调整了差异确定性成分的通常方式第页1t吨投影的(不可行)残差x个t吨在第一个第页滞后时间Δt吨,并让S公司11=1T型第页+1T型第页1t吨第页1t吨'。然后,使用αα'=0,

B类T型'S公司10α=B类T型'S公司10负极S公司11βα'α
=γ¯'1T型第页+1T型第页1t吨第页0t吨负极αβ'第页1t吨'α.

现在,啜饮t吨T型||(第页0t吨负极αβ第页1t吨)负极εt吨||=第页(T型负极0.5)自从第页0t吨αβ第页1t吨已针对确定性进行了调整以同样的方式; 请参阅引理2.1的证明。使用的参数约翰森(1995)第148页第二次显示,如下所示o个第页(1)项不影响B类T型'S公司10α,

B类T型'S公司10α=1T型t吨=第页+1T型γ¯'x个˜t吨负极1εt吨'α+o个第页(1),

哪里εt吨=εt吨负极t吨=第页+1T型εt吨Δ(f)(t吨)t吨=第页+1T型Δ(f)(t吨)Δ(f)(t吨)负极1Δ(f)(t吨),其中Δ(f)(t吨)=(f)(t吨)负极(f)(t吨负极1)因此,1T型t吨=第页+1T型γ¯x个˜t吨负极1εt吨可以写成两个词的差,

1T型t吨=第页+1T型γ¯x个˜t吨负极1εt吨

1T型1.5t吨=第页+1T型γ¯x个~t吨负极1T型D类T型负极1Δ(f)(t吨)1T型t吨=2T型T型D类T型负极1Δ(f)(t吨)T型D类T型负极1Δ(f)(t吨)负极1
×1T型t吨=2T型T型D类T型负极1Δ(f)(t吨)εt吨.

对于第一项,我们利用了Kurtz和Protter(1991),其条件与εt吨从假设2(特别是的鞅差性质γ¯x个~t吨负极1εt吨和等式[17])导致

1T型t吨=第页+1T型γ¯x个˜t吨负极1εt吨d日01B类˜τd日V(V).

对于第二项,使用等式[17]、CMT和引理1

[19]01B类˜τ01d日V(V)第页τ˙(第页)01τ˙(第页)τ˙(第页)d日第页负极1τ˙()d日

根据需要。

6.注意B类T型S公司11β是一个递归调整的I(1)过程与一个递归调节的I(0)过程的样本叉积矩。行为的推导γ¯第页1t吨示例2表明,I(0)过程的递归调整只有一个第页(1)对…的影响B类T型S公司11β,和引理2.2可用于获得所需的数量级。

命题3的证明

仅对迹统计量证明了该命题;最大特征值检验统计量的结果使用了相同的参数。

定理11.1证明的论据约翰森(1995)以一对一的方式应用于我们的情况:唯一的区别在于样本叉积矩矩阵的行为S公司00,S公司01、和S公司11在递归调整下(更准确地说,S公司00,βS公司11β,βS公司10,T型负极1B类T型S公司11B类T型,B类T型S公司10α、和B类T型S公司11). 所需结果在引理3中导出;与的引理10.3相比约翰森(1995).

根据导致方程式(11.16)的论点约翰森(1995),最大的第页特征值λˆ1,,λˆ第页被视为概率收敛于行列式方程的解

λΣββ负极Σβ0Σ00负极1Σ0β=0,

和其余的K(K)负极第页特征值以概率收敛到0。参考文献第159页的论点,以及引理3的初步结果,建立了e(电子)c(c)τ统计,从而完成证明。

确认

作者感谢克里斯托夫·汉克(Christoph Hanck)、赫尔穆特·吕特科波尔(Helmut Lütkepohl)、米查利斯·斯塔马托吉安尼斯(Michalis Stamatogiannis)和卡斯滕·特伦克勒(Carsten Trenkler)提出的非常有益的意见和建议。

工具书类

阿格伦,N。、和M。尤塞利乌斯.2012. “协整秩和初始条件的检验.”实证经济学42():66791.10.1007/s00181-010-0442-z在谷歌学者中搜索

拜耳公司,C、。、和C、。汉克.2013. “组合非整合测试.”时间序列分析杂志34(1):8395.10.1111/j.1467-9892.2012.00814.x号在谷歌学者中搜索

,.2002. “具有横截面相关性面板的非线性IV单位根检验.”计量经济学杂志110(2):26192.10.1016/S0304-4076(02)00095-7在谷歌学者中搜索

,十、。、和邮政编码:。菲利普斯.2009. “半参数协整秩选择.”计量经济学杂志12(供应S1):第83节104.10.1111/j.1368-423X.2008.00270.x在谷歌学者中搜索

德梅特雷斯库,M。,H。吕特凯波尔、和第页。赛科宁.2009. “具有不确定确定性趋势项的向量自回归过程的协整秩检验.”计量经济学杂志12():41435.10.1111/j.1368-423X.2009.00297.x在谷歌学者中搜索

多尔尼克,J.A.公司。,D.F.公司。亨德利、和B。尼尔森.1998. “协整模型中的推断:英国M1重新审视.”经济调查杂志12(5):53372.10.1111/1467-6419.00067在谷歌学者中搜索

埃利奥特,G.公司。、和英国。米勒.2006. “最小化初始条件对单位根测试的影响.”计量经济学杂志135:285310.2016年10月10日/j.jeconom.2005.07.024在谷歌学者中搜索

埃利奥特,G.公司。,T·J。罗森伯格、和J.H。股票.1996. “自回归单位根的有效检验.”计量经济学64(4):81336.10.2307/2171846在谷歌学者中搜索

哈维,D.I.公司。,S.J.公司。莱伯恩、和A.M.R.公司。泰勒.2009. “实践中的单位根检验:处理趋势和初始条件的不确定性.”计量经济学理论25():587636.10.1017/S026646660809018X号在谷歌学者中搜索

约翰森,S公司.1995.协整向量自回归模型中基于似然的推断.牛津:牛津大学出版社.10.1093/0198774508.001.0001在谷歌学者中搜索

库尔茨,T.G.公司。、和第页。普罗特.1991. “随机积分和随机微分方程的弱极限定理.”概率年报19():103570.10.1214/aop/1176990334在谷歌学者中搜索

库津,V(V).2005. “递归去意义与确定性季节性.”统计与概率信件72():195204.2016年10月10日/j.spl.2004.11.025在谷歌学者中搜索

莱伯恩,美国。,总高度。基姆、和第页。纽黑德.2005. “对Dickey–Fuller测试的一些更有力的修改的检验.”时间序列分析杂志26():35569.10.1111/j.1467-9892004.00406.x在谷歌学者中搜索

吕特凯波尔,H(H).2005.多时间序列分析新导论.柏林:施普林格.10.1007/978-3-540-27752-1在谷歌学者中搜索

吕特凯波尔,H。、和第页。赛科宁.2000. “具有时间趋势的VAR过程协整秩的检验.”计量经济学杂志95(1):17798.10.1016/S0304-4076(99)00035-4在谷歌学者中搜索

吕特凯波尔,H。,第页。赛科宁、和C、。特伦克勒.2001. “VAR过程协整秩的最大特征值与迹检验.”计量经济学杂志4(2):8.10.1111/1368-423X.00068在谷歌学者中搜索

吕特凯波尔,H。,第页。赛科宁、和C、。特伦克勒.2004. “未知时间水平位移VAR过程协整秩的检验.”计量经济学72(2):64762.10.1111/j.1468-0262.2004.00505.x号在谷歌学者中搜索

,M。,H。,M.H.先生。Cho公司、和J。.2013. “初始观测对递归去均值和去趋势法单位根检验的影响.”经济学快报120(2):19599.2016年10月10日/j.econlet.2013.03.033在谷歌学者中搜索

米勒,英国。、和G.公司。埃利奥特.2003. “单位根和初始条件的测试.”计量经济学71(4):126986.10.1111/1468-0262.00447在谷歌学者中搜索

尼尔森,B。、和答:。拉贝克.2000. “协整分析中的相似性问题.”牛津经济统计公报62(1):522.10.1111/1468-0084.00157在谷歌学者中搜索

菲利普斯,邮政编码:B.1987. “单位根时间序列回归.”计量经济学55(2):277301.10.2307/1913237在谷歌学者中搜索

菲利普斯,邮政编码:B.1988. “基于鞅逼近的样本协方差矩阵到随机积分的弱收敛性.”计量经济学理论4():52833.10.1017/S026646660001344X号在谷歌学者中搜索

菲利普斯,邮政编码:。、和序号。杜劳夫.1986. “具有集成过程的多时间序列回归.”经济研究综述53(4):47395.10.2307/2297602在谷歌学者中搜索

波斯基特,D.S.公司。、和H。吕特凯波尔.1995. “向量自回归模型中协整关系数的一致性估计.“在计量经济学理论与实践。Hans Schneewei的节日,编辑人R。加拉塔H。库申霍夫,87100.柏林:物理-Verlag.在谷歌学者中搜索

罗德里格斯岛,下午.2006. “递归趋势调整单位根检验的性质.”经济学快报91():41319.2016年10月10日/j.econlet.2005.12.026在谷歌学者中搜索

罗德里格斯岛,下午.2013. “递归调整、单位根检验和结构断裂.”时间序列分析杂志34(1):6282.10.1111/j.1467-9892.2012.00813.x号在谷歌学者中搜索

赛科宁,第页。、和R。卢克科宁.1997. “无穷阶向量自回归过程中的协整检验.”计量经济学杂志81(1):93126.10.1016/S0304-4076(97)00036-5在谷歌学者中搜索

赛科宁,第页。、和H。吕特凯波尔.1999. “VAR过程协整秩的局部似然比检验.”计量经济学理论15(1):5078.10.1017/S026646669915103X号在谷歌学者中搜索

赛科宁,第页。、和H。吕特凯波尔.2000年. “带截距VAR过程协整秩的检验.”计量经济学理论16():373406.10.1017/S0266466600163042在谷歌学者中搜索

赛科宁,第页。、和H。吕特凯波尔.2000亿. “向量自回归过程协整秩检验前的趋势调整.”时间序列分析杂志21(4):43556.10.1111/1467-9892.00192在谷歌学者中搜索

小腿,D.W.公司。、和学士学位。所以.2001. “单位根检验的递推平均平差.”时间序列分析杂志22(5):595612.10.1111/1467-9892.00243在谷歌学者中搜索

所以,学士学位。、和D.W.公司。小腿.1999. “时间序列推断中的递归平均调整.”统计与概率信件43(1):6573.10.1016/S0167-7152(98)00247-8在谷歌学者中搜索

南部,D。,邮政编码:。菲利普斯、和C.Y.公司。.2005. “HAC估计中的预白化偏差.”牛津经济与统计公报67(4):51746.10.1111/j.1468-0084.2005.00130.x在谷歌学者中搜索

泰勒,上午/下午.2002. “季节性和非季节性时间序列基于回归的递归平均调整单位根检验.”商业与经济统计杂志20(2):26981.10.1198/073500102317352001在谷歌学者中搜索

托达,H·Y.1994. “存在线性趋势时协整秩似然比检验的有限样本性质.”经济学与统计学综述76(1):6679.10.2307/2109827在谷歌学者中搜索

特伦克勒,C类.2003. “一组新的确定项事先调整的系统协整检验临界值.”经济学公告(11):19.在谷歌学者中搜索

齐沃特,E.公司。、和D.W.K.博士。安德鲁斯.1992. “关于大崩盘、油价冲击和单位根假设的进一步证据.”商业与经济统计杂志10():25170.在谷歌学者中搜索

在线发布:2015-5-21
印刷出版:2015-7-1

©2015 De Gruyter版权所有

于2014年6月12日从下载https://www.degruyter.com/document/doi/10.1515/jtse-2013-0005.html
滚动到顶部按钮