Parameter Estimation for Periodically Stationary Time Series

Abstract. The innovations algorithm can be used to obtain parameter estimates for periodically stationary time series models. In this paper, we compute the asymptotic distribution for these estimates in the case, where the innovations have a finite fourth moment. These asymptotic results are useful to determine which model parameters are significant. In the process, we also develop asymptotics for the Yule–Walker estimates.     

Maximum likelihood estimation for nearly non‐stationary stable autoregressive processes

The maximum likelihood estimate (MLE) of the autoregressive coefficient of a near‐unit root autoregressive process Y_t = ρ_nY_t?1 + ?_t with α‐stable noise {?_t} is studied in this paper. Herein ρ_n = 1 ? γ/n, γ ≥ 0 is a constant, Y₀ is a fixed random variable and ε_t is an α‐stable random variable with characteristic function φ(t,θ) for some parameter θ. It is shown that when 0 < α < 1 or α > 1 and E?₁ = 0, the limit distribution of the MLE of ρ_n and θ are mixtures of a stable process and Gaussian processes. On the other hand, when α > 1 and E?₁ ≠ 0, the limit distribution of the MLE of ρ_n and θ are normal. A Monte Carlo simulation reveals that the MLE performs better than the usual least squares procedures, particularly for the case when the tail index α is less than 1.     

Stationary bootstrapping for non‐parametric estimator of nonlinear autoregressive model

We consider stationary bootstrap approximation of the non‐parametric kernel estimator in a general kth‐order nonlinear autoregressive model under the conditions ensuring that the nonlinear autoregressive process is a geometrically Harris ergodic stationary Markov process. We show that the stationary bootstrap procedure properly estimates the distribution of the non‐parametric kernel estimator. A simulation study is provided to illustrate the theory and to construct confidence intervals, which compares the proposed method favorably with some other bootstrap methods.     

Non‐stationary autoregressive processes with infinite variance

Consider an AR(p) process , where {?_t} is a sequence of i.i.d. random variables lying in the domain of attraction of a stable law with index 0<α<2. This time series {Y_t} is said to be a non‐stationary AR(p) process if at least one of its characteristic roots lies on the unit circle. The limit distribution of the least squares estimator (LSE) of for {Y_t} with infinite variance innovation {?_t} is established in this paper. In particular, by virtue of the result of Kurtz and Protter (1991) of stochastic integrals, it is shown that the limit distribution of the LSE is a functional of integrated stable process. Simulations for the estimator of β and its limit distribution are also given.     

Wavelet‐Based Tests for Comparing Two Time Series with Unequal Lengths

Test procedures for assessing whether two stationary and independent time series with unequal lengths have the same spectral density (or same auto‐covariance function) are investigated. A new test statistic is proposed based on the wavelet transform. It relies on empirical wavelet coefficients of the logarithm of two spectral densities' ratio. Under the null hypothesis that two spectral densities are the same, the asymptotic normal distribution of the empirical wavelet coeffcients is derived. Furthermore, these empirical wavelet coefficients are asymptotically uncorrelated. A test statistic is proposed based on these results. The performance of the new test statistic is compared to several recent test statistics, with respect to their exact levels and powers. Simulation studies show that our proposed test is very comparable to the current test statistics in most cases. The main advantage of our proposed test statistic is that it is constructed very simply and is easy to implement.     

Robust Wilcoxon‐Type Estimation of Change‐Point Location Under Short‐Range Dependence

We introduce a robust estimator of the location parameter for the change‐point in the mean based on Wilcoxon statistic and establish its consistency for L₁ near‐epoch dependent processes. It is shown that the consistency rate depends on the magnitude of the change. A simulation study is performed to evaluate the finite sample properties of the Wilcoxon‐type estimator under Gaussianity as well as under heavy‐tailed distributions and disturbances by outliers, and to compare it with a CUSUM‐type estimator. It shows that the Wilcoxon‐type estimator is equivalent to the CUSUM‐type estimator under Gaussianity but outperforms it in the presence of heavy tails or outliers in the data.     

On the Spectral Density of the Wavelet Coefficients of Long‐Memory Time Series with Application to the Log‐Regression Estimation of the Memory Parameter

