ADF检验中滞后长度的选择.docx

1、ADF检验中滞后长度的选择ADF 检验 中滞后长度 的选择基于ARIMA(0,1,q)过程的模拟证据摘要】在进行 ADF检验时如何确定一个最优的滞后长度一直是研究者们关注的冋题。最近 的研究表明，不同的滞后长度选择方法对 ADF检验的统计推断影响很大。本文在已有研究的基础上，模拟了更为一般的 ARIMA(0,1,q)过程，分析了在不同的数据生成过程、检验式以及样本容量下，各种滞后长度选择方法对 ADF检验功效和实际检验水平的影响，最后认为修正的信息准则通常具有较合理的实际检验水平，而从一般到特殊法具有更为稳健的 ADF检验性质。关键词ADF检验滞后长度信息准则修正的信息准则从一般到特殊法Abs

2、tract: The optimal lag len gth in estimati ng Augme nted Dickey-Fuller statistics have bee n concen trated on for years. Previous research in dicated that differe nt leg len gth select ion models affect a lot on the statistical inference of ADF test. Based on all the researches available, this paper

3、 simulates a more gen eral ARIMA(0,1,q) process and an alyzes the in flue nce of lag len gth select ion criterio ns to the size and power of the ADF test with differe nt data gen erat ing processes, ADF regressi ons, and sample sizes. Fin ally, it is proved that the Modified In formati on Criteria a

4、lways shows a more proper size and the General to Special Criteria has more robust properties in ADF test.Keywords: ADF test Lag Len gth In formatio n Criteria Modified In formatio n Criteria Gen eral to Specific一、弓I 言随着时间序列非平稳问题的提出， 单位根检验目前已经成为宏观数据建模前首先要进行的工作。为此，Dickey和Fuller(1979, 1981)1提出了著名的 ADF

5、检验，并推导了当时间序列 yt是ARIMA(p,1,0)过程且满足检验式中滞后差分项长度 k 时ADF检验统计量的极限分布。然而，在实际运用 ADF检验时，真实的p是不知道的，因此需要研究者自己确定 k。总的来说滞后长度的选择方法主要分为两类。一类是经验法 (rule of thumb)。这种方法是研究者任意选择k,或将k表示为样本容量的函数。另外一类就是根据数据来选择 k。这种方法主要有 Akaike(1973)信息准则(Akaike In formation Criteria，以下简写为 AIC)、Schwarz(1978) 信息准则(Schwarz In formation Criter

6、ia，以下简写为 SIC)、Ha nnan 和 Qui nn (1979)信息准则 (Ha nnan and Qui nn In formation Criteria，以下简写为 HQIC)、从一般至 U特殊法贝 U (Ge neral to Special Criteria，以下简写为 GSC)、从特殊到一般法则(Special to Ge neral Criteria，以下简写 为SGC)等。此外，在后来的研究中， Weber(1998)又提出了非自相关法则 (No AutocorrelationCriteria)，即从一个比较简化的模型开始，逐渐增加滞后差分项直到残差不能拒绝非自相关 的

7、原假设。2001年他又提出了一种考虑滞后长度 k在特定区间kmin, kmax内的从特殊到一般法，该方法运用了一系列 F检验，确定的最优滞后长度是使得比其大的直到 kmax的所有滞后差分项对应参数的联合检验均不显著的最小的 k。然而很多学者都指出，ADF检验的结论对滞后长度 k的选择非常敏感。Phillips和Perron(1988)模拟发现当真实数据生成过程为随机游走时， 随着检验式中差分项滞后长度的增加， 会导致ADF检验的功效和水平都降低。另外， Schwert(1989)、Agiakloglou和Newbold(1992)以及 Harris(1992)等也指出不同的滞后长度选择方法对

8、ADF检验的实际水平和功效有明显影响。这就引发了关于不同方法确定滞后长度是否以及如何影响 ADF统计量极限分布的讨论。其实早在ADF检验提出不久，Said和Dickey(1984)就证明了对阶数未知的 ARMA过程检验单位根时，只要检验式中的滞后长度 k满足一定的上界条件和下界条件，仍可以用 ADF统计量来检验原过程中单位根的存在。紧接着， Lewis和Reinsel(1985)提出了一个与 Said和Dickey(1984)下界条件等价的条件，并证明当满足该下界条件和 Said和Dickey(1984)上界条件时检验式中滞后差分项的参数估计量具有一致性和渐近正态性。 Hannan和Deist

