机器学习-04-降维

降维

我们知道，解决过拟合的问题除了正则化和添加数据之外，降维就是最好的方法。降维的思路来源于维度灾难的问题，我们知道 \(n\) 维球的体积为：

\[CR^n \]

那么在球体积与边长为 \(2R\) 的超立方体比值为：

\[\lim\limits_{n\rightarrow0}\frac{CR^n}{2^nR^n}=0 \]

这就是所谓的维度灾难，在高维数据中，主要样本都分布在立方体的边缘，所以数据集更加稀疏。(个人认为这个例子不好，其实本质上讲就是空间大小随着维度是指数级增长的，但是我们的数据集顶多只能线性增长，这就导致过拟合的可能性大大增加。)

降维的算法分为：

1. 直接降维，特征选择
2. 线性降维，PCA，MDS等
3. 分线性，流形包括 Isomap，LLE 等

为了方便，我们首先将协方差矩阵（数据集）写成中心化的形式：

\[\begin{align*}S&=\frac{1}{N}\sum\limits_{i=1}^N(x_i-\overline{x})(x_i-\overline{x})^T\nonumber\\ &=\frac{1}{N}(x_1-\overline{x},x_2-\overline{x},\cdots,x_N-\overline{x})(x_1-\overline{x},x_2-\overline{x},\cdots,x_N-\overline{x})^T\nonumber\\ &=\frac{1}{N}(X^T-\frac{1}{N}X^T\mathbb{I}_{N1}\mathbb{I}_{N1}^T)(X^T-\frac{1}{N}X^T\mathbb{I}_{N1}\mathbb{I}_{N1}^T)^T\nonumber\\ &=\frac{1}{N}X^T(E_N-\frac{1}{N}\mathbb{I}_{N1}\mathbb{I}_{1N})(E_N-\frac{1}{N}\mathbb{I}_{N1}\mathbb{I}_{1N})^TX\nonumber\\ &=\frac{1}{N}X^TH_NH_N^TX\nonumber\\ &=\frac{1}{N}X^TH_NH_NX=\frac{1}{N}X^THX \end{align*} \]

这个式子利用了中心矩阵 \( H\)的对称性，这也是一个投影矩阵。

线性降维-主成分分析 PCA

损失函数

主成分分析中，我们的基本想法是将所有数据投影到一个子空间中，从而达到降维的目标，为了寻找这个子空间，我们基本想法是：

1. 所有数据在子空间中更为分散
2. 损失的信息最小，即：在补空间的分量少

其实一句话说白了，就是做一个基变换，使得样本在新的基上的投影的方差越大越好。

原来的数据很有可能各个维度之间是相关的，于是我们希望找到一组 \(p\) 个新的线性无关的单位基 \(u_i\)，降维就是取其中的 \(q\) 个基。于是对于一个样本 \(x_i\)，经过这个坐标变换后：

\[\hat{x_i}=\sum\limits_{i=1}^p(u_i^Tx_i)u_i=\sum\limits_{i=1}^q(u_i^Tx_i)u_i+\sum\limits_{i=q+1}^p(u_i^Tx_i)u_i \]

对于数据集来说，我们首先将其中心化然后再取上面的式子的第一项，并使用其系数的平方平均作为损失函数并最大化：

\[\begin{align*}J&=\frac{1}{N}\sum\limits_{i=1}^N\sum\limits_{j=1}^q((x_i-\overline{x})^Tu_j)^2\nonumber\\ &=\sum\limits_{j=1}^qu_j^TSu_j\ ,\ s.t.\ u_j^Tu_j=1 \end{align*} \]

由于每个基都是线性无关的，于是每一个 \(u_j\) 的求解可以分别进行，使用拉格朗日乘子法：

\[\mathop{argmax}_{u_j}L(u_j,\lambda)=\mathop{argmax}_{u_j}\space u_j^TSu_j+\lambda(1-u_j^Tu_j) \]

于是：

\[Su_j=\lambda u_j \]