Abstract. In recent years, methods to estimate the memory parameter using wavelet analysis have gained popularity in many areas of science. Despite its widespread use, a rigorous semi‐parametric asymptotic theory, comparable with the one developed for Fourier methods, is still lacking. In this article, we adapt to the wavelet setting, the classical semi‐parametric framework introduced by Robinson and his co‐authors for estimating the memory parameter of a (possibly) non‐stationary process. Our results apply to a class of wavelets with bounded supports, which include but are not limited to Daubechies wavelets. We derive an explicit expression of the spectral density of the wavelet coefficients and show that it can be approximated, at large scales, by the spectral density of the continuous‐time wavelet coefficients of fractional Brownian motion. We derive an explicit bound for the difference between the spectral densities. As an application, we obtain minimax upper bounds for the log‐scale regression estimator of the memory parameter for a Gaussian process and we derive an explicit expression of its asymptotic variance.     

A p‐Order signed integer‐valued autoregressive (SINAR(p)) model

In this article, we propose an extension of integer‐valued autoregressive INAR models. Using a signed version of the thinning operator, we define a larger class of ‐valued processes, called SINAR, which can have positive as well as negative correlations. Using a Markov chain method, conditions for stationarity and the existence of moments are investigated. In particular, it is shown that the autocorrelation function of any real‐valued AR process can be recovered with a SINAR process, which improves INAR modeling.     

Wasserstein autoregressive models for density time series

Data consisting of time-indexed distributions of cross-sectional or intraday returns have been extensively studied in finance, and provide one example in which the data atoms consist of serially dependent probability distributions. Motivated by such data, we propose an autoregressive model for density time series by exploiting the tangent space structure on the space of distributions that is induced by the Wasserstein metric. The densities themselves are not assumed to have any specific parametric form, leading to flexible forecasting of future unobserved densities. The main estimation targets in the order-p Wasserstein autoregressive model are Wasserstein autocorrelations and the vector-valued autoregressive parameter. We propose suitable estimators and establish their asymptotic normality, which is verified in a simulation study. The new order-p Wasserstein autoregressive model leads to a prediction algorithm, which includes a data driven order selection procedure. Its performance is compared to existing prediction procedures via application to four financial return data sets, where a variety of metrics are used to quantify forecasting accuracy. For most metrics, the proposed model outperforms existing methods in two of the data sets, while the best empirical performance in the other two data sets is attained by existing methods based on functional transformations of the densities.     

A negative binomial integer‐valued GARCH model

This article discusses the modelling of integer‐valued time series with overdispersion and potential extreme observations. For the problem, a negative binomial INGARCH model, a generalization of the Poisson INGARCH model, is proposed and stationarity conditions are given as well as the autocorrelation function. For estimation, we present three approaches with the focus on the maximum likelihood approach. Some results from numerical studies are presented and indicate that the proposed methodology performs better than the Poisson and double Poisson model‐based methods.     

Asymptotic self‐similarity and wavelet estimation for long‐range dependent fractional autoregressive integrated moving average time series with stable innovations

Abstract. Methods for parameter estimation in the presence of long‐range dependence and heavy tails are scarce. Fractional autoregressive integrated moving average (FARIMA) time series for positive values of the fractional differencing exponent d can be used to model long‐range dependence in the case of heavy‐tailed distributions. In this paper, we focus on the estimation of the Hurst parameter H = d + 1/α for long‐range dependent FARIMA time series with symmetric α‐stable (1 < α < 2) innovations. We establish the consistency and the asymptotic normality of two types of wavelet estimators of the parameter H. We do so by exploiting the fact that the integrated series is asymptotically self‐similar with parameter H. When the parameter α is known, we also obtain consistent and asymptotically normal estimators for the fractional differencing exponent d = H ? 1/α. Our results hold for a larger class of causal linear processes with stable symmetric innovations. As the wavelet‐based estimation method used here is semi‐parametric, it allows for a more robust treatment of long‐range dependent data than parametric methods.     

A test for second‐order stationarity of a time series based on the discrete Fourier transform

We consider a zero mean discrete time series, and define its discrete Fourier transform (DFT) at the canonical frequencies. It can be shown that the DFT is asymptotically uncorrelated at the canonical frequencies if and only if the time series is second‐order stationary. Exploiting this important property, we construct a Portmanteau type test statistic for testing stationarity of the time series. It is shown that under the null of stationarity, the test statistic has approximately a chi‐square distribution. To examine the power of the test statistic, the asymptotic distribution under the locally stationary alternative is established. It is shown to be a generalized non‐central chi‐square, where the non‐centrality parameter measures the deviation from stationarity. The test is illustrated with simulations, where is it shown to have good power.     

Limit theory for a general class of GARCH models with just barely infinite variance

