MEGA做为进化界常用的软件,虽然我自己没用过,还是需要支持一下的,它的输出有三种,newick, nexus和tabular,Newick自然不在话下,NEXUS同样也是,当然它也是非标准的NEXUS,跟BEAST的输出比较像,treeio也支持。最后就是tabular了,它实质上是把树的信息导出为表格,方便用户在excel中查看。而我要做的,就是把表格中树相关的信息,变成一个phylo对象,而相关的数据信息,整合在一起变成一个treedata对象。
假设有以下一个表格:
用最新的read.mega_tabular函数读进来是这样子滴:
一个表格,树和数据都有了。
当然假如你喜欢的话,树和数据是分分钟可以变成表格的,只需要使用tidytree包。
有treeio的解析,当然这些信息就少不了可以使用ggtree进行可视化。
三个包的联合,照顾到进化树的解析、输出、操作和可视化的方方面面,恰巧都是我写的。
最搞笑的是MEGA输出newick, nexus和tabular三种格式的树文件,然后自己的输出只能读newick一个。我也是很服气的。
领取专属 10元无门槛券
私享最新 技术干货