Asymptotic Results for Spatial ARMA Models

Causal quadrantal-type spatial ARMA(p, q) models with independent and identically distributed innovations are considered. In order to select the orders (p, q) of these models and estimate their autoregressive parameters, estimators of the autoregressive coefficients, derived from the extended Yule–Walker equations are defined. Consistency and asymptotic normality are obtained for these estimators. Then, spatial ARMA model identification is considered and simulation study is given.     

Order identification for gaussian moving averages using the covariation

In the traditional Box-Jenkins procedure for fitting ARMA time series models to data, the first step is order identification. The sample autocorrelation function can be used to identify pure moving average behavior. In this paper we consider using the autocovariation function identify the order of a univariate Gaussian time series. Simulation evidence indicates the suggested method may be a superior order identification tool when at least 100 observations are taken.     

基于Boosting-ARMA的碳价预测

自回归滑动平均(ARMA)模型是最流行的预测模型之一,而模型选择却是使用ARMA进行预测的难点,尤其是当真实模型的阶数较高时,因此提出Boosting-ARMA预测算法,利用Boosting算法进行最优子集ARMA寻找,自动且高效地完成ARMA模型的识别。模拟实验显示,Boosting-ARMA优于其他方法,用新算法预测碳价实证分析发现,Boosting-ARMA算法可以获得较高的碳价预测准确性并且方便快捷。     

Theory of Bilinear Time Series Models

ABSTRACT

In this paper, we prove some theoretic properties of bilinear time series models which are extension of ARMA models. The sufficient conditions for asymptotic stationarity and ivertibility of some types of bilinear models are derived. The structural theory of discussed bilinear models is similar to that of ARMA models. For illustration, a bilinear model has been fitted to the Wolfer sunspot numbers and a substantial reduction in sum of squared residuals is obtained as comparing with Box-Jenkins ARMA model.     

Calculating the autocovariances and the likelihood for periodic V ARMA models

This paper concerns the autocovariance calculation and likelihood evaluation for periodic vector ARMA models (PV ARMA). Based on a state space representation of PV ARMA models, we derive an algorithm for computing the PV ARMA autocovariances. The proposed method computes the autocovariances for distinct seasons separately, thereby facilitating efficient calculation for models with a large period. As a result, the obtained autocovariance calculation procedure is exploited in a periodic Chandrasekhar-type filter to evaluate the exact likelihood for Gaussian PV ARMA series. Empirical evidence shows the superiority of the periodic Chandrasekhar algorithm for likelihood evaluation over the Kalman-based one.     

The effect of temporal aggregation on the estimation accuracy of ARMA models

ABSTRACT

Autoregressive Moving Average (ARMA) time series model fitting is a procedure often based on aggregate data, where parameter estimation plays a key role. Therefore, we analyze the effect of temporal aggregation on the accuracy of parameter estimation of mixed ARMA and MA models. We derive the expressions required to compute the parameter values of the aggregate models as functions of the basic model parameters in order to compare their estimation accuracy. To this end, a simulation experiment shows that aggregation causes a severe accuracy loss that increases with the order of aggregation, leading to poor accuracy.     

On the pade tablé and its relationship to the r and sarrays and arm a modeling

The Box-Jenkins method is a popular and important technique for modeling and forecasting of time series. Unfortunately the problem of determining the appropriate ARMA forecasting model (or indeed if an ARMA model holds) is a major drawback to the use of the Box-Jenkins methodology. Gray et al. (1978) and Woodward and Gray (1979) have proposed methods of estimating p and qin ARMA modeling based on the R and Sarrays that circumvent some of these modeling difficulties.

In this paper we generalize the R and S arrays by showing a relationship to Padé approximunts and then show that these arrays have a much wider application than in just determining model order. Particular non-ARMA models can be identified as well. This includes certain processes that consist of deterministic functions plus ARMA noise, indeed we believe that the combined R and S arrays are the best overall tool so fur developed for the identification of general 2nd order (not just stationary) time scries models.     

ARMA MODELS REALIZATION AND IMPULSE RESPONSES

A parallel is made between the role played by covariances in the determination of auto-regressive models and the role played by impulse responses in the determination of ARMA models.

