在R语言中合并两棵系统发育树并保留根边的问题求助

问题原因 #

• 直接使用bind.tree(treeA, treeB)会合并两棵树的根节点，导致原有的根边被丢弃——合并后仅保留一个根节点，无法同时容纳两条根边。
• 你之前构建的anc_tree可能不是合法的phylo对象（比如结构不完整、为空），才会触发incorrect number of dimensions错误。

解决方案步骤 #

1. 构建合法的基础祖先树
   先创建一个包含根节点和两个子分支的基础树，作为连接两棵目标树的载体，同时可按需设置根边长度：

   library(ape)
   # 创建带两个临时末端节点的基础树，可自定义分支长度
   anc_tree <- read.tree(text = "(temp1:0.2, temp2:0.2);")
   # 给基础树添加根边（如果需要保留整体的根边）
   anc_tree$root.edge <- 0.3
   

2. 分别绑定两棵目标树
   将treeA和treeB依次绑定到基础树的两个临时节点上，这样两棵树的根边会被保留为连接基础节点与各自树根的分支：

   # 绑定treeA到临时节点temp1
   node_temp1 <- which(anc_tree$tip.label == "temp1")
   combined_tree <- bind.tree(anc_tree, treeA, where = node_temp1)
   
   # 绑定treeB到临时节点temp2（注意此时节点编号已更新，需重新获取）
   node_temp2 <- which(combined_tree$tip.label == "temp2")
   final_tree <- bind.tree(combined_tree, treeB, where = node_temp2)
   
   # 移除临时节点的标签（可选，根据需求保留或删除）
   final_tree$tip.label <- final_tree$tip.label[-c(node_temp1, node_temp2)]
   

3. 验证结果
   用plot(final_tree)查看树结构，确认两棵树的根边（原treeA和treeB的根边）以及基础树的根边都被保留。

关键注意事项 #

• 确保treeA和treeB都是带根边的有根树（可通过is.rooted(treeA)检查）。
• 基础树必须是结构完整的phylo对象，不能是空树或格式错误的树，否则会触发维度错误。

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