如何手动创建树状图（或“hclust”）对象？ （在 R 中）

Posted 2023-03-12

【中文标题】如何手动创建树状图（或“hclust”）对象？ （在 R 中）

【英文标题】：How do I manually create a dendrogram (or "hclust") object ? (in R)

【发布时间】：2011-01-19 15:14:24

【问题描述】：

我有一个树状图作为图像提供给我。由于它不是很大，我可以“手动”将它构造成一个 R 对象。

所以我的问题是，当我只有树状图图像时，如何手动创建树状图（或“hclust”）对象？

我看到有一个名为“as.dendrogram”的函数，但我找不到如何使用它的示例。

（p.s：这篇文章是针对我来自here的问题）

非常感谢， 塔尔

【参考方案1】：

我认为您最好创建一个hclust 对象，然后使用as.dendrogram 将其转换为树状图，然后尝试直接创建树状图。查看?hclust 帮助页面以了解hclust 对象的元素的含义。

下面是一个简单的例子，有四个叶子 A、B、C 和 D，首先是 A-B，然后是 C-D，最后是 AB-CD：

a <- list()  # initialize empty object
# define merging pattern: 
#    negative numbers are leaves, 
#    positive are merged clusters (defined by row number in $merge)
a$merge <- matrix(c(-1, -2,
                    -3, -4,
                     1,  2), nc=2, byrow=TRUE ) 
a$height <- c(1, 1.5, 3)    # define merge heights
a$order <- 1:4              # order of leaves(trivial if hand-entered)
a$labels <- LETTERS[1:4]    # labels of leaves
class(a) <- "hclust"        # make it an hclust object
plot(a)                     # look at the result   

#convert to a dendrogram object if needed
ad <- as.dendrogram(a)

【参考方案2】：

完全手动...不建议用于大型数据集。

tree<-list();
attributes(tree)<-list(members=18,height=3);
class(tree)<-"dendrogram";

tree[[1]]<-list();
attributes(tree[[1]])<-list(members=4,height=2,edgetext="1");
 tree[[1]][[1]]<-list();
 attributes(tree[[1]][[1]])<-list(members=2,height=1,edgetext="H");
  tree[[1]][[1]][[1]]<-list();
  attributes(tree[[1]][[1]][[1]])<-list(members=1,height=0,edgetext="H",label="1HH",leaf=TRUE);
  tree[[1]][[1]][[2]]<-list();
  attributes(tree[[1]][[1]][[2]])<-list(members=1,height=0,edgetext="T",label="1HT",leaf=TRUE);
 tree[[1]][[2]]<-list();
 attributes(tree[[1]][[2]])<-list(members=2,height=1,edgetext="T");
  tree[[1]][[2]][[1]]<-list();
  attributes(tree[[1]][[2]][[1]])<-list(members=1,height=0,edgetext="H",label="1TH",leaf=TRUE);
  tree[[1]][[2]][[2]]<-list();
  attributes(tree[[1]][[2]][[2]])<-list(members=1,height=0,edgetext="T",label="1TT",leaf=TRUE);
tree[[2]]<-list();
attributes(tree[[2]])<-list(members=2,height=2,edgetext="2");
 tree[[2]][[1]]<-list();
 attributes(tree[[2]][[1]])<-list(members=1,height=0,edgetext="H",label="2H",leaf=TRUE);
 tree[[2]][[2]]<-list();
 attributes(tree[[2]][[2]])<-list(members=1,height=0,edgetext="T",label="2T",leaf=TRUE);
tree[[3]]<-list();
attributes(tree[[3]])<-list(members=4,height=2,edgetext="3");
 tree[[3]][[1]]<-list();
 attributes(tree[[3]][[1]])<-list(members=2,height=1,edgetext="H");
  tree[[3]][[1]][[1]]<-list();
  attributes(tree[[3]][[1]][[1]])<-list(members=1,height=0,edgetext="H",label="3HH",leaf=TRUE);
  tree[[3]][[1]][[2]]<-list();
  attributes(tree[[3]][[1]][[2]])<-list(members=1,height=0,edgetext="T",label="3HT",leaf=TRUE);
 tree[[3]][[2]]<-list();
 attributes(tree[[3]][[2]])<-list(members=2,height=1,edgetext="T");
  tree[[3]][[2]][[1]]<-list();
  attributes(tree[[3]][[2]][[1]])<-list(members=1,height=0,edgetext="H",label="3TH",leaf=TRUE);
  tree[[3]][[2]][[2]]<-list();
  attributes(tree[[3]][[2]][[2]])<-list(members=1,height=0,edgetext="T",label="3TT",leaf=TRUE);
tree[[4]]<-list();
attributes(tree[[4]])<-list(members=2,height=2,edgetext="4");
 tree[[4]][[1]]<-list();
 attributes(tree[[4]][[1]])<-list(members=1,height=0,edgetext="H",label="4H",leaf=TRUE);
 tree[[4]][[2]]<-list();
 attributes(tree[[4]][[2]])<-list(members=1,height=0,edgetext="T",label="4T",leaf=TRUE);
tree[[5]]<-list();
attributes(tree[[5]])<-list(members=4,height=2,edgetext="5");
 tree[[5]][[1]]<-list();
 attributes(tree[[5]][[1]])<-list(members=2,height=1,edgetext="H");
  tree[[5]][[1]][[1]]<-list();
  attributes(tree[[5]][[1]][[1]])<-list(members=1,height=0,edgetext="H",label="5HH",leaf=TRUE);
  tree[[5]][[1]][[2]]<-list();
  attributes(tree[[5]][[1]][[2]])<-list(members=1,height=0,edgetext="T",label="5HT",leaf=TRUE);
 tree[[5]][[2]]<-list();
 attributes(tree[[5]][[2]])<-list(members=2,height=1,edgetext="T");
  tree[[5]][[2]][[1]]<-list();
  attributes(tree[[5]][[2]][[1]])<-list(members=1,height=0,edgetext="H",label="5TH",leaf=TRUE);
  tree[[5]][[2]][[2]]<-list();
  attributes(tree[[5]][[2]][[2]])<-list(members=1,height=0,edgetext="T",label="5TT",leaf=TRUE);
tree[[6]]<-list();
attributes(tree[[6]])<-list(members=2,height=2,edgetext="6");
 tree[[6]][[1]]<-list();
 attributes(tree[[6]][[1]])<-list(members=1,height=0,edgetext="H",label="6H",leaf=TRUE);
 tree[[6]][[2]]<-list();
 attributes(tree[[6]][[2]])<-list(members=1,height=0,edgetext="T",label="6T",leaf=TRUE);


windows(width=3,rescale="fixed");
par(ps=8);
plot(rev(tree),center=TRUE,horiz=TRUE);

【讨论】：

谢谢。虽然正如你所说 - 不适合胆小的人:)