Davide La Vecchia

Infinitesimal Robustness for Diffusions

We develop infinitesimally robust statistical procedures for the general diffusion processes. We first prove the existence and uniqueness of the times-series influence function of conditionally unbiased M-estimators for ergodic and stationary diffusions, under weak conditions on the (martingale) estimating function used. We then characterize the robustness of M-estimators for diffusions and derive a class of conditionally unbiased optimal robust estimators. To compute these estimators, we propose a general algorithm, which exploits approximation methods for diffusions in the computation of the robust estimating function. Monte Carlo simulation shows a good performance of our robust estimators and an application to the robust estimation of the exchange rate dynamics within a target zone illustrates the methodology in a real-data application.

Higher-Order Infinitesimal Robustness

Using the von Mises expansion, we study the higher-order infinitesimal robustness of a general M-functional and characterize its second-order properties. We show that second-order robustness is equivalent to the boundedness of both the estimator’s estimating function and its derivative with respect to the parameter. It implies, at the same time, (i) variance robustness and (ii) robustness of higher-order saddlepoint approximations to the estimator’s finite sample density. The proposed construction of second-order robust M-estimators is fairly general and potentially useful in a variety of relevant settings. Besides the theoretical contributions, we discuss the main computational issues and provide an algorithm for the implementation of second-order robust M-estimators. Finally, we illustrate our theory by Monte Carlo simulation and in a real-data estimation of the maximal losses of Nikkei 225 index returns. Our findings indicate that second-order robust estimators can improve on other widely applied robust estimators, in terms of efficiency and robustness, for moderate to small sample sizes and in the presence of deviations from ideal parametric models. Supplementary materials for this article are available online.

Robust heart rate variability analysis by generalized entropy minimization

Typical heart rate variability (HRV) times series are cluttered with outliers generated by measurement errors, artifacts and ectopic beats. Robust estimation is an important tool in HRV analysis, since it allows clinicians to detect arrhythmia and other anomalous patterns by reducing the impact of outliers. A robust estimator for a flexible class of time series models is proposed and its empirical performance in the context of HRV data analysis is studied. The methodology entails the minimization of a pseudo-likelihood criterion function based on a generalized measure of information. The resulting estimating functions are typically re-descending, which enable reliable detection of anomalous HRV patterns and stable estimates in the presence of outliers. The infinitesimal robustness and the stability properties of the new method are illustrated through numerical simulations and two case studies from the Massachusetts Institute of Technology and Bostons Beth Israel Hospital data, an important benchmark data set in HRV analysis.