Improved functional portmanteau tests

Functional time series is a popular method of forecasting in functional data analysis. The Box-Jenkins methodology for model building, with the aim of forecasting, includes three iterative steps of model identification, parameter estimation and diagnostic checking. Portmanteau tests are one of the most popular diagnostic checking tools. In particular, they are applied to find if the residuals of the fitted model are white noise. Gabrys and Kokoszka [Portmanteau test of independence for functional observations. J Am Stat Assoc. 2007;102(480):1338–1348.] proposed a portmanteau test of independence for functional observation based on Box and Pierce's statistic. Their statistic is too sensitive to the lag value, specially when the sample size is small. Here, two modifications of Gabrys and Kokoszka statistic are presented, which have superior properties in small samples. The efficiency of the modified statistics is demonstrated through a simulation study.     

On consistent testing for serial correlation in seasonal time series models

The author considers serial correlation testing in seasonal time series models. He proposes a test statistic based on a spectral approach. Many tests of this type rely on kernel-based spectral density estimators that assign larger weights to low order lags than to high ones. Under seasonality, however, large autocorrelations may occur at seasonal lags that classical kernel estimators cannot take into account. The author thus proposes a test statistic that relies on the spectral density estimator of Shin (2004), whose weighting scheme is more adapted to this context. The distribution of his test statistic is derived under the null hypothesis and he studies its behaviour under fixed and local alternatives. He establishes the consistency of the test under a general fixed alternative. He also makes recommendations for the choice of the smoothing parameters. His simulation results suggest that his test is more powerful against seasonality than alternative procedures based on classical weighting schemes. He illustrates his procedure with monthly statistics on employment among young Americans.     

Testing time series for interpolability and whiteness

We propose a test to decide if a time series is represented by its linear interpolator better than by its mean value. The same test can be employed to decide if a time series has to be considered white noise. The test is based on a new estimate of the index of linear determinism (Battaglia, 1983, Inverse autocovariances and a measure of linear determinism for a stationary process, J. Time Series Anal. 4, 79-87) and its asymptotic distribution is derived. Comparison with the popular Ljung-Box portmanteau test has been performed based on both asymptotic power and a simulation experiment. The new test     

Portmanteau tests for linearity of stationary time series

This article considers the problem of testing for linearity of stationary time series. Portmanteau tests are discussed which are based on generalized correlations of residuals from a linear model (that is, autocorrelations and cross-correlations of different powers of the residuals). The finite-sample properties of the tests are assessed by means of Monte Carlo experiments. The tests are applied to 100 time series of stock returns.     

Non parametric portmanteau tests for detecting non linearities in high dimensions

ABSTRACT

We propose two non parametric portmanteau test statistics for serial dependence in high dimensions using the correlation integral. One test depends on a cutoff threshold value, while the other test is freed of this dependence. Although these tests may each be viewed as variants of the classical Brock, Dechert, and Scheinkman (BDS) test statistic, they avoid some of the major weaknesses of this test. We establish consistency and asymptotic normality of both portmanteau tests. Using Monte Carlo simulations, we investigate the small sample properties of the tests for a variety of data generating processes with normally and uniformly distributed innovations. We show that asymptotic theory provides accurate inference in finite samples and for relatively high dimensions. This is followed by a power comparison with the BDS test, and with several rank-based extensions of the BDS tests that have recently been proposed in the literature. Two real data examples are provided to illustrate the use of the test procedure.     

Modeling nonlinear time series with local mixtures of generalized linear models

The authors consider a novel class of nonlinear time series models based on local mixtures of regressions of exponential family models, where the covariates include functions of lags of the dependent variable. They give conditions to guarantee consistency of the maximum likelihood estimator for correctly specified models, with stationary and nonstationary predictors. They show that consistency of the maximum likelihood estimator still holds under model misspecification. They also provide probabilistic results for the proposed model when the vector of predictors contains only lags of transformations of the modeled time series. They illustrate the consistency of the maximum likelihood estimator and the probabilistic properties via Monte Carlo simulations. Finally, they present an application using real data.     

On testing for causality in variance between two multivariate time series

