Testing for a change in persistence in the presence of non-stationary volatility

In this paper we consider tests for the null of (trend-) stationarity against the alternative of a change in persistence at some (known or unknown) point in the observed sample, either from I(0)$I\left(0\right)$ to I(1)$I\left(1\right)$ behaviour or vice versa, of, inter alia, [Kim, J., 2000. Detection of change in persistence of a linear time series. Journal of Econometrics 95, 97–116]. We show that in circumstances where the innovation process displays non-stationary unconditional volatility of a very general form, which includes single and multiple volatility breaks as special cases, the ratio-based statistics used to test for persistence change do not have pivotal limiting null distributions. Numerical evidence suggests that this can cause severe over-sizing in the tests. In practice it may therefore be hard to discriminate between persistence change processes and processes with constant persistence but which display time-varying unconditional volatility. We solve the identified inference problem by proposing wild bootstrap-based implementations of the tests. Monte Carlo evidence suggests that the bootstrap tests perform well in finite samples. An empirical illustration using US price inflation data is provided.     

Beyond panel unit root tests: Using multiple testing to determine the nonstationarity properties of individual series in a panel

Most panel unit root tests are designed to test the joint null hypothesis of a unit root for each individual series in a panel. After a rejection, it will often be of interest to identify which series can be deemed to be stationary and which series can be deemed nonstationary. Researchers will sometimes carry out this classification on the basis of n$n$ individual (univariate) unit root tests based on some ad hoc significance level. In this paper, we suggest and demonstrate how to use the false discovery rate (FDR)$\left(FDR\right)$ in evaluating I(1)/I(0)$I\left(1\right)/I\left(0\right)$ classifications.     

Ratio-based estimators for a change point in persistence

We study estimation of the date of change in persistence, from I(0)$I\left(0\right)$ to I(1)$I\left(1\right)$ or vice versa. Contrary to statements in the original papers, our analytical results establish that the ratio-based break point estimators of Kim [Kim, J.Y., 2000. Detection of change in persistence of a linear time series. Journal of Econometrics 95, 97–116], Kim et al. [Kim, J.Y., Belaire-Franch, J., Badillo Amador, R., 2002. Corringendum to “Detection of change in persistence of a linear time series”. Journal of Econometrics 109, 389–392] and Busetti and Taylor [Busetti, F., Taylor, A.M.R., 2004. Tests of stationarity against a change in persistence. Journal of Econometrics 123, 33–66] are inconsistent when a mean (or other deterministic component) is estimated for the process. In such cases, the estimators converge to random variables with upper bound given by the true break date when persistence changes from I(0)$I\left(0\right)$ to I(1)$I\left(1\right)$. A Monte Carlo study confirms the large sample downward bias and also finds substantial biases in moderate sized samples, partly due to properties at the end points of the search interval.     

Estimating deterministic trends with an integrated or stationary noise component

We propose a test for the slope of a trend function when it is a priori unknown whether the series is trend-stationary or contains an autoregressive unit root. The procedure is based on a Feasible Quasi Generalized Least Squares method from an AR(1) specification with parameter α$\alpha$, the sum of the autoregressive coefficients. The estimate of α$\alpha$ is the OLS estimate obtained from an autoregression applied to detrended data and is truncated to take a value 1 whenever the estimate is in a T^−δ$T^{-\delta }$ neighborhood of 1. This makes the estimate “super-efficient” when α=1$\alpha =1$ and implies that inference on the slope parameter can be performed using the standard Normal distribution whether α=1$\alpha =1$ or |α|<1$\left|\alpha \right|<1$. Theoretical arguments and simulation evidence show that δ=1/2$\delta =1/2$ is the appropriate choice. Simulations show that our procedure has better size and power properties than the tests proposed by [Bunzel, H., Vogelsang, T.J., 2005. Powerful trend function tests that are robust to strong serial correlation with an application to the Prebish–Singer hypothesis. Journal of Business and Economic Statistics 23, 381–394] and [Harvey, D.I., Leybourne, S.J., Taylor, A.M.R., 2007. A simple, robust and powerful test of the trend hypothesis. Journal of Econometrics 141, 1302–1330].     

Testing functional inequalities

This paper develops tests for inequality constraints of nonparametric regression functions. The test statistics involve a one-sided version of _Lp$L_p$-type functionals of kernel estimators (1≤p<∞)$(1\le p<\infty )$. Drawing on the approach of Poissonization, this paper establishes that the tests are asymptotically distribution free, admitting asymptotic normal approximation. In particular, the tests using the standard normal critical values have asymptotically correct size and are consistent against general fixed alternatives. Furthermore, we establish conditions under which the tests have nontrivial local power against Pitman local alternatives. Some results from Monte Carlo simulations are presented.     

Panel unit root tests in the presence of a multifactor error structure

This paper extends the cross-sectionally augmented panel unit root test (CIPS$CIPS$) proposed by Pesaran (2007) to the case of a multifactor error structure, and proposes a new panel unit root test based on a simple average of cross-sectionally augmented Sargan–Bhargava statistics (CSB$CSB$). The basic idea is to exploit information regarding the m$m$ unobserved factors that are shared by k$k$ observed time series in addition to the series under consideration. Initially, we develop the tests assuming that m⁰$m^0$, the true number of factors, is known and show that the limit distribution of the tests does not depend on any nuisance parameters, so long as k≥m⁰−1$k\ge m^0-1$. Small sample properties of the tests are investigated by Monte Carlo experiments and are shown to be satisfactory. Particularly, the proposed CIPS$CIPS$ and CSB$CSB$ tests have the correct size for all   combinations of the cross section (N$N$) and time series (T$T$) dimensions considered. The power of both tests rises with N$N$ and T$T$, although the CSB$CSB$ test performs better than the CIPS$CIPS$ test for smaller sample sizes. The various testing procedures are illustrated with empirical applications to real interest rates and real equity prices across countries.     

