Cluster.R
From Organic Design wiki
clusterPlot <- function(x, dmethod = "euclidian", hmethod="complete", ...) {
oldpar <- par(mfrow=c(2,1)) xDist <- dist(x, method=dmethod)
- 1) Tree dendrogram
xClust <- hclust(xDist) plclust(xDist)
- 2) Image plot...
image(x = 1:nrow(x), y = 1:ncol(x), z = x[xClust$order, ]) par(oldpar) invisible()
}
- Test
x <- matrix(rnorm(40), nc=4) clusterPlot(x)