Illustrate:
hc <- hclust(dist(crs$dataset[1:17,c(2,7,9:10,12)]), "ave") plot(hc) This looks like 4 clusters. km <- kmeans(crs$dataset[1:17,c(2,7,9:10,12)], 4) km$cluster