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

首先看下官方示例:

library(vegan)

data(BCI)

sp1

plot(sp1, ci.type="poly", col="blue", lwd=2, ci.lty=0, ci.col="lightblue")

boxplot(sp1, col="yellow", add=TRUE, pch="+")

出来的结果图如下:

那么这幅图表明了什么含义呢?

首先看下输入数据

> head(BCI[, 1:3])

Abarema.macradenia Vachellia.melanoceras Acalypha.diversifolia

1 0 0 0

2 0 0 0

3 0 0 0

4 0 0 0

5 0 0 0

6 0 0 0

我们简单的看一下BCI这个数据,它的每一行代表了一个样本,不同样本采样的地点不同,每一列是1个物种的丰度

最终的物种累计曲线中,横坐标是样本个数,纵坐标是发现的物种个数,随着样本个数的增加,发现的物种个数也不断增加;

其实,物种累计曲线反应的就是抽样个数对物种多样性的影响;可以看到,当抽样个数较少时,发现的物种并不全面,并不能表征整个群落结构,随着抽样个数的上升,发现的物种数越来越多,也更能表征这个群落结构;

在实际分析中,什么样的情况表明我们的采样量足够了呢,主要看曲线的末端,如果曲线末端部分还呈现 急剧上升的趋势,表明抽样量不足;增加样本量,还能继续发现新的物种;当曲线末端上升趋势趋于平缓时,则表明采样量足够,

vegan稀释曲线 基因丰度_R语言 vegan包计算物种累计曲线相关推荐

  1. vegan稀释曲线 基因丰度_蒙古沙冬青及其伴生植物AM真菌物种多样性

    蒙古沙冬青(Ammopiptanthus mongolicus)隶属豆科沙冬青属, 是西北荒漠生境中唯一常绿阔叶灌木, 耐干旱.抗逆性强, 在保持水土和防治荒漠化方面作用显著[.与蒙古沙冬青相伴而生的 ...

  2. vegan稀释曲线 基因丰度_基于OTU的稀释曲线(Rarefaction curves) + ggplot2

    1. 简介 稀释曲线(Rarefaction curves)是从样品中随机抽取一定测序量的数据(序列条数),统计它们所对应的OTUs种类(代表物种),并以抽取的测序数据量与对应的代表OTUs来构建曲线 ...

  3. 科研文献|粪便污染可以解释人为影响环境中抗生素耐药基因丰度

    ​研究介绍 对耐抗生素细菌对人类健康构成的威胁的日益关注,也使注意力转向了这个问题的环境层面.处理后的污水排放导致抗生素耐药菌.耐药基因和抗生素残留释放到环境中.然而,尚不清楚污水和受污水影响的环境中 ...

  4. 宏基因组实战6. 不比对快速估计基因丰度Salmon

    前情提要 如果您在学习本教程中存在困难,可能因为缺少背景知识,建议先阅读本系统前期文章 宏基因组分析理论教程 微生物组入门圣经+宏基因组分析实操课程 1背景知识-Shell入门与本地blast实战 2 ...

  5. 在线作图|微生物多样性分析——物种累计曲线

    物种累计曲线 物种累计曲线( species accumulation curves)用于衡量和预测群落中物种丰富度随样本量扩大而增加的幅度,在生物多样性和群落调查中,被广泛用于判断样本量是否充分的并 ...

  6. c语言程序计算华氏度,C语言:一个计算摄氏度华氏度和绝对温度三种数转换的程序...

    C语言:一个计算摄氏度华氏度和绝对温度三种数转换的程序 答案:2  信息版本:手机版 解决时间 2020-03-31 21:49 已解决 2020-03-31 06:26 初学C语言 有道题大概是 要 ...

  7. R语言raster包计算多个栅格图像平均值、标准差的方法

      本文介绍基于R语言中的raster包,批量读取多张栅格图像,对多个栅格图像计算平均值.标准差,并将所得新的栅格结果图像保存的方法.   在文章R语言raster包批量读取单一或大量栅格图像(htt ...

  8. lda主题模型困惑度_r语言lda主题模型代码 lda主题模型 案例分析

    r语言做灰色预测模型代码为什么会出错? 我正在MATLAB下编写一个Grey system的函数库, 我想其他的网友也做过类似的工作吧. 灰色模型的MATLAB求解代码应该很容易找到吧.如果不想用MA ...

  9. 语言模拟蒲丰问题_R语言小数定律的保险业应用:泊松分布模拟索赔次数

    原文链接: 拓端数据科技 / Welcome to tecdat​tecdat.cn 在保险业中,由于分散投资,通常会在合法的大型投资组合中提及大数定律.在一定时期内,损失"可预测" ...

最新文章

  1. AI智能手机会是什么样?至少有这十个功能
  2. Linux下的kill函数的用法
  3. glance was not installed properly
  4. 44行代码AC_卡片换位(DFS变形题 视频讲解 )
  5. java 自动装拆箱
  6. 【广州】Web前端工程师
  7. php文件操作(上传文件)1
  8. chrome浏览器开发者模式header请求参数“Request Payload“找不到的问题
  9. android sqlcipher 加密,Android 数据库加密 SQLCipher使用方法
  10. 一名亚马逊软件工程师的辛酸苦楚
  11. 3D数学之四元组应用及实现
  12. Tencent笔试题收集
  13. 小规模纳税人和一般纳税人的区别
  14. 通达信标记符号_通达信49个图标,高手指教,通达信软件里的股票标记符号,能增加吗...
  15. 新手网站建设优化,这些网站为你提供数之不尽的免费素材!(3)
  16. java table 增加行_使用POI给word中的表格增加行
  17. SVN报错:can't open file db/txn-current-lock:permission denied 解决方法
  18. 作为3-5年的iOS开发者:你为什么迟迟进不去大厂呢?
  19. 这样软件开发人才级别的划分你同意吗
  20. 一篇文章,带你走进Java

热门文章

  1. java set 接口_java笔记四:Set接口
  2. 初级篇第三期:初识UI
  3. linux7ip路由,Centos7源地址访问路由(双IP双网关配置)
  4. pytorch 入门学习使用逻辑斯蒂做二分类-6
  5. c语言,求素数个数,关于求素数个数的话题
  6. html云文件系统,一种HTML5云文件系统
  7. springboot通过各种不同类型参数获取容器中的bean工具类
  8. PHP中利用Ffmpeg获得flv视频缩略图和播放时间
  9. Spring中bean属性注入方式总结:构造方法、setter注入、p命令空间注入、SpEL注入、集合注入
  10. linux 动态内存分配,具体来说,fork()如何处理Linux中malloc()动态分配的内存?