Hi there, I use cluster.stats to test which cluster is better, for example: library(fpc) data(xclara) a3 <- clara(xclara, 3) a4 <- clara(xclara, 4) d <- dist(xclara) val1 <- cluster.stats(d, a3$clust) val2 <- cluster.stats(d, a4$clust) barplot(c(val1$avg, val2$avg)) Which value in the output of cluster.stats can tell us 3 cluster is better than 4? I also try to test "kmeans" method with "clara" method by using cluster.stats. The same problem shows up. Thank you! [[alternative HTML version deleted]]