之前发表了一篇推文给出了如何拆分OTU表格从而得到样本中丰富和稀有物种的方法,前两天在群里有人提出了问题,说是拆分后的表格加回去之后和原始的表格并不一致,这才发现之前的代码存在一定的漏洞。
之前的代码在计算条件稀有物种和条件丰富物种的时候忽略的“在部分样本”中这一条件,从而导致条件稀有物种中含有稀有物种、条件丰富物种中含有丰富物种。
本次更新的代码中对这一问题进行了处理,目前的代码应该是没有问题了,我试验了一下,拆分后得到的6个子表格合并之后与原始的表格是一致的了。
物种丰度划分
我第一次看到物种丰度划分的方法是在北京大学温东辉老师于2016年发表在《FEMS Microbiology Ecology》上的文章,文章是针对杭州湾沉积物的一项研究,提出了按照物种丰度对物种进行划分,从而找到影响微生物群落的关键物种的方法,感兴趣的朋友可以点击文末的“阅读原文”获取该文献。
物种划分方法
在这篇文章中将物种按照其在样本中的检出丰度分为了5类。