SENSITIVITY ANALYSIS FOR MIRROR-STRATIFIABLE CONVEX FUNCTIONS

被引:14
作者
Fadili, Jalal [1 ]
Malick, Jerome [2 ,3 ]
Peyre, Gabriel [4 ,5 ]
机构
[1] Normandie Univ, GREYC, CNRS, ENSICAEN, Caen, France
[2] CNRS, St Martin Dheres, France
[3] LJK, St Martin Dheres, France
[4] ENS Paris, CNRS, Paris, France
[5] ENS Paris, DMA, Paris, France
基金
欧洲研究理事会;
关键词
convex analysis; inverse problems; sensitivity; active sets; first-order splitting algorithms; applications in imaging and machine learning; LOCAL LINEAR CONVERGENCE; MODEL SELECTION; PARTIAL SMOOTHNESS; CONSISTENCY; REGRESSION; SHRINKAGE; ALGORITHM;
D O I
10.1137/17M113825X
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
This paper provides a set of sensitivity analysis and activity identification results for a class of convex functions with a strong geometric structure that we have called "mirror-stratifiable." These functions are such that there is a bijection between a primal and a dual stratification of the space into partitioning sets, called strata. This pairing is crucial for tracking the strata that are identifiable by solutions of parametrized optimization problems or by iterates of optimization algorithms. This class of functions encompasses all regularizers routinely used in signal and image processing, machine learning, and statistics. We show that this "mirror-stratifiable" structure enjoys a nice sensitivity theory, allowing us to study stability of solutions of optimization problems with respect to small perturbations, as well as activity identification of first-order proximal splitting-type algorithms. Existing results in the literature typically assume that, under a nondegeneracy condition, the active set associated with a minimizer is stable with respect to small perturbations and is identified in finite time by optimization schemes. In contrast, our results do not require a nondegeneracy assumption: in consequence, the optimal active set is not necessarily stable anymore, but we are able to track precisely the set of identifiable strata. We show that these results have crucial implications when solving challenging ill-posed inverse problems via regularization, a typical scenario in which the nondegeneracy condition is not fulfilled. Our theoretical results, illustrated by numerical simulations, allow us to characterize the instability behaviour of the regularized solutions by locating the set of all low-dimensional strata that can be potentially identified by these solutions.
引用
收藏
页码:2975 / 3000
页数:26
相关论文
共 60 条
[1]  
[Anonymous], SIAM P APPL MATH
[2]  
[Anonymous], 1998, GRUNDLEHREN MATH WIS
[3]  
[Anonymous], 2002, THESIS STANFORD U
[4]  
[Anonymous], 2009, WAVELET TOUR SIGNAL
[5]  
[Anonymous], 2014, CVX MATLAB SOFTWARE
[6]  
[Anonymous], 2006, Journal of the Royal Statistical Society, Series B
[7]  
[Anonymous], FDN TRENDS OPTIM, DOI DOI 10.1561/2400000003
[8]  
[Anonymous], PRINCET MATH SER
[9]  
[Anonymous], 2016, GREYC CNRS UMR 6072
[10]  
Bach FR, 2008, J MACH LEARN RES, V9, P1019