Auto-regressive models are known to maximize the Burg-entropy under covariance constraints. Auto-regressive-moving-average models give the maximum of the Burg-entropy among processes sharing the same covariances and impulse responses up to a certain lag. Such models are constructed by iterative or algebraic methods under the different constraints.

A new recursive method of identification of the order of an ARMA model is also developed, based on the generalized reflection coefficients.     

A NOTE ON ARMA PARAMETERISATION*

This paper considers the relationship between ARMA parameterisations of models for y(t) and Ay(t), where A is invertible and y(t) is a vector time series (t = 0,±1,…). An ARMA model for the transformed series Ay(t) may have fewer parameters than a model for y(t). This paper shows that such a saving is illusory because the apparently saved parameters are exactly balanced by the number of new parameters appearing in A.     

Identification of composite (∑+II) arma models by relatively simpler models

In the real world situations, many time series are aggregates of two or more time series. An aggregation may take place due to an addition or the product or both of two or more time series. We are often interested in the study of the properties of aggregates which are, in turn, dependent on the properties of the constituent series. Motivated by this problem, the authors study in this paper the properties of models generated by the operator (Σ+II) on autoregressive-moving-average (ARMA) processes of orders (p_i,q_i), i = l→n . A few practical examples where such models have been used are given in the introduction and an illustrative numerical example is discussed at the end of the paper.     

Asymptotic Inefficiency of Mean-Correction on Parameter Estimation for a Periodic First-Order Autoregressive Model

A common practice in time series analysis is to fit a centered model to the mean-corrected data set. For stationary autoregressive moving-average (ARMA) processes, as far as the parameter estimation is concerned, fitting an ARMA model without intercepts to the mean-corrected series is asymptotically equivalent to fitting an ARMA model with intercepts to the observed series. We show that, related to the parameter least squares estimation of periodic ARMA models, the second approach can be arbitrarily more efficient than the mean-corrected counterpart. This property is illustrated by means of a periodic first-order autoregressive model. The asymptotic variance of the estimators for both approaches is derived. Moreover, empirical experiments based on simulations investigate the finite sample properties of the estimators.     

Estimation in periodic restricted EXPAR(1) models

ABSTRACT

This article is devoted to study the problem of estimation in the periodic restricted exponential autoregressive EXPAR(1) models. The estimation procedure that is used is the least-square method. Simulation studies are carried out in order to check the asymptotic properties. An application to monthly flow data for the Fraser River in British Columbia is included.     

On the identification of ARMA echelon-form models

An identification procedure for multivariate autoregressive moving average (ARMA) echelon-form models is proposed. It is based on the study of the linear dependence between rows of the Hankel matrix of serial correlations. To that end, we define a statistical test for checking the linear dependence between vectors of serial correlations. It is shown that the test statistic t?_n considered is distributed asymptotically as a finite linear combination of independent chi-square random variables with one degree of freedom under the null hypothesis, whereas under the alternative hypothesis, t?_N/N converges in probability to a positive constant. These results allow us, in particular, to compute the asymptotic probability of making a specification error with the proposed procedure. Links to other methods based on the application of canonical analysis are discussed. A simulation experiment was done in order to study the performance of the procedure. It is seen that the graphical representation of t?_N, as a function of N, can be very useful in identifying the dynamic structure of ARMA models. Furthermore, for the model considered, the proposed identification procedure performs very well for series of 100 observations or more and reasonably well with short series of 50 observations.     

Value-at-risk forecasting based on Gaussian mixture ARMA–GARCH model

In this paper, we develop a new forecasting algorithm for value-at-risk (VaR) based on ARMA–GARCH (autoregressive moving average–generalized autoregressive conditional heteroskedastic) models whose innovations follow a Gaussian mixture distribution. For the parameter estimation, we employ the conditional least squares and quasi-maximum-likelihood estimator (QMLE) for ARMA and GARCH parameters, respectively. In particular, Gaussian mixture parameters are estimated based on the residuals obtained from the QMLE of GARCH parameters. Our algorithm provides a handy methodology, spending much less time in calculation than the existing resampling and bias-correction method developed in Hartz et al. [Accurate value-at-risk forecasting based on the normal-GARCH model, Comput. Stat. Data Anal. 50 (2006), pp. 3032–3052]. Through a simulation study and a real-data analysis, it is shown that our method provides an accurate VaR prediction.     

