Markov regime switching in mean and in fractional integration parameter

We propose a specific general Markov-regime switching estimation both in the long memory parameter d and the mean of a time series. We employ Viterbi algorithm that combines the Viterbi procedures in two state Markov-switching parameter estimation. It is well-known that existence of mean break and long memory in time series can be easily confused with each other in most cases. Thus, we aim at observing the deviation and interaction of mean and d estimates for different cases. A Monte Carlo experiment reveals that the finite sample performance of the proposed algorithm for a simple mixture model of Markov-switching mean and d changes with respect to the fractional integrating parameters and the mean values for the two regimes.     

On the usability of the fluctuation test statistic to identify multiple cointegration break points

The fluctuation test suggested by Hansen and Johansen [Some tests for parameter constancy in cointegrated VAR models, Econometrics J. 2 (1999), pp. 306–333] intends to distinguish between the presence of zero and one break in cointegration relations. In this article, we provide evidence by Monte Carlo simulations that it also serves as a graphical device to detect even multiple break locations. It suffices to consider a simplified and easy-to-implement version of the original fluctuation test. Its break detection performance depends on the sign of change in cointegration parameters and the break height. The sign issue can be approached successfully by a backward application of the test statistic. If breaks are observable, the break locations are detected at the true location on average. We apply the graphical procedure to assess the cointegration of bond yields of Spain, Italy and Portugal with German yields for the period 1995–2013 which is surprisingly supported by the trace test. However, the recursive cointegration approach shows that a stable relationship with German yields is only present for sub-periods between the introduction of the Euro and the global financial crisis which is in line with expectations. The statistical robustness of these results is supported by a forward and backward application of the cointegration breakdown test by Andrews and Kim [Tests for cointegration breakdown over a short time period, J. Bus. Econom. Stat. 24 (2006), pp. 379–394].     

Goodness-of-fit tests for exponentiality based on a loss-of-memory type functional equation

Three goodness-of-fit tests for exponentiality based on the functional equation characterization 1?F(2x)=[1?F(x)]² for every x?0 are proposed and shown to compare well to several popular tests against common alternative cdf's. Small-sample critical values for α=0.10,0.05 are developed for the two superior test statistics and the asymptotic null-distributions are characterized.     

THE BEHAVIOR OF HEGY TESTS FOR QUARTERLY TIME SERIES WITH SEASONAL MEAN SHIFTS

ABSTRACT

This paper studies the behavior of the HEGY statistics for quarterly data, for seasonal autoregressive unit roots, when the analyzed time series is deterministic seasonal stationary but exhibits a change in the seasonal pattern. We analyze also the HEGY test for the nonseasonal unit root. the data generation process being trend stationary too. Our results show that when the break magnitudes are finite, the HEGY test statistics are not asymptotically biased toward the nonrejection of the seasonal and nonseasonal unit root hypotheses. However, the finite sample power properties may be substantially affected, the behavior of the tests depending on the type of the break.     

A test for fractional cointegration using the sieve bootstrap

The concept of fractional cointegration (Cheung and Lai in J Bus Econ Stat 11:103–112, 1993) has been introduced to generalize traditional cointegration (Engle and Granger in Econometrica 55:251–276, 1987) to the long memory framework. In this work we propose a test for fractional cointegration with the sieve bootstrap and compare by simulations the performance of our proposal with other semiparametric methods existing in literature: the three steps technique of Marinucci and Robinson (J Econom 105:225–247, 2001) and the procedure to determine the fractional cointegration rank of Robinson and Yajima (J Econom 106:217–241, 2002).     

Test of hypotheses in panel data models when the regressor and disturbances are possibly non-stationary

This paper considers the problem of hypothesis testing in a simple panel data regression model with random individual effects and serially correlated disturbances. Following Baltagi et al. (Econom. J. 11:554–572, 2008), we allow for the possibility of non-stationarity in the regressor and/or the disturbance term. While Baltagi et al. (Econom. J. 11:554–572, 2008) focus on the asymptotic properties and distributions of the standard panel data estimators, this paper focuses on testing of hypotheses in this setting. One important finding is that unlike the time-series case, one does not necessarily need to rely on the “super-efficient” type AR estimator by Perron and Yabu (J. Econom. 151:56–69, 2009) to make an inference in the panel data. In fact, we show that the simple t-ratio always converges to the standard normal distribution, regardless of whether the disturbances and/or the regressor are stationary.     

