Sparsity and Some Basics of L1 Regularization

Sparsity 是当今机器学习领域中的一个重要话题。John Lafferty 和 Larry Wasserman 在 2006 年的一篇评论中提到：

Some current challenges … are high dimensional data, sparsity, semi-supervised learning, the relation between computation and risk, and structured prediction.John Lafferty and Larry Wasserman. Challenges in statistical machine learning. Statistica Sinica. Volume 16, Number 2, pp. 307-323, 2006.

Sparsity 的最重要的“客户”大概要属 high dimensional data 了吧。现在的机器学习问题中，具有非常高维度的数据随处可见。例如，在文档或图片分类中常用的 bag of words 模型里，如果词典的大小是一百万，那么每个文档将由一百万维的向量来表示。高维度带来的的一个问题就是计算量：在一百万维的空间中，即使计算向量的内积这样的基本操作也会是非常费力的。不过，如果向量是稀疏的的话（事实上在 bag of words 模型中文档向量通常都是非常稀疏的），例如两个向量分别只有和个非零元素，那么计算内积可以只使用次乘法完成。因此稀疏性对于解决高维度数据的计算量问题是非常有效的。

当然高维度带来的问题不止是在计算量上。例如在许多生物相关的问题中，数据的维度非常高，但是由于收集数据需要昂贵的实验，因此可用的训练数据却相当少，这样的问题通常称为“small , large problem”——我们一般用表示数据点的个数，用表示变量的个数，即数据维度。当的时候，不做任何其他假设或者限制的话，学习问题基本上是没法进行的。因为如果用上所有变量的话，越大，通常会导致模型越复杂，但是反过来有很小，于是就会出现很严重的 overfitting 问题。例如，最简单的线性回归模型：

使用 square loss 来进行学习的话，就变成最小化如下的问题

这里是数据矩阵，而是由标签组成的列向量。该问题具有解析解

eq: 1 »

然而，如果的话，矩阵将会不是满秩的，而这个解也没法算出来。或者更确切地说，将会有无穷多个解。也就是说，我们的数据不足以确定一个解，如果我们从所有可行解里随机选一个的话，很可能并不是真正好的解，总而言之，我们 overfitting 了。

解决 overfitting 最常用的办法就是 regularization ，例如著名的 ridge regression 就是添加一个 regularizer ：

直观地来看，添加这个 regularizer 会使得模型的解偏向于 norm 较小的。从凸优化的角度来说，最小化上面这个等价于如下问题：

其中是和一一对应的是个常数。也就是说，我们通过限制的 norm 的大小实现了对模型空间的限制，从而在一定程度上（取决于的大小）避免了 overfitting 。不过 ridge regression 并不具有产生稀疏解的能力，得到的系数仍然需要数据中的所有特征才能计算预测结果，从计算量上来说并没有得到改观。

不过，特别是在像生物或者医学等通常需要和人交互的领域，稀疏的解除了计算量上的好处之外，更重要的是更具有“可解释性”。比如说，一个病如果依赖于 5 个变量的话，将会更易于医生理解、描述和总结规律，但是如果依赖于 5000 个变量的话，基本上就超出人肉可处理的范围了。

在这里引入稀疏性的方法是用 regularization 代替 regularization，得到如下的目标函数

eq: 2 »

该问题通常被称为 LASSO (least absolute shrinkage and selection operator) 。LASSO 仍然是一个 convex optimization 问题，不过不再具有解解析解。它的优良性质是能产生稀疏性，导致中许多项变成零。

可是，为什么它能产生稀疏性呢？这也是一直让我挺感兴趣的一个问题，事实上在之前申请学校的时候一次电话面试中我也被问到了这个问题。我当时的回答是背后的理论我并不是很清楚，但是我知道一个直观上的理解。下面我们就先来看一下这个直观上的理解。

首先，很 ridge regression 类似，上面形式的 LASSO 问题也等价于如下形式：

也就是说，我们将模型空间限制在的一个 -ball 中。为了便于可视化，我们考虑两维的情况，在平面上可以画出目标函数的等高线，而约束条件则成为平面上半径为的一个 norm ball 。等高线与 norm ball 首次相交的地方就是最优解。如图 (fig: 1) 所示：

fig: 1 »

-ball meets quadratic function. -ball has corners. It’s very likely that the meet-point is at one of the corners.
-ball meets quadratic function. -ball has no corner. It is very unlikely that the meet-point is on any of axes."

可以看到，-ball 与 -ball 的不同就在于他在和每个座标轴相交的地方都有“角”出现，而目标函数的测地线除非位置摆得非常好，大部分时候都会在角的地方相交。注意到在角的位置为产生稀疏性，例如图中的相交点就有，而更高维的时候（想象一下三维的 -ball 是什么样的？）除了角点以外，还有很多边的轮廓也是既有很大的概率成为第一次相交的地方，又会产生稀疏性。

相比之下，-ball 就没有这样的性质，因为没有角，所以第一次相交的地方出现在具有稀疏性的位置的概率就变得非常小了。这就从直观上来解释了为什么 regularization 能产生稀疏性，而 regularization 不行的原因了。

不过，如果只限于 intuitive 的解释的话，就不那么好玩了，但是背后完整的理论又不是那么容易能够搞清楚的，既然这次的标题是 Basics ，我们就先来看一个简单的特殊情况好了。

接下来我们考虑 orthonormal design 的情况：。然后看看 LASSO 的解具体是什么样子。注意 orthonormal design 实际上是要求特征之间相互正交。这可以通过对数据进行 PCA 以及模长 normalize 来实现。

注意到 LASSO 的目标函数 (eq: 2) 是 convex 的，根据 KKT 条件，在最优解的地方要求 gradient 。不过这里有一点小问题： -norm 不是光滑的，不存在 gradient ，所以我们需要用一点 subgradient 的东西。

def: 1 »

定义 subgradient; subdifferential

对于在维欧氏空间中的凸开子集上定义的实值函数，一个向量维向量称为在一点处的 subgradient ，如果对于任意，满足

由在点处的所有 subgradient 所组成的集合称为处的 subdifferential ，记为。

注意 subgradient 和 subdifferential 只是对凸函数定义的。例如一维的情况，，在处的 subdifferential 就是这个区间（集合）。注意在的 gradient 存在的点，subdifferential 将是由 gradient 构成的一个单点集合。这样就将 gradient 的概念加以推广了。这个推广有一个很好的性质。