Estimation of the parameters of symmetric stable ARMA and ARMA–GARCH models

In this article, we first propose the modified Hannan–Rissanen Method for estimating the parameters of autoregressive moving average (ARMA) process with symmetric stable noise and symmetric stable generalized autoregressive conditional heteroskedastic (GARCH) noise. Next, we propose the modified empirical characteristic function method for the estimation of GARCH parameters with symmetric stable noise. Further, we show the efficiency, accuracy and simplicity of our methods with Monte-Carlo simulation. Finally, we apply our proposed methods to model the financial data.     

Analysis of autoregressive models with symmetric stable innovations

This paper develops algorithms for fitting autoregressive models with symmetric stable innovations using auto-covariation function. A recursive algorithm is proposed for generalized Yule-Walker estimation of autoregressive coefficients and partial auto-covariation function. It also introduces a new information criterion, useful for consistent order selection. Applications of the proposed methods are illustrated using observations simulated from autoregressive models with symmetric stable innovations as well as by analysing a set of real data.     

The specification of vector autoregressive moving average models

In this paper we propose a new identification method based on the residual white noise autoregressive criterion (Pukkila et al., 1990) to select the order of VARMA structures. Results from extensive simulation experiments based on different model structures with varying number of observations and number of component series are used to demonstrate the performance of this new procedure. We also use economic and business data to compare the model structures selected by this order selection method with those identified in other published studies.     

Testing for residual correlation of any order in the autoregressive process

We are interested in the implications of a linearly autocorrelated driven noise on the asymptotic behavior of the usual least-squares estimator in a stable autoregressive process. We show that the least-squares estimator is not consistent and we suggest a sharp analysis of its almost sure limiting value as well as its asymptotic normality. We also establish the almost sure convergence and the asymptotic normality of the estimated serial correlation parameter of the driven noise. Then, we derive a statistical procedure enabling to test for correlation of any order in the residuals of an autoregressive modelling, giving clearly better results than the commonly used portmanteau tests of Ljung–Box and Box–Pierce, and appearing to outperform the Breusch–Godfrey procedure on small-sized samples.     

Normal tests for unit roots based on instrumental variable estimators

For estimating unit roots of autoregressive processes, we introduce a new instrumental variable (IV) method which discounts large values of regressors corresponding to the unit roots. Based on the IV estimator, we propose new unit root tests whose limiting null distributions are standard normal. Observation at time t is adjusted for mean recursively by the sample mean of observations up to the time t. The powers of the proposed tests are better than those of the Dickey–Fuller tests and are comparable to those of the tests based on the weighted symmetric estimator, which are known to have the best power against stationary alternatives.     

Indirect estimation of randomized generalized autoregressive conditional heteroskedastic models

The class of generalized autoregressive conditional heteroskedastic (GARCH) models can be used to describe the volatility with less parameters than autoregressive conditional heteroskedastic (ARCH)-type models, their distributions are heavy-tailed, with time-dependent conditional variance, and are able to model clustering of volatility. Despite all these facts, the way that GARCH models are built imposes limits on the heaviness of the tails of their unconditional distribution. The class of randomized generalized autoregressive conditional heteroskedastic (R-GARCH) models includes the ARCH and GARCH models allowing the use of stable innovations. Estimation methods and empirical analysis of R-GARCH models are the focus of this work. We present the indirect inference method to estimate the R-GARCH models, some simulations and an empirical application.     

Denoising ozone concentration measurements with BAMS filtering

We propose a method for filtering self-similar geophysical signals infected by an autoregressive noise using a combination of non-decimated wavelet transform and a Bayesian model. In the application part, we consider separating the instrumentation noise from high frequency ozone concentration measurements sampled in the atmospheric boundary layer. The elicitation of priors needed to specify the statistical model in this application is guided by the well-known Kolmogorov K41-theory, which describes the statistical structure of turbulent high frequency scalar concentration fluctuations.     

Nonlinear semiparametric autoregressive model with finite mixtures of scale mixtures of skew normal innovations

