Difference between revisions of "Cluster.R"
From Organic Design wiki
m |
m (another version) |
||
| Line 1: | Line 1: | ||
| − | clusterPlot <- function(x, dmethod = "euclidian", hmethod="complete", ...) { | + | clusterPlot <- function(x, dmethod = "euclidian", hmethod = "complete", ...) { |
| − | oldpar <- par(mfrow=c(2,1)) | + | oldpar <- par(mfrow=c(2,1)) |
| − | xDist <- dist(x, method=dmethod) | + | xDist <- dist(x, method = dmethod) |
| − | + | xClust <- hclust(xDist, method = hmethod) | |
| − | |||
| − | xClust <- hclust(xDist, method=hmethod) | ||
plclust(xClust, sub="") | plclust(xClust, sub="") | ||
| − | + | image(list(x = 1:nrow(x), y = 1:ncol(x), z = x[clust$order,]), | |
| − | image(x = 1:nrow(x), y = 1:ncol(x), z = x[ | + | col=heat.colors(nrow(x)), ...) |
par(oldpar) | par(oldpar) | ||
invisible() | invisible() | ||
} | } | ||
| + | |||
| + | clusterPlot(x) | ||
# Test | # Test | ||
Revision as of 03:24, 16 November 2006
clusterPlot <- function(x, dmethod = "euclidian", hmethod = "complete", ...) {
oldpar <- par(mfrow=c(2,1))
xDist <- dist(x, method = dmethod)
xClust <- hclust(xDist, method = hmethod)
plclust(xClust, sub="")
image(list(x = 1:nrow(x), y = 1:ncol(x), z = x[clust$order,]),
col=heat.colors(nrow(x)), ...)
par(oldpar)
invisible()
}
clusterPlot(x)
- Test
x <- matrix(rnorm(40), nc=4) colnames(x) <- LETTERS[1:4] rownames(x) <- letters[1:10] clusterPlot(x)
