曲线聚类_【第 45 期】如何用聚类模型 (kmeans) 做数据分析？

Dear 丶 2023-01-06 14:54 232阅读 0赞

![9ee0c38459d363c3038eab2f54433ba3.png][]

正文开始~

k-means 属于无监督学习算法，无监督算法的内涵是观察无标签数据集自动发现隐藏结构和层次，在无标签数据中寻找隐藏规律。   聚类模型在数据分析当中的应用： **既可以作为一个单独过程，用于寻找数据内在规律，也可以作为分类等其他分析任务的前置探索** 。 例如我们想探寻我们产品站内都有哪些社交行为群体，刚开始拍脑门想可能并不会很容易，这时候可以根据用户属性、行为对用户进行聚类，根据结果将每个簇定义为一类社交群体，基于这些类训练后续的分类模型，给用户打标签后进行个性化推荐、运营。 ![b756815eabba9aa1692bd762c4face42.png][]   **k-means算法与距离** K-means聚类的目标，是将n个观测数据点按照一定标准划分到k个聚类中，数据点根据相似度划分。 每一个聚类有一个质心，质心是对聚类中所有点的位置求平均值得到的点。 每个观测点属于距离它最近的质心所代表的聚类。   模型最终会选择n个观测点到所属聚类质心 **距离平方和(损失函数)最小**的聚类方式作为模型输出。K-means聚类分析中，特征变量需要是 **数值变量，**以便于计算距离。   我们使用距离来测量两个样本的相似性，距离的实质是他将两个具有多维特征数据的样本的比较映射成一个数字，可以通过这个数字的大小来衡量距离。   几个常见距离计算方法：

*  欧几里得距离-直线距离，不适合高维度数据，对某一维度大数值差异更加敏感；
 *  曼哈顿距离-也叫出租车距离，用来标明两个点在标准坐标系上的绝对轴距总和，只计算水平或垂直距离，对某一维度大数值差异不敏感；
 *  Hamming距离-可用来测量含有分类值的向量之间的距离；
 *  余弦距离-通过计算两个向量的夹角余弦值来评估相似度，适用于结果与数据中特征的取值大小无关，而与不同特征之间比值有关的案例；

k-means的实质是每次都把质心移动到群内所有点的‘means’上，不是建立在距离这个基础上，而是建立在最小化方差和的基础上， **方差恰好是欧几里得距离平方，如果采用其他距离但依然去最小化方差和，会导致整个算法无法收敛，所以k-means使用欧几里得方法。** ** **

### **![8a6ad531745419b49df65bce71a6cab8.png][]** ###

### **k-means算法原理** ###

第一步：数据归一化、离群点处理后，随机选择k个聚类质心，k的选择下一节详细讲； 第二步：所有数据点关联划分到离自己最近的质心，并以此为基础划分聚类； 第三步：将质点移动到当前划分聚类包含所有数据点的中心(means)； 重复第二步、第三步n次，直到所有点到其所属聚类质心的距离平方和最小；  ![7a8f4309f0dcde1165eddacf3ebb5789.png][] ![81bb25e74f89e741e7edcb693c6d4822.png][] **多次随机：** 选择聚几类，则开始随机选择几个聚类质心。 基于不同随机初始质心(第一步的质心)的尝试： 尝试1：第一步→第二步→第三步 … (重复二、三步~10次)→聚类方式1 尝试2：第一步→第二步→第三步 … (重复二、三步~10次)→聚类方式2 … 在所有尝试结果中，选择所有点到其所属聚类质心的距离平方和(方差和)最小的聚类方式；

### **![88493aa1cbd9f20cfbe9e7f0208ce39c.png][]** ###

**k值选择方法**K值的选择是k-means算法很重要的一步，K值选择方法有肘部法则、拍脑袋法、gap statistic法、轮廓系数等，本篇主要介绍肘部法则及gap statistic两种常用方法。   **1、肘部法则**  我们可以尝试不同的K值，并将不同K值所对应的损失函数画成折线，横轴为K的取值，纵轴为距离平方和定义的损失函数，距离平方和随着K的变大而减小。当K=5时，存在一个拐点，像人的肘部一样，当k<5时，曲线极速下降，当K>5时，曲线趋于平稳，拐点即为K的最佳取值。 ![95ac3b85070e05ea140ae9b083c7329a.png][] **2、间隔统计量(gap statistic)**  手肘法则是强依赖经验的方法，而gap statistic方法则不强依赖经验，只需要找到最大gap statistic 所对应的K即可，我们继续使用上面的损失函数，当分为K组时，对应的损失函数为Dk，则gap statistic定义为： ![c2f88025b8f144dd6e42c44f11f92459.png][] 这里的E(logDk)是logDk的期望，通过蒙特卡洛模拟产生，我们在样本所在的区域内按照均匀分布随机产生和原始样本数一样多的随机样本，并对这些随机样本做k-means，得到对应的损失函数Dk，重复多次即可得出E(logDk)的近似值，从而最终可以计算gap statistic。而gap statistic取得最大值所对应的K就是最佳的分组数。如下图所示，K=3时，Gap(K)的取值最大，所以3为K的最佳取值。 ![e38b829a4060769f32fba1a63bafb371.png][]

## **![27c5e51889f70af98ec76c8fcd98ad34.png][]****k-means数据分析实战案例** ##

