Linux系统下TASSEL5.0的安装及使用

· · 753 次点击 · · 开始浏览    
这是一个创建于 的文章,其中的信息可能已经有所发展或是发生改变。

之前做SSR标记和简化测序都是用的桌面化的TASSEL进行运算,对于做关联分析的初学者非常简单,十分好入门。但是关于命令行的操作就要麻烦一些,但是运算起来不像GAPIT那样消耗大量的内存,下面做一些简单的示例。


安装

Java JDK 8.0的安装

#下载
wget http://download.oracle.com/otn-pub/java/jdk/8u112-b15/jdk-8u112-linux-x64.rpm
#安装
rpm -ivhc jdk-8u112-linux-x64.rpm
#测试安装
java -vsrsion
#查看java版本
[root@lilibei java]# java -version
java version "1.8.0_112"
Java(TM) SE Runtime Environment (build 1.8.0_112-b15)
Java HotSpot(TM) 64-Bit Server VM (build 25.112-b15, mixed mode)
  • 更改环境变量
vim + /etc/profile

按照如下添加

JAVA_HOME=/usr/java/jdk1.8.0
JRE_HOME=/usr/java/jdk1.8.0/jre
PATH=$PATH:$JAVA_HOME/bin:$JRE_HOME/bin
CLASSPATH=.:$JAVA_HOME/lib/dt.jar:$JAVA_HOME/lib/tools.jar:$JRE_HOME/lib
export JAVA_HOME JRE_HOME PATH CLASSPATH
#编辑完后执行脚本
source /etc/profile
[root@cricaas java]# echo $PATH
/usr/lib64/qt-3.3/bin:/usr/local/sbin:/usr/local/bin:/sbin:/bin:/usr/sbin:/usr/bin:/usr/local/git/bin:/home/lilibei/bin:/root/bin:/usr/local/git/bin:/home/lilibei/bin:/usr/java/jdk1.8.0/bin:/usr/java/jdk1.8.0/jre/bin
[root@cricaas java]# echo $JAVA_HOME
/usr/java/jdk1.8.0
[root@cricaas java]#

TASSEL安装

TASSEL的安装很简单,按照如下即可

#github clone
git clone https://bitbucket.org/tasseladmin/tassel-5-standalone.git
cd tassel-5-standalone

使用

示例文件均存放于

/home/lilibei/tassel-5-standalone/TASSELTutorialData/data/

计算kinship矩阵

#适合小数据集,重测序数据需要更改内存配置参数
./run_pipeline.pl -fork1 -h TASSELTutorialData/data/mdp_genotype.hmp.txt -ck -export kinship.txt

LD分析

#绘制LD连锁图,很难看,可以使用hapview或者LDheatmap
./run_pipeline.pl -fork1 -h TASSELTutorialData/data/mdp_genotype.hmp.txt -ld -ldd png -o chr_5000sites_ld.png
#计算LD
./run_pipeline.pl -fork1 -h TASSELTutorialData/data/mdp_genotype.hmp.txt -ld -export LD_5000sites

提取染色体

#提取每条染色体
./run_pipeline.pl -fork1 -h TASSELTutorialData/data/mdp_genotype.hmp.txt -separate  -export chromosome
#提取1和3号染色体
./run_pipeline.pl -fork1 -h TASSELTutorialData/data/mdp_genotype.hmp.txt -separate 1,3  -export chromosome

MLM分析

./run_pipeline.pl -fork1 -h TASSELTutorialData/data/mdp_genotype.hmp.txt -filterAlign -filterAlignMinFreq 0.05 -fork2 -r TASSELTutorialData/data/mdp_traits.txt -fork3 -q TASSELTutorialData/data/mdp_population_structure.txt -excludeLastTrait -fork4 -k TASSELTutorialData/data/mdp_kinship.txt -combine5 -input1 -input2 -input3 -intersect -combine6 -input5 -input4 -mlm -export lilibei -runfork1 -runfork2 -runfork3 -runfork4

参考资料:

TASSEL5.0命令行官方使用手册
https://my.oschina.net/golang/blog/210306
https://bytebucket.org/tasseladmin/tassel-5-source/wiki/docs/Tassel_Pipeline_Tutorial20110511.pdf


有疑问加站长微信联系(非本文作者)

本文来自:简书

感谢作者:

查看原文:Linux系统下TASSEL5.0的安装及使用

入群交流(和以上内容无关):加入Go大咖交流群,或添加微信:liuxiaoyan-s 备注:入群;或加QQ群:692541889

753 次点击  
加入收藏 微博
上一篇:01Gin源码解读
下一篇:Golang:切片
暂无回复
添加一条新回复 (您需要 登录 后才能回复 没有账号 ?)
  • 请尽量让自己的回复能够对别人有帮助
  • 支持 Markdown 格式, **粗体**、~~删除线~~、`单行代码`
  • 支持 @ 本站用户;支持表情(输入 : 提示),见 Emoji cheat sheet
  • 图片支持拖拽、截图粘贴等方式上传