9、ler(1988)2则提出了各信息准则确定一个平稳可逆的 ARMA过程滞后长度的若干性质。随后，Ng和Perron(1995)明确解答了哪些滞后长度选择方法满足这些上界与下界条件，以 及运用它们确定滞后长度如何影响 ADF检验统计量极限分布的问题。首先，该文讨论了检验式中滞后长度 k不满足Said和Dickey(1984)或Lewis和Reinsel(1985)下界条件对 ADF检 验统计量极限分布的影响。他们认为这时仍渐近服从标准 DF分布，同时滞后差分项的参数估计量仍具有一致性，但其向真值收敛的速度要小于 (T为样本容量，下同)。接着，Ng和Perron(1995)将滞后长度的选择准则与上

10、述极限分布条件相比较， 证明了在ADF检验中，利用各信息准则确定的滞后长度时不满足下界条件， 但统计量仍服从标准 DF分布。而当运用GSC时，如果我们确定的滞后长度最大值满足上界条件和 Lewis和Reinsel(1985)下界条件，则滞后差分项的参数估计量具有一致性和渐近正态性，可以用 t统计量、F统计量和 Wald统计量检验其显著性。最后通过模拟重点讨论了当数据生成过程为 ARIMA(0,1,1)时各方法确定的滞后长度以及对 ADF检验功效和实际检验水平的影响。类似地，Hall(1994)还从一个纯自相关过程入手，给出了当真实数据生成过程是一个 ARIMA(p,1,0)过程时，ADF统计量

11、服从DF分布应满足的假设条件。并讨论了不同滞后长度 选择准则对ADF统计量极限分布的影响。他认为当运用 AIC、SIC、HQIC以及GSC确定滞后长度时，满足上述条件，因此 ADF统计量仍服从标准 DF分布，而运用SGC时不能满足上述条件，从而 ADF统计量的极限分布发生变化，不再服从标准 DF分布。最后对于不同的ARIMA(p,1,0)过程，模拟了基于各种准则的 ADF检验功效与实际检验水平。此外，随着研究的不断深入，学者们又从一些新的角度对滞后长度选择的问题进行了探讨。比如 Ng 和 Perron(2001)将 Elliott、Rothenberg、和 Stock(1996)3以及 Duf

12、our 和 King(1991)4 提出的局部GLS退势法与Perron和Ng(1996)5提出的修正的单位根检验统计量相结合， 提出了一系列MGLS统计量来检验单位根。在这种检验中，他们首度运用了一系列修正的信 息准则(Modified In formation Criteria ，以下简写为 MIC)来确定滞后长度， 并给出了其局部渐近性质。MIC与一般信息准则的本质区别就在于它考虑到检验式中一阶滞后项参数估计量 的偏差与滞后长度是高度相关的，进而通过加入一个包含一阶滞后项参数估计量的修正项对 信息准则拟和不足的问题进行了一定的校正。 Ng和Perron(2005)又重点探讨了在运用各种信

