wangbingcai2017的个人博客分享 http://blog.sciencenet.cn/u/wangbingcai2017

博文

关联分析TASSEL软件管道命令

已有 5985 次阅读 2017-4-29 16:32 |个人分类:工作|系统分类:科研笔记| Pipeline, Tassel

WIN7下TASSEL5软件在窗口下可以进行多任务处理,但是不能够自动化,TASSEL5软件另外提供了管道命令(pipeline),可以进行批处理来减少重复性劳动。管道命令介绍参照附件Tassel_Pipeline_Tutorial20160330.pdf



典型的MLM(混合线性模型)分析管道命令如下:

       perl run_pipeline.pl

       -fork1 -h genotype.hmp -filterAlign -filterAlignMinFreq 0.05        注:导入基因型数据并过滤

       -fork2 -r trait.txt                                         注:导入表型数据

       -fork3 -r pop_structure.txt -excludeLastTrait               注:导入群体结构数据

       -fork4 -k kinship.txt -combine5 -input1 -input2 -input3 -intersect -combine6 -input5 -input4 -mlm -mlmVarCompEst P3D -mlmCompressionLevel None -export result       注:导入kinship矩阵,合并表型、基因型和群体结构,设定MLM参数

不同的Plugin以灰白背景区分。管道命令的具体用法见附件Tassel5PipelineCLI.pdf



在此基础上结合批处理实现对33个基因型数据的MLM分析:

       for /l %%i in (1,1,33) do

       perl run_pipeline.pl

       -fork1 -h ./hmp/%%i.hmp -filterAlign -filterAlignMinFreq 0.05

       -fork2 -r trait.txt

       -fork3 -r pop_structure -excludeLastTrait

       -fork4 -k kinship -combine5 -input1 -input2 -input3 -intersect -combine6 -input5 -input4 -mlm -mlmVarCompEst P3D -mlmCompressionLevel None -export result%%i

将以上脚本保存为bat格式,放在TASSEL5的安装目录里,其它数据也放在安装目录中。

Tassel_Pipeline_Tutorial20160330.pdf

Tassel5PipelineCLI.pdf




https://blog.sciencenet.cn/blog-3323464-1051904.html


收藏 IP: 220.249.99.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (1 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-5-22 03:40

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部