Abstract:High-dimensional linear regression under heavy-tailed noise or outlier corruption is challenging, both computationally and statistically. Convex approaches have been proven statistically optimal but suffer from high computational costs, especially since the robust loss functions are usually non-smooth. More recently, computationally fast non-convex approaches via sub-gradient descent are proposed, which, unfortunately, fail to deliver a statistically consistent estimator even under sub-Gaussian noise. In this paper, we introduce a projected sub-gradient descent algorithm for both the sparse linear regression and low-rank linear regression problems. The algorithm is not only computationally efficient with linear convergence but also statistically optimal, be the noise Gaussian or heavy-tailed with a finite 1 + epsilon moment. The convergence theory is established for a general framework and its specific applications to absolute loss, Huber loss and quantile loss are investigated. Compared with existing non-convex methods, ours reveals a surprising phenomenon of two-phase convergence. In phase one, the algorithm behaves as in typical non-smooth optimization that requires gradually decaying stepsizes. However, phase one only delivers a statistically sub-optimal estimator, which is already observed in the existing literature. Interestingly, during phase two, the algorithm converges linearly as if minimizing a smooth and strongly convex objective function, and thus a constant stepsize suffices. Underlying the phase-two convergence is the smoothing effect of random noise to the non-smooth robust losses in an area close but not too close to the truth. Numerical simulations confirm our theoretical discovery and showcase the superiority of our algorithm over prior methods.

Outlier-resilient model fitting via percentile losses: Methods for general and convex residuals

A highly robust estimator through partially likelihood function modeling and its application in computer vision

Percentile-Based Residuals for Model Assessment

Outlier Detection by Energy Minimization in Quantized Residual Preference Space for Geometric Model Fitting

GORFLM: Globally Optimal Robust Fitting for Linear Model

A Robust Learning Approach for Regression Models Based on Distributionally Robust Optimization.

A Distributionally Robust Optimization Approach for Outlier Detection *

Robust Optimization Using a Methodology Based on Cross Entropy Methods

Sparse and Robust Linear Regression: An Optimization Algorithm and Its Statistical Properties

A robust mixed‐effects parametric quantile regression model for continuous proportions: quantifying the constraints to vitality in cushion plants

A redescending M-estimator approach for outlier-resilient modeling

Double Probability Integral Transform Residuals for Regression Models with Discrete Outcomes

Computationally Efficient and Statistically Optimal Robust High-Dimensional Linear Regression

Globally Optimal Linear Model Fitting with Unit-Norm Constraint

A Robust Learning Algorithm for Regression Models Using Distributionally Robust Optimization under the Wasserstein Metric

Model Misfit Minimization

Outlier Detection Using Robust Optimization with Uncertainty Sets Constructed from Risk Measures.

Risk-Based Robust Statistical Learning by Stochastic Difference-of-Convex Value-Function Optimization

Revisiting Robust Model Fitting Using Truncated Loss

Bootstrap estimation of the proportion of outliers in robust regression

Model-Robust Designs for Quantile Regression