step1  原材料:

1. OTU表

2. OTU.fasta

step2  利用qiime转换为GreenGene ID替换原来的OTU_num

#qiime 的 OTU 序列注释过程
assign_taxonomy.py -i otu.fasta -r gg_13_5_otus/rep_set/97_otus.fasta -t gg_13_5_otus/taxonomy/97_otu_taxonomy.txt --uclust_max_accepts 1 -o gg97 

#根据具体的注释细节,对应 greengene id(这里用了一个手写的 R 脚本实现转换,在网盘附件中)
grep 'H' gg97/otu_tax_assignments.log > gg97/otu_tax_assignments2.log
Rscript replace_id.r gg97/otu_tax_assignments2.log otu_table.txt otu_table2.txt

#转为 biom 格式
sed -i '1i\# Constructed from biom file' otu_table2.txt
biom convert -i otu_table2.txt -o otu_table2.biom --table-type="OTU table" --to-json
 

step3  可选校正16S拷贝数

#校正 16

对16S数据使用PICRUSt预测输出KEGG和COG表相关推荐

  1. 16S的细菌群落功能预测工具PICRUSt2学习

    2018年推出了全新版本的PICRUSt,即PICRUSt2(https://github.com/picrust/picrust2) PICRUSt2 (Phylogenetic Investiga ...

  2. 基于16S的细菌群落功能预测工具Tax4Fun2

    基于16S的细菌群落功能预测工具Tax4Fun2 回想第一次接触微生物组分析时是2017年,那会儿主流的两个基于16S的细菌群落功能预测工具就是PICRUSt和Tax4Fun,均可以获得类似于宏基因组 ...

  3. ​易生信-宏基因组积微学术论坛:基于大数据整合准确预测土壤的枯萎病发生...

    博彩众家之长,积微成就突破.为促进我国宏基因组研究领域的学术交流和技术分享,推动微生物组领域的发展,"宏基因组"公众号联合国内外优秀人才组织"易生信-宏基因组 积微学术论 ...

  4. 混合时空图卷积网络:利用导航数据改进交通预测效果 | KDD 2020

    作者 | 高德机器学习团队 出品 | AI科技大本营(ID:rgznai100) 引言 时空预测(spatio-temporal forecasting)在天气预报.运输规划等领域有着重要的应用价值. ...

  5. DL之LSTM之MvP:基于TF利用LSTM基于DIY时间训练csv文件数据预测后100个数据(多值预测)状态

    DL之LSTM之MvP:基于TF利用LSTM基于DIY时间训练csv文件数据预测后100个数据(多值预测)状态 目录 数据集csv文件内容 输出结果 设计思路 训练记录全过程 数据集csv文件内容 输 ...

  6. bp神经网络预测未来五年数据_基于小波神经网络的数据中心KPI预测

    随着软件和微服务的发展,智能运维越来越受到人们的重视.在大量的运维数据里,最不可忽视的就是各种关键性能指标数据(Key Performance Indicators,KPI),它们在数学上都可以被表达 ...

  7. (四)比特币时间序列数据的AI预测

    目录 介绍 创建LSTM NN回归器 训练模型 测试模型 下一步 下载异常检测 - 17.9 MB 介绍 本系列文章将指导您完成使用AI开发功能齐全的时间序列预测器和异常检测器应用程序所需的步骤.我们 ...

  8. 支持向量机识别数字集(数据采集+模型训练+预测输出)

    SVM简单实现图像中的数字分类 支持向量机 使用SVM识别一个物体 必须了解的两个参数 SVM的识别流程 代码实现 图像采集 保存采集的图像 数据集处理 训练数据 预测数字 支持向量机 支持向量机(S ...

  9. 如何用神经网络预测数据,神经网络做预测的原理

    神经网络算法原理? 神经网络预测学习样本中的驾驶行为特征. 如图显示了某个驾驶场景的行驶路径深度学习训练,通过神经网络可以学习驾驶人的行为,并根据当前获取的环境信息决策行驶轨迹,进而可以控制车辆的转向 ...

最新文章

  1. 2天训练出15亿参数大模型,国产开源项目力克英伟达Megatron-LM,来自LAMB作者团队...
  2. 砂.随笔.二十五.如果你是氧气
  3. PHP之旅4 php 超全局变量
  4. Linux 下的推迟执行
  5. python open函数 创建变量文件_python的open函数怎么用
  6. git修改commit注释_【Slog】Git之多人同feature的同分支开发
  7. Google、Facebook等不可错过的 AI 开源工具,你想要的都在这里!
  8. 剑指offer(C++)-JZ76:删除链表中重复的结点(数据结构-链表)
  9. 服务器开机信号,机柜服务器BMC开关机控制系统及方法
  10. [转]coolfire黑客入门教程系列之(六)
  11. Visio如何画箭头
  12. 三星R463无线网卡驱动,声卡驱动,FN快捷键驱动,Easy_Display_Manager
  13. 计算机点击应用时总是跳到顶端,电脑总是跳出当前应用程序
  14. dll修复精灵无法链接服务器,教你一键dll修复精灵怎么使用
  15. Mysql给一个大表加一列_MySQL 大表添加一列的实现
  16. 立方体三维之立方体平行投影
  17. AcWing 1127 香甜的黄油
  18. Linux clamav杀毒
  19. ​无人机视频直播推流方案
  20. SDH系统实现PCM信号传输以及时分复用的simulink仿真

热门文章

  1. qtableview选中第一行时表头会变色_亲民的“网红潜水表”精工鱼罐头 Prospex,你会选择吗?...
  2. 更深的技术探索,更多的场景实践M etaCon元宇宙技术大会圆满召开
  3. 三光(可见光、红外光、激光)云台产品调研
  4. linux交换分区的命令,LInux下的交换分区以及相关查看命令
  5. Android中ExpandableListView控件的用法详解
  6. WebRTC 实现P2P音视频通话——搭建信令服务器
  7. 李宁卖咖啡,意欲何为?
  8. mysql怎么定位错误信息_如何快速定位MySQL 的错误日志(Error Log)?
  9. 华脉智联校园可视化智慧安全感知处置系统解决方案
  10. 【数据结构+算法】浙传OJ Contest 2290:13信息1 Java 6 问题 F: Hexagon Game