|
|
最近在做16S高通量结果分析,主要是利用qiime来做,目前遇到瓶颈了,不知道为什么我在做beta diversity里面的make_2d_plots.py这步中,出来的结果没有标签lable(即哪种颜色哪种形状代表什么样品),不知道该怎么调试,,有哪位大神可以告知一二。。。运行的命令主要是以下几个:
single_rarefaction.py -i otu_table_taxonomy.biom -o div_beta//otu_table_taxonomy_even11723.biom -d 11723
make_prefs_file.py -m typical_Map.txt -o div_beta//prefs_out.txt --mapping_headers_to_use Description
beta_diversity.py -i div_beta//otu_table_taxonomy_even11723.biom -o div_beta/ --metrics weighted_unifrac -t rep_phylo.tre
mv div_beta//weighted_unifrac_otu_table_even3703.txt div_beta//weighted_unifrac_dm.txt
principal_coordinates.py -i div_beta//weighted_unifrac_dm.txt -o div_beta//weighted_unifrac_pc.txt
make_2d_plots.py -i div_beta//weighted_unifrac_pc.txt -m typical_Map.txt -p div_beta//prefs_out.txt -o 2d_plots
我自己有时想想,是不是prefs_out.txt的文档出错了,但对于我这个初学者而言,有点难,,跪求大神告知
![qiime求助]()
w74h4833202_1479816598_854.jpg |
» 猜你喜欢
极限拉扯,如何抉择
已经有5人回复
课题组会帮你通过考核么?
已经有6人回复
OpenAI又攻克数学难题了,离AI做有机合成还有多远?
已经有7人回复
售一区SCI文章T0P,我:8O.551.O54,科目全,可十急
已经有9人回复
售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急
已经有8人回复
售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急
已经有9人回复
售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急
已经有7人回复
售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急
已经有10人回复
售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急
已经有9人回复
|