TASSEL 关联分析软件 起步教程

目录

这里下载的是win的64为系统,截图如下:

01

安装成功后,打开菜单如下:

02

截图如下:

03

打开data,load,选择Make Best Guess

04

选择几个示例数据:

05

打开后的数据如下

06

里面包括系谱数据、性状数据和基因型数据(snp)。

3.1清洗数据

选中mdp_trait,

07

然后选择:Data中的Transform Phenotype,

08

可以对数据进行转化、标准化等操作,注意,要先对数据进行选择,然后再进行操作:

09

也可以对缺失值的数据进行删除,点击imput,Numerical impute,就会生成没有缺失值的数据,这只是缺失值的不同替换方法。

10

11

3.2主成分分析

主成分分析(PCA)是一种统计方法,它可以将相互关联的变量转化为独立的主成分(PC),第一种成分包含最多的组分,其它依次降低。另一个主成分的作用可以用标记的主成分来代表群体结构。这种方法比最大似然法节省时间。因为大部分的分子标记都是字符,需要先将其转化为数值,然后再进行主成分分析,一般将纯合的标记用0代替,另一个纯合子用2代替,杂合的用1代替。PCA要求变量不能有缺失值,因此,在进行主成分分析时,需要对数据进行清洗,去除缺失值。

去掉频率小于0.05的标记,可以选择Data,选择Site,然后在最小频率的框中键入0.05,然后选择Removeminor SNP status,然后点击Filter,进行过滤,模型如下:

12

选择PCA,然后选择5个主成分(默认项),点击确定,就会生成结果,模型如下:

13

14

结果如下:

15

3.3用基因标记估计系谱

利用主成分分析可以判断群体的结构特征,但是如果利用系谱信息,这种结果会更加准确。可以用基因型数据生成系谱信息,首先选中基因型数据,点击Analysis,选择Kinship

16

结果如下:

17

3.4用一般线性模型分析GLM

下面我们用GLM模型来分析示例数据,mdp_genotype.hmp.txt是snp数据,里面有3093个标记,281个玉米自交系,另一个文件是mdp_population_structure.txt,里面是282个玉米自交系的群体结构,还有一个是mdp_traits.txt,里面是282玉米自交系的表型数据。

首先对基因型数据进行过滤,去掉频率小于0.05的,最小的数目是150,点击过滤,生成过滤后的基因型数据:

18

然后对数据进行个过滤,选择开花期dpoll这个性状,

19

进行协变量选择,即选择群体结构的文件,这里我们去掉Q3,数据如下:

20

合并数据,将这三个过滤好的数据,选中进行合并,点击Data Intersect Join,

21

数据合并如下:

22

然后选中合并后的数据,用analysis,GLM来进行分析

23

24

运行结果如下:

25

QQ图:

26

P-value值:

27

3.5用混合线性模型分析

混合模型需要添加系谱矩阵

28

点击run

29

结果:

30

相关图形:

31

32

转载本文请联系原作者获取授权,同时请注明本文来自邓飞科学网博客。

链接地址:http://blog.sciencenet.cn/blog-2577109-962946.html

上一篇:联合方差分析的误差方差齐次性检测 R语言程序

下一篇:农业上对不同品种照相,按照品种名称批量改变相片名称的方法

