A likelihood ratio test for stationarity of rating transitions

We study the time-stationarity of rating transitions, modelled by a time-continuous discrete-state Markov process and derive a likelihood ratio test. For multiple Markov processes from a multiplicative intensity model, maximum likelihood parameter estimates can be written as martingale transform of the processes, counting transitions between the rating states, so that the profile partial likelihood ratio is asymptotically χ²${\chi }^2$-distributed. An application to an internal rating data set reveals highly significant instationarity.     

On the properties of the coefficient of determination in regression models with infinite variance variables

We examine the asymptotic properties of the coefficient of determination, R²$R^2$, in models with α-stable$\alpha \text{-stable}$   random variables. If the regressor and error term share the same index of stability α<2$\alpha <2$, we show that the R²$R^2$  statistic does not converge to a constant but has a nondegenerate distribution on the entire [0,1]$[0,1]$ interval. We provide closed-form expressions for the cumulative distribution function and probability density function of this limit random variable, and we show that the density function is unbounded at 0 and 1. If the indices of stability of the regressor and error term are unequal, we show that the coefficient of determination converges in probability to either 0 or 1, depending on which variable has the smaller index of stability, irrespective of the value of the slope coefficient. In an empirical application, we revisit the Fama and MacBeth (1973) two-stage regression and demonstrate that in the infinite-variance case the R²$R^2$  statistic of the second-stage regression converges to 0 in probability even if the slope coefficient is nonzero. We deduce that a small value of the R²$R^2$  statistic should not, in itself, be used to reject the usefulness of a regression model.     

Quasi ML estimation of the panel AR(1) model with arbitrary initial conditions

In this paper we show that the Quasi ML estimation method yields consistent Random and Fixed Effects estimators for the autoregression parameter ρ$\rho$ in the panel AR(1) model with arbitrary initial conditions and possibly time-series heteroskedasticity even when the error components are drawn from heterogeneous distributions. We investigate both analytically and by means of Monte Carlo simulations the properties of the QML estimators for ρ$\rho$. The RE(Q)MLE for ρ$\rho$ is asymptotically at least as robust to individual heterogeneity and, when the data are i.i.d. and normal, at least as efficient as the FE(Q)MLE for ρ$\rho$. Furthermore, the QML estimators for ρ$\rho$ only suffer from a ‘weak moment conditions’ problem when ρ$\rho$ is close to one if the cross-sectional average of the variances of the errors is (almost) constant over time, e.g. under time-series homoskedasticity. However, in this case the QML estimators for ρ$\rho$ are still consistent when ρ$\rho$ is local to or equal to one although they converge to a non-normal possibly asymmetric distribution at a rate that is lower than N^1/2$N^{1/2}$ but at least N^1/4$N^{1/4}$. Finally, we study the finite sample properties of two types of estimators for the standard errors of the QML estimators for ρ$\rho$, and the bounds of QML based confidence intervals for ρ$\rho$.     

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.     

Estimating quadratic variation consistently in the presence of endogenous and diurnal measurement error

We propose an econometric model that captures the effects of market microstructure on a latent price process. In particular, we allow for correlation between the measurement error and the return process and we allow the measurement error process to have a diurnal heteroskedasticity. We propose a modification of the TSRV estimator of quadratic variation. We show that this estimator is consistent, with a rate of convergence that depends on the size of the measurement error, but is no worse than n^−1/6$n^{-1/6}$. We investigate in simulation experiments the finite sample performance of various proposed implementations.     

Edgeworth expansions and normalizing transforms for inequality measures

Finite sample distributions of studentized inequality measures differ substantially from their asymptotic normal distribution in terms of location and skewness. We study these aspects formally by deriving the second-order expansion of the first and third cumulant of the studentized inequality measure. We state distribution-free expressions for the bias and skewness coefficients. In the second part we improve over first-order theory by deriving Edgeworth expansions and normalizing transforms. These normalizing transforms are designed to eliminate the second-order term in the distributional expansion of the studentized transform and converge to the Gaussian limit at rate O(n⁻¹)$O\left(n^{-1}\right)$. This leads to improved confidence intervals and applying a subsequent bootstrap leads to a further improvement to order O(n^−3/2)$O\left(n^{-3/2}\right)$. We illustrate our procedure with an application to regional inequality measurement in Côte d’Ivoire.     

Distribution-free tests of stochastic monotonicity

