育种数据分析之放飞自我分享 http://blog.sciencenet.cn/u/yijiaobai 关注:生物统计,数量遗传,混合线性模型,生物信息,R,Perl,Python,GWAS,GS相关方法,文章及代码

博文

如何计算遗传相关(How to)

已有 9197 次阅读 2018-5-7 17:51 |个人分类:数量遗传学|系统分类:科研笔记

Q1

植物或林木育种中,我的数据类型包括基因型(品种),地点(环境),重复(或者无重复),还有一系列性状,比如株高、产量等等,我怎么计算性状间的遗传相关,并且给出显著性?

示例数据

数据描述:location为地点,Genotype为基因型,Rep为重复,DWKPlant, DWKEar, TW, KDM, TKW为性状,想要计算这5个性状之间,两两的遗传相关和表型相关,并给出显著性
图片.png

常见误区

直接计算表型值的相关系数r,r这里并不是遗传相关和表型相关,而只是性状间的相关系数

图片.png

正确做法

首先要运行多性状模型,计算性状间的遗传协方差组分(Gcov12)和遗传方差组分(G11,G22)

图中,黄色的表示方差组分,绿色的表示协方差组分,蓝色的是计算的遗传相关:
计算公式:
图片.png

图片.png

GenStat界面

数据

图片.png

模型

图片.png

asreml-r代码

mod <- asreml(cbind(DWKEar, DWKPlant, KDM, TKW, TW) ~ trait + trait:(location + Rep), random=~ us(trait):Genotype, rcov =~ units:us(trait), data=dat)



https://blog.sciencenet.cn/blog-2577109-1112805.html

上一篇:R bug
下一篇:混合线性模型中BLUE值 VS BLUP值
收藏 IP: 125.37.212.*| 热度|

1 杨正瓴

发表评论 评论 (0 个评论)

数据加载中...
扫一扫,分享此博文

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

GMT+8, 2024-4-26 22:27

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部