本篇内容主要讲解“怎么用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包计算物种累计曲线”有了更深的了解，不妨来实际操作一番吧！这里是亿速云网站，更多相关内容可以进入相关频道进行查询，关注我们，继续学习！
免责声明：本站发布的内容（图片、视频和文字）以原创、转载和分享为主，文章观点不代表本网站立场，如果涉及侵权请联系站长邮箱：is@yisu.com进行举报，并提供相关证据，一经查实，将立刻删除涉嫌侵权内容。
