1. MOTIVATION
In the past decade or so, there has been much interest in testing for the presence of a unit root in panel data. Many researchers have proposed statistics to test the hypothesis of a common unit autoregressive root. Recent surveys by Baltagi and Kao (2000), Choi (2004), Hurlin and Mignon (2004), and Breitung and Pesaran (2005) provide an overview of these developments.
In this context, Breitung (2000) proposed a t-ratio-type test statistic for testing a panel unit root. Through numerical analysis, he claimed that his test has “nice” power properties within a certain local neighborhood of unity. The present paper investigates analytically the asymptotic power properties of Breitung's test and clarifies some of the analytic results in Breitung (2000). Specifically, we show that the limiting distribution of the Breitung test is the same under the null and the O(n−1/2T−1) local alternatives considered by Breitung, so that the test has trivial power in such narrow neighborhoods. We provide expressions for the local asymptotic power of this test for wider n−1/4T−1 local departures from the null, discuss comparative results for other procedures such as point optimal tests, and study the accuracy of the asymptotic approximations in finite samples.
2. BREITUNG'S TEST STATISTIC AND CLAIMED POWER
Suppose that panel data yit is generated by the following simple components model:
where the unobserved error term xit follows
Our main interest is in testing the presence of a unit root in all cross-sectional units, namely,
For this, we assume that εit ∼ iid(0,σ2) with E(εit4) < ∞, the initial observations xi0 are independent and identically distributed (i.i.d.) across i with E(xi04) < ∞ and independent of εit for all t ≥ 1 and i.1
These restrictions are made for simplicity in the following analysis and can be relaxed to cover more general cases.
Notice that this testing problem is invariant to the following linear transformation: yit* = yit + μi* + βi*t. To construct a test that is invariant to the transformation, Breitung (2000) suggested the use of the following transformed data:
for t = 1,…,T − 1, where st2 = (T − t)/(T − t + 1), and
for t = 2,…,T.2
Note that (Δyit)* and yit-1* correspond to the terms in equations (16) and (17), respectively, of Breitung (2000).
where
is a consistent estimator of σ2. This is a t-ratio statistic from a pooled regression of the transformed data (Δyit)* on yit−1*. Define
. Under the null hypothesis, a direct calculation (given in the Appendix) shows that for 2 ≤ t ≤ T − 1
Hence, by central limit theory (see, e.g., Phillips and Moon, 1999), it is possible to show that
as n,T → ∞. On the other hand,
so that
as n,T → ∞. Therefore, under the null hypothesis, it follows that as n,T → ∞,
Breitung's test for H0 with size α rejects H0 if
, where zα is the (1 − α) quantile of the N(0,1) distribution.
To analyze the local power of the test based on
, Breitung (2000) considered the following local parameterization:
and he claimed (see Breitung, 2000, Thm. 5) that the test BnT has asymptotically significant local power against the local alternative
To verify the claim, because
under the local alternative (2), we can express the transformed variables as follows:
Note that Ait and Cit do not depend on c. Then, the preceding formulas lead to the following decomposition of
:
LEMMA 1. The following conditions hold under the local alternative ρi = 1 − (c/n1/2T) as n,T → ∞.
It follows from Lemma 1 that
. Together with the limit of
, under the local alternative in (2), we therefore obtain
and have the following result.
THEOREM 2. Under the local alternative ρi = 1 − (c/n1/2T), assumed in Breitung (2000), the Breitung test statistic
has asymptotic power equal to the size of the test.
The preceding result does not imply a mistake in Breitung (2000). Specifically, Breitung stated that the power of his test depends on a quantity that in our notation is
. However, the preceding results demonstrate that this expectation has zero limit. Instead of an analytical evaluation of the limit, Breitung relied on a numerical approximation to the expectation that suggested nonnegligible power against alternatives of the form (2). Our results indicate that the numerical approximation is poor and that consequently the implications regarding power in neighborhoods of the form (2) are misleading.
3. ASYMPTOTIC LOCAL POWER OF BREITUNG'S TEST
The main problem with Theorem 2 is that the local neighborhood in (2) shrinks too fast to unity and so the asymptotic power is trivial. To correct the result, this section considers local neighborhoods that shrink to unity more slowly and derives an asymptotic local power function for the Breitung test.
First, consider the local-to-unity region defined by
where the ci may be defined as the realizations of a sequence of i.i.d. random variables whose support lies in an interval of the form [−Mlc,Muc] for some 0 < Mlc < ∞ and 0 ≤ Muc < ∞. Let ci be independent of εjs for all i,j,s. Under these conditions, the numerator of Breitung's test statistic can be decomposed as
where the components Ait, Bit, Cit−1, and Dit−1 are defined in the previous section.
LEMMA 3. Under the local alternative (4), the following conditions hold as n,T → ∞.
Using Lemma 3, the limit distribution under the local alternative hypothesis (4) is
from which we deduce the following theorem.
THEOREM 4. Against the local alternative in (4), the asymptotic local power of Breitung's test is
, where Φ( ) is the standard normal cumulative distribution function.
Remarks.
(1) Contrary to the Breitung (2000) claim, nontrivial local power is defined in neighborhoods that shrink toward the null hypothesis at the rate 1/n1/4T. This is the same rate worked out by Moon, Perron, and Phillips (2005) in defining the power envelope for panel unit root testing in the context of incidental trends for models in this form.
(2) In contrast to Breitung (2000), the preceding result is obtained against heterogeneous alternatives. The test therefore has significant power against this type of hypothesis despite pooling.
(3) The power of Breitung's test depends on the second moments of the local-to-unity parameters. Thus, for a given mean autoregressive parameter, the more heterogeneous the alternatives are, the easier they are to detect.
(4) Moon et al. (2005) derive the power envelope for the preceding testing problem, suggest a common point optimal (CPO) test, and discuss the local asymptotic power of other tests such as those proposed by Levin, Lin, and Chu (2002), Ploberger and Phillips (2002), and Moon and Phillips (2004). According to those results, the test based on BnT is more powerful than the Levin et al. and Moon and Phillips tests but less powerful than the Ploberger and Phillips and CPO tests.
A small-scale simulation was conducted to assess the accuracy of these asymptotic results in finite samples. We use the following data generating process:
The heterogeneous trend coefficients are taken to be iid N(0,1). We assume that the error term is independent in both time and cross-section dimensions with a Gaussian distribution and identical variances. We consider four values for n (10, 25, 100, and 250) and three values for T (50, 100, and 250). All tests are carried out at the 5% significance level, and the number of replications is set at 10,000.
The autoregressive parameters are generated according to (4). We consider the following nine distributions for the local-to-unity parameters:
(0) ci = 0 ∀i (size),
(1) ci ∼ iidU [0,2],
(2) ci ∼ iidU [0,4],
(3) ci ∼ iidU [0,8],
(4) ci ∼ iidχ2(1),
(5) ci ∼ iidχ2(2),
(6) ci ∼ iidχ2(4),
(7) ci = 1 ∀i,
(8) ci = 2 ∀i.
These distributions enable us to examine performance of the tests as the mass of the distribution of the localizing parameters moves away from the null hypothesis. We can also look at the effect of homogeneous versus heterogeneous alternatives (cases (1) and (4) versus (7) and cases (2) and (5) versus (8)) together with the role of the higher order moments of the distribution. For instance, case (7) has the same mean as cases (1) and (4) but smaller higher order moments than the other two cases. The same situation arises for cases (2), (5), and (8) and for cases (3) and (6). Note that the alternatives with χ2 distributions do not fit our asymptotic framework because they have unbounded support.
Table 1 presents the results. The second column provides the size and power predicted by our asymptotic theory using the moments of ci. The other columns in the table report the size and size-adjusted power of the tests for the various combinations of n and T. If asymptotic theory were a reliable guide to finite-sample behavior, all columns in the table would be very close.
Size and size-adjusted power of Breitung's test
Overall, we see that the test performs much more closely to the asymptotic theory as both n and T increase. The test underrejects for large n relative to T. Power is usually below what is predicted by asymptotic theory, especially for the more distant alternatives, but the discrepancy diminishes with increases in either n or T or both. Finally, experiments where the local-to-unity parameters have a fatter tailed distribution tend to have higher power as predicted. Thus, for a given mean autoregressive parameter, more heterogeneous alternatives are easier to detect (despite the pooling approach used in the test). There also seems no sign that the unboundedness of the χ2 distributions affects the validity of the asymptotic theory.
APPENDIX: Proofs
Proof of Lemma 1.
Part (i).
Part (i) follows because the distribution of
is identical to that of
under the null because neither Ait nor Cit−1 depends on c.
Part (ii).
Under the given assumptions, by a weak law of large numbers,
Recall the definition
. For 2 ≤ t ≤ T − 1, we have
Also,
Combining these, we have the required result for part (ii) because
Parts (iii) and (iv).
These follow by a weak law of large numbers. █
Proof of Lemma 3. In this proof we use the notation
, where xi0* = xi0. Also, for notational convenience, we write εi0 = xi0 = xi0*. Then, by definition, for t ≥ 1,
Parts (i) and (iv).
They follow in the same fashion as parts (i) and (iv) of Lemma 1.