13、息准则时，可用观测值个数 (即调整的样本容量 卜计算均方误差时的自由度、以及计算惩罚因子(penalty factor)时使用的观测值个数对滞后长度选择的影响。结果表明在有限样本下 AIC与SIC选择的滞后长度对上述三个因素非常敏感。综上所述，已有的研究主要集中在对 ARIMA(p,1,0)和ARIMA(0,1,1)过程进行单位根检验时，各方法确定的滞后长度以及相应的单位根检验的功效与实际水平上。 而对ARIMA(0,1,q)即含有单位根的高阶移动平均过程的研究则比较少。另外，也鲜见 MIC与其他方法比较的 相关研究。针对这些问题，本文对 Hall(1994) , Ng和Perron(1995

14、, 2001)的方法和结论进行扩展，在接下来的部分中用蒙特卡罗模拟的方法在有限样本下研究一个更一般的 ARIMA(0,1,q)过程，对模拟结果中不同滞后期选择方法尤其是 MIC的优劣进行比较，以期找到一种能应用在更一般的数据生成过程中，并使 ADF检验推断更真实可靠的滞后长度选择方法。最后一部分是对全文的 总结，并提出了一些滞后项选择及 ADF检验中需要注意的问题。二、模拟结果根据Hall(1994)，Ng和Perron(1995, 2001)文章中的结论，运用信息准则和 GSC确定滞后长度时，ADF统计量仍服从标准 DF分布。其中运用GSC时滞后差分项以的速度收敛于真值， 从而使ADF检验有

15、一个更优的有限样本性质。 MIC是对通常信息准则的修正。因此本文选取AIC、SIC、MAIC、MSIC以及GSC五种方法来确定 ADF检验式中的滞后长度。重点考 察小样本下当误差项为高阶移动平均过程时基于各准则的 ADF检验功效和实际检验水平的特征，以及MIC与其他方法相比对 ADF检验统计推断的影响和滞后长度选择的异同。 各方法确定滞后长度的原理如下：首先，AIC与SIC具有相似的形式，选择的滞后长度 k满足使(1)式的值最小。其中 AIC准则中CT=2，SIC准则中CT=logT，表示估计方程的误差均方，它往往随着滞后长度的增加 而下降。是ADF检验式中的解释变量个数，它等于滞后差分项个数

16、 k加上常数项以及时间趋势项，会随滞后长度的增加而变大，代表了对过度拟和的惩罚。因此选择 k使(1)最小意味着在较少参数和较小的残差平方和之间做出选择。(1)另外，Ng和Perron(2001)提出了一系列的修正的信息准则即 MIC。其选择的滞后长度是使得目标方程(2)的值最小的k，依据CT的表达式不同 MIC又分别称为MAIC与MSIC。它与一般的信息准则的不同就是增加了一个修正因子，其表达式为：其中是ADF检验式中一阶滞后项的参数估计量。 Ng和Perron(2001)证明会随着ADF检验式中滞后差分项个数 k的增加而减小，尤其当数据生成过程的移动平均部分含有负根时， 这种减小更加明显，因

17、此可以有效地校正一般信息准则拟和不足的问题。GSC则是在ADF检验式中选取r=j+m个滞后差分项，并通过对最后 m个参数(i=1,m)的显著性进行联合检验来完成的，其中 j 0, jmax。该检验的 Wald形式为：其中(5)它代表所有解释变量的方差协方差矩阵，是中右下方 mxm阶的块矩阵。代表该检验式回归函数的误差均方，其中代表回归式的残差。检验规则为：j从最大的取值jmax开始，依次降低其取值直到 (4)式表示的统计量显著。该统计量服从自由度为 m的x2分布。基于显著性水平 a滞后长度k的取值为k = j +1 ,当是统计量所有值中第一个大于临界值的值时。 k = 0 ,当统计量所有值均小

18、于临界值时。为了考察误差项为高阶移动平均过程时 ADF检验中滞后长度的选择问题，我们对形如 (7)式的数据生成过程共10种情况运用上述五种方法选择滞后长度继而进行 ADF检验。其中L是滞后因子，ut是白噪声，yO=O。0 4=0.0;10 种数据生成过程如下： 0 1=0.8, 0 2=0.0, 0 3=0.0, 4=0.5;, 0 2=0.0, 0 3=0.0, B 1=0.5, 0 2=0.0, 0 3=0.0, 0 4=001；=0.8, 0 2=0.0, 0 3=04=0.0; 0 1=0.8, 0 2=0.5,0 3=0.0, 0 4=0.00 1=0.5, 0 2=0.3, 0 3

19、=0.0, 0 4=152,; 0 2=0.3, 0 3=0.0,( 4=0.0;3=0.2,0 1=).8, 0 2=0.5, 0 3=0.0,( 0=0.08, 0-0.5, 0 3=0.0, 0 4=0901 =0.5, 0 2=0.3, 00 4=0.1这10种情况描述了误差项移动平均部分的根在个数、大小、正负等方面的不同情 形。ADF检验的原假设 H0: p =；备择假设H1: p 1ADF检验式如下：(b)为考察不同情形下 ADF检验的功效和实际检验水平，我们对每种数据生成过程分别取 3 =10.95、0.85，用Rats6.2模拟样本容量T=100时基于两种检验式(a)和(b)的

20、上述检验过程以及 T=250时基于检验式(a)的上述检验过程。对每种情况重复 10000次，计算ADF统计量小于临界值的概率，同时记录每次选择的滞后长度， 最后计算滞后长度的均值和标准差。 当真实数据生成过程是单位根过程即 3 = 1时，ADF统计量小于临界值的概率就是犯弃真错误的概率，即实际检验水平。而当真实数据生成过程为平稳过程即 3 1时，ADF统计量小于临界值的概率则是1-犯取伪错误的概率，即检验功效。这里运用 GSC确定滞后长度时取 m=1 ,即计算单个参数的t统计量，显著性水平取 5%，各准则的最大滞后长度取 kmax = 206。根据模拟结果，我们重点比较了各方法在滞后长度选择及

21、其相应的 ADF检验功效和实际检验水平方面的异同，并考察了误差项为高阶移动平均时的各种数据生成过程、 不同检验式以及样本容量对滞后长度选择及 ADF检验统计推断的影响，从而对 Hall(1994) , Ng和Perro n(1995)的结论做了一定的补充。1.不同准则的比较滞后长度及其标准差方面(见表1、3、5)，总得来说，GSC、AIC与MAIC选择的滞后长度 通常要高于SIC与MSIC，前者往往倾向于过度拟和。具体来说，当 DGP中移动平均部分的根为正时，AIC与MAIC、SIC与MSIC选择的平均滞后长度很接近。当移动平均部分含 有负根时，修正的信息准则往往比相应的原准则选择更大的平均滞

22、后长度， 并且这一差距随着根更接近于-1或负根个数增加而更加明显，同时滞后长度的标准差也随之增加。表1滞后长度的均值和标准差检验式 (a) T=1003 ARIMA(p,d,q) 0 1 0 2 0 3 0 4 AIC SIC MAIC MSIC GSC1.0 ARIMA(0,1,1) 0.8 0.0 0.0 0.0 5.54 3.06 2.79 1.18 5.29 2.94 2.75 1.23 8.14 5.381.0 ARIMA(0,1,1) 0.5 0.0 0.0 0.0 2.99 2.72 1.46 0.71 2.81 2.44 1.48 0.82 6.61 6.181.0 ARIMA

23、(0,1,1) -0.5 0.0 0.0 0.0 2.85 2.68 1.34 0.75 3.01 2.55 1.77 1.08 6.71 6.301.0 ARIMA(0,1,1) -0.8 0.0 0.0 0.0 4.62 3.28 2.03 1.34 6.39 3.60 4.50 2.56 7.65 5.731.0 ARIMA(0,1,2) 0.8 0.5 0.0 0.0 7.39 3.14 4.34 1.30 7.06 2.97 4.26 1.41 9.20 4.771.0 ARIMA(0,1,2) 0.5 0.3 0.0 0.0 3.79 2.70 2.18 0.81 3.62 2.4

24、1 2.16 0.85 7.05 5.891.0 ARIMA(0,1,2) -0.5 0.3 0.0 0.0 2.15 2.78 0.44 0.80 2.39 2.62 0.96 1.17 6.25 6.531.0 ARIMA(0,1,2) -0.8 0.5 0.0 0.0 3.96 3.19 1.31 1.33 5.43 3.39 3.47 2.15 7.29 5.871.0 ARIMA(0,1,2) -0.8 -0.5 0.0 0.0 6.38 3.56 3.22 2.02 8.47 3.75 6.66 3.38 8.71 5.141.0 ARIMA(0,1,4) 0.5 0.3 0.2

25、0.1 5.21 3.63 1.63 1.32 4.91 3.43 1.82 1.44 8.30 5.440.95 ARIMA(1,0,1) 0.8 0.0 0.0 0.0 5.57 3.14 2.77 1.17 5.33 3.17 2.76 1.37 8.09 5.360.95 ARIMA(1,0,1) 0.5 0.0 0.0 0.0 2.87 2.57 1.44 0.73 3.06 2.75 1.51 1.08 6.59 6.200.95 ARIMA(1,0,1) -0.5 0.0 0.0 0.0 2.55 2.61 1.08 0.76 3.94 3.20 2.68 1.69 6.29 6

26、.310.95 ARIMA(1,0,1) -0.8 0.0 0.0 0.0 2.80 3.11 0.75 1.06 8.54 4.54 7.17 3.82 6.74 6.290.95 ARIMA(1,0,2) 0.8 0.5 0.0 0.0 7.31 3.09 4.32 1.32 7.11 3.23 4.22 1.62 9.10 4.800.95 ARIMA(1,0,2) 0.5 0.3 0.0 0.0 3.74 2.55 2.18 0.85 3.78 2.83 2.24 0.96 7.08 5.950.95 ARIMA(1,0,2) -0.5 0.3 0.0 0.0 1.76 2.58 0.

27、27 0.67 3.24 3.13 1.92 1.57 6.11 6.660.95 ARIMA(1,0,2) -0.8 0.5 0.0 0.0 2.72 3.03 0.57 1.01 6.96 4.56 5.22 3.39 6.56 6.210.95 ARIMA(1,0,2) -0.8 -0.5 0.0 0.0 3.22 3.69 0.56 1.27 11.11 4.96 10.38 4.77 7.56 5.930.95 ARIMA(1,0,4) 0.5 0.3 0.2 0.1 5.06 3.53 1.61 1.31 5.14 3.58 2.41 1.72 8.17 5.440.85 ARIM

28、A(1,0,1) 0.8 0.0 0.0 0.0 5.45 3.03 2.73 1.18 5.62 3.96 2.73 1.96 8.00 5.360.85 ARIMA(1,0,1) 0.5 0.0 0.0 0.0 2.83 2.59 1.36 0.72 3.61 3.57 1.15 1.66 6.55 6.190.85 ARIMA(1,0,1) -0.5 0.0 0.0 0.0 2.02 2.54 0.58 0.73 5.81 4.57 4.41 3.29 6.17 6.550.85 ARIMA(1,0,1) -0.8 0.0 0.0 0.0 1.10 2.49 0.08 0.35 8.72

29、 6.24 8.03 6.06 5.74 6.760.85 ARIMA(1,0,2) 0.8 0.5 0.0 0.0 7.25 3.03 4.28 1.30 7.36 3.86 4.07 2.21 9.08 4.750.85 ARIMA(1,0,2) 0.5 0.3 0.0 0.0 3.75 2.64 2.11 0.82 4.17 3.61 2.32 1.61 7.01 5.940.85 ARIMA(1,0,2) -0.5 0.3 0.0 0.0 1.36 2.61 0.12 0.46 4.73 4.35 3.25 2.92 5.93 6.740.85 ARIMA(1,0,2) -0.8 0.

30、5 0.0 0.0 2.02 2.79 0.42 0.71 6.61 6.25 5.40 5.65 6.07 6.570.85 ARIMA(1,0,2) -0.8 -0.5 0.0 0.0 1.63 2.61 0.32 0.58 6.90 7.46 6.75 7.37 6.25 6.630.85 ARIMA(1,0,4) 0.5 0.3 0.2 0.1 4.98 3.52 1.54 1.27 5.56 4.04 2.84 2.22 8.14 5.44注：表中所列和分别代表模拟 10000次时各准则确定的(1- 3 L)yt滞后长度的均值及其标准差。对同一数据生成过程而言，基于各种信息准则的实际

31、检验水平由小到大分别为： MAIC、MSIC、AIC、SIC。尤其是当DGP的移动平均部分含有较大负根时， AIC与SIC的检验尺度扭曲非常严重，在样本容量为 100时甚至达到了 50%以上(见表2)。而这时MAIC犯第一类错误的概率都能保持在 10%以下。此时基于GSC的检验尺度扭曲要小于一般的信息准则，但仍大于修正的信息准则，在其他情况下通常介于 AIC与SIC之间。从检验功效来看，基于SIC的检验功效最高，在移动平均部分中含有较大负根时其检验功效非常接近于 1。其次是基于AIC与GSC的检验功效，且前者的略高。基于修正的信息准则的检验功效最低，尤 其是在样本容量为 100时更加明显(见表

32、2)。表2 ADF检验的功效和实际检验水平检验式 (a) T=1003 ARIMA(p,d,q) 0 1 0 2 0 3 0 4 AIC SIC MAIC MSIC GSC1.0 ARIMA(0,1,1) 0.8 0.0 0.0 0.0 0.0566 0.0569 0.0239 0.0211 0.05961.0 ARIMA(0,1,1) 0.5 0.0 0.0 0.0 0.0611 0.0674 0.0313 0.0248 0.06291.0 ARIMA(0,1,1) -0.5 0.0 0.0 0.0 0.1008 0.1580 0.0516 0.0676 0.09811.0 ARIMA(0,1,1) -0.8 0.0 0.0 0.0 0.2436 0.4257 0.