用户工具

站点工具


card双圈图

简介

meta分析流程中,会遇到客户需求按分组绘制CARD双圈图。

此处整理补充了meta CARD分组双圈图分析思路。

功能

CARD分组双圈图分析。

数据准备

1. unigene丰度文件准备:

双圈图的作图文件为Unigenes.absolute.total.tax.xls和Unigenes.protein.rgi.txt,整体性思路是按照分组将2个表格处理成分组表格, 首先处理Unigenes.absolute.total.tax.xls将同一个分组的样本列提取出来,然后使用脚本处理该文件的行,脚本如下:

perl /TJPROJ5/META_ASS/16s/chenjiawei/script/META/CARD/fenzu-circos/shaixuan.pl Unigenes.absolute.total.tax.xls

此处挑选出分组unigene id形成一个list,命名为shaixuan.list。

2. rgi预测文件准备:

根据挑选出的unigene id list,即shaixuan.list对rgi文件进行筛选,选出该分组的rgi预测结果。

perl shaixuan2.pl Unigenes.protein.rgi.txt shaixuan.list Unigenes.protein.rgi-new.txt

数据分析

1. 生成分组双圈图各层级绘图文件:

perl /PUBLIC/software/MICRO/share/MetaGenome_pipeline/MetaGenome_pipeline_V5.1/lib/05.Function/lib/CARD/lib/mechanism.tax.pl --mech /PUBLIC/software/MICRO/share/MetaGenome_pipeline/MetaGenome_pipeline_V5.1/lib/05.Function/lib/CARD/lib/../database/ARO_Mechnism.list --anno Unigenes.protein.rgi.txt --table Unigenes.absolute.total.tax.xls  --outdir Mech_tax --prefix meta.tax

2. 绘制双圈图:

perl /PUBLIC/software/MICRO/share/MetaGenome_pipeline/MetaGenome_pipeline_V5.1/lib/05.Function/lib/CARD/lib/circos_overview4.0.pl  --table Mech_tax/meta.tax.p.xls

完整脚本路径:

/TJPROJ5/META_ASS/16s/chenjiawei/script/META/CARD/fenzu-circos
card双圈图.txt · 最后更改: 2023/01/05 03:10 由 chenjiawei