Johannes Radinger
2014-Jul-09 12:06 UTC
Cutting hierarchical cluster tree at specific height fails
Hi, I'd like to cut a hierachical cluster tree calculated with hclust at a specific height. However ever get following error message: "Error in cutree(hc, h = 60) : the 'height' component of 'tree' is not sorted (increasingly)" Here is a working example to show that when specifing a height in cutree() the code fails. In contrast, specifying the number of clusters in cutree() works. What is the exact problem and how can I solve it? x <- c(rnorm(100,50,10),rnorm(100,200,25),rnorm(100,80,15)) y <- c(rnorm(100,50,10),rnorm(100,200,25),rnorm(100,150,25)) df <- data.frame(x,y) plot(df) hc <- hclust(dist(df,method = "euclidean"), method="centroid") plot(hc) df$memb <- cutree(hc, h = 60) # this does not work df$memb <- cutree(hc, k = 3) # this works! plot(df$x,df$y,col=df$memb) Thank you for your hints! Best regards, Johannes [[alternative HTML version deleted]]