Implicit differentiation of Lasso-type models for hyperparameter optimization
This addresses the difficulty of setting regularization parameters for Lasso-type estimators, which is crucial in practice but often incremental as it builds on existing bi-level optimization approaches.
The paper tackles the challenge of hyperparameter optimization for Lasso-type models by introducing an efficient implicit differentiation algorithm that avoids matrix inversion and leverages solution sparsity, demonstrating superior performance in optimizing held-out data error or SURE compared to standard methods.
Setting regularization parameters for Lasso-type estimators is notoriously difficult, though crucial in practice. The most popular hyperparameter optimization approach is grid-search using held-out validation data. Grid-search however requires to choose a predefined grid for each parameter, which scales exponentially in the number of parameters. Another approach is to cast hyperparameter optimization as a bi-level optimization problem, one can solve by gradient descent. The key challenge for these methods is the estimation of the gradient with respect to the hyperparameters. Computing this gradient via forward or backward automatic differentiation is possible yet usually suffers from high memory consumption. Alternatively implicit differentiation typically involves solving a linear system which can be prohibitive and numerically unstable in high dimension. In addition, implicit differentiation usually assumes smooth loss functions, which is not the case for Lasso-type problems. This work introduces an efficient implicit differentiation algorithm, without matrix inversion, tailored for Lasso-type problems. Our approach scales to high-dimensional data by leveraging the sparsity of the solutions. Experiments demonstrate that the proposed method outperforms a large number of standard methods to optimize the error on held-out data, or the Stein Unbiased Risk Estimator (SURE).