可见，我们需要的基就是协方差矩阵的特征向量。损失函数最大那就取特征值的前 \(q\) 个最大值和其对应的特征向量。

下面看其损失的信息最少这个条件，同样适用系数的平方平均作为损失函数，并最小化：

\[\begin{align*}J&=\frac{1}{N}\sum\limits_{i=1}^N\sum\limits_{j=q+1}^p((x_i-\overline{x})^Tu_j)^2\nonumber\\ &=\sum\limits_{j=q+1}^pu_j^TSu_j\ ,\ s.t.\ u_j^Tu_j=1 \end{align*} \]

同样的：

\[\mathop{argmin}_{u_j}L(u_j,\lambda)=\mathop{argmin}_{u_j}\space u_j^TSu_j+\lambda(1-u_j^Tu_j) \]

损失函数最小取在本征值剩下的个最小的几个值。数据集的协方差矩阵可以写成 \(S=U\Lambda U^T\)，直接对这个表达式当然可以得到本征矢。

SVD 与 PCoA

下面使用实际训练时常常使用的 SVD 直接求得这个 \(q\) 个本征矢。

对中心化后的数据集进行奇异值分解：

\[HX=U\Sigma V^T,U^TU=E_N,V^TV=E_p,\Sigma:N\times p \]

于是：

\[S=\frac{1}{N}X^THX=\frac{1}{N}X^TH^THX=\frac{1}{N}V\Sigma^T\Sigma V^T \]

因此，我们直接对中心化后的数据集进行 SVD，就可以得到特征值和特征向量 \(V\)，在新坐标系中的坐标就是：

\[HX\cdot V \]

由上面的推导，我们也可以得到另一种方法 PCoA 主坐标分析，定义并进行特征值分解：

\[T=HXX^TH=U\Sigma\Sigma^TU^T \]

由于：

\[TU\Sigma=U\Sigma(\Sigma^T\Sigma) \]

于是可以直接得到坐标。这两种方法都可以得到主成分，但是由于方差矩阵是 \(p\times p\) 的，而 \(T\) 是 \(N\times N\) 的，所以对样本量较少的时候可以采用 PCoA的方法。

p-PCA

下面从概率的角度对 PCA 进行分析，概率方法也叫 p-PCA。我们使用线性模型，类似之前 LDA，我们选定一个方向，对原数据 \(x\in\mathbb{R}^p\) ，降维后的数据为 \(z\in\mathbb{R}^q,q<p\)。降维通过一个矩阵变换（投影）进行：

\[\begin{align*} z&\sim\mathcal{N}(\mathbb{O}_{q1},\mathbb{I}_{qq})\\ x&=Wz+\mu+\varepsilon\\ \varepsilon&\sim\mathcal{N}(0,\sigma^2\mathbb{I}_{pp}) \end{align*} \]

对于这个模型，我么可以使用EM算法进行学习，在进行推断的时候需要求得 \(p(z|x)\)，推断的求解过程和线性高斯模型类似。

\[\begin{align*} &p(z|x)=\frac{p(x|z)p(z)}{p(x)}\\ &\mathbb{E}[x]=\mathbb{E}[Wz+\mu+\varepsilon]=\mu\\ &Var[x]=WW^T+\sigma^2\mathbb{I}_{pp}\\ \Longrightarrow p(z|x)=\mathcal{N}(W^T(WW^T+&\sigma^2\mathbb{I})^{-1}(x-\mu),\mathbb{I}-W^T(WW^T+\sigma^2\mathbb{I})^{-1}W) \end{align*} \]

小结

降维是解决维度灾难和过拟合的重要方法，除了直接的特征选择外，我们还可以采用算法的途径对特征进行筛选，线性的降维方法以 PCA 为代表，在 PCA 中，我们只要直接对数据矩阵进行中心化然后求奇异值分解或者对数据的协方差矩阵进行分解就可以得到其主要维度。非线性学习的方法如流形学习将投影面从平面改为超曲面。

转载自：https://www.yuque.com/bystander-wg876/yc5f72?
对原文略作修正