THE CONVERGENCE RATE OF SEQUENTIAL FIXED-WIDTH CONFIDENCE INTERVALS FOR REGULAR FUNCTIONALS

Let X₁X₂,.be i.i.d. random variables and let U_n= (n r)^-1S?^(n,r) h (X_i1,., X_ir,) be a U-statistic with EU_n= v, v unknown. Assume that g(X₁) =E[h(X₁,.,X_r) - v |X₁]has a strictly positive variance s?². Further, let a be such that φ(a) - φ(-a) =α for fixed α, 0 < α < 1, where φ is the standard normal d.f., and let S²_n be the Jackknife estimator of n Var U_n. Consider the stopping times N(d)= min {n: S²_n: + n^-1≤²a^-2},d > 0, and a confidence interval for v of length 2d,of the form I_n,d= [U_n,-d, Un + d]. We assume that Var U_n is unknown, and hence, no fixed sample size method is available for finding a confidence interval for v of prescribed width 2d and prescribed coverage probability α Turning to a sequential procedure, let I_N(d),d be a sequence of sequential confidence intervals for v. The asymptotic consistency of this procedure, i.e. lim_{d → 0}P(v ∈ I_N(d),d)=α follows from Sproule (1969). In this paper, the rate at which |P(v ∈ I_N(d),d) converges to α is investigated. We obtain that |P(v ∈ I_N(d),d) - α| = 0 (d^{1/2-(1+k)/2(1+m)}), d → 0, where K = max {0,4 - m}, under the condition that E|h(X₁, X_r)|^m < ∞m > 2. This improves and extends recent results of Ghosh & DasGupta (1980) and Mukhopadhyay (1981).     

Infant mortality rates: time trends and fractional integration

This paper examines the existence of time trends in the infant mortality rates in a number of countries in the twentieth century. We test for the presence of deterministic trends by adopting a linear model for the log-transformed data. Instead of assuming that the error term is a stationary I(0), or alternatively, a non-stationary I(1) process, we allow for the possibility of fractional integration and hence for a much greater degree of flexibility in the dynamic specification of the series. Indeed, once the linear trend is removed, all series appear to be I(d) with 0<d<1, implying long-range dependence. As expected, the time trend coefficients are significantly negative, although of a different magnitude from those obtained assuming integer orders of differentiation.     

Improved inference for a linear mixed-effects model when the subpopulation effects are clustered

We provide Bayesian methodology to relax the assumption that all subpopulation effects in a linear mixed-effects model have, after adjustment for covariates, a common mean. We expand the model specification by assuming that the m subpopulation effects are allowed to cluster into d groups where the value of d, 1?d?m, and the composition of the d groups are unknown, a priori. Specifically, for each partition of the m effects into d groups we only assume that the subpopulation effects in each group are exchangeable and are independent across the groups. We show that failure to take account of this clustering, as with the customary method, will lead to serious errors in inference about the variances and subpopulation effects, but the proposed, expanded, model leads to appropriate inferences. The efficacy of the proposed method is evaluated by contrasting it with both the customary method and use of a Dirichlet process prior. We use data from small area estimation to illustrate our method.     

On tail parameter estimation in certain point process models

Let μ(ds, dx) denote Poisson random measure with intensity dsG(dx) on (0, ∞) × (0, ∞), for a measure G(dx) with tails varying regularly at ∞. We deal with estimation of index of regular variation α and weight parameter ξ if the point process is observed in certain windows K_n = [0, T_n] × [Y_n, ∞), where Y_n → ∞ as n → ∞. In particular, we look at asymptotic behaviour of the Hill estimator for α. In certain submodels, better estimators are available; they converge at higher speed and have a strong optimality property. This is deduced from the parametric case G(dx) = ξαx^−α−1 dx via a neighbourhood argument in terms of Hellinger distances.     

Monitoring Change in Persistence Against the Null of Difference-Stationarity in Infinite Variance Observations

In this article, we propose a moving kernel-weighted variance ratio statistic to monitor persistence change in infinite variance observations. We focus on I(1) to I(0) persistence change for sequences in the domain of attraction of a stable law and local-to-finite variance sequences. The null distribution of the monitoring statistic and its consistency are proved. In particular, a bootstrap procedure is proposed to determine the critical values for the derived asymptotic distribution depends on unknown tail index. The small sample performances of proposed monitoring procedure are illustrated by both simulation and application to a high frequency financial data.     

