y35 •Εz

 

Cross-Sectional Effects of Common and Heterogeneous Regressors on Asymptotic Properties of Panel Autoregressive Unit Root Tests

 

Katsuto Tanaka
Faculty of Economics, Gakushuin University, Tokyo, Japan

 

 

The present paper deals with nonstationary panel autoregressive iARj models, and examines cross-sectional effects of regressors on the asymptotic properties of panel unit root tests for the ARi1j coefficient. We consider various types of common and heterogeneous regressors and compute limiting local powers of tests as T ¨‡ for each N, where T and N are the time and cross section dimensions, respectively. Dealing with tests based on the ordinary least squares estimator iOLSEj and the generalized LSEiGLSEj, we examine how common and heterogeneous regressors affect the tests as N becomes large. It is shown that the existence of common regressors does not affect the tests asymptotically as N ¨‡. This means that the power of the tests remains the same even if the model contains common regressors. We further derive the limiting power envelopes of the most powerful invariant iMPIj tests, which yields the conclusion that the GLSE-based tests are asymptotically efficient, unlike the time series case.

 

Keywords Asymptotically efficient test, Common regressor, Cross-sectional effect, Heterogeneous regressor, Moment generating function, Numerical integration, Panel unit root tests.

 

Address correspondence to Katsuto Tanaka:

 

Faculty of Economics
Gakushuin University
Mejiro, Toshima-ku, Tokyo 171-8588
Japan
e-mail:

 

1. INTRODUCTION

 

Nonstationary panel AR models were extensively discussed in Moon and Perron i2008j and Moon, Perron, and Phillipsi2007j, where the former deals with the case of heterogeneous intercepts, whereas the latter discusses the case of heterogeneous trends. In these papers, the limiting local powers of various panel AR unit root tests are computed as T and N jointly tend to ‡ under the local alternative that shrinks to the null at the rate of 1/iTNƒΘj, where T is the time series dimension and N is the cross section y36 •Εz dimension with 0 ƒ ƒΘ ƒ 1.
@@Unlike the above works, the present paper examines the effect of the cross section dimension N on the unit root tests as T ¨‡. This may be useful when T is bigger than N and it is desirable to see the intermediate situation rather than the final situation as both T and N go to ‡. We consider four types of regressors:i1j a common intercept and trend,i2j heterogeneous intercepts and a common trend,i3j a common intercept and heterogeneous trends, i4j heterogeneous intercepts and trends. For these models we conduct panel unit root tests based on the OLSE and GLSE of the AR coefficient, and some other tests based on these residuals. To see the cross-sectional effect, we compute limiting local powers of these unit root tests as T ¨ ‡ for each intermediate N under the AR coefficient close to unity in the order of 1/T. It is theoretically and graphically shown that, as N becomes large, the existence of common regressors does not affect the asymptotic properties of these tests, although that of heterogeneous regressors does affect. This fact was also partly observed in panel AR models discussed in Breitung i2000j and Moon et al.i2007j. We give more detailed analysis of this fact for each intermediate value of N. We also derive the limiting powers of these tests and envelopes of the most powerful invariant iMPIj tests as N ¨‡, utilizing the joint moment generating functions im.g.f.sj associated with the test statistics obtained in Nabeya and Tanakai1990j, Tanakai1996, Chap. 7j, and Tanakai2017, Chap. 10j.
@@The outline of the paper is as follows. In Section 2 we present panel AR models to be dealt with in this paper. In Section 3 we compute limiting local powers of various unit root tests. In Section 3.1 we deal with OLSE-based tests, followed by GLSE-based tests in Section 3.2. The limiting power envelopes are derived in Section 3.3, and it is found that the GLSE-based tests are asymptotically efficient, unlike the time series case. The effect of temporal or cross-sectional dependence of the error term on the tests is discussed in Section 3.4. Section 4 concludes the paper. Proofs of theorems are provided in the Appendix.

 

2. PANEL AR MODELS

 

The panel AR models to be discussed in this paper are the following types:

where i refers to cross section, whereas t refers to time series. The processoƒΕit p is defined for all models by