We propose data generating structures which can be represented as the nonlinear autoregressive models with single and finite mixtures of scale mixtures of skew normal innovations. This class of models covers symmetric/asymmetric and light/heavy-tailed distributions, so provide a useful generalization of the symmetrical nonlinear autoregressive models. As semiparametric and nonparametric curve estimation are the approaches for exploring the structure of a nonlinear time series data set, in this article the semiparametric estimator for estimating the nonlinear function of the model is investigated based on the conditional least square method and nonparametric kernel approach. Also, an Expectation–Maximization-type algorithm to perform the maximum likelihood (ML) inference of unknown parameters of the model is proposed. Furthermore, some strong and weak consistency of the semiparametric estimator in this class of models are presented. Finally, to illustrate the usefulness of the proposed model, some simulation studies and an application to real data set are considered.     

Inference for linear and nonlinear stable error processes via estimating functions

This paper describes an estimating function approach for parameter estimation in linear and nonlinear times series models with infinite variance stable errors. Joint estimates of location and scale parameters are derived for classes of autoregressive (AR) models and random coefficient autoregressive (RCA) models with stable errors, as well as for AR models with stable autoregressive conditionally heteroscedastic (ARCH) errors. Fast, on-line, recursive parametric estimation for the location parameter based on estimating functions is discussed using simulation studies. A real financial time series is also discussed in some detail.     

LS estimation of periodic autoregressive models with non-Gaussian errors: a simulation study

In this article, the least squares (LS) estimates of the parameters of periodic autoregressive (PAR) models are investigated for various distributions of error terms via Monte-Carlo simulation. Beside the Gaussian distribution, this study covers the exponential, gamma, student-t, and Cauchy distributions. The estimates are compared for various distributions via bias and MSE criterion. The effect of other factors are also examined as the non-constancy of model orders, the non-constancy of the variances of seasonal white noise, the period length, and the length of the time series. The simulation results indicate that this method is in general robust for the estimation of AR parameters with respect to the distribution of error terms and other factors. However, the estimates of those parameters were, in some cases, noticeably poor for Cauchy distribution. It is also noticed that the variances of estimates of white noise variances are highly affected by the degree of skewness of the distribution of error terms.     

A new estimator for the unit root

We compare the ordinary least squares, weighted symmetric, modified weighted symmetric (MWS), maximum likelihood, and our new modification for least squares (MLS) estimator for first-order autoregressive in the case of unit root using Monte Carlo method. The Monte Carlo study sheds some light on how well the estimators and the predictors perform on different samples sizes. We found that MLS estimator is less biased and has less mean squared error (MSE) than any other estimators, and MWS predictor error performs well, in the sense of MSE, than any other predictors’ methods. The sample percentiles for the distribution of the τ statistic for the first, second, and third periods in the future, for alternative estimators, are reported to know if it agrees with those of normal distribution or not.     

L p –Estimators as Estimates of a Parameter of Location for a Sharp–pointed Symmetric Density

We study the asymptotics of L _p estimators, p > 0, over a sample having a symmetric density with a sharp–point at the centre of symmetry of the distribution. The rates of convergence of the L _p estimators in this situation depend on p and on the shape of the density. To obtain some of the limit distributions, we present new results in the asymptotics of M–estimators. We extend the delta method to the case when the Euclidean norm of the conveniently normalized M–estimators converge to a power of the Euclidean norm of a (possibly Gaussian) stable distribution.     

Some results on unilateral ARMA lattice processes

This paper presents several linked results on unilateral autoregressive moving average processes on a rectangular lattice. It is shown that axially symmetric two-dimensional quadrant processes must be separable. Exact forms for the inverse variance matrix are obtained in some cases, which allow exact Gaussian maximum likelihood estimation and simulation. It is shown that generating functions can be used for extrapolation. The herringbone simulation method is discussed.     

On the estimation problem of periodic autoregressive time series: symmetric and asymmetric innovations

Periodic autoregressive (PAR) models with symmetric innovations are widely used on time series analysis, whereas its asymmetric counterpart inference remains a challenge, because of a number of problems related to the existing computational methods. In this paper, we use an interesting relationship between periodic autoregressive and vector autoregressive (VAR) models to study maximum likelihood and Bayesian approaches to the inference of a PAR model with normal and skew-normal innovations, where different kinds of estimation methods for the unknown parameters are examined. Several technical difficulties which are usually complicated to handle are reported. Results are compared with the existing classical solutions and the practical implementations of the proposed algorithms are illustrated via comprehensive simulation studies. The methods developed in the study are applied and illustrate a real-time series. The Bayes factor is also used to compare the multivariate normal model versus the multivariate skew-normal model.     

