CMStatistics 2016: Start Registration
View Submission - CMStatistics
B0589
Title: Least trimmed squares for functional principal component analysis Authors:  Holger Cevallos Valdiviezo - Ghent University (Belgium)
Matias Salibian-Barrera - The University of British Columbia (Canada)
Stefan Van Aelst - University of Leuven (Belgium) [presenting]
Abstract: Classical functional principal component analysis can yield erroneous approximations in presence of outliers. To reduce the influence of atypical data we propose two methods based on trimming: a multivariate least trimmed squares (LTS) estimator and a componentwise variant. The multivariate LTS minimizes the least squares criterion over subsets of curves. The componentwise version minimizes the sum of univariate LTS scale estimators in each of the components. In general the curves can be considered as realizations of a random element on a separable Hilbert space. For a fixed dimension $q$, we then aim to robustly estimate the $q$-dimensional linear subspace that gives the best approximation to the functional data. Our estimators use smoothing to first represent irregularly spaced curves in a high-dimensional space and then calculates the LTS solution on these multivariate data. The solution of the multivariate data is subsequently mapped back onto the Hilbert space. Poorly fitted observations can then be flagged as outliers. A simulation study and real data applications show that our estimators yield competitive results, both in identifying outliers and approximating regular data when compared to other existing methods.