16s之稀释曲线绘制

稀释曲线是利用已测得16S rDNA序列中已知的各种OTU的相对比例,来计算抽取n个(n小于测得reads序列总数)reads时出现OTU数量的期望值,然后根据一组n值(一般为一组小于总序列数的等差数列)与其相对应的OTU数量的期望值做出曲线来

通过qiime 的alpha_rarefaction.py (计算稀释丰度)获取各个类型的矩阵表达表

## 1

$ alpha_rarefaction.py -i otu_tax_table.biom -f -p alpha_params.txt -o result -m sample_info.txt

## 2
$ head alpha_params.txt
alpha_diversity:metrics observed_species,shannon,chao1,simpson

## 3
$ head sample_info.txt
#sampleid    group
S1      1
S2      2
S3      3
T1      1
T2      2
T3      3


## 参考format
http://qiime.org/documentation/file_formats.html#id-to-taxonomy-map
http://qiime.org/documentation/qiime_parameters_files.html

## 将上述的得到的结果用做绘制R图的输入文件

setwd("D:/16S/16sRscript")
library("ggplot2")
library("reshape2")
library("tidyr")


inputfile = c("observed_species.txt")
output
  • 1
    点赞
  • 6
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值