This article proposes a nonparametric test of monotonicity for conditional distributions and its moments. Unlike previous proposals, our method does not require smooth estimation of the derivatives of nonparametric curves. Distinguishing features of our approach are that critical values are pivotal under the null in finite samples and that the test is invariant to any monotonic continuous transformation of the explanatory variable. The test statistic is the sup-norm of the difference between the empirical copula function and its least concave majorant with respect to the explanatory variable coordinate. The resulting test is able to detect local alternatives converging to the null at the parametric rate n^−1/2$n^{-1/2}$, with n$n$ the sample size. The finite sample performance of the test is examined by means of a Monte Carlo experiment and an application to testing intergenerational income mobility.     

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.     

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.     

Semiparametric binary regression models under shape constraints with an application to Indian schooling data

We consider estimation of the regression function in a semiparametric binary regression model defined through an appropriate link function (with emphasis on the logistic link) using likelihood-ratio based inversion. The dichotomous response variable Δ$\Delta$ is influenced by a set of covariates that can be partitioned as (X,Z)$(X,Z)$ where Z$Z$ (real valued) is the covariate of primary interest and X$X$ (vector valued) denotes a set of control variables. For any fixed X$X$, the conditional probability of the event of interest (Δ=1$\Delta =1$) is assumed to be a non-decreasing function of Z$Z$. The effect of the control variables is captured by a regression parameter β$\beta$. We show that the baseline conditional probability function (corresponding to X=0$X=0$) can be estimated by isotonic regression procedures and develop a likelihood ratio based method for constructing asymptotic confidence intervals for the conditional probability function (the regression function) that avoids the need to estimate nuisance parameters. Interestingly enough, the calibration of the likelihood ratio based confidence sets for the regression function no longer involves the usual χ²${\chi }^2$ quantiles, but those of the distribution of a new random variable that can be characterized as a functional of convex minorants of Brownian motion with quadratic drift. Confidence sets for the regression parameter β$\beta$ can however be constructed using asymptotically χ²${\chi }^2$ likelihood ratio statistics. The finite sample performance of the methods are assessed via a simulation study. The techniques of the paper are applied to data sets on primary school attendance among children belonging to different socio-economic groups in rural India.     

Asymptotic distribution of factor augmented estimators for panel regression

In this paper we derive an asymptotic theory for linear panel regression augmented with estimated common factors. We give conditions under which the estimated factors can be used in place of the latent factors in the regression equation. For the principal components estimate of the factor space it is shown that these conditions are satisfied when T/N→0$T/N\to 0$ and N/T³→0$N/T^3\to 0$ under regularity. Monte Carlo studies verify the asymptotic theory.     

Robust adaptive rate-optimal testing for the white noise hypothesis

A new test is proposed for the weak white noise null hypothesis. The test is based on a new automatic selection of the order for a Box–Pierce (1970) test statistic or the test statistic of Hong (1996). The heteroskedasticity and autocorrelation-consistent (HAC) critical values from Lee (2007) are used, allowing for estimation of the error term. The data-driven order selection is tailored to detect a new class of alternatives with autocorrelation coefficients which can be o(n^−1/2)$o\left(n^{-1/2}\right)$ provided there are sufficiently many of such coefficients. A simulation experiment illustrates the good statistical properties of the test both under the weak white noise null and the alternative.     

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.     

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].     

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.     

Generalized spectral tests for the martingale difference hypothesis

This article proposes a test for the martingale difference hypothesis (MDH) using dependence measures related to the characteristic function. The MDH typically has been tested using the sample autocorrelations or in the spectral domain using the periodogram. Tests based on these statistics are inconsistent against uncorrelated non-martingales processes. Here, we generalize the spectral test of Durlauf (1991) for testing the MDH taking into account linear and nonlinear dependence. Our test considers dependence at all lags and is consistent against general pairwise nonparametric Pitman's local alternatives converging at the parametric rate n^-1/2,$n^{-1/2},$ with n$n$ the sample size. Furthermore, with our methodology there is no need to choose a lag order, to smooth the data or to formulate a parametric alternative. Our approach could be extended to specification testing of the conditional mean of possibly nonlinear models. The asymptotic null distribution of our test depends on the data generating process, so a bootstrap procedure is proposed and theoretically justified. Our bootstrap test is robust to higher order dependence, in particular to conditional heteroskedasticity. A Monte Carlo study examines the finite sample performance of our test and shows that it is more powerful than some competing tests. Finally, an application to the S&P 500 stock index and exchange rates highlights the merits of our approach.