A comparison of estimation methods in non-stationary ARFIMA processes

This paper reports an extensive Monte Carlo simulation study based on six estimators for the long memory fractional parameter when the time series is non-stationary, i.e., ARFIMA(p, d, q) process for d?>?0.5. Parametric and semiparametric methods are compared. In addition, the effect of the parameter estimation is investigated for small and large sample sizes and non-Gaussian error innovations. The methodology is applied to a well known data set, the so-called UK short interest rates.     

Admissible minimax model selection rule

This paper considers the problem of choosing one between the simple model N(0, I_d) and the full model N(0 I_d) based on the observation X from N(θ I_d) where X, θεR^d, 0 is the null vector in R^d and I_d is the d×d identity matrix. It is shown that the selection rule which chooses the full model if |x| > a_o , for some a₀ > 0 and the simple model otherwise is an admissible minimax model selection rule relative to a loss function which takes into account both inaccuracy and complexity.     

A generalized ARFIMA process with Markov-switching fractional differencing parameter

We propose a general class of Markov-switching-ARFIMA (MS-ARFIMA) processes in order to combine strands of long memory and Markov-switching literature. Although the coverage of this class of models is broad, we show that these models can be easily estimated with the Durbin–Levinson–Viterbi algorithm proposed. This algorithm combines the Durbin–Levinson and Viterbi procedures. A Monte Carlo experiment reveals that the finite sample performance of the proposed algorithm for a simple mixture model of Markov-switching mean and ARFIMA(1, d, 1) process is satisfactory. We apply the MS-ARFIMA models to the US real interest rates and the Nile river level data, respectively. The results are all highly consistent with the conjectures made or empirical results found in the literature. Particularly, we confirm the conjecture in Beran and Terrin [J. Beran and N. Terrin, Testing for a change of the long-memory parameter. Biometrika 83 (1996), pp. 627–638.] that the observations 1 to about 100 of the Nile river data seem to be more independent than the subsequent observations, and the value of differencing parameter is lower for the first 100 observations than for the subsequent data.     

Minimal 2-coverings of a finite affine space based on GF(2)

Let EG(m, 2) denote the m-dimensional finite Euclidean space (or geometry) based on GF(2), the finite field with elements 0 and 1. Let T be a set of points in this space, then T is said to form a q-covering (where q is an integer satisfying 1?q?m) of EG(m, 2) if and only if T has a nonempty intersection with every (m-q)-flat of EG(m, 2). This problem first arose in the statistical context of factorial search designs where it is known to have very important and wide ranging applications. Evidently, it is also useful to study this from the purely combinatorial point of view. In this paper, certain fundamental studies have been made for the case when q=2. Let N denote the size of the set T. Given N, we study the maximal value of m.     

The group inverse of finite homogeneous QBD processes

Generalized inverses of I?P, where P is a stochastic matrix, play an important role in the theory of Markov chains. In particular, the group inverse (I?P)^# has a probabilistic interpretation and is well suited for algorithmic implementation. We determine (I?P)^# for finite homogeneous quasi-birth-and-death (QBD) processes by exploiting both the structure of the process and the probabilistic properties of the group inverse.     

The Estimation and Testing of the Cointegration Order Based on the Frequency Domain

ABSTRACT

This article proposes a method to estimate the degree of cointegration in bivariate series and suggests a test statistic for testing noncointegration based on the determinant of the spectral density matrix for the frequencies close to zero. In the study, series are assumed to be I(d), 0 < d ? 1, with parameter d supposed to be known. In this context, the order of integration of the error series is I(d ? b), b ∈ [0, d]. Besides, the determinant of the spectral density matrix for the dth difference series is a power function of b. The proposed estimator for b is obtained here performing a regression of logged determinant on a set of logged Fourier frequencies. Under the null hypothesis of noncointegration, the expressions for the bias and variance of the estimator were derived and its consistency property was also obtained. The asymptotic normality of the estimator, under Gaussian and non-Gaussian innovations, was also established. A Monte Carlo study was performed and showed that the suggested test possesses correct size and good power for moderate sample sizes, when compared with other proposals in the literature. An advantage of the method proposed here, over the standard methods, is that it allows to know the order of integration of the error series without estimating a regression equation. An application was conducted to exemplify the method in a real context.     