**案例背景：** O2O平台为了更好地为线下店面服务，增加一个增值服务，即利用自己拥有的地理位置数据为线下店面选址，数据如下： ![5ee20f7a58e3e2c0f81200483e71b806.png][]   每一条数据是一个兴趣点(POI - Point of Interest)的特征，具体指的是以这个位置为中心的500米半径圆里，各类商家/设施的数量，数据中我们隐藏掉了每个POI的具体名称、坐标、类型。选址的商家将试图从这些位置中选择一个作为下一个店面的位置。 ** ** 商家想知道这40个潜在店面位置之间是否有显著的差异。我们可以将所有POI按照相似程度，划分成几个类别？   **步骤：** ![45610176e70d5f4ad55a26ffaffed8a1.png][]

*  数据准备：数据获取、数据清洗、数据变换等步骤，重点是针对分析目的，进行特征选择以及特征标准化；
 *  数据建模：使用k-means算法进行数据建模；
 *  后续分析：聚类模型的特征描述分析，基于业务问题，进一步分析；

**1、读取数据** ![1755eab7a6c551a69ef764733906779f.png][] **2、特征选取** ![42ccbbfbf530c6b6260ccd3d668d1fb2.png][] **3、** **标准化** ![46d81d8fead9320e5ff39faab6b4ee95.png][] **4、** **k值选择** ![af2dadedc89b90ab3664d7808f8524cc.png][] **5、模型建立** ![fa74e0d7b8453661cfde2204e0b59734.png][] **5、** **聚类分析，对每一聚类进行进一步分析和描述** ![c0281bf977001ec7fe7ff8d458044364.png][] 根据描述分析，生成poi聚类结果画像，利用聚类，对人群/商品/行为标签进行精细划分，作为进一步商业决策的基础； ![a2929c6375930e71623f4eb33955595d.png][]

--------------------

k-means是面试常考算法，其算法原理简单，且参数很少，也是业界的常用算法。但是其缺点也很明显，例如其对离群点较敏感、容易产生局部最优(多次随机来解决)、需要提前确定k值、聚类结果依赖初始位置选择等。在接下来的文章中，会继续为大家介绍其他的常用聚类算法，例如层次聚类、密度聚类(dbscan)，敬请期待吧～

(完)

--------------------

### 上期文章： ###

**【第 44 期】知乎部署系统演进**

### 相关文章： ###

▪【第 36 期】如何用决策树模型做数据分析？▪【第 35 期】如何用线性回归模型做数据分析？

![5ad7ba4b60e696ca596cd7d033228024.png][]

[9ee0c38459d363c3038eab2f54433ba3.png]: /images/20221119/4a25abfdf8b24283b26fc28b804725ee.png
[b756815eabba9aa1692bd762c4face42.png]: /images/20221119/f00590c58c1645029a133a8572ab7287.png
[8a6ad531745419b49df65bce71a6cab8.png]: /images/20221119/698b3e334f6d475384d4aa7d5e770497.png
[7a8f4309f0dcde1165eddacf3ebb5789.png]: https://img-blog.csdnimg.cn/img_convert/7a8f4309f0dcde1165eddacf3ebb5789.png
[81bb25e74f89e741e7edcb693c6d4822.png]: https://img-blog.csdnimg.cn/img_convert/81bb25e74f89e741e7edcb693c6d4822.png
[88493aa1cbd9f20cfbe9e7f0208ce39c.png]: https://img-blog.csdnimg.cn/img_convert/88493aa1cbd9f20cfbe9e7f0208ce39c.png
[95ac3b85070e05ea140ae9b083c7329a.png]: https://img-blog.csdnimg.cn/img_convert/95ac3b85070e05ea140ae9b083c7329a.png
[c2f88025b8f144dd6e42c44f11f92459.png]: https://img-blog.csdnimg.cn/img_convert/c2f88025b8f144dd6e42c44f11f92459.png
[e38b829a4060769f32fba1a63bafb371.png]: https://img-blog.csdnimg.cn/img_convert/e38b829a4060769f32fba1a63bafb371.png
[27c5e51889f70af98ec76c8fcd98ad34.png]: https://img-blog.csdnimg.cn/img_convert/27c5e51889f70af98ec76c8fcd98ad34.png
[5ee20f7a58e3e2c0f81200483e71b806.png]: https://img-blog.csdnimg.cn/img_convert/5ee20f7a58e3e2c0f81200483e71b806.png
[45610176e70d5f4ad55a26ffaffed8a1.png]: https://img-blog.csdnimg.cn/img_convert/45610176e70d5f4ad55a26ffaffed8a1.png
[1755eab7a6c551a69ef764733906779f.png]: https://img-blog.csdnimg.cn/img_convert/1755eab7a6c551a69ef764733906779f.png
[42ccbbfbf530c6b6260ccd3d668d1fb2.png]: https://img-blog.csdnimg.cn/img_convert/42ccbbfbf530c6b6260ccd3d668d1fb2.png
[46d81d8fead9320e5ff39faab6b4ee95.png]: https://img-blog.csdnimg.cn/img_convert/46d81d8fead9320e5ff39faab6b4ee95.png
[af2dadedc89b90ab3664d7808f8524cc.png]: https://img-blog.csdnimg.cn/img_convert/af2dadedc89b90ab3664d7808f8524cc.png
[fa74e0d7b8453661cfde2204e0b59734.png]: https://img-blog.csdnimg.cn/img_convert/fa74e0d7b8453661cfde2204e0b59734.png
[c0281bf977001ec7fe7ff8d458044364.png]: https://img-blog.csdnimg.cn/img_convert/c0281bf977001ec7fe7ff8d458044364.png
[a2929c6375930e71623f4eb33955595d.png]: https://img-blog.csdnimg.cn/img_convert/a2929c6375930e71623f4eb33955595d.png
[5ad7ba4b60e696ca596cd7d033228024.png]: https://img-blog.csdnimg.cn/img_convert/5ad7ba4b60e696ca596cd7d033228024.png