We present a new algorithmic framework for grouped variable selection that is based on discrete mathematical optimization. While there exist several appealing approaches based on convex relaxations and nonconvex heuristics, we focus on optimal solutions for the $\ell_0$-regularized formulation, a problem that is relatively unexplored due to computational challenges. Our methodology covers both high-dimensional linear regression and nonparametric sparse additive modeling with smooth components. Our algorithmic framework consists of approximate and exact algorithms. The approximate algorithms are based on coordinate descent and local search, with runtimes comparable to popular sparse learning algorithms. Our exact algorithm is based on a standalone branch-and-bound (BnB) framework, which can solve the associated mixed integer programming (MIP) problem to certified optimality. By exploiting the problem structure, our custom BnB algorithm can solve to optimality problem instances with $5 \times 10^6$ features and $10^3$ observations in minutes to hours -- over $1000$ times larger than what is currently possible using state-of-the-art commercial MIP solvers. We also explore statistical properties of the $\ell_0$-based estimators. We demonstrate, theoretically and empirically, that our proposed estimators have an edge over popular group-sparse estimators in terms of statistical performance in various regimes. We provide an open-source implementation of our proposed framework.
翻译:我们为基于离散数学优化的分组变量选择提出了一个新的算法框架。 虽然存在一些基于康韦克斯放松和非康维克斯超光度的吸引性方法,但我们侧重于以美元=0.0美元正规化配方的最佳解决方案,由于计算方面的挑战,这一问题相对没有探讨。我们的方法涵盖高维线回归和非参数性稀释添加模型,并包含平滑的成分。我们的算法框架包括近度和精确算法。近似算法基于协调血统和本地搜索,运行时间与流行的稀释算法相仿。我们精确的算法基于独立分支和约束(BnB)框架,它能够解决相关的混合整流方案(MIP)优化问题。通过利用问题结构,我们的定制的BnB算法可以解决最佳性问题,5美元=6美元特征和每小时10美元的拟议观察。我们提出的开放性算法系统比目前可能使用的最新商业MIP解算法的运行时间要大1 000多倍。我们还探索了我们拟议中以美元为核心的统计学底框架的统计属性。