使用ape包绘制进化树报错:tip不足且对象为多棵树
解决ape包绘制10k Trees进化树的报错问题
你的问题核心是:从10k Trees下载的Nexus文件包含100棵树的集合(multiPhylo对象),而plot.phylo()默认只能处理单棵树(phylo对象),直接传入集合就会触发"少于2个末端"的错误。
快速解决步骤:
先确认对象类型,执行:
class(primatetree)输出应该是
"multiPhylo",证明这是多树集合。绘制单棵树:
提取集合里的某一棵树(比如第一棵)再绘图,执行:plot.phylo(primatetree[[1]])这样就能正常显示单棵树的结构,也会输出树的末端数、节点数信息。
批量绘制所有树:
如果需要看全部100棵树,用循环批量生成:# 循环遍历每棵树,每棵树开新窗口绘制 for (tree_idx in 1:length(primatetree)) { dev.new() plot.phylo(primatetree[[tree_idx]], main = paste("第", tree_idx, "棵树")) }或者直接调用ape的多树绘图函数,自动排版所有树:
plot(primatetree)(注意:100棵树排版会比较密集,适合小窗口查看概览)
补充说明:
从10k Trees下载的Nexus文件默认包含多棵抽样树,加载后会自动生成multiPhylo对象,这是正常情况。如果只需要单棵树,加载后通过[[索引]]提取即可,不需要修改下载的文件。
内容的提问来源于stack exchange,提问作者Rob
相关产品推荐
相关产品推荐

