Cluster.R

From Organic Design wiki
Revision as of 23:17, 7 December 2006 by Sven (talk | contribs) (Add alt naming functionality)

clusterPlot <- function(x, rowNames = NULL, dmethod = "euclidian", hmethod="complete", ...) { 

 oldpar <- par(mfrow=c(2,1))  

 if(is.character(rowNames) & (length(rowNames) == nrow(x))) {
   rownames(x) <- rowNames
 }
 xDist <- dist(x, method=dmethod)
  1. 1) Tree dendrogram
 xClust <- hclust(xDist, method=hmethod)
 plclust(xClust, sub="")
  1. 2) Image plot...
 image(x = 1:nrow(x), y = 1:ncol(x), z = x[xClust$order, ], col=heat.colors(nrow(x)), ...)
 par(oldpar)
 invisible()

}

  1. Test

x <- matrix(rnorm(40), nc=4) colnames(x) <- LETTERS[1:4] rownames(x) <- letters[1:10] clusterPlot(x) clusterPlot(x, rowNames = LETTERS[1:10])

Retrieved from "https://organicdesign.nz/wiki/index.php?title=Cluster.R&oldid=42889"