【自用、更新】humann2代码整理

本文主要介绍了如何使用humann2进行差异分析,并详细解释了如何通过调整参数进行分组依据设定。此外,还涉及到了如何将结果转换为不同注释格式,如酶、蛋白结构域、GO和反应等,以及如何添加注释文件。
摘要由CSDN通过智能技术生成

1.差异分析

# input是输出的数据,--focal-metadatum --last-metadatum这俩后面的参数是找差异的分组依据

# 激活conda
conda activate humann2
# input是输出的数据,--focal-metadatum --last-metadatum这俩后面的参数是找差异的分组依据
humann2_associate --input pathabundance_acq.pcl  --focal-metadatum Group --focal-type categorical --last-metadatum Group  --output associate_acq.txt
# 生成条形图 这里我前面写了循环,实际可以不写
humann2_barplot --sort sum metadata --input pathabundance.pcl --focal-feature ${id%:*} --focal-metadatum Group --last-metadatum Group --output barplot_${id%:*}.pdf

 2.ko注释

# umann2_regroup_table 修改-g参数,转换为eggNOG(-g uniref90_eggnog),酶(-g uniref90_level4ec),pfam蛋白结构域(-g uniref90_pfam),go(-g uniref90_go)和rxn(-g uniref90_rxn),再合并
#添加相应注释map_eggnog_name.txt.gz,map_go_name.txt.gz……

# humann2转换 数据库那里可以换成别的
humann2_regroup_table -i A.tsv  -g uniref90_ko -o B.tsv
# 解压缩数据库
zcat utility_mapping/map_ko_
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值