|
|
最近在做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 |
» 猜你喜欢
售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急
已经有3人回复
售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急
已经有4人回复
售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急
已经有3人回复
售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急
已经有3人回复
售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急
已经有3人回复
售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急
已经有3人回复
售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急
已经有3人回复
基金正文30页指的是报告正文还是整个申请书
已经有4人回复
今年春晚有几个节目很不错,点赞!
已经有6人回复
球磨粉体时遇到了大的问题,请指教!
已经有15人回复
|