ng机器学习视频笔记（十二） ——PCA实现样本特征降维

ng机器学习视频笔记（十二）

——PCA实现样本特征降维

（转载请附上本文链接——linhxx）

 一、概述

         所谓降维（dimensionality reduction），即降低样本的特征的数量，例如样本有10个特征值，要降维成5个特征值，即通过一些方法，把样本的10个特征值映射换算成5个特征值。

         因此，降维是对输入的样本数据进行处理的，并没有对预测、分类的结果进行处理。

         降维的最常用的方法叫做主成分分析（PCA，principal component analysis）。最常用的业务场景是数据压缩、数据可视化。该方法只考虑样本的特征，不考虑样本的结果标签，因此PCA是一种无监督学习算法。

二、PCA基础

1、目的

         PCA的目的，是找到一个低维的超平面，让每个样本投影到该平面的距离的总和最短。通常样本投影之前还要进行均值化和归一化。

2、降维

         当需要使用PCA，将样本数据从二维降到一维，即可以理解成在平面上，找到一个向量，使得平面上的所有点到该向量的垂直距离的总和最短，如下图所示：

         三维降维到二维如下图所示：

         其中，每个样本到该向量的垂直距离，称为投影误差（projection error），u⁽ⁱ⁾表示第i个样本的投影误差，正值表示和向量同向，负值表示反向。下图的红线表示向量，蓝色的线表示投影误差。

         上图画出了红线和粉线，粉色的即错误的pca的结果，可以看出所有点到这个粉线的投影误差都非常大这个就是不正确的pca。而红色的线，相比之下，所有点到其的投影误差就非常小了。

3、定义

         对PCA更通用的定义：

n维特征，要降维到k维，泽需要找出k个向量u⁽¹⁾ ,u⁽²⁾ ,…u^(k)，让所有样本到由这些向量组成的超平面的距离之和最小。如下图所示：

4、PCA与线性回归区别

         上面PCA的例子，看起来非常像线性回归，然而实际上，PCA并不是线性回归。

         线性回归，其二维图像的含义是，对于1个特征x，输出结果是y，即线性回归的纵轴是输出的标签y。因此，线性回归最终找到的线，其目的是让所有的样本的点距离这条线的纵向距离的和最小。

         相比之下，PCA的图，是两个特征组成的图，x1和x2没有关联，只不过是一个样本的两个特征。其拟合的线，目的是使每个样本到这个线的垂直距离（即最短距离）的和最小。如下图所示，左边为线性回归的图，右边为PCA的图：

三、PCA计算过程

1、数据预处理

         1）把各特征进行归一化，令这些特征在一个数量级中。

         2）计算所有样本，每个特征值的均值。

         3）把每个样本中的每个特征值，换成原特征值减去该特征值的均值。

2、计算所有样本的协方差，生成协方差矩阵Σ（注意这里的Σ是矩阵的标记，而不是求和的标记），Σ是一个n*n维的矩阵。

3、计算Σ的特征值和特征向量。

4、将特征值按照从大到小的顺序排序，选择其中最大的k个，然后将其对应的k个特征向量分别作为列向量组成特征向量矩阵。

5、将样本点投影到选取的特征向量上。

6、说明

         对于octave语言，求出Σ后，调用函数[U,S,V]=svd(Sigma)，其中的U，是一系列由向量u⁽ⁱ⁾组成的矩阵，前面k个组成的矩阵即为PCA所需要的矩阵，这里把其称为U_reduce，注意这个U是一个正交矩阵，即U^T=U^-1，这个性质后面会用到。

         最终投影的结果Z=U_reduce^T*X_approx，这里的X不考虑x₀=1这一项。

         过程如下图所示：

四、选择合适的k

         k表示的是降维后的维度，诸如数据压缩等场景，没有明确的降维的维度要求，则需要有一个概念来衡量降维的效果。这里引入一个公式，可以衡量降维的效果。

         降维的效果 = 平均投影误差 / 总方差，如下图所示：

         当计算结果为0.01，表示样本的99%的方差性会被保留，这个值越大表示样本损失的特征越少。这里根据需要，可以设定不同的结果，例如要求0.05，则是保留95%的方差性。

         因此，判断k值的方法如下：

         1）k初始值设置为1。

2）计算Z、U_reduce、X_approx(即预处理后的输入样本)。

         3）按上面的公式，计算结果是否大于要求（例如要求是0.01）。

         4）如果不是，则k+1，然后重复2、3步骤，直到满足条件。

         这样计算量很大，在octave中，有简化算法，上面公式[U,S,V]=svd(Sigma)中的S，即特征值排序后组成的对角矩阵，则公式有简化算法，如下图所示：

五、重建特征值

         重建（reconstruction）特征值，即把被压缩的特征值还原回原来的维度。根据式子Z=U_reduce^T*X_approx，因此X_approx=( U_reduce^T)^-1*Z。

另外，由于上面提到，U是一个正交矩阵，其转置矩阵和逆矩阵的结果是一样的，因此公式可以写成X_approx=U_reduce*Z。

考虑到X_approx≈X，因此重建特征值后，结果即为X≈X_approx= U_reduce*Z，如下图所示：

六、PCA使用建议

1、使用

         在监督学习中，假设输入的样本是(x⁽¹⁾,y⁽¹⁾), (x⁽²⁾,y⁽²⁾)…(x^(m),y^(m))，现考虑使用logistic回归，用PCA来简化计算，方法如下：

         1）提取出x⁽¹⁾, x⁽²⁾,…x^(m)，即不考虑y。

         2）用PCA把其映射成z⁽¹⁾, z⁽²⁾,…z^(m)。

         3）拼接回y，形成(z ⁽¹⁾,y⁽¹⁾), (z ⁽²⁾,y⁽²⁾)…(z ^(m),y^(m))。

         4）进行参数计算，计算出θ，得到公式h_θ(z)=1/(1+e^-θTz)。

         5）对于新输入的x，映射成z后，带入上面的h(z)即可。

2、反例

         1）不宜用于解决过拟合

         PCA算法，可以减少特征值的数量，因此对于过拟合的缓解是有作用的。但是，考虑到还有更优的解决过拟合的方式——正则化，因此不要用PCA来解决过拟合。

         其中主要的问题，在于PCA的压缩过程，会丢失一些样本的特性，而正则化不会丢失太多的样本特性。

         2）不能无脑使用PCA

         PCA是有助于加速机器学习的学习过程，但是不能无脑的使用。建议所有的算法，都先用全量特征参与学习，只有发现真的速度太慢、占用太多的内存或硬盘，才考虑使用PCA。

3、实际应用

         PCA主要用于数据压缩和数据可视化。

         1）数据压缩

         数据压缩，主要目的是减少内存和硬盘的使用，并且加速机器学习的速度。对于此目的，需要用上面介绍的方式，设定一个压缩后数据丢失的容忍值，然后计算出最合适的最小的降维维度k。

         2）数据可视化

         数据可视化，通常只能在二维、三维，因此此时直接将数据降维到二维或者三维即可。

——written by linhxx

更多最新文章，欢迎关注微信公众号“决胜机器学习”，或扫描右边二维码。