微生物组测序数据相对丰富度转化

（一）不同分类水平注释结果相对丰富度转化

在微生物组测序中我们会从taxa_bar.qvz文件中下载到不同水平（level1,level2,level3……）的注释结果，一般我们用这些不同注释结果做微生物组的组成的堆积图，但是这些并不是相对丰富度，因此我们需要对其进行转换。

以level-6属水平注释结果转化为例：


rm(list = ls())
#读入文件
all.data<-read.csv("level-6.csv",row.name=1,header=T)
dim(all.data)
#相对丰富度转化
data.prop<-all.data/rowSums(all.data)
#文件保存
write.csv(data.prop,"genus.csv")

（二）特征表ASVs水平相对丰富度转化

在某些情况我们也需要对特征表中的ASVs的测序结果也转换为相对丰富度，特征表与不同分类水平的注释结果表略有不同，因此在转化中也有略微不同，^…^.

#读入文件
all.data<-read.csv("ASVs特征表.csv",row.name=1,header=T)
#转置
df<-t(all.data)
#相对丰富度转化
data.prop<-df/rowSums(df)
#由于数据量大，转化完直接写出后打开不方便，因此对数据进行转置
ASV<-t(data.prop)
write.csv(ASV,"ASVs特征表_相对丰度.csv")

微生物组测序数据相对丰富度转化相关推荐

QIIME2 分析微生物组测序数据
#启动qiime2 source activate qiime2-2020.2 #创建工作目录并进入工作目录 wd=/mnt/d/Ju/ mkdir -p ${wd} cd ${wd} (一)自建库数 ...
内蒙古农大孙志宏教授证实超深度混合宏基因组测序能够对人类肠道微生物组中的低丰度物种进行基因组和功能表征...
导读人类肠道微生物群中已经发现了大量微生物基因组,但由于目前大多数研究中使用的测序深度相对较浅,在个体水平上了解低丰度物种的作用仍具有挑战.为了提高基因组的组装性能,本研究采用了Illumina H ...
MITRE：利用微生物组时间序列数据推断与宿主状态变化相关的特征
文章目录 MITRE:利用微生物组时间序列数据推断与宿主状态变化相关的特征热心肠日报写在前面主要结果图1. MITRE通过对时间序列的微生物群落的系统发育信息进行学习,同宿主的状态改变联系起来 ...
GB：徐健/王师/黄适合作开发微生物组测序新方法2bRAD-M
2022年1月26日,中国科学院青岛生物能源与过程研究所徐健研究团队联合中国海洋大学王师教授团队在Genome Biology发表了题为 Species-resolved sequencing of ...
积微论坛报告视频+PPT：用微生物组时序数据重现生物膜装配动态过程
点击阅读原文观看视频!后台回复"王金锋"获得视频和讲稿下载链接. 用微生物组时序数据重现生物膜装配动态过程主讲人介绍王金锋,博士.2011年入职中国科学院北京生命科学研究院工作 ...
积微论坛--用微生物组时序数据重现生物膜装配动态过程PPT对应的讲解和提问
生物信息学习的正确姿势 NGS系列文章包括NGS基础.高颜值在线绘图和分析.转录组分析 (Nature重磅综述|关于RNA-seq你想知道的全在这).ChIP-seq分析 (ChIP-seq基本分析流 ...
Nature Methods：Rob Knight发布Striped UniFrac算法轻松分析微生物组大数据
Striped UniFrac微生物组大尺度分析算法 Striped UniFrac: enabling microbiome analysis at unprecedented scale Natu ...
生物工程学报-微生物组测序与分析专刊-邀请函
<生物工程学报> <微生物组测序与分析专刊> 邀请函经过近一年的准备,<微生物组测序与分析专家共识>将在<生物工程学报>发表,并且同时以"微 ...
微生物16S测序数据的正确打开方式
16S rRNA基因测序(也称16S rDNA测序)是最常用的菌群多样性分析的手段.对于新手,如果收到一份不讲"人话"的16S测序分析报告,很快就会被各种生态学术语.各种指数.各种 ...

微生物组测序数据相对丰富度转化

微生物组测序数据相对丰富度转化相关推荐

最新文章

热门文章