为了账号安全,请及时绑定邮箱和手机立即绑定

Qiime1-4.准备参数文件

标签:
Java

由于在我们使用qiime1进行分析的过程中会使用到一些包含了多个脚本的命令,为了更改不同脚本的参数,我们会使用参数文件(Parameters file),包括更改统计方法,绘图时某一组的颜色等等。

例子

针对一些步骤,我们不想使用默认参数,所以我们需要生成一个参数文件修改参数。下面的代码会生成一个16S_pickotu_param.txt的参数文件。

## This will output a text file to whichever your directory you're currently working in## so be sure you know where you are saving the file!# Enable reverse strand matching for uclust so that it double checks both sequence possibilities.echo "pick_otus:enable_rev_strand_match True" >> 16S_pickotu_param.txt# Change to program to RDP for assigning taxnomy of the OTU's instead of uclust echo "assign_taxonomy:assignment_method rdp" >> 16S_pickotu_param.txt# Change to confidence to 0.5 in the RDP classiferecho "assign_taxonomy:confidence 0.5" >> 16S_pickotu_param.txt

最后的txt文件如下:

webp

具体的格式:
涉及的脚本名称+:+参数全称+空格+具体参数
assign_taxonomy:assignment_method rdp为例
可以在qiime1的官网上搜索到assign_taxonomy.py的具体参数(http://qiime.org/scripts/assign_taxonomy.html

webp

assign_taxonomy.py

-m, --assignment_method
Taxon assignment method, must be one of rdp, blast, rtax, mothur, uclust, sortmerna [default: uclust]

发现其中有一个参数时设置分类器,默认的事uclust,我们想使用rdp分类器。那就可以在参数文件中输入这个脚本的名称,然后加上冒号,接着输入这个参数的全称+空格+rdp。



作者:jlyq617
链接:https://www.jianshu.com/p/d839d68c5867


点击查看更多内容
TA 点赞

若觉得本文不错,就分享一下吧!

评论

作者其他优质文章

正在加载中
  • 推荐
  • 评论
  • 收藏
  • 共同学习,写下你的评论
感谢您的支持,我会继续努力的~
扫码打赏,你说多少就多少
赞赏金额会直接到老师账户
支付方式
打开微信扫一扫,即可进行扫码打赏哦
今天注册有机会得

100积分直接送

付费专栏免费学

大额优惠券免费领

立即参与 放弃机会
意见反馈 帮助中心 APP下载
官方微信

举报

0/150
提交
取消