Testing for Structural Changes in the Presence of Long Memory
Walter Krämer* and Philipp SibbertsenFachbereich Statistik, Universität Dortmund, Germany
Abstract
We derive the limiting null distributions of the standard and OLS-based CUSUM- tests for a structural change of the coefficients of a linear regression model in the context of long-memory disturbances. We show that both tests behave fundamentally different in a long-memory environment, as compared to short memory, and that long memory is easily mistaken for structural change when standard critical values are employed.
Key words: long memory; structural change; CUSUM- tests JEL classification: C12; C22; C52
1. Introduction
It is by now well known that long memory and structural change are easily confused [Lobato and Savin (1998), Engle and Smith (1999), Granger and Hyung (1999), Diebold and Inoue (2001), Davidson and Sibbertsen (2002), and many oth-ers]. It is therefore of interest to know about the stochastic properties of procedures for detecting and measuring long memory when there is only structural change and for detecting and measuring structural change when there is only long memory.
While the former problem has attracted considerable attention, there has been rather little work on the latter. Nunes et al. (1995), Kuan and Hsu (1998), and Hsu (2001) show that conventional procedures for detecting and dating structural changes tend to find spurious breaks, usually in the middle of the sample, when there is in fact only long memory [Bai 1998)]. As concerns testing, Hidalgo and Robinson (1996) propose a test for structural change which is robust to long mem-ory but assumes a single change point at a known point in time, while Wright (1998) considers the CUSUM-test where no such assumptions are made. Otherwise the lit-erature on testing for structural change in the presence of long memory is rather thin.
Below we follow White (1998) and consider the behaviour of the standard and the OLS-based CUSUM-tests, whose limiting distributions are well understood in the context of various regressor-sequences and i.i.d. or short memory disturbances
Received January 2, 2002, accepted November 25, 2002.
* Correspondence to: Fachbereich Statistik, Universität Dortmund, D-44221 Dortmund, Germany. Re-search was supported by Deutsche Forschungsgemeinschaft. We are grateful to A. Zeileis for expert com-putational assistance and to an anonymous referee for constructive criticism and comments.
[Krämer et al. (1988) and Ploberger and Krämer (1992, 1996)]. As shown by Wright (1998) for the OLS-based CUSUM-test and the special case of polynomial regres-sors, these limiting distributions are not robust to departures from short memory—in fact, the OLS-based CUSUM-test has an asymptotic size of unity. This is so because the cumulated sums of the OLS-residuals, although still converging to some well-defined limit process, need a larger normalizing factor for doing so.
The present paper allows for more general regressor sequences and covers the conventional CUSUM-test based on recursive residuals as well. The intuition behind looking at the conventional CUSUM-test was that recursive residuals, other than OLS-residuals, might not inherit the long-memory characteristics of the true distur-bances which lead to the exaggerated size of the OLS-based CUSUM-test. However, this is not the case. As is shown below, the null distribution of both the standard and the OLS-based CUSUM-tests tends to infinity, for the same reason already found by Wright (1998): cumulated sums of both types of residuals tend to infinity much faster than in the case of short memory. This is a rather negative result which con-firms related theorems from the structural-change-vs-long-memory-literature: Simi-lar to structural change being mistaken for long memory, long memory is likewise easily mistaken for structural change [see also Krämer et al. (2001)]. While there do exist solutions for some special cases, such as Künsch's (1986) procedure to dis-criminate between long memory and monotonic trends, or Teverovsky and Taqqu's (1997) suggestion based on empirical variances, a general treatment of this problem is still missing and it remains an open problem to efficiently discriminate between the two.
2. Two Unpleasant Theorems
We consider the standard linear regression model t
t t β x ε
y c ,
t = 1, 2,…, T (1) with non-stochastic, fixed regressors x and stationary, zero-mean Gaussian dis-t turbances ε . We assume that t
¦
T o f t t c x T 1 1 and (2)¦
T co t t t Q x x T 1 (finite,nonsingular). 1 (3) These are standard assumptions in linear regression large sample asymptotics; they exclude trending data, which require separate treatment and proofs which differ from the ones below.We are concerned with testing the model (1) against the alternative of unspeci-fied structural change in the regression coefficients β . We consider first the OLS-based CUSUM-test as proposed by Ploberger and Krämer (1992). This test
rejects the null hypothesis of no structural change for large values of TS := sup ( ) 1 0λ CT λ
,
where (4) > @¦
Tλ t t ε T λ T σ e C 1 , 1 2 1 ˆ : ) ( (5)and where et: ytxt'βˆare the OLS-residuals from (1), ˆ ( ) ( )
1 1 ' 1 t T t t T t xtxt xy β ¦ ¦
is the OLS estimator, and σˆ: ( e2/T)12
t
¦
.The limiting null distribution of TS is well known for white noise and short memory disturbances. Our first theorem extends these results to stationary long-memory disturbances, where the ε follow a stationary ARFIMA(p,d,q) t process, i.e., ( 1 B) εt is a stationary ARMA(p,q) process, and where stationarity requires that 0 < d < 21 . For such processes, it is well known that autocorrelations decay only hyperbolically according to
, ) ( ) ( d k t tε L k k ε E kof (6)
where L(k) is slowly varying at infinity.
In the following, “” denotes weak convergence and “oP ” denotes con-vergence in probability.
Theorem 1: In the regression model (1) with disturbance as in (6), we have
), ( ) ( ) (λ B λ cQ 1ξ λ C Td T d c (7)
where )Bd(λ is fractional Brownian Motion with self-similarity parameter d and ) , 0 ( ~ ) (λ N λσ2Q ξ ε .
Remark: Fractional Brownian motion is defined as a Gaussian process with sta-tionary increments and variance 2( ) 2d1
d λ kλ
EB , where k is a positive constant. It is self-similar with parameter d. This means that for every factor c, cdBd( λc ) has the same distribution as the original process Bd(λ). The standard Brownian mo-tion is self-similar with parameter 1/2. Fracmo-tional Brownian momo-tion has the follow-ing integral representation which can be found, for example, in Tsay and Chung (2000): ³ ³ f u * * * 0 0 )) ( ] ) ( ) [( ) ( ) ( ( ) 2 1 ( ) 1 ( ) 1 ( ) 2 1 ( ) (λ dd dd λ sdB s λ s λ dB s Bd λ d d ,
where B(s) denotes standard Brownian motion and *(s) denotes the Gamma func-tion.
}, ) ( { ˆ ) ( [ ] 1 ] [ 1 1 1 1 1 2 1
¦
¦
¦
¦
Tλ c c t λ T t T t t t T t t x t t ε T λ T σ ε x xx xε C (8)where
> @
Tλ is the largest integer smaller than Tλ . This implies that ε λ T t T t t T x t t d d T dC λ T zTλ T x xx xε σ T ( ) { [ ] [ ] ( ) }/ˆ 1 1 1 1 2 1 2 1¦
¦
¦
c c,
(9)where zt zt1εt, z0 0 and where the limit in distribution of this expression is
to be derived. To see this limit, note that ), ( ] [ 2 / 1 Z σ B λ Td Tλ ε d (10)
¦
] o [ 1 1 Tλ t t c λ x T , (11) , ) 1 ( [ ] 1 1 1 1 o c¦
xx λ Q T λ T t t t and (12) ). ( ] [ 1 2 1 λ ξ ε x T Tλ t t t d¦
(13) The last convergence follows from a limit theorem of Taqqu (1975), where it is shown that ¦> @Tλt t
d ε
T 1/2 1 converges to a Gaussian random variable. Because we
assume non-stochastic regressors x , this result can be also applied in our situation. t The covariance structure of the Gaussian random variable ξ(λ)follows by direct algebra and mean zero follows because ε was assumed to be a mean zero random t variable.
As to the variance, it is easily seen that
¦
T¦
o t T t ε P t ε eT εT p σ σ 1 1 2 2 2 2 0 (1) . ˆ (14)Therefore, collecting together the limits in (10) through (14), the limiting relation-ship (7) follows.
From (7), it is immediately seen that TS oP f under
0
H , so the OLS-based CUSUM-test is extremely non-robust to long-memory disturbances in the sense that long memory is easily mistaken for structural change when conven-tional critical values are employed. The convergence to infinity of the test statistics becomes faster for larger values of d. Our Monte Carlo experiments in section 3 be-low will show that this leads to quite dramatic errors of the first kind even for rather small samples and small values of the long-memory parameter d.
Next we consider the standard CUSUM-test based on recursive residuals
t t t t t y xfβ e ) 1 ( ˆ ~ c , where
) 1 ( ) 1 ( 1 1 ) 1 ( ) 1 ( ( ) ˆ t X t cXt X t cyt β , and (15) 2 1 1 ) 1 ( ) 1 ( ) ) ( 1 ( t t t t t x X X x f c c t = K+1, K+2,…,T (16)
and where K is the number of regressors and the superscript t−1 means that only observations 1, 2,…, t−1 are used. The standard CUSUM-test rejects for large val-ues of ), 2 1 /( ) ( sup 1 0 λ W λ λ ST T (17) where
¦
[ ( )] 1 1 2 1 ~ ˆ : ) ( K λ T K K t t ε T λ T σ e W . (18)Our next result shows that this process likewise converges to a well-defined limit process only after additional normalization:
Theorem 2: In the regression model (1) with disturbances as in (6), we have
), ( ) (λ B λ W Td T o d (19)
where again Bd(λ) is fractional Brownian Motion with self-similarity parameter d.
Proof : Following Krämer et al. (1988), we write WT(λ) as
. ) ˆ ( 1 ) ( [ ( )] ( 1) 1 )] ( [ 1 t t K T λ K K t t K T λ K K t T ε β β x T λ W c
¦
¦
(20)Let Qj : 1T
¦
ij1xixi' . First we show that. ) ln )(ln ( ] ) [( ) ˆ ( max 2 1 2 1 1 ) ( β y x β x Q o T T β T t K d p t K j j j j j t c d d
¦
(21)Let .St:
¦
tj 1(yjxcjβ)xj By the law of the iterated logarithm for the sums of long-memory Gaussian random variables [see Taqqu (1977)] we obtain some bounds for the maximum of S . We have for some slowly varying function L(T) t), 1 ( )] ( ) ( ) / 2 [( max 1 ) 1 2( 12) 21 2 1 ( 2 ) 2 1 ( p d t O T ln ln T L T S T td d d d (22)
so (22) follows directly from Lemma 3.1 of Jureckova and Sen (1984).
Thus we can replace the recursive least squares estimator in (21) by the ap-proximation in (22). Combining (21) and (22) gives
¦ c ¦ ¦ c ] [ 1 1 2 1 ) 2 1 ( ] [ 1 ) 1 ( ) 1 ( ) ( lnln ) , ˆ ( 1 Tλ k t t j d j j ij λ T t t t t c y xβ oT T T β x y T (23) where °¯ ° ® ! c j i j i j i x Q x c i i j ij 0 1 1 ) 1 ( (24)
[see also Sibbertsen (2000)]. In view of a result by Sen (1984), we can now see that the variance of the expression (23) converges to 1. We have
) 1 ( 1 2 i O c i j ij ¦ d . (25)
So from theorem 5.1 of Taqqu (1975), one can see as in Theorem 1 that TdωT(λ) converges to fractional Brownian motion. The statement of the theorem is then an immediate consequence of (23).
Theorem 2 shows that the null distribution of the standard CUSUM-test tends to infinity as well, so the standard CUSUM-test has likewise an asymptotic size of unity. Next we check via Monte Carlo how serious this problem is in finite samples.
3. Some Finite Sample Monte Carlo Evidence
To see how far the empirical rejection rates overstate the nominal rejection rates, we generate data according to (1), with disturbances generated via the R statis-tical software package, which is available free of charge from the internet under http://cran.r-project.org. For the generation of the long-memory disturbances, the routine fracdiff was applied.
Figure 1 below gives the empirical rejection rates using 1000 runs and standard critical values from the i.i.d. disturbance case for the OLS-based CUSUM-test with a nominal significance level of α= 5 when the disturbances are in fact AR-FIMA(0,d,0). It confirms our theoretical results: rejection rates increase with d and sample size and produce misleading evidence even for small d and T. For instance, for a sample size of T = 200 and d = 0, 2, the true rejection probability is roughly 50% when in fact the test is done under the assumption that it is only 5%.
Figure 2 gives the corresponding empirical rejection rates for the standard CUSUM-test. Not surprisingly, and confirming the intuition which has led us to con-sider the standard CUSUM-test in the first place, the empirical size is not as far off the mark as for the OLS-based CUSUM-test, but the test is misleading here as well. For instance, for a sample size of T = 200 and d = 0, 2, we have a true rejection rate of about 20%, which is still four times as large as the nominal rejection rate of 5%.
Fig. 1.
Fig. 2.
Summing up, the paper shows that long-memory disturbances can easily induce the appearance of a structural break, confirming Nunes et al. (1995), Kuan and Hsu (1998), and Hsu (2001), where it was shown that the least squares change point estimator tends to identify a spurious change when there is only long memory.
References
Bai, J., (1998), “A Note on Spurious Break,” Econometric Theory, 14, 663-669. Davidson, J. and P. Sibbertsen, (2002), “Generating Schemes of Long-Memory
Processes,” Technical Report, 46/2002, SFB 475, Universität Dortmund. Diebold, F. S. and A. Inoue, (2001), “Long Memory and Regime Switching,”
Engle, R. F. and A. Smith, (1999), “Stochastic Permanent Breaks,” Review of Eco-nomics and Statistics, 81, 553-574.
Granger, C. W. J. and N. Hyung, (1999), “Occasional Structural Breaks and Long Memory,” UCSD Discussion Paper, 99-114.
Hidalgo, J. and P. M. Robinson, (1996), “Testing for Structural Change in a Long-Memory Environment,” Journal of Econometrics, 70,159-174.
Hsu, C.-C., (2001), “Change Point Estimation in Regressions with I(d)-Variables,” Economics Letters, 70, 147-155.
Jureckova, J. and P. K. Sen, (1984), “On Adaptive Scale-Equivariant M-Estimators in Linear Models,” Statistics and Decision, 2, 31-46.
Künsch, H., (1986), “Discrimination between Monotonic Trends and Long Range Dependence,” Journal of Applied Probability, 23, 1025-1030.
Krämer, W., W. Ploberger, and R. Alt, (1988), “Testing for Structural Change in Dynamic Models,” Econometrica, 56, 1355-1369.
Krämer, W., P. Sibbertsen, and C. Kleiber, (2001), “Long Memory vs. Structural Change in Financial Time Series,” Allgemeines Statistisches Archiv, 86, 83-96. Kuan, C.-M. and C.-C. Hsu, (1998), “Change Point Estimation of Fractionally
Inte-gral Processes,” Journal of Time Series Analysis, 19, 552-575.
Lobato, I. N. and N. E. Savin, (1998), “Real and Spurious Long-Memory Properties of Stock Market Data (with Discussion and Reply),” Journal of Business and Economic Statistics, 16, 261-283.
Nunes, L. C., C.-M. Kuan, and P. Newbold, (1995), “Spurious Break,” Econometric Theory, 11, 736-749.
Ploberger, W. and W. Krämer, (1992), “The CUSUM-Test with OLS Residuals,” Econometrica, 60, 271-286.
Ploberger, W. and W. Krämer, (1996), “A Trend-Resistent Test for Structural Change Based on OLS Residuals,” Journal of Econometrics, 70, 175-185. Sibbertsen, P., (2000), “Robust CUSUM-M-Test in the Presence of Long Memory
Disturbances,” Technical Report, 19/2000, SFB 475, Universität Dortmund. Taqqu, M. S., (1975), “Weak Convergence of Fractional Brownian Motion to the
Rosenblatt Process,” Zeitschrift Wahrscheinlichkeitstheorie und verwandte Ge-biete, 31, 287-302.
Taqqu, M. S., (1977), “Law of the Iterated Logarithm for Sums of Non-Linear Func-tions of Gaussian Variables that Exhibits a Long Range Dependence,”
Zeitschrift Wahrscheinlichkeitstheorie und verwandte Gebiete, 40, 203-238. Teverovsky, V. and M. Taqqu, (1997), “Testing for Long-Range Dependence in the
Presence of Shifting Means or a Slowly Declining Trend using a Vari-ance-Type Estimator,” Journal of Time Series Analysis, 18, 279-304. Tsay, W. J. and C. F. Chung, (2000), “The Spurious Regression of Fractionally
In-tegrated Processes,” Journal of Econometrics, 96, 155-182.
Wright, J., (1998), “Testing for a Structural Break at Unknown Date with Long-Memory Disturbances,” Journal of Time Series Analysis, 19, 369-379.