y37 •Εz

where it is assumed thatoƒΕit p starts from ƒΕi0 0 for each i, and is driven byoƒΓit p. We initially assumeoƒΓit p` i.i.d.i0, ƒΠ2 j for simplicity of presentation. The case of temporal or cross-sectional dependence will be discussed in Section 3.4.
Model A is the most restricted model with common intercept and trend. Model B has heterogeneous intercepts, whereas Model C has heterogeneous trends. Model D is the most unrestricted model with heterogeneous intercepts and trends. Note that these four models coincide with each other when N1.
@@For the above models we consider the panel AR unit root test

where we assume that, under H1, ƒΟi takes the following form:

with c „ 0 and 0 ƒ ƒΘ ƒ 1. This is a simple extension of the time series unit root test. A more general alternative allows the true value of ƒΟi to be different among cross sections. Moon and Perroni2008j and Moon, Perron, and Phillips i2007j assume such an alternative, but we maintain i7j to simplify subsequent discussions.
@@Under the above setting we shall explore asymptotic properties of various unit root tests. For this purpose we define the Ornstein-UhlenbeckiO-Uj process by

where r Έ m0, 1n and oWi irjp is the standard Brownian motion independent of oWkirjpii‚kj so that Y1 irj, . . . , YN irj are i.i.d. for any r Έm0, 1n.

 

3. LIMITING POWERS AND POWER ENVELOPES

 

We first compute the limiting local powers of various unit root tests for Models A through D. In Section 3.1 we deal with OLSE-based tests, followed by GLSE-based tests in Section 3.2. The limiting power envelopes of the MPI tests are derived in Section 3.3. The effect of temporal or cross-sectional dependence of the error term is discussed in Section 3.4.

 

3.1. OLSE-Based Tests

The present test was earlier considered in Moon et al. i2007j, and Moon and Perron i2008j. The limiting local power was also computed in these works as both T and N go to ‡ under a more general setting. Here we examine the cross-sectional effect of regressors as T ¨‡ for each N.

y38 •Εz

Let itiM jbe the OLS residual obtained from Model M iMA, B, C, Dj. Then we compute the estimator iM j of ƒΟi ƒΟ under H0 by

where

@The following theorem describes the asymptotic distribution of iM j as T ¨ ‡ for each N, the proof of which is given in the Appendix.

 

Theorem 1. As T ¨‡ with N fixed under ƒΟi 1|cN / T , the asymptotic distribution of iM j in Model M iMA, B, C, Dj follows

where

y39 •Εz

@Some remarks follow.

 

iaj When N1, that is, in the time series case, the distribution of QNiM j in i12j reduces to Q1iD jU1iD j/V1iD j for all M. Note also that U1iA j/V1iA j corresponds to the popular near-unit root distribution associated with the time series model xtƒΟxt|1{ƒΓt with ƒΟ1|c/T .

 

ibj As N becomes large, it holds that

where the distribution of this last quantity is obtained from Model M without common regressors, which means that the effect of common regressors fades away as N becomes large.

 

icj We can deal with some other variations of the above models, for which we can also consider the statistics T iiEj|1j. For example, we can show that, as T ¨‡ with N fixed under ƒΟi 1| cN /T ,

Thus we also conclude that, for these models, the existence of common regressors does not affect the asymptotic behavior of the OLSE-based tests as N ¨‡, which was also described inibj.

 

idj Ui iC j and Vi iC j behave differently from the other quantities, which may be because Ui iC j/Vi iC j results from the restricted regression without intercept yit ƒΐi t{ƒΕit . It can also be shown that Ui iM j and Vi iM j are uncorrelated under ƒΟ1 for MB, D, but are correlated for MA, C. In fact, it holds that CoviUi iA j, Vi iA jj1/3 and CoviUi iC j, Vi iC jj1/175 when cN0. These can be computed easily from the joint moment generating functionim.g.f.j described below.

 

To compute the distribution of QNiM j in i12j for each N, we use the joint m.g.f.miM jix, y j of Ui iM j and Vi iM j defined by

y40 •Εz

where, by putting ƒΚ , we havemTanakai2017, Chap. 10jn

Then the distribution of QNiM j can be computed by using Imhoffs formulamImhof i1961jn

Numerical computation like Simpsonfs formula can be used to compute i14j by taking care of the computation of square roots of complex-valued quantitiesmTanakai1996jn.
@@Figure 1 draws the probability densities of QNiB j and QNiD j for various values of N under H0icN0jto examine the cross-sectional effect of N. As Theorem 2 below indicates, these distributions converge to |3 and |15/2, respectively, as N becomes large. Note that Q1iB jQ1iD j. The distributions QNiB j for N „1 are shifted from Q1iD j, whereas QNiD j for N „ 1 are just the convolution since P iQNiD jz jPiizVi iD j|Ui iD jj≥0j, as is seen fromi14j. The general feature of QNiA j and QNiC j are the same as QNiB j, although those densities are not presented here.

 

@@We next compute limiting powers of the tests based on QNiM j as N ¨ ‡ under ƒΟ1|cN with cN y41 •Εz c/NƒΘ. We need to find the limiting distribution of normalized QNiM j by suitably choosing ƒΘ. For this purpose, let us put

The joint m.g.f. miM jix, y j of Ui iM j and Vi iM j shown above can be used to compute these moments using the Taylor expansion, as is shown in the Appendix. We have, by the week law of large numbersiWLLNjand the central limit theoremiCLTj,

 

 

y42 •Εz

where

Then it is recognized that, for the asymptotic distribution of normalized QNiM j to be nondegenerate, cNOi1/j when a1b0|a0b1‚0, and cNOi1/N 1/4j when a1b0|a0b10 and a2b0|a0b2‚0. It is shown in the Appendix that cNOi1/j for QNiA j and QNiB j , whereas cNOi1/N 1/4j for QNiC j and QNiD j, and we have the the following theorem.

 

Theorem 2. The limiting powers of the tests based on QNiM jiMA, B, C, Dj under ƒΟ1| c /iNƒΘT j at the 100ƒΑ% level are given as follows:

where ƒ³iEj is the distribution function of Ni0, 1j, and zƒΑ is the 100ƒΑ% point of Ni0, 1j, whereas ƒΘ1/2 for Models A and B, and ƒΘ1/4 for Models C and D.

 

@@It follows that the OLSE-based unit root tests in Models A and B have nontrivial powers in a N |1/2T |1 neighborhood of unity, whereas the powers for Models C and D are nontrivial in a N |1/4T |1 neighborhood of unity. It is also seen that the limiting power decreases as the model complexity increases.
Figure 2 shows powers of the QNiB j - and QNiD j -tests against ccN NƒΘΈm0, 20n at the 5% level for N1, 10, 100, ‡, where the powers for N ƒ ‡ are obtained fromi14j by putting z at the 5% point of the null distribution of QNiM j, whereas those for N‡ are obtained from Theorem 2. It seems that the powers for N100 are still not well approximated by the limiting powers. This is particularly true of Model D. The powers of the QNiB j -test are higher than those of the QNiD j -test, which is also evident from Theorem 2. This means that the existence of heterogeneous trends decreases the power.

 

@@In the next subsection we consider the GLSE-based tests, which will be shown to be better than the

y43 •Εz

 

OLSE-based tests.

 

3.2.GLSE-Based Tests

Let us express Model M iMA, B, C, Dj as y X iM jƒΑiM j{ƒΕ, where X iM j and ƒΑiM j are the regression matrix and parameter vector in Model M, respectively, whereas

y44 •Εz

Then we define the GLS residual by iM jy | X iM jiM j, where

Here is the Kronecker product and C is the T ~ T lower triangular matrix with is, t j-th element being 1 for st and 0 otherwise. The GLSE iM j of ƒΟ can be computed followingi9j with iM j replaced by iM j.
@@The following theorem describes the asymptotic distribution of iM j as T ¨ ‡ for each N, the proof of which is given in the Appendix.

 

Theorem 3. As T ¨‡ with N fixed under ƒΟ1| cN /T , the asymptotic distribution of iM j for Model M iM  A, B, C, Dj follows

where

@@It is noticed that the distributional structure of the GLSE-based statistics RNiM j remains the same as that of the OLSE-based statistics QNiM j. It is also seen that RNiA j coincides with RNiB j. The same is true of RNiC j and RNiD j, and these properties are also shared in the time series casemTanakai1996jn. The densities of RNiA ji  RNiB jj under H0 are drawn at the top of Figure 3 for N1, 10, 30, whereas those of RNiC jiRNiD jj at the bottom for N1, 10, 50. The former densities are seen to be shifted from the latter as N becomes large. Both RNiA j and RNiB j converge to 0, whereas both RNiC j and RNiD j converge to |3, as is described in Theorem 4 below.

 

@@We next consider limiting powers of the tests based on RNiM j as N ¨ ‡, which is described in the following theorem, the proof of which is given in the Appendix.

 

Theorem 4. The limiting powers of the tests based on RNiM jiMA, B, C, Dj as N ¨‡ under ƒΟ1|c /iNƒΘT j at the 100ƒΑ% level are given as follows:

y45 •Εz

where ƒΘ1/2 for Models A and B, and ƒΘ1/4 for Models C and D.

 

@@It follows from Theorem 4 that RNiM j converges to 0 for M  A, B, whereas it converges to |3 for MC, D, as was mentioned before. It is also noticed from Theorems 4 and 2 that the GLSE-based tests are better than the OLSE-based tests in Models B, C, D, although those are the same in Model A. The top of Figure 4 shows powers of RNiA jiRNiB jj-tests at the 5% level for N1, 10, 50, ‡, whereas the bottom of Figure 4 those of the RNiC ji RNiD jj-tests. It is seen that, for Models A and B, the powers for N50 are reasonably well approximated by the limiting powers, whereas, for Models C and D, the aprroximation is still not good enough for N50. It is seen that the former powers are higher than the latter, as is anticipated from Theorem 4. This means that the existence of heterogeneous trends decreases the power, as in the OLSE-based tests.

 

 

y46 •Εz

 

3.3. Limiting Power Envelopes

In previous subsections, dealing with Models A through D, we considered panel unit root tests based on OLSE and GLSE, for which the limiting local powers were computed and power comparisons were made among those tests, examining the cross-sectional effects. In this subsection we derive the power envelopes, from which the performance of these tests can be evaluated. The idea was earlier developed in the time series context by Elliott et al. i1996j, and was extended to the nonstationary panel data by Moon et al. i2007j. Here we derive the power envelopes for Models A through D, paying attention to the cross-sectional effects.
@@Let us consider the testing problem

where ƒΖN ƒΖ/NƒΘ with ƒΖ being a known positive constant. We assume that the true value of ƒΟ under H1 is y47 •Εz given by ƒΟicj1|cN /T with cNc/NƒΘ. AssumingoƒΓit p` NIDi0, ƒΠ2j, the Neyman-Pearson lemma tells us that the test which rejects H0 for small values of

