如何在R语言中基于指定结构的矩阵绘制树形图？

当然有合适的R包能搞定这个需求！结合你给出的矩阵结构（每行对应一个叶子节点的完整层级路径），我推荐用data.tree包来实现，它对这种路径型树形数据的处理特别顺手，下面一步步给你演示：

第一步：数据准备与格式转换 #

首先把你的矩阵读入R，然后将每行的层级值拼接成data.tree能识别的路径字符串（用/分隔层级）：

# 先定义你的矩阵M
M <- matrix(c(1,1,3,2,
              2,2,1,1,
              3,2,3,1,
              2,2,2,2,
              1,1,1,1,
              3,2,3,2,
              1,1,3,1,
              2,1,1,1),
            nrow = 8, byrow = TRUE,
            dimnames = list(paste0("Row", 1:8), paste0("Level", 1:4)))

# 加载data.tree包
library(data.tree)

# 将每行的层级值拼接成路径字符串，比如"1/1/3/2"
paths <- apply(M, 1, function(x) paste(x, collapse = "/"))
# 直接转换成树形结构对象
tree <- as.Node(data.frame(pathString = paths))

第二步：可视化树形图 #

快速绘图（data.tree自带函数） #

用data.tree的内置plot()函数就能快速生成树形图，还能调整节点大小、字体大小：

plot(tree, nodeSize = 10, fontSize = 12)

定制化绘图（配合DiagrammeR） #

如果想要更灵活的样式（比如节点颜色、形状、布局方向），可以搭配DiagrammeR包，它能生成更美观的矢量图：

library(DiagrammeR)

# 把data.tree对象转换成DiagrammeR的图结构
diagram <- ToDiagrammeRGraph(tree)

# 自定义样式：纵向布局、圆形节点、统一配色
diagram %>%
  add_global_graph_attrs(attr = "rankdir", value = "TB", attr_type = "graph") %>%
  add_global_graph_attrs(attr = "shape", value = "circle", attr_type = "node") %>%
  add_global_graph_attrs(attr = "color", value = "#617a89", attr_type = "node") %>%
  render_graph()

其他可选工具 #

• ggtree：如果你习惯用ggplot2的语法，可以把树形结构转换成phylo对象后用ggtree绘制，适合需要叠加其他ggplot图层的场景，但转换步骤稍多。
• igraph：可以先构建父子节点的边列表，再用igraph绘制，适合需要同时分析树形图结构属性的场景。

小提示 #

你的矩阵中第1、3层级有3种节点值，其余层级是2种，data.tree会自动识别共享的父节点，不需要手动去重；如果需要给不同层级的节点设置不同样式，可以用tree$Do()方法批量添加属性，再在可视化时调用。

内容的提问来源于stack exchange，提问作者Mark