Verifying the existence of a relationship between two multivariate time series represents an important consideration. In this article, the procedure developed by Cheung and Ng [A causality-in-variance test and its application to financial market prices, J. Econom. 72 (1996), pp. 33–48] designed to test causality in variance for univariate time series is generalized in several directions. A first approach proposes test statistics based on residual cross-covariance matrices of squared (standardized) residuals and cross products of (standardized) residuals. In a second approach, transformed residuals are defined for each residual vector time series, and test statistics are constructed based on the cross-correlations of these transformed residuals. Test statistics at individual lags and portmanteau-type test statistics are developed. Conditions are given under which the new test statistics converge in distribution towards chi-square distributions. The proposed methodology can be used to determine the directions of causality in variance, and appropriate test statistics are presented. Monte Carlo simulation results show that the new test statistics offer satisfactory empirical properties. An application with two bivariate financial time series illustrates the methods.     

Priors and component structures in autoregressive time series models

New approaches to prior specification and structuring in autoregressive time series models are introduced and developed. We focus on defining classes of prior distributions for parameters and latent variables related to latent components of an autoregressive model for an observed time series. These new priors naturally permit the incorporation of both qualitative and quantitative prior information about the number and relative importance of physically meaningful components that represent low frequency trends, quasi-periodic subprocesses and high frequency residual noise components of observed series. The class of priors also naturally incorporates uncertainty about model order and hence leads in posterior analysis to model order assessment and resulting posterior and predictive inferences that incorporate full uncertainties about model order as well as model parameters. Analysis also formally incorporates uncertainty and leads to inferences about unknown initial values of the time series, as it does for predictions of future values. Posterior analysis involves easily implemented iterative simulation methods, developed and described here. One motivating field of application is climatology, where the evaluation of latent structure, especially quasi-periodic structure, is of critical importance in connection with issues of global climatic variability. We explore the analysis of data from the southern oscillation index, one of several series that has been central in recent high profile debates in the atmospheric sciences about recent apparent trends in climatic indicators.     

Simple diagnostic tools for inverstigating linear trends in time series

This paper presents a simple diagnostic tool for time series. Based on a coefficient α that veries between 1 and 0, the tool measures the approximation of a time series to an arithmetic progression (i.e., a linear function of time). The proposed α is based on the ratio of the average squared second difference to the average squared first difference of the ginven series. As such, α reduces to the Von Neumann ratio η of the series of first differences, namely, α = 1-η/4. For an arithmetic progression α = 1, and deviations therefrom cause it to decrease. Unlike the correlation coefficient (between the entries and the indics), α is sensitive to local, or piecewise, linearity. Here α is evaluated for an assortment of simple time series models such as random walk, AR(1) and MA(1). Large-sample distribution yields a number of commonly used stochastic models including non-normal process. For most standard deterministic and stochastic models, α stabilizes as n approaches infinity, and provides a statistic that is capable of distinguishing between many different standard random and deterministic models. A further measure τ, which together with α distinguisches between random walks and deterministic trend plus i.i.d., is also suggested. Some examples based on empirical data are also studied.     

Bootstrap tests for variance components in generalized linear mixed models

In many applications of generalized linear mixed models to clustered correlated or longitudinal data, often we are interested in testing whether a random effects variance component is zero. The usual asymptotic mixture of chi‐square distributions of the score statistic for testing constrained variance components does not necessarily hold. In this article, the author proposes and explores a parametric bootstrap test that appears to be valid based on its estimated level of significance under the null hypothesis. Results from a simulation study indicate that the bootstrap test has a level much closer to the nominal one while the asymptotic test is conservative, and is more powerful than the usual asymptotic score test based on a mixture of chi‐squares. The proposed bootstrap test is illustrated using two sets of real‐life data obtained from clinical trials. The Canadian Journal of Statistics © 2009 Statistical Society of Canada     

Nonlinearity tests in time series analysis

Summary It is widely recognized that the class of ARIMA models may fail to capture fully the dynamics of real phenomena since these are often characterized by strong nonlinear components. Thus, it is important that any preliminary analysis (or evaluation of model adequacy) includes a check on the linearity of the generating process. The paper reviews recent developments in the theory of testing nonlinearity in time series analysis.     

Structural changes in autoregressive models for binary time series

We study autoregressive models for binary time series with possible changes in their parameters. A procedure for detection and testing of a single change is suggested. The limiting behavior of the test statistic is derived. The performance of the test is analyzed under the null hypothesis as well as under different alternatives via a simulation study. Application of the method to a real data set on US recession is provided as an illustration.     

