[R] how can I reorder a dendrogram?

Thomas Walter Thomas.Walter at embl.de
Thu Mar 20 11:39:13 CET 2008


Oups ... sorry, actually I found the problem.

The problem lies not with the reordering but with the construction of 
the dendrogram: if a hierarchical method is used, there are never more 
than 2 branches for a node. Therefore, the reordering does not give the 
expected result. I will have to try something different ...

Thanks anyway.

Thomas.


Thomas Walter wrote:

>Hi!
>
>I am trying to reorder a dendrogram via reorder.dendrogram. However, I 
>observed some problems with this, and I will illustrate them with an 
>example.
>
>Take the following clustering problem:
>
>datamatrix <- matrix(c(2,2,2.5,2,1.5,2,2,1.5,2,2.5, 
>6,2,6.5,2,5.5,2,6,1.5,6,2.5, 4,4,4.5,4,3.5,4,4,3.5,4,4.5), ncol=2, 
>byrow=TRUE)
>distmatrix <- dist(datamatrix, method="manhattan")
>hc <- hclust(distmatrix, method="single")
>dendro <- as.dendrogram(hc)
>
>The datamatrix contains three equidistant (for manhattan distance) 
>clusters, each of which contains 5 points.
>Now, I want to impose an order:
>
>weights <- c(2.0, 2.0, 2.0, 2.0, 2.0, 6.0, 6.0, 6.0, 6.0, 6.0, 4.0, 4.0, 
>4.0, 4.0, 4.0)
>ddd <- reorder(dendro, weights, agglo.FUN=mean)
>
>but if you compare the order of ddd with dendro, you see no change:
>
>unlist(ddd)
> [1] 15 14 13 11 12  5  4  3  1  2 10  9  8  6  7
>
>unlist(dendro)
> [1] 15 14 13 11 12  5  4  3  1  2 10  9  8  6  7
>
>I would have expected something like:
> 5  4  3  1  2 15 14 13 11 12 10  9  8  6  7
>
>or something of the sort. (I still do not know, if the order should be 
>ascending or descending, but in the obtained result, it is neither nor). 
>I do not see, where my mistake is ...
>
>Thanks for your advice!
>
>Thomas.
>
>______________________________________________
>R-help at r-project.org mailing list
>https://stat.ethz.ch/mailman/listinfo/r-help
>PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
>and provide commented, minimal, self-contained, reproducible code.
>
>  
>


-- 
-------------------------
Thomas Walter

EMBL Heidelberg
Meyerhofstrasse 1
69117 Heidelberg
Germany

Tel: +49 (0)6221 387-8857



More information about the R-help mailing list