We consider a discrete optimization formulation for learning sparse classifiers, where the outcome depends upon a linear combination of a small subset of features. Recent work has shown that mixed integer programming (MIP) can be used to solve (to optimality) $\ell_0$-regularized regression problems at scales much larger than what was conventionally considered possible. Despite their usefulness, MIP-based global optimization approaches are significantly slower compared to the relatively mature algorithms for $\ell_1$-regularization and heuristics for nonconvex regularized problems. We aim to bridge this gap in computation times by developing new MIP-based algorithms for $\ell_0$-regularized classification. We propose two classes of scalable algorithms: an exact algorithm that can handle $p\approx 50,000$ features in a few minutes, and approximate algorithms that can address instances with $p\approx 10^6$ in times comparable to the fast $\ell_1$-based algorithms. Our exact algorithm is based on the novel idea of \textsl{integrality generation}, which solves the original problem (with $p$ binary variables) via a sequence of mixed integer programs that involve a small number of binary variables. Our approximate algorithms are based on coordinate descent and local combinatorial search. In addition, we present new estimation error bounds for a class of $\ell_0$-regularized estimators. Experiments on real and synthetic data demonstrate that our approach leads to models with considerably improved statistical performance (especially, variable selection) when compared to competing methods.
翻译:我们考虑为学习分散的分类器制定离散的优化配方,其结果取决于一小组特征的线性组合。最近的工作表明,混合整数编程(MIP)可以用来解决(最优化)$\ell_0美元(美元)常规回归问题,其规模远大于常规认为可行的范围。尽管其有用性,但基于MIP的全球优化方法比相对成熟的计算法要慢得多,因为美元1美元(美元1美元)的正规化和超常性(美元1美元)的常规化问题。我们的目标是通过为美元=ell_0美元常规分类开发基于MIP的新算法来弥合计算时间上的这一差距。我们提出了两类可缩略微缩缩缩缩缩缩缩缩算法:精确算法可以在几分钟内处理$p\approx 50,000美元特征,以及近似于美元=1美元(美元)的快速算法。我们精确的算法基于新改进的计算方法,即基于基于基于正缩缩缩缩缩缩缩缩缩缩算法的计算方法,我们用原始变数进行原始变校正数据。