data("Psych24")
## create a dist object and also get rid of the one negative entry in the
## correlation matrix
d <- as.dist(1 - abs(Psych24))
pimage(d)
## do hclust as in Ling (1973)
hc <- hclust(d, method = "complete")
plot(hc)
pimage(d, hc)
## use seriation
order <- seriate(d, method = "tsp")
#order <- seriate(d, method = "tsp", control = list(method = "concorde"))
pimage(d, order)
Run the code above in your browser using DataLab