数据分析系列 之python语言中的聚类分析

本文详细介绍了K-means算法的基本原理与实现步骤,并通过具体案例演示了如何使用Python中的scipy和scikit-learn库来实现K-means算法。此外,还探讨了模型的选择与评估方法。

1 基础算法
(1) K-means算法:对于给定的样本集,按照样本之间的距离大小,将样本集划分为K个簇。让簇内的点尽量紧密的连在一起,而让簇间的距离尽量的大。
(2) K-means算法是局部最优解,初始聚类中心一般是随机选择,有可能运行两次的结果稍有不同。
(3) 距离公式常采用欧式距离和余弦相似度公式,前者越小代表距离越小,后者越大代表越相似。

2 算法实现

import numpy as np
from scipy.cluster.vq import * vq,kmeans,whiten
listA=[86.0,72.0,94.0,83.0]
listB=[89.0,85.0,97.0,93.0]
listC=[86.0,76.0,96.0,82.0]
data=np.array([listA,listB,listC])
whiten=whiten(data)
#kmeans功能是对数据进行聚类,2代表类别,返回结果是元组
#A,_表示只需要取元组的第一个内容
centroids,_=kmeans(whiten,2)
#vq是矢量量化函数
result,_=vq(whiten,centroids)
print(result)

3 scikit learn工具包解决方法

import numpy as np
from sklearn.cluster import KMeans
from scipy.spatial.distance import cdist
import matplotlib.pyplot as plt
from sklearn import datasets
listA=[86.0,72.0,94.0,83.0]
listB=[89.0,85.0,97.0,93.0]
listC=[86.0,76.0,96.0,82.0]
X=np.array([listA,listB,listC])
#fit是一个训练过程
kmeans=KMeans(n_clusters=2).fit(X)
#predict是一个预测过程
pred=kmeans.predict(X)
print(pred)

#注意聚类和分类的区别
#直观理解起来就是:聚类是无标签的,分类是有标签的
#以下为分类的示例,利用SVM算法
from sklearn import datasets
from sklearn import svm
clf=svm.SVC(gamma=0.001,C=100.)
digits=datasets.load_digits()
#fit方法学习,data[:-1]从位置0到位置-1之前的数
clf.fit(digit.data[:-1],digits.target[:-1])
#predict方法预测,data[-1]最后一个位置的数
result=clf.predict(digit.data[-1])
print(result)

4 模型的选择和评估
“肘”方法:绘制出K值和SSE(误差平方和)的对应关系,找到趋于平缓的位置K值即为最佳聚类数。

参考资料:
https://www.icourse163.org/learn/NJU-1001571005?tid=1463102441&from=study#/learn/content?type=detail&id=1240380198&sm=1 用python玩转数据

本程序是在python中完成,基于sklearn.cluster中的k-means聚类包来实现数据的聚类,对于里面使用的数据格式如下:(注意更改程序中的相关参数) 138 0 124 1 127 2 129 3 119 4 127 5 124 6 120 7 123 8 147 9 188 10 212 11 229 12 240 13 240 14 241 15 240 16 242 17 174 18 130 19 132 20 119 21 48 22 37 23 49 0 42 1 34 2 26 3 20 4 21 5 23 6 13 7 19 8 18 9 36 10 25 11 20 12 19 13 19 14 5 15 29 16 22 17 13 18 46 19 15 20 8 21 33 22 41 23 69 0 56 1 49 2 40 3 52 4 62 5 54 6 32 7 38 8 44 9 55 10 70 11 74 12 105 13 107 14 56 15 55 16 65 17 100 18 195 19 136 20 87 21 64 22 77 23 61 0 53 1 47 2 33 3 34 4 28 5 41 6 40 7 38 8 33 9 26 10 31 11 31 12 13 13 17 14 17 15 25 16 17 17 17 18 14 19 16 20 17 21 29 22 44 23 37 0 32 1 34 2 26 3 23 4 25 5 25 6 27 7 30 8 25 9 17 10 12 11 12 12 12 13 7 14 6 15 6 16 12 17 12 18 39 19 34 20 32 21 34 22 35 23 33 0 57 1 81 2 77 3 68 4 61 5 60 6 56 7 67 8 102 9 89 10 62 11 57 12 57 13 64 14 62 15 69 16 81 17 77 18 64 19 62 20 79 21 75 22 57 23 73 0 88 1 75 2 70 3 77 4 73 5 72 6 76 7 76 8 74 9 98 10 90 11 90 12 85 13 79 14 79 15 88 16 88 17 81 18 84 19 89 20 79 21 68 22 55 23 63 0 62 1 58 2 58 3 56 4 60 5 56 6 56 7 58 8 56 9 65 10 61 11 60 12 60 13 61 14 65 15 55 16 56 17 61 18 64 19 69 20 83 21 87 22 84 23 41 0 35 1 38 2 45 3 44 4 49 5 55 6 47 7 47 8 29 9 14 10 12 11 4 12 10 13 9 14 7 15 7 16 11 17 12 18 14 19 22 20 29 21 23 22 33 23 34 0 38 1 38 2 37 3 37 4 34 5 24 6 47 7 70 8 41 9 6 10 23 11 4 12 15 13 3 14 28 15 17 16 31 17 39 18 42 19 54 20 47 21 68 22
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值