微生物群落分析思考:Simpson(D)指数的多种形式及其与Niche width (B)的关系

导师在帮忙修改论文时,给出了关于微生物群落Simpson指数及其与Niche width关系的批注,大意是说Simpson指数有多种形式,要合理选取,并且Simpson指数与Niche width给出的信息相重合,同时出现的话略显冗余,因此略作学习小计。

总的来说,Simpson (D)有1-D、1/D、-log(1-D)几种形式。

一般利用R语言的vegan包计算出来的Simpson是我们常说的D,计算方法也就是1 - eq?%5Csum pi^2,简单理解就是1-OTU相对丰度的和

#使用vegan包计算多样性指数
Simpson <- diversity(df, index = "simpson", MARGIN = 2, base =  exp(1))

但是这个值与群落多样性呈现出一个相反的关系,因此,后来出现了1-D和1/D,可以根据自己的研究选取,在方法部分表明即可。

另外,还有一种呈现形式是-log(1-D),能够让Simpson和Shannon指数具有相同的量表,简单理解就是两者单位统一了。可以参考文献Haegeman et al. 2013Robust estimation of microbial diversity in theory and in practice | The ISME Journal (nature.com)

然后就是Simpson与niche width的关系,利用'spaa'计算Levins niche width (B)的公式是1 / eq?%5Csum pi^2,其实也就是B=1/(1-D),给出的信息确实是相重叠的。

欢迎交流指正。

 

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值