tassel软件使用linux,科学网—TASSEL 关联分析软件 起步教程 - 邓飞的博文相关推荐

  1. tassel软件使用linux,Tassel 命令行版本使用笔记

    1. 常见的基因型数据格式 1.1 Hapmap 数据格式 1 2 3 4 5 6 7 8 9 10 11 12 rs# alleles chrom pos strand assembly# cent ...

  2. gpt efi win7 linux,科学网—UEFI+GPTSSD+HDD 成功安装win10+Centos linux7 过程 - 陈长云的博文...

    主机为:Dell Precision T7910台式工作站  SSD + HDD 双硬盘 1.必备工具 ① Disk Genious ② Win10系统(光盘和U盘,我用的是U盘),网上可以下载 cn ...

  3. 生信软件c语言,科学网—[转载]没有docker我真的不想动这样的生信软件 - 张成岗的博文...

    没有docker我真的不想动这样的生信软件 2020-03-26阅读 2620 C语言源代码需要编译的软件 最开始开发者都是C语言流派, 所以标准的源代码安装三部曲即可,即使 configure+ma ...

  4. docker安装gamit_科学网—Ubuntu系统GAMIT/GLOBK程序安装 - 陈超的博文

    最近开始学习GAMIT,网上资料还是蛮多的,但是感觉都是东拼西凑的,一点都不系统,一点不适合初学者.安装教程也是乱七八糟的,下面把我的安装过程分享一下:(我也是参考的网上一个教程,很久以前下载的,地址 ...

  5. pfamscan 的使用_科学网—[转载]InterProScan的使用教程 - 黄顺谋的博文

    InterPro是一个数据库,其提供蛋白序列的功能分析并归纳为一个个蛋白家族,同时还预测了presence of domains和important sites.为了将蛋白分类,InterPro使用先 ...

  6. 乐符识别matlab,科学网—[原][Matlab][04] Midi音乐键盘 - 王楠的博文

    以前在网上读到一个用matlab播canon音乐的源码感觉很有意思,但声音还不够好. 于是想做个自己的音乐播放器,有键盘,有对应的钢琴声音,读取乐谱就能播.如下. (1)按键的音高与频率 从1到下一个 ...

  7. 信息与计算机,科学网—信息与计算机(1) - 姜咏江的博文

    信息与计算机(1)信息的概念姜咏江 当今科学研究中使用最多的概念是什么?信息!然而信息却是那样地让人们困惑.有句话说:"把复杂的问题说简单了,那是学问:将简单的问题说复杂了,那是蒙人!&qu ...

  8. matlab音乐键盘模拟,科学网—[原][Matlab][04] Midi音乐键盘 - 王楠的博文

    以前在网上读到一个用matlab播canon音乐的源码感觉很有意思,但声音还不够好. 于是想做个自己的音乐播放器,有键盘,有对应的钢琴声音,读取乐谱就能播.如下. (1)按键的音高与频率 从1到下一个 ...

  9. 主成分分析碎石图_科学网—主成分分析之R篇 - 林元震的博文

    R中psych包可以进行主成分分析,其分析的步骤为: (1) 判断主成分的个数: (2) 提取主成分: (3) 获取主成分得分: (4) 列出主成分方程,解释主成分意义. [例子] 测定了20株杨树树 ...

  10. meep php,科学网—centos5.4下meep配置 - 郑改革的博文

    centos5.4下meep配置 已有 4453 次阅读 2010-5-20 19:19 |个人分类:未分类|系统分类:科研笔记| meep 1. install BLAS sudo yum inst ...

最新文章

  1. qtdll在linux系统运行,在QT下编写带DLL的程序
  2. 图解排序算法之3种简单排序(选择,冒泡,直接插入)
  3. VsCode 配置 C/C++ 开发环境,真的很简单
  4. jQuery处理点击父级checkbox所有子级checkbox都选中,取消选中所有子级checkbox都取消...
  5. 如何查看sqlserver日志的方法
  6. 上位机和下位机的概念,理解如何实现PC从PLC中读取数据?
  7. 初试Node —— node.js的安装
  8. php 禁用外部实体,php – Doctrine 2 – 从实体外部禁用PrePersist
  9. 【java】java 随机数 Random ThreadLocalRandom SecureRandom
  10. w ndows10电脑配置看哪里,win10对电脑配置要求及如何查看电脑配置
  11. php中的rtrim_PHP函数rtrim()使用中的怪异现象
  12. java——jui的应用
  13. 米兔机器人恐龙拼图手册_米兔积木机器人(多形态组合、自平衡系统、App联动)...
  14. 【代码审计】51 TP5框架、无框架 变量覆盖反序列化
  15. 群晖服务器显示灯,【群晖 DS119j 单盘位NAS 网络存储服务器使用总结】功能|配置|机身|指示灯_摘要频道_什么值得买...
  16. 用计算机弹小猪佩奇,【绘本时间】小猪佩奇之弹弹球
  17. kafka连接mysql数据库,Kafka-connect-jdbc-source连接mysql数据库实战
  18. Windows终端美化
  19. mysql的gis_MySQL的GIS功能
  20. 深度解析用户画像标签体系构建方法

热门文章

  1. 手机 人人网android 2.2,人人网客户端安卓版
  2. 64 The HTTP header line [ssl-client-cert:12330x11] does not conform to RFC 7230 and has been ignored
  3. 程序员也要学英语——动词相关合集
  4. jdk1.8新特性:函数式接口、方法引用、函数式编程、常用函数式接口
  5. Youtube——如何将视频中的英文字幕转换成中文字幕
  6. 面试官问你还有什么要问我的吗?
  7. 三角函数π/2转化_浅谈三角函数导数综合题(1)
  8. java 间接引用_Java中的四种引用方式的区别
  9. 毕业设计——基于STM32的智能家具系统(语音识别控制、步进电机、舵机)
  10. 中标麒麟Neokylin7桌面版安装指南——基于VirtualBox虚拟机