vegan稀释曲线 基因丰度_R语言 vegan包计算物种累计曲线
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包计算物种累计曲线相关推荐
- vegan稀释曲线 基因丰度_蒙古沙冬青及其伴生植物AM真菌物种多样性
蒙古沙冬青(Ammopiptanthus mongolicus)隶属豆科沙冬青属, 是西北荒漠生境中唯一常绿阔叶灌木, 耐干旱.抗逆性强, 在保持水土和防治荒漠化方面作用显著[.与蒙古沙冬青相伴而生的 ...
- vegan稀释曲线 基因丰度_基于OTU的稀释曲线(Rarefaction curves) + ggplot2
1. 简介 稀释曲线(Rarefaction curves)是从样品中随机抽取一定测序量的数据(序列条数),统计它们所对应的OTUs种类(代表物种),并以抽取的测序数据量与对应的代表OTUs来构建曲线 ...
- 科研文献|粪便污染可以解释人为影响环境中抗生素耐药基因丰度
研究介绍 对耐抗生素细菌对人类健康构成的威胁的日益关注,也使注意力转向了这个问题的环境层面.处理后的污水排放导致抗生素耐药菌.耐药基因和抗生素残留释放到环境中.然而,尚不清楚污水和受污水影响的环境中 ...
- 宏基因组实战6. 不比对快速估计基因丰度Salmon
前情提要 如果您在学习本教程中存在困难,可能因为缺少背景知识,建议先阅读本系统前期文章 宏基因组分析理论教程 微生物组入门圣经+宏基因组分析实操课程 1背景知识-Shell入门与本地blast实战 2 ...
- 在线作图|微生物多样性分析——物种累计曲线
物种累计曲线 物种累计曲线( species accumulation curves)用于衡量和预测群落中物种丰富度随样本量扩大而增加的幅度,在生物多样性和群落调查中,被广泛用于判断样本量是否充分的并 ...
- c语言程序计算华氏度,C语言:一个计算摄氏度华氏度和绝对温度三种数转换的程序...
C语言:一个计算摄氏度华氏度和绝对温度三种数转换的程序 答案:2 信息版本:手机版 解决时间 2020-03-31 21:49 已解决 2020-03-31 06:26 初学C语言 有道题大概是 要 ...
- R语言raster包计算多个栅格图像平均值、标准差的方法
本文介绍基于R语言中的raster包,批量读取多张栅格图像,对多个栅格图像计算平均值.标准差,并将所得新的栅格结果图像保存的方法. 在文章R语言raster包批量读取单一或大量栅格图像(htt ...
- lda主题模型困惑度_r语言lda主题模型代码 lda主题模型 案例分析
r语言做灰色预测模型代码为什么会出错? 我正在MATLAB下编写一个Grey system的函数库, 我想其他的网友也做过类似的工作吧. 灰色模型的MATLAB求解代码应该很容易找到吧.如果不想用MA ...
- 语言模拟蒲丰问题_R语言小数定律的保险业应用:泊松分布模拟索赔次数
原文链接: 拓端数据科技 / Welcome to tecdattecdat.cn 在保险业中,由于分散投资,通常会在合法的大型投资组合中提及大数定律.在一定时期内,损失"可预测" ...
最新文章
- AI智能手机会是什么样?至少有这十个功能
- Linux下的kill函数的用法
- glance was not installed properly
- 44行代码AC_卡片换位(DFS变形题 视频讲解 )
- java 自动装拆箱
- 【广州】Web前端工程师
- php文件操作(上传文件)1
- chrome浏览器开发者模式header请求参数“Request Payload“找不到的问题
- android sqlcipher 加密,Android 数据库加密 SQLCipher使用方法
- 一名亚马逊软件工程师的辛酸苦楚
- 3D数学之四元组应用及实现
- Tencent笔试题收集
- 小规模纳税人和一般纳税人的区别
- 通达信标记符号_通达信49个图标,高手指教,通达信软件里的股票标记符号,能增加吗...
- 新手网站建设优化,这些网站为你提供数之不尽的免费素材!(3)
- java table 增加行_使用POI给word中的表格增加行
- SVN报错:can't open file db/txn-current-lock:permission denied 解决方法
- 作为3-5年的iOS开发者:你为什么迟迟进不去大厂呢?
- 这样软件开发人才级别的划分你同意吗
- 一篇文章,带你走进Java
热门文章
- java set 接口_java笔记四:Set接口
- 初级篇第三期:初识UI
- linux7ip路由,Centos7源地址访问路由(双IP双网关配置)
- pytorch 入门学习使用逻辑斯蒂做二分类-6
- c语言,求素数个数,关于求素数个数的话题
- html云文件系统,一种HTML5云文件系统
- springboot通过各种不同类型参数获取容器中的bean工具类
- PHP中利用Ffmpeg获得flv视频缩略图和播放时间
- Spring中bean属性注入方式总结:构造方法、setter注入、p命令空间注入、SpEL注入、集合注入
- linux 动态内存分配,具体来说,fork()如何处理Linux中malloc()动态分配的内存?