source("clus1.s") # two uneven clusters n <- 10 x <- c(rnorm(70*n,-2,1.5),rnorm(30*n,3,0.5)) break.kmeans(x,2) break.hclust(x,2) # big cluster vs. small cluster with same variance x <- c(rnorm(70*n,-3,1),rnorm(30*n,3,1)) break.kmeans(x,6) # three clusters x <- c(rnorm(70*n,-2,0.5),rnorm(100*n,2.5,1.5),rnorm(50*n,7,0.1)) break.kmeans(x,3) break.hclust(x,3) hc <- hclust(dist(x),method="ward") plot(hc) # in S: # hc <- hclust(dist(x),method="ave") # plclust(hc) plot.hclust.breaks(hc,x)