In this article, limit theory is established for a general class of generalized autoregressive conditional heteroskedasticity models given by ?_t = σ_tη_t and σ_t = f (σ_t?1, σ_t?2,…, σ_t?p, ?_t?1, ?_t?2,…, ?_t?q), when {?_t} is a process with just barely infinite variance, that is, {?_t} is a process with infinite variance but in the domain of normal attraction. In particular, we show that under some regular conditions, converges weakly to a Gaussian process. Applications of the asymptotic results to statistical inference, such as unit root test and sample autocorrelation, are also investigated. The obtained result fills in a gap between the classical infinite variance and finite variance in the literature. Further, when applying our limiting result to Dickey–Fuller (DF) test in a unit root model with integrated generalized autoregressive conditional heteroskedasticity (IGARCH) errors, it just confirms the simulation result of Kourogenis and Pittis (2008) that the DF statistics with IGARCH errors converges in law to the standard DF distribution.     

Central limit theorems for nonparametric estimators with real‐time random variables

In this article, asymptotic theories for nonparametric methods are studied when they are applied to real‐time data. In particular, we derive central limit theorems for nonparametric density and regression estimators. For this we formally introduce a sequence of real‐time random variables indexed by a parameter related to fine gridding of time domain (or fine discretization). Our results show that the impact of fine gridding is greater in the density estimation case in the sense that strong dependence due to fine gridding severely affects the major strength of nonparametric density estimator (or its data‐adaptive property). In addition, we discuss some issues about nonparametric regression model with fine gridding of time domain.     

A Robbins–Monro Algorithm for Non‐Parametric Estimation of NAR Process with Markov Switching: Consistency

We approach the problem of non‐parametric estimation for autoregressive Markov switching processes. In this context, the Nadaraya–Watson‐type regression functions estimator is interpreted as a solution of a local weighted least‐square problem, which does not admit a closed‐form solution in the case of hidden Markov switching. We introduce a non‐parametric recursive algorithm to approximate the estimator. Our algorithm restores the missing data by means of a Monte Carlo step and estimates the regression function via a Robbins–Monro step. We prove that non‐parametric autoregressive models with Markov switching are identifiable when the hidden Markov process has a finite state space. Consistency of the estimator is proved using the strong α‐mixing property of the model. Finally, we present some simulations illustrating the performances of our non‐parametric estimation procedure.     

Optimal Rate of Convergence for Empirical Quantiles and Distribution Functions for Time Series

Given a stationary sequence , we are interested in the rate of convergence in the central limit theorem of the empirical quantiles and the empirical distribution function. Under a general notion of weak dependence, we show a Berry–Esseen result with optimal rate n^?1/2. The setup includes many prominent time series models, such as functions of ARMA or (augmented) GARCH processes. In this context, optimal Berry–Esseen rates for empirical quantiles appear to be novel.     

Statistical tests for a single change in mean against long‐range dependence

Statistical tests are introduced for distinguishing between short‐range dependent time series with a single change in mean, and long‐range dependent time series, with the former making the null hypothesis. The tests are based on estimation of the self‐similarity parameter after removing the change in mean from the series. The focus is on the GPH (Geweke and Porter‐Hudak, 1983) and local Whittle estimation methods in the spectral domain. Theoretical properties of the resulting estimators are established when testing for a single change in mean, and small sample properties of the tests are examined in simulations. The introduced tests improve on the BHKS ( Berkes et al., 2006 ) test which is the only other available test for the considered problem. It is argued that the BHKS test has a low power against long‐range dependence alternatives and that this happens because the BHKS test statistic involves estimation of the long‐run variance. The BHKS test could be improved readily by considering its R/S‐like regression version which estimates the self‐similarity parameter and which does not involve the long‐run variance. Yet better alternatives are to use more powerful estimation methods (such as GPH or local Whittle) and lead to the tests introduced here.     

Semiparametric Estimation in Time‐Series Regression with Long‐Range Dependence

Abstract. We consider semiparametric estimation in time‐series regression in the presence of long‐range dependence in both the errors and the stochastic regressors. A central limit theorem is established for a class of semiparametric frequency domain‐weighted least squares estimates, which includes both narrow‐band ordinary least squares and narrow‐band generalized least squares as special cases. The estimates are semiparametric in the sense that focus is on the neighbourhood of the origin, and only periodogram ordinates in a degenerating band around the origin are used. This setting differs from earlier studies on time‐series regression with long‐range dependence, where a fully parametric approach has been employed. The generalized least squares estimate is infeasible when the degree of long‐range dependence is unknown and must be estimated in an initial step. In that case, we show that a feasible estimate which has the same asymptotic properties as the infeasible estimate, exists. By Monte Carlo simulation, we evaluate the finite‐sample performance of the generalized least squares estimate and the feasible estimate.