继续浏览精彩内容
慕课网APP
程序员的梦工厂
打开
继续
感谢您的支持,我会继续努力的
赞赏金额会直接到老师账户
将二维码发送给自己后长按识别
微信支付
支付宝支付

Qiime1-4.准备参数文件

慕虎7371278
关注TA
已关注
手记 1259
粉丝 203
获赞 873

由于在我们使用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


打开App,阅读手记
0人推荐
发表评论
随时随地看视频慕课网APP