is MPI, where itiM ji0j and itiM ji1j are the GLS residuals obtained from Model M under H0 and H1, respectively. The residual itiM ji0j is the same as the GLS residual dealt with in the last subsection, that is, iM ji0jiMj, whereas iM ji1jy | X iM jƒΑiM ji1j, where

with ƒΆiƒΖjIN CiƒΟiƒΖjjC'iƒΟiƒΖjj. Here CiƒΟiƒΖjj is the T ~ T lower triangular matrix withis, T j-th element being ƒΟ|s-t|iƒΖj for st and 0 otherwise. The test based on SNTiM jiƒΖj with fixed ƒΖ is called the point optimal invariantiPOIj testmKingi1987jn.
@@The following theorem gives the weak convergence of SNTiM jiƒΖj as T ¨‡ for each N, the proof of which is given in the Appendix.

 

Theorem 5. As T ¨‡ under ƒΟ1| cN /T for each N, the MPI test statisticSNTiM jiƒΖj ini18j follows

where

with ƒΒN 1{ƒΖN {ƒΖ2N .

 

@@It is seen that the expression for SNiM jiƒΖj ini19j is of a similar nature to QNiM j ini12j and RNiMj in i16j. It is also noticed that the distribution of SNiM jiƒΖj depends on ƒΖ that is the value under H1. Thus the MPI test based on SNiM jiƒΖj is not uniformly best, but we can modify SNiM jiƒΖj so that the distribution of the modified statistic does not depend on ƒΖ as N ¨‡. Then we can compute the limiting power of the test based on a modified statistic which yields the limiting power envelope of all the invariant tests for Model M. The following theorem gives such statistics and the power envelopes.

y48 •Εz

Theorem 6. The limiting powers of the tests based on the MPI statistics SNiM jiƒΖj ini19j at the 100ƒΑ% level as N ¨‡ under ƒΖN ƒΖ/NƒΘ and cNc/NƒΘ are given by

where ƒΘ1/2 for Models A and B, and ƒΘ1/4 for Models C and D.

 

@@The limiting powers of the modified tests give the power envelope of all the invariant tests. Comparing Theorem 6 with Theorem 4 it is seen that the power functions of the GLSE-based tests coincide with the power envelopes. Thus the GLSE-based tests are asymptotically efficient, unlike the time series case. This is a merit of panel tests as N ¨‡.
@@There are some other tests that are asymptotically efficient. Here we take up two such tests. Define

The test that rejects H0 for KNTiM j small is locally best invariantiLBIj, although the test is inapplicable to Models C and D because KNTiM j ί 0 for M  C, D, whereas the test that rejects H0 for LNTiM j small is LBI and unbiasediLBIUj for MC,D mTanakai2017, Chap. 10jn. We have, as T ¨‡ for each N,

Since it can be shown that, for MA, B,

we have iKNiM j|1j Λ Ni|c, 2j by putting cNc /, which implies that the KNiM j-tests for MA, B are asymptotically efficient. It is evident that the LNiM j-tests are asymptotically efficient for MC, D.

y49 •Εz

Note that the LBI and LBIU tests in the time series caseiN1j are asymptotically inefficientmTanaka i1996, Chap. 9jn. We also note in passing that, if the GLS residual in i23j is replaced by the OLS residual, the resulting statistic is essentially the Durbin-Watson statistic and the corresponding test is asymptotically inefficient.

 

3.4. Effect of Temporal or Cross-Sectional Dependence

Here we consider the situation where there exists temporal or cross-sectional dependence of the error term oƒΓit p in i5j and examine the effect of such dependence on the test statistics obtained in previous subsections.
@@Let us first consider temporal dependence. For this purpose we assume

where ƒΣi iLj1{ƒΣi 1L{ƒΣi 2L 2{E E E with L being the lag-operator. The distributional properties of the statistics T iiM j|1j in i12j and T iiM j|1j in i16j are affected by this relaxation. In fact, it can be shown that, as T ¨‡ for each N,

where Ui iM j and Vi iM j are defined in i12j, and Wi iM j and Xi iM j are defined ini16j, whereas ƒΙi  is the ratio of the short-run to long-run variances of oƒΓit p given byi1j. The above statistics depend on the short-run and long-run variances of the error term that characterize temporal dependence.
@@We next consider cross-sectional dependence, for which we assume that

@It then follows that

y50 •Εz

Here iiM j, iiM j, iiM j, and iiM j replace Ui iM j, Vi iM j, Wi iM j, and Xi iM j, respectively, with Yi irj replaced by iirj, where

and oiirjp is the standard Brownian motion with

The test statistics depend on the covariances ƒΠik of the error term that characterize cross-sectional dependence.
@@It is recognized from the above observations that, to use the asymptotic results obtained in previous subsections, we need to modify the statistics to make them independent of nuisance parameters. This remains to be done.

 

4. CONCLUDING REMARKS

 

Under a simple setting, we have presented a unified approach to deriving the limiting local powers of panel AR unit root tests, paying attention to the cross-sectional effect of N. For this purpose it is necessary to compute moments up to the second order of the limiting statistic in the time series direction. We found it easier to use its m.g.f., unlike in the literature. It happened that the tests that were not powerful in the time series case become more powerful in the panel case. It was also found that the existence of a common intercept and/or a common trend does not affect the asymptotic behavior of the tests. This holds for not only the tests based on OLS and GLS residuals, but also power envelopes.
@@The present approach can be applied to unit root tests for other types of panel models such as panel moving average models or panel error components models. Some simple extensions are found in Tanaka i2017, Chap. 10j. For these models the panel LBI or LBIU tests can be used and the corresponding statistics have a distributional structure similar to the panel AR unit root tests discussed in this paper. Details are reported in Tanakai2018j.

 

5. APPENDIX: PROOFS OF THEOREMS

 

Proof of Thereom 1: We first deal with Model D. Given the OLSEs i  and i of ƒΏi  and ƒΐi , respectively, the OLS residual is

y51 •Εz

The continuous mapping theoremiCMTj yields, as T ¨‡ with N fixed,

Thus the relationi12j is proved for Model D by the CMT.
@@We next deal with Model B, for which the OLSEs and of ƒΏ and ƒΐ are given by

where IN is the identity matrix of order N, is the Kronecker product, and

y52 •Εz

We then have

where ƒΕtiƒΕ1t , . . . , ƒΕNt j'. It follows that T iiB j|1jUitiB j/VitiB j, where it holds that, as T ¨‡ with N fixed,

Let us put Y irjiY1irj, . . . , YN irjj'. Following Nabeya i2000j, consider Z irjHY irj, where H  is the N ~N orthogonal matrix with the first row being i'N / so that oZirjpoYirjp and Y'irjiNY'irjH'HiN  Z1irj. Then it holds that

y53 •Εz

which shows thati12j holds for MB by the CMT. We can provei12j for MA and C similarly.

 

Proof of Theorem 2: Let us consider Model B. The joint m.g.f of Ui iB j and Vi iB j is given in the text, which can be used to compute first and second moments of these quantities as cN ¨ 0 by employing computerized algebra and the Taylor expansion. Here we show how to compute EiUi iB jj and VariUi iB jj. It follows fromi13j that

where

with gi0jecN , a1isinhcNj/cN and a2coshcN. Then we have

Here it holds that

Then, using the expansions for e-cN , e-2cN , a1, and a2 as cN ¨ 0,

y54 •Εz

we can obtain moments of Ui iB j. We now have

It follows from the WLLN that QNiB j ¨ |3 in probability. Thus we consider

The m.g.f. of Ui iB j{3Vi iB j can be obtained from miB jix, y j by putting y3x, which yields

Putting cNc/ , we now have

which proves the second relation in the theorem. We note in passing that

so that it holds that

are uncorrelated.
@@For the other models the first two moments of Ui iM j and Vi iM j as cN ¨ 0 can be obtained from their joint m.g.f.s given in Tanakai2017, Chap. 10j. More specifically we have

y55 •Εz

These moments yield the limiting powers for MA, C, D shown in the theorem, which establishes Theorem 2. We note in passing that, when cN0,

so that Ui iC j and Vi iC j are correlated, whereas Ui iD j and Vi iD j are uncorrelated.

 

Proof of Theorem 3: We first deal with Model D. Defining the T ~T  lower triangular matrix C with C is, t j1 for st and C is, t j0 for s ƒ t, the GLSEs of ƒΏi  and ƒΐi  are given by

so that we have

Consider T i|1j Uit iD j/Vit iD j. It holds that, as T ¨‡ with N  fixed,

y56 •Εz

which provesi16j for MD.
@@We next consider Model B. The GLSEs of ƒΏ =iƒΏ1, . . . , ƒΏN j' and ƒΐ are given by

which yields

Then, applying the orthogonal transformation ZirjHYirj used in the proof of Theorem 1, it holds that, as T ¨‡ with N  fixed,

Noting thatoZ irjpoY irjp, the relationi16j is proved for MB. We can provei16j for MA and C similarly, which establishes Theorem 3.

y57 •Εz

Proof of Theorem 4: It follows from Theorem 3 that, as N  becomes large,

where Wi iM j and Xi iM j are defined in i16j. The limiting distribution of RNiM j for MA, B is the same as that of QNiA j in i12j, which proves the first relation. For Models C and D, Wi iM j|1/2, whereas the m.g.f. of Xi iM j is given in Tanakai2017, Chap. 10j and

which yields RNiM j ¨|3 in probability. Then we consider

where

Putting cNc/N 1/4, we now have

which proves the second relation, and Theorem 4 has been established.

 

Proof of Theorem 5: The relationi19j for Model D was proved in Tanakai2017, Chap. 10j. Let us consider Model B, for which it holds that

where

It holds that

y58 •Εz

which gives

Thus it holds that, as T ¨‡ with N  fixed,

where

Applying the orthogonal transformation ZirjHYirj, we obtain

where

y59 •Εz

Noting that oZ irjpoY irjp and iitiB ji0j|i,t|1iB ji0jj2/T ¨ ƒΠ2 in probability, we prove i19j for MB after some manipulations. We can prove i19j for MA and C similarly, which establishes Theorem 5.

 

Proof of Theorem 6: For Models A and B we have SNiA jiƒΖj  Zi iA jiƒΖj/N{opi1j, where Zi iA jiƒΖjis defined in i19j and the m.g.f. of iZi iA jiƒΖj{ƒΖN j/ƒΖ2N  is given in Tanakai2017, Chap. 10j as

where ƒΚ  Differentiation of miA jixj gives us

It follows from the CLT that, as N ¨‡ under ƒΖN = ƒΖ/N 1/2 and cN = c/N 1/2,

which yieldsi20j.
@@For Models C and D we have SNiC jiƒΖjZi iC jiƒΖj/N{opi1j, where Zi iC jiƒΖj is defined ini19j and the m.g.f. of iZi iC jiƒΖj{ƒΖN j/ƒΖ2N  is given in Tanakai2017, Chap. 10j as

where ƒΚ  . Differentiation of miC jixj gives us

It follows from the CLT that, as N ¨‡ under ƒΖN =ƒΖ/N 1/4 and cN =c/N 1/4,

y60 •Εz

which provesi21j, and Theorem 6 has been established.

 

REFERENCES

Breitung, J. i2000j. The local power of some unit root tests for panel data. Advances in Econometrics 15:161-178.

 

Elliott, G., T. Rothenberg, and J. Stock i1996j. Efficient tests for an autoregressive unit root. Econometrica 64:813-836.

 

Imhof, J.P. i1961j. Computing the distribution of quadratic forms in normal variables. Biometrika 48:419-426.

 

King, M.L. i1987j. Towards a theory of point optimal testing. Econometric Reviews 6:169-218.

 

Moon, H.R. and B. Perron i2008j. Asymptotic local power of pooled t-ratio tests for unit roots in panels with fixed effects. Econometrics Journal 11:80-104.

 

Moon, H.R., B. Perron, and P.C.B. Phillips i2007j. Incidental trends and the power of panel unit root tests. Journal of Econometrics 141:416-459.

 

Nabeya, S. i2000j. Asymptotic distributions for unit root test statistics in nearly integrated seasonal autoregressive models. Econometric Theory 16:200-230.

 

Nabeya, S. and K. Tanaka i1990j. A general approach to the limiting distribution for estimators in time series regression with nonstable autoregressive errors. Econometrica 58:145-163.

 

Tanaka, K.i1996j. Time Series Analysis: Nonstationary and Noninvertible Distribution Theory. New York: John Wiley & Sons.

 

Tanaka, K. i2017j. Time Series Analysis: Nonstationary and Noninvertible Distribution Theory. 2nd ed. New York: John Wiley & Sons.

 

Tanaka, K. i2018j. Computing limiting local powers and power envelopes of panel MA unit root tests and stationarity tests. to appear in Econometric Theory.