R语言无监督学习:PCA主成分分析可视化

作者: 拓端tecdat | 来源:发表于2020-03-24 14:23 被阅读0次

原文链接:http://tecdat.cn/?p=9839

总览

在监督学习中,我们通常可以访问n个 观测值的p个 特征 集,并在相同观测值上测得的Y。

无监督学习是一组没有相关的变量Y的方法。在这里,我们重点介绍两种技术…

主成分分析:用于数据可视化或在其他监督学习方法之前进行预处理的工具。

聚类:发现数据中未知组的方法。

无监督学习的挑战

通常,无监督学习比主观学习更具挑战性,因为它更具主观性。分析没有简单的目标,例如预测响应。无监督学习通常用作探索性数据分析的一部分。此外,由于没有普遍接受的交叉验证或验证方法,因此很难评估获得的结果的准确性。简而言之,除了简单的直觉或手头上的过程的理论知识外,我们无法真正在无人监督的情况下检查工作。但是,无监督方法有许多用途:

通过识别患者亚组来了解癌症行为。

网站(尤其是电子商务)通常会根据您之前的活动尝试向您推荐产品。

Netflix电影推荐。

主成分分析

当出现大量相关变量时,主要成分使我们能够将集合概括为较少数量的代表变量,这些变量共同解释了原始集合中的大多数可变性。

主成分分析(PCA)是指计算主成分的过程,以及随后在理解数据中使用这些成分的过程。PCA还可以用作数据可视化的工具。

什么是主要成分

假设我们希望通过对一组p个 特征的测量值来可视化n个观测值,以用于探索性数据分析的一部分。具体来说,我们希望找到一种数据的低维表示形式,该表示形式可以捕获尽可能多的信息。PCA提供了一种执行此操作的方法。PCA会寻求少量尽可能有趣的维度,其中有趣的概念通过观察值在整个维度上的变化量来度量。

我们还可以通过利用主要组件来衡量丢失了多少信息。为此,我们可以计算每个主成分解释的方差比例(PVE)。通常最好将其解释为累积图,以便我们可以可视化每个成分的PVE和所解释的总方差。一

确定要使用的主成分数

总的来说,我们希望使用最少数量的主成分来充分理解数据。可以说,做到这一点的最好方法是在scree图中可视化数据,我们将在后面演示。它只是累积PVE的图。与我们选择其他学习技术的最佳调整参数的方式类似,查看百分比变化何时下降,这样,添加主要成分并不会真正增加大量的方差。我们可以结合一些对数据的理解来使用这种技术。

大多数统计方法都可以适应于使用主成分作为预测变量,这有时会导致噪声较小。

可视化

我们执行PCA。

states <- rownames(USArrests) states

need-to-insert-img

## [1] "Alabama" "Alaska" "Arizona" "Arkansas" ## [5] "California" "Colorado" "Connecticut" "Delaware" ## [9] "Florida" "Georgia" "Hawaii" "Idaho" ## [13] "Illinois" "Indiana" "Iowa" "Kansas" ## [17] "Kentucky" "Louisiana" "Maine" "Maryland" ## [21] "Massachusetts" "Michigan" "Minnesota" "Mississippi" ## [25] "Missouri" "Montana" "Nebraska" "Nevada" ## [29] "New Hampshire" "New Jersey" "New Mexico" "New York" ## [33] "North Carolina" "North Dakota" "Ohio" "Oklahoma" ## [37] "Oregon" "Pennsylvania" "Rhode Island" "South Carolina" ## [41] "South Dakota" "Tennessee" "Texas" "Utah" ## [45] "Vermont" "Virginia" "Washington" "West Virginia" ## [49] "Wisconsin" "Wyoming"

need-to-insert-img

数据集的列包含四个变量。

names(USArrests)

need-to-insert-img

## [1] "Murder" "Assault" "UrbanPop" "Rape"

need-to-insert-img

让我们来探讨一下数据。

kable(summary(USArrests))

need-to-insert-img

need-to-insert-img

我们可以看到数据具有不同的均值和方差。此外,这些变量是在完全不同的尺度上测量的。例如UrbanPop,以百分比为单位,每10万个人测量次数。如果我们不对数据进行标准化,那就麻烦了。

执行PCA 提供主成分载荷。

need-to-insert-img

我们已经可以确定每个主成分所代表的内容。例如,第一个部分似乎解释了与犯罪有关的信息与城市人口之间的差异。这也是第一个组成部分,从直观上来说,这是最大的差异。第二部分肯定解释了城市环境的影响,第三和第四部分显示了其他犯罪的区别。

我们可以绘制第一个主成分的图。

Biplot

need-to-insert-img

在这里我们可以看到很多信息。首先查看轴,轴上的PC1 x和轴上的PC2y。箭头显示了它们如何在两个维度上移动。黑色状态显示每个状态在PC方向上如何变化。例如,加利福尼亚州既有高犯罪率,又是城市人口最多的国家之一。

该$sdev属性输出每个组件的标准偏差。每个分量解释的方差可以通过对这些平方进行平方来计算:

## [1] 2.4802 0.9898 0.3566 0.1734

need-to-insert-img

然后,为了计算每个主成分解释的方差比例,我们先将其除以总方差。

## [1] 0.62006 0.24744 0.08914 0.04336

need-to-insert-img

在这里,我们看到第一PC解释了大约62%的数据,第二PC解释了大约24%的数据。我们还可以绘制此信息。

碎石图

par(mfrow=c(1,2)) plot(pve, xlab='Principal Component', ylab='Proportion of Variance Explained', ylim=c(0,1), type='b') plot(cumsum(pve), xlab='Principal Component', ylab='Cumuative Proportion of Variance Explained', ylim=c(0,1), type='b')

need-to-insert-img

need-to-insert-img

相关文章

网友评论

    本文标题:R语言无监督学习:PCA主成分分析可视化

    本文链接:https://www.haomeiwen.com/subject/fhpoyhtx.html