Bounded Influence Propagation {\tau}-Estimation: A New Robust Method for ARMA Model Estimation

TL;DR

This paper introduces a robust and efficient estimator for ARMA models, called the BIP-τ-estimator, which effectively handles outliers and impulsive noise, with proven statistical properties and practical algorithms.

Contribution

The paper proposes the BIP-τ-estimator for ARMA models, incorporating an auxiliary model to prevent outlier propagation and providing algorithms for robust parameter estimation.

Findings

The BIP-τ-estimator is strongly consistent and asymptotically normal.

Explicit influence function derived for AR(1) with additive outliers.

Numerical experiments show superior performance over existing methods.

Abstract

A new robust and statistically efficient estimator for ARMA models called the bounded influence propagation (BIP) {\tau}-estimator is proposed. The estimator incorporates an auxiliary model, which prevents the propagation of outliers. Strong consistency and asymptotic normality of the estimator for ARMA models that are driven by independently and identically distributed (iid) innovations with symmetric distributions are established. To analyze the infinitesimal effect of outliers on the estimator, the influence function is derived and computed explicitly for an AR(1) model with additive outliers. To obtain estimates for the AR(p) model, a robust Durbin-Levinson type and a forward-backward algorithm are proposed. An iterative algorithm to robustly obtain ARMA(p,q) parameter estimates is also presented. The problem of finding a robust initialization is addressed, which for orders p+q>2 is a non-trivial matter. Numerical experiments are conducted to compare the finite sample performance of the proposed estimator to existing robust methodologies for different types of outliers both in terms of average and of worst-case performance, as measured by the maximum bias curve. To illustrate the practical applicability of the proposed estimator, a real-data example of outlier cleaning for R-R interval plots derived from electrocardiographic (ECG) data is considered. The proposed estimator is not limited to biomedical applications, but is also useful in any real-world problem whose observations can be modeled as an ARMA process disturbed by outliers or impulsive noise.