有人知道我如何(a)获取TnT以打印出新树中的内部节点,或者我如何复制TnT命名过程,最好使用python(ete2 / ete3)?
TIA
答案 0 :(得分:1)
只是偶然发现了您的帖子。
TnT不允许使用内部节点ID或分支长度
这在两个方面都是不正确的; TNT确实允许这些。在做出这样的声明之前,您应该确保已经仔细阅读并理解了程序的文档。
对于newick格式的分支长度(在树n上),请结合使用“ ttag”命令和“ blength”命令,然后是“ export”:
ttag- ;; ttag =; blength * n;导出>文件名;
ttag命令允许操作/串联树标签。例如。如果您将树标记设置为打开并打印显示节点编号的树(使用“ naked-”命令),则可以使用“ tsave * filename”(打开树文件)将标记保存到新树中。通过“保存*”(保存标签):
naked-; ttag=; tplot n; tsav * con ; sav*;
所有这些都在TNT的在线帮助中进行了说明。
由于TNT还允许处理(通过脚本编写)分支长度,节点号等的值,因此您可能可以完成TNT本身所需的许多工作,而无需使用Python等。
对于内部节点编号,是的,TNT使用其自己的内部规则来确保对于相同进化枝,两个相同的树始终具有相同的编号。规则是:将根节点编号为T(分类单元数),然后从第一个分类单元开始向下编号节点(从T + 1开始,然后是T + 2,T + n,直到到达根)。移至下一个分类单元,并开始逐步向下进行根编号,直到找到根或已经为其分配了编号的节点为止。这样,您可以在需要的任何地方重现内部节点编号。