我正在使用identify.hclust在R中手动剪切从hclust创建的树状图。该函数的默认返回是每个组中观察值的in。我需要这些信息,但我也需要知道这个团队的高度。有什么办法可以做到这一点吗?非常感谢!
可重现数据:
set.seed(1)
dat = rnorm(100,0,1)
hca = hclust(dist(dat))
plot(hca, hang=-1, sub="", xlab="", labels=F)
heightsAndIDs = identify(hca) #Gives only IDs例如,我使用identify在以下高度剪切树状图,并希望获得分支的合并高度:
segments(3,2,8, col="red")
segments(15,1,18, col="green")
segments(20,1,24,col="blue")
segments(38,1.5,45,col="purple")
segments(75, 1.5, 82,col="cyan")https://stackoverflow.com/questions/41288886
复制相似问题