picrust2功能预测-从qiime2安装到数据分析

更新时间:2023-05-08 20:32:39 阅读: 评论:0

如果使⽤以前的版本不会出现这样的问题。
有个⼩插曲:不知道怎么解决这个问题的时候,⽆论输⼊什么数据都显⽰出这个结果,⼼态差点就崩了,然后花了点钱请淘宝上做⽣信分析的⼈,⽤我的数据帮我跑了⼀遍,可以出结果。最后观察版本之间有什么差别的时候就发现最新版本优化了,然后淘宝上那个⽣信技术厉害的客服⽤最新版本也出不了结果,他就去官⽹发了⼀个issue(就是上⾯图⽚⾥的),开发者很快就给回复了,解释了原因。
问题3:使⽤varch跑出的结果使⽤out.fa和out_table.biom两个⽂件可以直接跑出预测结果,但使⽤qiime跑出的结果中两个⽂件夹的out名称不⼀致,要删掉后⾯的其他字符可采⽤以下⽅法。
修改fastq⽂件OTU 名字`
awk '{print $1 }' FASTA_IN > FASTA_OUT
其中FASTA_IN和FASTA_OUT是输⼊和输出FASTA⽂件的名称。
$2:⼀⾏⼀⾏的读取指定的⽂件, 以空格作为分隔符,打印第⼆个字段
⽐如有这样⼀个⽂件
a1 b1 c1 d1
a2 b2 c2 d2
执⾏的结果是,输出:
b1
b2
问题4:会出现以下报错,⼤概率是电脑内存不够,picrust2官⽹有⼀个针对这个错误的issue。
picrust2_pipeline.py –s out.fa -i out_table.biom \
-o picrust2_out_pipeline -p 4
-p 4是使⽤4线程的意思
输出的结果包括EC,KO,PATHWAY的丰度信息。
但是但是结果⽂件中的结果⽐较简洁,⽐如只输⼊EC1.1.1,⽽不标注对应的酶名字,这时候需要加⼯加⼯。根据需要酌情添加描述:add_descriptions.py -i EC_metagenome_out/pred_metagenome_unstrat.tsv \
-m EC  -o EC_metagenome_out/pred_metagenome_unstrat_descrip.tsv
add_descriptions.py -i KO_metagenome_out/pred_metagenome_unstrat.tsv \
-m KO  -o KO_metagenome_out/pred_metagenome_unstrat_descrip.tsv
add_descriptions.py -i COG_metagenome_out/pred_metagenome_unstrat.tsv \
-m COG  -o COG_metagenome_out/pred_metagenome_unstrat_descrip.tsv
add_descriptions.py -i PFAM_metagenome_out/pred_metagenome_unstrat.tsv \
-m PFAM  -o PFAM_metagenome_out/pred_metagenome_unstrat_descrip.tsv
add_descriptions.py -i TIGRFAM_metagenome_out/pred_metagenome_unstrat.tsv\
-m TIGRFAM  -o TIGRFAM_metagenome_out/pred_metagenome_unstrat_descrip.tsv
add_descriptions.py -i pathways_out/path_abun_unstrat.tsv
-m METACYC -o pathways_out/path_abun_unstrat_descrip.tsv
结果输出中没有KEGG结果,可通过以下命令⽣成KEGG pathway 丰度表和注释结果
pathway_pipeline.py -i KO_metagenome_out/pred_metagenome_unstrat.tsv \
-o KEGG_pathways_out --no_regroup --map \
/home/qiime2/miniconda/envs/picrust/lib/python3.6/site-package/
picrust2/default_files/pathway_mapfiles/KEGG_pathways_to_KO.tsv
这⾥KEGG_pathways_to_KO.tsv需要你在你的虚拟机上搜索后输⼊完整并且正确的⽂件路径。
# 添加功能描述
add_descriptions.py -i KEGG_pathways_out/path_abun_ \
--custom_map_table /home/qiime2/miniconda/envs/picrust/lib/python3.6/site-
package/picrust2/default_files/description_mapfiles/KEGG_pathways_ \
-o KEGG_pathways_out/path_abun_unstrat_
这⾥/KEGG_pathways_需要你在你的虚拟机上搜索后输⼊完整并且正确的⽂件路径。

本文发布于:2023-05-08 20:32:39,感谢您对本站的认可!

本文链接:https://www.wtabcd.cn/fanwen/fan/78/559020.html

版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系,我们将在24小时内删除。

标签:结果   版本   问题   输出   时候   打印   出现   解释
相关文章
留言与评论(共有 0 条评论)
   
验证码:
推荐文章
排行榜
Copyright ©2019-2022 Comsenz Inc.Powered by © 专利检索| 网站地图