如何基于使用R的子项的标签来标记树形图中的每个节点

时间:2012-11-20 22:14:33

标签: r dendrogram hierarchical-clustering

我在R中有一个树形图,其中每个叶子都有一个值。我喜欢通过对其子节点的值求和来定义每个节点的值。我熟悉dendrapply,但是我不知道如何在函数中访问节点的子节点以及如何递归地编写函数。

这是开头的代码:

library("stats")
library("fastcluster")
library("cluster")
D = rbind( + c(1,1,1,1,1), 
 + c(1,2,1,1,1),
 + c(2,2,2,2,2), 
 + c(3,4,5,6,9)

)
dnd = as.dendrogram(hclust.vector(D))

apply_text <<- function(n) {
   if (!is.leaf(n)) {

      attr(n, "edgetext") <- add the value of the branches
   }
   if (is.leaf(n)) {
      attr(n, "edgetext") <- 1
   }
   n
}

tmp <- dendrapply(dnd, apply_text)
plot(tmp)

1 个答案:

答案 0 :(得分:0)

这可能是一个答案,然而,它正在重新实现dendrapply。

apply_text <<- function(n){
  if (!is.leaf(n)) {
    cutversion = cut(n, h = attributes(n)$height)
    leftLabel = attr(apply_text(cutversion$lower[[1]]), "edgetext")  
    rightLabel= attr(apply_text(cutversion$lower[[2]]), "edgetext")
    attr(n, "edgetext") = as.numeric(as.character(leftLabel)) + as.numeric(as.character(rightLabel)) 
     }
  if(is.leaf(n)) {
    attr(n,"edgetext") <- 1
  }
    n
}

tmp <- dendrapply(dnd, apply_text)

是否有人知道如何删除标签上的多边形?其他人似乎也要求将其删除。有什么进展吗?