[R] hclust function

[Philippe Hupé]

I am trying to use hclust function from mva package but I can't 
understand what the variable height stands for when I use the ward 
criterion. Here is a small example :

library(mva)
x <- c(1,2,3,4)
cluster <- hclust(dist(x), method="ward")
cluster$height

cluster$height = 1 1 3 but it is supposed to be 0.5 0.5 4

the distance between clusters is Ni*Nj/(Ni+Nj)*||Gi-Gj||^2 and 
cluster$height does not match with this distance.

Does someone has any idea ?

Thanks in advance.

Philippe Hupé

-- 

--------------------------------------------------

Philippe Hupé
Institut Curie - Equipe Bioinformatique
26, rue d'Ulm - 75005 PARIS France
+33 (0)1 42 34 65 29

Philippe.Hupe at curie.fr <mailto:Philippe.Hupe at curie.fr>