A Comparison of Unit-Root Test Criteria

During the past 15 years, the ordinary least squares estimator and the corresponding pivotal statistic have been widely used for testing the unit-root hypothesis in autoregressive processes. Recently, several new criteria, based on maximum likelihood estimators and weighted symmetric estimators, have been proposed. In this article, we describe several different test criteria. Results from a Monte Carlo study that compares the power of the different criteria indicate that the new tests are more powerful against the stationary alternative. Of the procedures studied, the weighted symmetric estimator and the unconditional maximum likelihood estimator provide the most powerful tests against the stationary alternative. As an illustration, the weekly series of one-month treasury-bill rates is analyzed.     

Robust image segmentation via Bayesian type criterion

Image segmentation plays an important role in image processing before image recognition or compression. Many segmentation solutions follow the information theoretic criteria and often have excellent results; however, they are not robust to reduce the noise effect in contaminated image data. To guarantee the optimal segmentation with possible noise, a robust Bayesian information criterion is proposed to segment a grayscale image and it is less sensitive to noise. The asymptotic properties are also studied. Monte Carlo numerical experiments along with a brain magnetic resonance image are conducted to evaluate the performance of the new method.     

A new lifetime distribution for a series-parallel system: properties,applications and estimations under progressive type-II censoring

A compound class of zero truncated Poisson and lifetime distributions is introduced. A specialization is paved to a new three-parameter distribution, called doubly Poisson-exponential distribution, which may represent the lifetime of units connected in a series-parallel system. The new distribution can be obtained by compounding two zero truncated Poisson distributions with an exponential distribution. Among its motivations is that its hazard rate function can take different shapes such as decreasing, increasing and upside-down bathtub depending on the values of its parameters. Several properties of the new distribution are discussed. Based on progressive type-II censoring, six estimation methods [maximum likelihood, moments, least squares, weighted least squares and Bayes (under linear-exponential and general entropy loss functions) estimations] are used to estimate the involved parameters. The performance of these methods is investigated through a simulation study. The Bayes estimates are obtained using Markov chain Monte Carlo algorithm. In addition, confidence intervals, symmetric credible intervals and highest posterior density credible intervals of the parameters are obtained. Finally, an application to a real data set is used to compare the new distribution with other five distributions.     

Maximum likelihood estimation for directional conditionally autoregressive models

A spatial process observed over a lattice or a set of irregular regions is usually modeled using a conditionally autoregressive (CAR) model. The neighborhoods within a CAR model are generally formed using only the inter-distances or boundaries between the regions. To accommodate directional spatial variation, a new class of spatial models is proposed using different weights given to neighbors in different directions. The proposed model generalizes the usual CAR model by accounting for spatial anisotropy. Maximum likelihood estimators are derived and shown to be consistent under some regularity conditions. Simulation studies are presented to evaluate the finite sample performance of the new model as compared to the CAR model. Finally, the method is illustrated using a data set on the crime rates of Columbus, OH and on the elevated blood lead levels of children under the age of 72 months observed in Virginia in the year of 2000.     

THE AUTOREGRESSIVE MOVING AVERAGE MODEL FOR SPATIAL ANALYSIS

A two dimensional autoregressive moving average spatial model is used to analyse spatial interaction. Maximum likelihood estimates of the unknown parameters are derived as the solution of a system of nonlinear equations, and are shown to be best asymptotic normal. One important computational procedure is discussed. The argument is extended to the general regression model with autoregressive moving average residuals. Explicit computational formulae are given.     

Maximum likelihood estimation for generalized conditionally autoregressive models of spatial data

Conditionally autoregressive (CAR) models are often used to analyze a spatial process observed over a lattice or a set of irregular regions. The neighborhoods within a CAR model are generally formed deterministically using the inter-distances or boundaries between the regions. To accommodate directional and inherent anisotropy variation, a new class of spatial models is proposed that adaptively determines neighbors based on a bivariate kernel using the distances and angles between the centroid of the regions. The newly proposed model generalizes the usual CAR model in a sense of accounting for adaptively determined weights. Maximum likelihood estimators are derived and simulation studies are presented for the sampling properties of the estimates on the new model, which is compared to the CAR model. Finally the method is illustrated using a data set on the elevated blood lead levels of children under the age of 72 months observed in Virginia in the year of 2000.