Все, что вам нужно, возвращается функцией силуэта.Просто запишите это и суммируйте так, как хотите.Вот пример использования встроенных данных радужной оболочки.
library(cluster)
Iris_KM3 = kmeans(iris[,1:4],3)
SIL = silhouette(Iris_KM3$cluster, dist(iris[,1:4]))
aggregate(SIL[,3], list(SIL[,1]), mean)
Group.1 x
1 1 0.07624005
2 2 0.49471909
3 3 0.62148628
Если вы запустите приведенный выше код, попробуйте набрать SIL
или str(SIL)
, чтобы увидеть, что дает вам функция.