Robust thin-plate splines for multivariate spatial smoothing

Abstract

A novel family of multivariate robust smoothers based on the thin-plate (Sobolev) penalty that is particularly suitable for the analysis of spatial data is proposed. The proposed family of estimators can be expediently computed even in high dimensions, is invariant with respect to rigid transformations of the coordinate axes and can be shown to possess optimal theoretical properties under mild assumptions. The competitive performance of the proposed thin-plate spline estimators relative to their non-robust counterpart is illustrated in a simulation study and a real data example involving two-dimensional geographical data on ozone concentration.