r - 从 R 中的切割树状图中提取标签成员/分类(即 : a cutree function for dendrogram)

我试图从 R 中的树状图中提取分类，我在某个高度有 cut。在 cutree 对象上使用 hclust 很容易做到这一点，但我不知道如何在 dendrogram 对象上执行此操作。

此外，我不能只使用来自原始 hclust 的集群，因为(令人沮丧)，来自 cutree 的类的编号与 cut 的类的编号不同。

hc <- hclust(dist(USArrests), "ave")

classification<-cutree(hc,h=70)

dend1 <- as.dendrogram(hc)
dend2 <- cut(dend1, h = 70)


str(dend2$lower[[1]]) #group 1 here is not the same as
classification[classification==1] #group 1 here

有没有办法让分类相互映射，或者从 dendrogram 对象(也许巧妙地使用 dendrapply ？)以更像 0x2518124313431 的格式提取较低的分支成员资格？

最佳答案

我建议您使用 dendextend 包中的 cutree 函数。它包括一个树状图方法(即: dendextend:::cutree.dendrogram )。

您可以从 its introductory vignette 了解有关该软件包的更多信息。

我应该补充一点，虽然您的函数( classify )很好，但使用来自 dendextend 的 cutree 有几个优点:

它还允许您使用特定的 k(簇数)，而不仅仅是 h(特定的高度)。

与你在 hclust 上从 cuttree 得到的结果一致(classify 不会)。

通常会更快。

以下是使用代码的示例:

# Toy data:
hc <- hclust(dist(USArrests), "ave")
dend1 <- as.dendrogram(hc)

# Get the package:
install.packages("dendextend")
library(dendextend)

# Get the package:
cutree(dend1,h=70) # it now works on a dendrogram
# It is like using:
dendextend:::cutree.dendrogram(dend1,h=70)

顺便说一下，在这个功能的基础上，dendextend允许用户做更多很酷的事情，比如基于切割树状图的颜色分支/标签:

dend1 <- color_branches(dend1, k = 4)
dend1 <- color_labels(dend1, k = 5)
plot(dend1)

r - 从 R 中的切割树状图中提取标签成员/分类(即 : a cutree function for dendrogram)-LMLPHP

最后，这里还有一些代码来演示我的其他观点:

# This would also work with k:
cutree(dend1,k=4)

# and would give identical result as cutree on hclust:
identical(cutree(hc,h=70)  , cutree(dend1,h=70)  )
   # TRUE

# But this is not the case for classify:
identical(classify(dend1,70)   , cutree(dend1,h=70)  )
   # FALSE


install.packages("microbenchmark")
require(microbenchmark)
microbenchmark(classify = classify(dend1,70),
               cutree = cutree(dend1,h=70)  )
#    Unit: milliseconds
#        expr      min       lq   median       uq       max neval
#    classify  9.70135  9.94604 10.25400 10.87552  80.82032   100
#      cutree 37.24264 37.97642 39.23095 43.21233 141.13880   100
# 4 times faster for this tree (it will be more for larger trees)

# Although (if to be exact about it) if I force cutree.dendrogram to not go through hclust (which can happen for "weird" trees), the speed will remain similar:
microbenchmark(classify = classify(dend1,70),
               cutree = cutree(dend1,h=70, try_cutree_hclust = FALSE)  )
# Unit: milliseconds
#        expr       min        lq    median       uq      max neval
#    classify  9.683433  9.819776  9.972077 10.48497 29.73285   100
#      cutree 10.275839 10.419181 10.540126 10.66863 16.54034   100

如果您正在考虑改进此功能的方法，请在此处进行修补:

https://github.com/talgalili/dendextend/blob/master/R/cutree.dendrogram.R

我希望你或其他人会发现这个答案很有帮助。

关于r - 从 R 中的切割树状图中提取标签成员/分类(即 : a cutree function for dendrogram)，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/25452472/