Least quantile regression via modern optimization

From MaRDI portal
Publication:482902

DOI10.1214/14-AOS1223zbMATH Open1302.62154arXiv1310.8625MaRDI QIDQ482902FDOQ482902

Rahul Mazumder, Dimitris Bertsimas

Publication date: 6 January 2015

Published in: The Annals of Statistics (Search for Journal in Brave)

Abstract: We address the Least Quantile of Squares (LQS) (and in particular the Least Median of Squares) regression problem using modern optimization methods. We propose a Mixed Integer Optimization (MIO) formulation of the LQS problem which allows us to find a provably global optimal solution for the LQS problem. Our MIO framework has the appealing characteristic that if we terminate the algorithm early, we obtain a solution with a guarantee on its sub-optimality. We also propose continuous optimization methods based on first-order subdifferential methods, sequential linear optimization and hybrid combinations of them to obtain near optimal solutions to the LQS problem. The MIO algorithm is found to benefit significantly from high quality solutions delivered by our continuous optimization based methods. We further show that the MIO approach leads to (a) an optimal solution for any dataset, where the data-points (yi,mathbfxi)'s are not necessarily in general position, (b) a simple proof of the breakdown point of the LQS objective value that holds for any dataset and (c) an extension to situations where there are polyhedral constraints on the regression coefficient vector. We report computational results with both synthetic and real-world datasets showing that the MIO algorithm with warm starts from the continuous optimization methods solve small (n=100) and medium (n=500) size problems to provable optimality in under two hours, and outperform all publicly available methods for large-scale (n=10,000) LQS problems.


Full work available at URL: https://arxiv.org/abs/1310.8625




Recommendations




Cites Work


Cited In (21)

Uses Software





This page was built for publication: Least quantile regression via modern optimization

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q482902)