Euler(p,q) Processes and Their Application to Non Stationary Time Series with Time Varying Frequencies

We introduce Euler(p, q) processes as an extension of the Euler(p) processes for purposes of obtaining more parsimonious models for non stationary processes whose periodic behavior changes approximately linearly in time. The discrete Euler(p, q) models are a class of multiplicative stationary (M-stationary) processes and basic properties are derived. The relationship between continuous and discrete mixed Euler processes is shown. Fundamental to the theory and application of Euler(p, q) processes is a dual relationship between discrete Euler(p, q) processes and ARMA processes, which is established. The usefulness of Euler(p, q) processes is examined by comparing spectral estimation with that obtained by existing methods using both simulated and real data.     

Forecasting Vector ARMA Processes With Systematically Missing Observations

The following two predictors are compared for time series with systematically missing observations: (a) A time series model is fitted to the full series X_t , and forecasts are based on this model, (b) A time series model is fitted to the series with systematically missing observations Y _τ, and forecasts are based on the resulting model. If the data generation processes are known vector autoregressive moving average (ARMA) processes, the first predictor is at least as efficient as the second one in a mean squared error sense. Conditions are given for the two predictors to be identical. If only the ARMA orders of the generation processes are known and the coefficients are estimated, or if the process orders and coefficients are estimated, the first predictor is again, in general, superior. There are, however, exceptions in which the second predictor, using seemingly less information, may be better. These results are discussed, using both asymptotic theory and small sample simulations. Some economic time series are used as illustrative examples.     

On Fan's adaptive Neyman tests for comparing two spectral densities

In this paper, we consider tests for assessing whether two stationary and independent time series have the same spectral densities (or same autocovariance functions). Both frequency domain and time domain test statistics for this purpose are reviewed. The adaptive Neyman tests are then introduced and their performances are investigated. Our tests are adaptive, that is, they are constructed completely by the data and do not involve any unknown smoothing parameters. Simulation studies show that our proposed tests are at least comparable to the current tests in most cases. Furthermore, our tests are much more powerful in some cases, such as against the long orders of autoregressive moving average (ARMA) models such as seasonal ARMA series.     

Empirical study of robust estimation methods for PAR models with application to the air quality area

Abstract

This paper compares three estimators for periodic autoregressive (PAR) models. The first is the classical periodic Yule-Walker estimator (YWE). The second is a robust version of YWE (RYWE) which uses the robust autocovariance function in the periodic Yule-Walker equations, and the third is the robust least squares estimator (RLSE) based on iterative least squares with robust versions of the original time series. The daily mean particulate matter concentration (PM10) data is used to illustrate the methodologies in a real application, that is, in the Air Quality area.     

Estimation in a bivariate integer-valued autoregressive process

ABSTRACT

A bivariate integer-valued autoregressive time series model is presented. The model structure is based on binomial thinning. The unconditional and conditional first and second moments are considered. Correlation structure of marginal processes is shown to be analogous to the ARMA(2, 1) model. Some estimation methods such as the Yule–Walker and conditional least squares are considered and the asymptotic distributions of the obtained estimators are derived. Comparison between bivariate model with binomial thinning and bivariate model with negative binomial thinning is given.     

Bayesian estimation for time-series regressions improved with exact likelihoods

We propose an estimation procedure for time-series regression models under the Bayesian inference framework. With the exact method of Wise [Wise, J. (1955). The autocorrelation function and spectral density function. Biometrika, 42, 151–159], an exact likelihood function can be obtained instead of the likelihood conditional on initial observations. The constraints on the parameter space arising from the stationarity conditions are handled by a reparametrization, which was not taken into consideration by Chib [Chib, S. (1993). Bayes regression with autoregressive errors: A Gibbs sampling approach. J. Econometrics, 58, 275–294] or Chib and Greenberg [Chib, S. and Greenberg, E. (1994). Bayes inference in regression model with ARMA(p, q) errors. J. Econometrics, 64, 183–206]. Simulation studies show that our method leads to better inferential results than their results.