Empirically feasible solutions and explicit dynamics for rational expectation models

Fairly general rational expectation (RE) models are proved to have linear vector autoregressive moving average (VARMA) models as their reduced forms, using Muth's method of undetermined coefficients (MUC). An advantage of using VARMAs is that RE estimation can benefit from well-developed theory of solutions of stochastic difference equations and computer packages for VARMA and state space (SS) models. For example, we report theoretical derivations of explicit dynamics associated with the RE structure for Muth's and Lucas-Sargent-Wallace's models, suggesting generalizations and new solutions. We illustrate an RE structure estimation by using energy prices and Fackler and Krieger's (J. Bus. Econom. Statist. 4 (1986), 71–80) study of five major macroeconomic variables.     

Double unit root tests for cross-sectionally dependent panel data

This paper proposes various double unit root tests for cross-sectionally dependent panel data. The cross-sectional correlation is handled by the projection method [P.C.B. Phillips and D. Sul, Dynamic panel estimation and homogeneity testing under cross section dependence, Econom. J. 6 (2003), pp. 217–259; H.R. Moon and B. Perron, Testing for a unit root in panels with dynamic factors, J. Econom. 122 (2004), pp. 81–126] or the subtraction method [J. Bai and S. Ng, A PANIC attack on unit roots and cointegration, Econometrica 72 (2004), pp. 1127–1177]. Pooling or averaging is applied to combine results from different panel units. Also, to estimate autoregressive parameters the ordinary least squares estimation [D.P. Hasza and W.A. Fuller, Estimation for autoregressive processes with unit roots, Ann. Stat. 7 (1979), pp. 1106–1120] or the symmetric estimation [D.L. Sen and D.A. Dickey, Symmetric test for second differencing in univariate time series, J. Bus. Econ. Stat. 5 (1987), pp. 463–473] are used, and to adjust mean functions the ordinary mean adjustment or the recursive mean adjustment are used. Combinations of different methods in defactoring to eliminate the cross-sectional dependency, integrating results from panel units, estimating the parameters, and adjusting mean functions yields various available tests for double unit roots in panel data. Simple asymptotic distributions of the proposed test statistics are derived, which can be used to find critical values of the test statistics.

We perform a Monte Carlo experiment to compare the performance of these tests and to suggest optimal tests for a given panel data. Application of the proposed tests to a real data, the yearly export panel data sets of several Latin–American countries for the past 50 years, illustrates the usefulness of the proposed tests for panel data, in that they reveal stronger evidence of double unit roots than the componentwise double unit root tests of Hasza and Fuller [Estimation for autoregressive processes with unit roots, Ann. Stat. 7 (1979), pp. 1106–1120] or Sen and Dickey [Symmetric test for second differencing in univariate time series, J. Bus. Econ. Stat. 5 (1987), pp. 463–473].     

On the Usefulness or Lack Thereof of Optimality Criteria for Structural Change Tests

Elliott and Müller (2006) considered the problem of testing for general types of parameter variations, including infrequent breaks. They developed a framework that yields optimal tests, in the sense that they nearly attain some local Gaussian power envelop. The main ingredient in their setup is that the variance of the process generating the changes in the parameters must go to zero at a fast rate. They recommended the so-called qL?L test, a partial sums type test based on the residuals obtained from the restricted model. We show that for breaks that are very small, its power is indeed higher than other tests, including the popular sup-Wald (SW) test. However, the differences are very minor. When the magnitude of change is moderate to large, the power of the test is very low in the context of a regression with lagged dependent variables or when a correction is applied to account for serial correlation in the errors. In many cases, the power goes to zero as the magnitude of change increases. The power of the SW test does not show this non-monotonicity and its power is far superior to the qL?L test when the break is not very small. We claim that the optimality of the qL?L test does not come from the properties of the test statistics but the criterion adopted, which is not useful to analyze structural change tests. Instead, we use fixed-break size asymptotic approximations to assess the relative efficiency or power of the two tests. When doing so, it is shown that the SW test indeed dominates the qL?L test and, in many cases, the latter has zero relative asymptotic efficiency.