Robust M-Estimators of Multivariate Location and Scatter

Let x₁,⋯, xₙ be a sample from an m-variate distribution which is spherically symmetric up to an affine transformation. This paper deals with the robust estimation of the location vector t and scatter matrix V by means of "M-estimators," defined as solutions of the system: ∑ᵢ u₁(dᵢ)(xᵢ - t) = 0 and n⁻¹∑ᵢ u₂(dᵢ²)(xᵢ - t)(xᵢ - t)' = V, where dᵢ² = (xᵢ - t)'V⁻¹(xᵢ - t). Existence and uniqueness of solutions of this system are proved under general assumptions about the functions u₁ and u₂. Then the estimators are shown to be consistent and asymptotically normal. The breakdown bound and the influence function are calculated, showing some weaknesses of the estimates for high dimensionality. An algorithm for the numerical calculation of the estimators is described. Finally, numerical values of asymptotic variances, and Monte Carlo small-sample results are exhibited.