温馨提示×

温馨提示×

您好,登录后才能下订单哦!

密码登录×
登录注册×
其他方式登录
点击 登录注册 即表示同意《亿速云用户服务条款》

怎么用R语言的vegan包计算物种累计曲线

发布时间:2022-03-21 10:42:50 来源:亿速云 阅读:424 作者:iii 栏目:开发技术

本篇内容主要讲解“怎么用R语言的vegan包计算物种累计曲线”,感兴趣的朋友不妨来看看。本文介绍的方法操作简单快捷,实用性强。下面就让小编来带大家学习“怎么用R语言的vegan包计算物种累计曲线”吧!

R语言 vegan包计算物种累计曲线

vegan 包是进行群落数据分析最常用的R包,其中的 specaccum 函数用来计算物种的累计曲线。

输入数据,otutable:列表示不同的物种,行表示不同的采样点,中间的数字代表物种丰度:

示例代码:
require(vegan)
otu <- read.table("otu_table.txt", header = T,
                       sep = "\t",row.names=1)
all <- specaccum(otu, method = "random")  
png(file="spe.png",h=2000,w=2000,res=300)
plot(all, ci.type = "poly", col = "blue", lwd = 2, ci.lty = 0,
    ci.col = "lightblue", main = "Species Accumulation Curves", xlab = "Number of samples",
    ylab = "Number of  species")
boxplot(all, col = "yellow", add = TRUE, pch = "+")
dev.off()

最终的物种累计曲线中,横坐标是样本个数,纵坐标是发现的物种个数,随着样本个数的增加,发现的物种个数也不断增加;物种累计曲线反应的就是抽样个数对物种多样性的影响;如果曲线末端部分呈现 急剧上升的趋势,表明抽样量不足;增加样本量,还能继续发现新的物种;当曲线末端上升趋势趋于平缓时,则表明采样量足够。

到此,相信大家对“怎么用R语言的vegan包计算物种累计曲线”有了更深的了解,不妨来实际操作一番吧!这里是亿速云网站,更多相关内容可以进入相关频道进行查询,关注我们,继续学习!

向AI问一下细节

免责声明:本站发布的内容(图片、视频和文字)以原创、转载和分享为主,文章观点不代表本网站立场,如果涉及侵权请联系站长邮箱:is@yisu.com进行举报,并提供相关证据,一经查实,将立刻删除涉嫌侵权内容。

AI