![](https://img-blog.csdnimg.cn/7f88ecb9a9fd411cbcd3ac9541cb3e3b.png?x-oss-process=image/resize,m_fixed,h_224,w_224)
微生物群落构建分析
文章平均质量分 87
运用并行计算加快生态零模型的模拟,节约统计分析时间。专栏以R语言为基础,实现Raup-Crick distance 、beta-NTI 和 null deviation of beta diversity的并行计算。
道阻且长1994
溯洄从之,道阻且长。溯游从之,宛在水中央。
展开
-
利用GPU并行计算beta-NTI,大幅减少群落构建计算时间
18个样本,抽平到8500条序列,4344个OTUs,计算beta-NTI共花费时间如下:> GPU(GTX1050):1分20秒> iCAMP包 的`bNTIn.p()` 函数4核并行:约16分钟如果更好的显卡,更大的数据量,节约的时间应该更加可观。原创 2023-05-15 18:06:51 · 2186 阅读 · 4 评论 -
检测OTU序列遗传发育信号的R实现
群落格局及其背后的构建机制是群落生态学研究的重点内容,Stegen et al. 在2011年使用了beta-NTI这一指标来推断群落构建,此后得到了微生物生态学家的广泛使用。使用beta-NTI推断群落构建有一个前提条,就是该基因必须在较小的遗传距离上具有遗传发育信号(Q1:为什么这是必须的呢?),即在较小的遗传距离下,遗传发育距离越大,其生态位距离越大(Q2:为什么是较小的遗传距离下有遗传距离发育信号即可,不需要整个的遗传距离呢?)。本文将介绍三种方法检测近距离下的遗传发育信号,并用R语言实现。原创 2022-04-01 20:17:36 · 1277 阅读 · 0 评论 -
冗余分析(RDA)中若包含生物学重复会怎样?
一般微生物测序实验会包含三个生物学重复,然后获得有重复的OTU table和环境因子的数据。许多文献在对OTUs进行冗余分析时都包含了重复,包含重复与不包含重复会有何种不同?是否会影响我们对分析结果的解读呢?下面我们以基于距离的冗余分析(db-RDA)作为例子,尝试分析一下,看看两者的区别。原创 2022-03-31 16:25:34 · 986 阅读 · 0 评论 -
Sloan中性群落模型(NCM)推断群落构建原理及其R实现
许多SCI文章都使用了一种叫NCM的方法来对微生物群落构建(中性理论 or 生态位理论?)进行预测,本文将对NCM的实现原理及其R实现做介绍。在介绍其原理的过程中会有一些数学公式,不要对这些数学公式产生抵触情绪,我会解释其每一项的生态学含义。与原理相比,其R实现则相对简单,阅读完后可将原理与R实现进行对比学习。原创 2022-03-23 15:48:16 · 10778 阅读 · 20 评论 -
R语言并行计算RC~bray-curtis~距离
;群落构建分析是微生物生态学分析的重要组成部分,成为目前文章发表的热点技术。之前我们介绍了计算beta-NTI(beta nearest taxon index)来进行群落构建分析。|beta-NTI| >2说明决定性过程主导,其中beta-NTI >2说明OTU的遗传距离发散,为生物交互作用主导,beta-NTI < -2则说明OUT的遗传距离收敛为环境选择主导。|beta-NTI|原创 2022-02-07 21:07:55 · 7807 阅读 · 47 评论 -
R语言并行计算 deviation of null beta diversity(beta多样性零偏差)
群落构建分析是微生物生态学分析的重要组成部分,成为目前文章发表的热点技术。之前我们介绍了计算beta-NTI(beta nearest taxon index)来进行群落构建分析。|beta-NTI| >2说明决定性过程主导,其中beta-NTI >2说明OTU的遗传距离发散,为生物交互作用主导,beta-NTI < -2则说明OUT的遗传距离收敛为环境选择主导。|beta-NTI|原创 2022-02-07 11:41:45 · 3620 阅读 · 0 评论 -
R语言并行计算beta-NTI值
R语言并行构建零模型,加快beta-NTI值的计算。原创 2021-07-06 11:33:40 · 14677 阅读 · 97 评论