聚类评估算法-轮廓系数(Silhouette Coefficient )

轮廓系数(Silhouette Coefficient),是聚类效果好坏的一种评价方式。最早由 Peter J. Rousseeuw 在 1986 提出。它结合内聚度和分离度两种因素。可以用来在相同原始数据的基础上用来评价不同算法、或者算法不同运行方式对聚类结果所产生的影响。

         方法:

                  1,计算样本i到同簇其他样本的平均距离ai。ai 越小,说明样本i越应该被聚类到该簇。将ai 称为样本i的簇内不相似度

                            簇C中所有样本的a i 均值称为簇C的簇不相似度。

                  2,计算样本i到其他某簇Cj 的所有样本的平均距离bij,称为样本i与簇Cj 的不相似度。定义为样本i的簇间不相似度:bi =min{bi1, bi2, ..., bik}

                            bi越大,说明样本i越不属于其他簇。

                   3,根据样本i的簇内不相似度a i 和簇间不相似度b i ,定义样本i的轮廓系数

                   4,判断:

                            si接近1,则说明样本i聚类合理;

                            si接近-1,则说明样本i更应该分类到另外的簇;

                            若si 近似为0,则说明样本i在两个簇的边界上。

 

         所有样本的s i 的均值称为聚类结果的轮廓系数,是该聚类是否合理、有效的度量。

  • 48
    点赞
  • 154
    收藏
    觉得还不错? 一键收藏
  • 5
    评论
好的,我理解了您的问题。您想知道如何使用 DIANA 聚类算法对鸢尾花数据集进行聚类,并将轮廓系数进行可视化展示。 首先,我们需要加载鸢尾花数据集和必要的库。可以使用以下代码进行加载: ``` from sklearn.datasets import load_iris from sklearn.cluster import AgglomerativeClustering from sklearn.metrics import silhouette_score import numpy as np import matplotlib.pyplot as plt # 加载鸢尾花数据集 iris = load_iris() X = iris.data y = iris.target ``` 接下来,我们可以使用 DIANA 聚类算法对数据进行聚类。DIANA 算法是一种层次聚类算法,可以通过以下代码进行使用: ``` # 使用 DIANA 算法进行聚类 diana = AgglomerativeClustering(n_clusters=3, linkage='single').fit(X) ``` 在上述代码中,我们将数据分为 3 类,并使用 DIANA 算法进行聚类。`linkage` 参数指定了距离计算的方式,这里使用 `single` 表示使用最近邻距离计算。 接下来,我们可以计算轮廓系数,并将其进行可视化展示。可以使用以下代码进行展示: ``` # 计算轮廓系数 silhouette_avg = silhouette_score(X, diana.labels_) sample_silhouette_values = silhouette_samples(X, diana.labels_) # 可视化轮廓系数 fig, ax = plt.subplots() y_lower = padding = 2 for i in range(3): ith_cluster_silhouette_values = sample_silhouette_values[diana.labels_ == i] ith_cluster_silhouette_values.sort() size_cluster_i = ith_cluster_silhouette_values.shape[0] y_upper = y_lower + size_cluster_i color = plt.cm.get_cmap("Spectral")(float(i) / 3) ax.fill_betweenx(np.arange(y_lower, y_upper), 0, ith_cluster_silhouette_values, facecolor=color, edgecolor=color, alpha=0.7) ax.text(-0.05, y_lower + 0.5 * size_cluster_i, str(i)) y_lower = y_upper + padding ax.axvline(x=silhouette_avg, color="red", linestyle="--") ax.set_yticks([]) ax.set_xlabel("Silhouette coefficient values") ax.set_ylabel("Cluster labels") plt.show() ``` 在上述代码中,我们使用 `silhouette_score` 函数计算了轮廓系数,使用 `silhouette_samples` 函数获取了每个数据点的轮廓系数,并将轮廓系数进行可视化展示。其中,每个条形图表示一个数据点,条形图的颜色表示所属的聚类,条形图的长度表示轮廓系数的大小。 最后,我们使用 `ax.axvline` 函数在轮廓系数图中添加了一个红色虚线,表示平均轮廓系数。 希望这个回答能够帮助到您!如果您有任何问题,请随时问我。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 5
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值