A complete asymptotic series for the autocovariance function of a long memory process

An infinite-order asymptotic expansion is given for the autocovariance function of a general stationary long-memory process with memory parameter d∈(−1/2,1/2)$d\in (-1/2,1/2)$. The class of spectral densities considered includes as a special case the stationary and invertible ARFIMA(p,d,q$p,d,q$) model. The leading term of the expansion is of the order O(1/k^1−2d)$O(1/k^{1-2d})$, where k$k$ is the autocovariance order, consistent with the well known power law decay for such processes, and is shown to be accurate to an error of O(1/k^3−2d)$O(1/k^{3-2d})$. The derivation uses Erdélyi’s [Erdélyi, A., 1956. Asymptotic Expansions. Dover Publications, Inc, New York] expansion for Fourier-type integrals when there are critical points at the boundaries of the range of integration - here the frequencies {0,2π}$\{0,2\pi \}$. Numerical evaluations show that the expansion is accurate even for small k$k$ in cases where the autocovariance sequence decays monotonically, and in other cases for moderate to large k$k$. The approximations are easy to compute across a variety of parameter values and models.     

A bootstrap theory for weakly integrated processes

This paper develops a bootstrap theory for models including autoregressive time series with roots approaching to unity as the sample size increases. In particular, we consider the processes with roots converging to unity with rates slower than n^-1$n^{-1}$. We call such processes weakly   integrated processes. It is established that the bootstrap relying on the estimated autoregressive model is generally consistent for the weakly integrated processes. Both the sample and bootstrap statistics of the weakly integrated processes are shown to yield the same normal asymptotics. Moreover, for the asymptotically pivotal statistics of the weakly integrated processes, the bootstrap is expected to provide an asymptotic refinement and give better approximations for the finite sample distributions than the first order asymptotic theory. For the weakly integrated processes, the magnitudes of potential refinements by the bootstrap are shown to be proportional to the rate at which the root of the underlying process converges to unity. The order of boostrap refinement can be as large as o(n^-1/2+?)$o(n^{-1/2+?})$ for any ?>0$?>0$. Our theory helps to explain the actual improvements observed by many practitioners, which are made by the use of the bootstrap in analyzing the models with roots close to unity.     

Nonstationary discrete choice: A corrigendum and addendum

We correct the limit theory presented in an earlier paper by Hu and Phillips [2004a. Nonstationary discrete choice. Journal of Econometrics 120, 103–138] for nonstationary time series discrete choice models with multiple choices and thresholds. The new limit theory shows that, in contrast to the binary choice model with nonstationary regressors and a zero threshold where there are dual rates of convergence (n^1/4$n^{1/4}$ and n^3/4$n^{3/4}$), all parameters including the thresholds converge at the rate n^3/4$n^{3/4}$. The presence of nonzero thresholds therefore materially affects rates of convergence. Dual rates of convergence reappear when stationary variables are present in the system. Some simulation evidence is provided, showing how the magnitude of the thresholds affects finite sample performance. A new finding is that predicted probabilities and marginal effect estimates have finite sample distributions that manifest a pile-up, or increasing density, towards the limits of the domain of definition.     

Adaptive estimation of the dynamics of a discrete time stochastic volatility model

This paper is concerned with the discrete time stochastic volatility model _Yi=exp(_Xi/2)_ηi$Y_i=exp(X_i/2){\eta }_i$, _Xi+1=b(_Xi)+σ(_Xi)_ξi+1$X_{i+1}=b\left(X_i\right)+\sigma \left(X_i\right){\xi }_{i+1}$, where only (_Yi)$\left(Y_i\right)$ is observed. The model is rewritten as a particular hidden model: _Zi=_Xi+_εi$Z_i=X_i+{\epsilon }_i$, _Xi+1=b(_Xi)+σ(_Xi)_ξi+1$X_{i+1}=b\left(X_i\right)+\sigma \left(X_i\right){\xi }_{i+1}$, where (_ξi)$\left({\xi }_i\right)$ and (_εi)$\left({\epsilon }_i\right)$ are independent sequences of i.i.d. noise. Moreover, the sequences (_Xi)$\left(X_i\right)$ and (_εi)$\left({\epsilon }_i\right)$ are independent and the distribution of ε$\epsilon$ is known. Then, our aim is to estimate the functions b$b$ and σ²${\sigma }^2$ when only observations _Z1,…,_Zn$Z_1,\dots ,Z_n$ are available. We propose to estimate bf$bf$ and (b²+σ²)f$(b^2+{\sigma }^2)f$ and study the integrated mean square error of projection estimators of these functions on automatically selected projection spaces. By ratio strategy, estimators of b$b$ and σ²${\sigma }^2$ are then deduced. The mean square risk of the resulting estimators are studied and their rates are discussed. Lastly, simulation experiments are provided: constants in the penalty functions defining the estimators are calibrated and the quality of the estimators is checked on several examples.     

Distribution-free bounds for serial correlation coefficients in heteroskedastic symmetric time series

We consider the problem of testing whether the observations _X1,…,_Xn$X_1,\dots ,X_n$ of a time series are independent with unspecified (possibly nonidentical) distributions symmetric about a common known median. Various bounds on the distributions of serial correlation coefficients are proposed: exponential bounds, Eaton-type bounds, Chebyshev bounds and Berry–Esséen–Zolotarev bounds. The bounds are exact in finite samples, distribution-free and easy to compute. The performance of the bounds is evaluated and compared with traditional serial dependence tests in a simulation experiment. The procedures proposed are applied to U.S. data on interest rates (commercial paper rate).