最近好几个人遇到了同样的问题,就是在将OTU/ASV和系统发育树对齐的时候,报错:
意思是OTU/ASV名字和系统发育树的节点全都不匹配,导致树上的节点全都被去掉了,树就变成了NULL。
这种问题的出现基本上都是进化树导致的。需要仔细查看进化树文件是否正确。在R中读入进化树:
查看一下树,正确的格式应该是这样的:
但是,如果你的数据是某公司分析的,他们返回的树文件可能是这样的:
OTU/ASV名字多了单引号。而单引号也会被当成树节点的一部分,因此就和OTU/ASV不一致,导致全不匹配。。。
解决的办法也很简单,把树节点单引号替换掉即可:
完整的代码如下: