Hi: I have calculated the Silhouette score and Dunn score after hierarchical clustering for 3 clusters: #Distance measure d <- dist(USArrests, method = "euclidean") #Hierarchical clustering hc <- hclust(dist(USArrests), "ave") #calculating silhouette value for 3 clusters sil<- silhouette(cutree(hc, k=3), d) #calculating Dunn index for 3 clusters clus <- cutree(hc, 3) dun <- dunn(d, clus) How can the best of the two score be obtained? Is there any package to automatically obtain the optimal (or best) of the Silhouette and the Dunn scores ? Regards: John