Если ваши данные выглядят примерно так:
num_clusters = 3
X, y = datasets.load_iris(return_X_y=True)
kmeans_model = KMeans(n_clusters=num_clusters, random_state=1).fit(X)
cluster_labels = kmeans_model.labels_
Вы можете использовать metrics.silhouette_samples
для вычисления коэффициентов силуэта для каждой выборки, а затем взять среднее значение для каждого кластера:
sample_silhouette_values = metrics.silhouette_samples(X, cluster_labels)
means_lst = []
for label in range(num_clusters):
means_lst.append(sample_silhouette_values[cluster_labels == label].mean())
print(means_lst)
[0.4173199215409322, 0.7981404884286224, 0.45110506043401194] # 1 mean for each of the 3 clusters