欢迎光临散文网 会员登陆 & 注册

R|rTASSEL进行基因组分析(一)

2023-03-29 23:37 作者:伏见锁骨上的印记_  | 我要投稿

TASSEL可用于评估性状关联,进化模式和连锁不平衡。本文介绍的rTASSEL是利用R作为前端,使用TASSEL的分析能力以及R数据处理和解析功能进行分析的R包。

1.安装rTASSEL

    方法一:

    方法二:

    *方法二安装出现报错

Failed to install 'rTASSEL' from GitHub

        *rTASSEL安装与运行需要rJava,使用命令行安装rJava(install.packages("rJava"))后继续执行上述命令。显示DONE(rTASSEL)表示安装成功。

rTASSEL安装成功界面

2.运行rTASSEL


正确运行rTASSEL

3.读取数据(演示使用rTASSEL中的数据进行)


基因型数据
表型数据
查看数据信息

*内置数据为279份样品,3093个SNP位点

4.计算亲缘关系

亲缘关系数据结果查看

5.计算混合线性模型(MLM)

数据结果查看

*结果包含三个文件

6.数据结果可视化绘制曼哈顿图(Manhattan)

曼哈顿图结果

7.数据结果可视化绘制连锁不平衡图(LD)

连锁不平衡图结果


预告:

1.绘制曼哈顿图有专门的R包,如CMplot,qqman。

2.数据前处理,过滤最小等位基因频率小于0.05的位点。

3.使用自己的数据进行分析。

……【后续想起来再补充】


文献参考:

Bradbury PJ, Zhang Z, Kroon DE, Casstevens TM, Ramdoss Y, Buckler ES. (2007)TASSEL: Software for association mapping of complex traits in diverse samples.Bioinformatics 23:2633-2635.

Monier et al., (2022). rTASSEL: An R interface to TASSEL for analyzing genomic diversity. Journal of Open Source Software, 7(76), 4530, https://doi.org/10.21105/joss.04530


【最复杂的代码,做最简单的事。】 

R|rTASSEL进行基因组分析(一)的评论 (共 条)

分享到微博请遵守国家法律