Matrix strategies for computing the least trimmed squares estimation of the general linear and SUR models

  • Authors:
  • Marc Hofmann;Erricos John Kontoghiorghes

  • Affiliations:
  • Wirtschaftswissenschaftliche Fakultät, Universität Basel, Switzerland;Department of Public and Business Administration, University of Cyprus, Cyprus and Department of Economics, Queen Mary, University of London, UK

  • Venue:
  • Computational Statistics & Data Analysis
  • Year:
  • 2010

Quantified Score

Hi-index 0.03

Visualization

Abstract

An algorithm for computing the exact least trimmed squares (LTS) estimator of the standard regression model has recently been proposed. The LTS algorithm is adapted to the general linear and seemingly unrelated regressions models with possible singular dispersion matrices. It searches through a regression tree to find the optimal estimates and has combinatorial complexity. The model is formulated as a generalized linear least squares problem. Efficient matrix techniques are employed to update the generalized residual sum of squares of a subset model. Specifically, the new algorithm utilizes previous computations to update a generalized QR decomposition by a single row. The sparse structure of the model is exploited. Theoretical measures of computational complexity are provided. Experimental results confirm the ability of the new algorithms to identify outlying observations.