A new multicollinearity diagnostic for generalized linear models

We propose a new collinearity diagnostic tool for generalized linear models. The new diagnostic tool is termed the weighted variance inflation factor (WVIF) behaving exactly the same as the traditional variance inflation factor in the context of regression diagnostic, given data matrix normalized. Compared to the use of condition number (CN), WVIF shows more reliable information on how severe the situation is, when data collinearity does exist. An alternative estimator, a by-product of the new diagnostic, outperforms the ridge estimator in the presence of data collinearity in both aspects of WVIF and CN. Evidences are given through analyzing various real-world numerical examples.     

Indirect inference for fractional time series models

In this paper we present an indirect estimation procedure for (ARFIMA) fractional time series models.The estimation method is based on an ‘incorrect’criterion which does not directly provide a consistent estimator of the parameters of interest,but leads to correct inference by using simulations.

The main steps are the following. First,we consider an auxiliary model which can be easily estimated.Specifically,we choose the finite lag Autoregressive model.Then, this is estimated on the observations and simulated values drawn from the ARFIMA model associated with a given value of the parameters of interest.Finally,the latter is calibrated in order to obtain close values of the two estimators of the auxiliary parameters.

In this article,we describe the estimation procedure and compare the performance of the indirect estimator with some alternative estimators based on the likelihood function by a Monte Carlo study.     

Efficient estimation for time series following generalized linear models

In this paper, we consider James–Stein shrinkage and pretest estimation methods for time series following generalized linear models when it is conjectured that some of the regression parameters may be restricted to a subspace. Efficient estimation strategies are developed when there are many covariates in the model and some of them are not statistically significant. Statistical properties of the pretest and shrinkage estimation methods including asymptotic distributional bias and risk are developed. We investigate the relative performances of shrinkage and pretest estimators with respect to the unrestricted maximum partial likelihood estimator (MPLE). We show that the shrinkage estimators have a lower relative mean squared error as compared to the unrestricted MPLE when the number of significant covariates exceeds two. Monte Carlo simulation experiments were conducted for different combinations of inactive covariates and the performance of each estimator was evaluated in terms of its mean squared error. The practical benefits of the proposed methods are illustrated using two real data sets.     

Zero-truncated compound Poisson integer-valued GARCH models for time series

Starting from the compound Poisson INGARCH models, we introduce in this paper a new family of integer-valued models suitable to describe count data without zeros that we name zero-truncated CP-INGARCH processes. For such class of models, a probabilistic study concerning moments existence, stationarity and ergodicity is developed. The conditional quasi-maximum likelihood method is introduced to consistently estimate the parameters of a wide zero-truncated compound Poisson subclass of models. The conditional maximum likelihood method is also used to estimate the parameters of ZTCP-INGARCH processes associated with well-specified conditional laws. A simulation study that compares some of those estimators and illustrates their finite distance behaviour as well as a real-data application conclude the paper.     

Change detection in linear regression with time series errors

A new test for detecting a change in linear regression parameters assuming a general weakly dependent error structure is given. It extends earlier methods based on cumulative sums assuming independent errors. The novelty is in the new standardization method and in smoothing when the time series is dominated by high frequencies. Simulations show the excellent performance of the test. Examples are taken from environmental applications. The algorithm is easy to implement. Testing for multiple changes can be done by segmentation. The Canadian Journal of Statistics 38:65–79; 2010 © 2009 Statistical Society of Canada     

Goodness-of-fit tests of generalized linear mixed models for repeated ordinal responses

Categorical longitudinal data are frequently applied in a variety of fields, and are commonly fitted by generalized linear mixed models (GLMMs) and generalized estimating equations models. The cumulative logit is one of the useful link functions to deal with the problem involving repeated ordinal responses. To check the adequacy of the GLMMs with cumulative logit link function, two goodness-of-fit tests constructed by the unweighted sum of squared model residuals using numerical integration and bootstrap resampling technique are proposed. The empirical type I error rates and powers of the proposed tests are examined by simulation studies. The ordinal longitudinal studies are utilized to illustrate the application of the two proposed tests.     

Change-in-mean problem for long memory time series models with applications

Test and estimation procedures for detecting a change in the mean are proposed in infinite moving average long memory time series models. The asymptotic properties of the test statistics and the change-point estimators are investigated. The method is illustrated through the analysis of real data sets from econometrics and climatology.