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

博文

R语言计算玉米测交试验的配合力

已有 467 次阅读 2020-9-15 20:52 |个人分类:农学统计|系统分类:科研笔记

1. 玉米测交试验

可以认为,玉米侧交试验,是一个NCII的试验,在玉米实际的育种中,由于其测用结合的特性,应用广泛。

试验设计如下:

骨干系:4个 待测系:7个 重复:2次 观测值:产量

「计算目的:」

  • 1 计算待测系的一般配合力(选择最优品系)
  • 2 计算待测系和骨干系的特殊配合力(选出最优组合)

2. 概念解释

「配合力:」 杂交育种成功的关键在于正确选配亲本,因为杂种的表现如何,不仅取决于自交系本身的农艺性状,而且决定于亲本间的相互关系,这种相互关系成为配合力。

「一般配合力(GCA)」 如果某一品系与许多品种、自交系杂交,其F1的平均产量成为该品种的一般配合力(GCA)。

「特殊配合力(SCA)」 将其中某个特定组合的表型与一般配合力之差数称为该品系的特殊配合力(SCA),它反映了两个亲本之间的特殊关系。

「配合力的作用」 育种过程中,可以选用一般配合力高的自交系组成轮回群体,以积累加性基因,提高一般配合力。选用特殊配合力高的亲本组配杂交组合,可以通过多基因位点的互作,获得较高的杂种优势效应。

刘纪麟《玉米育种学》 第二版,中国农业出版社 P94

3. 计算流程

  • 1 读取数据
  • 2 查看数据结构,转化为因子
  • 3 构建模型
  • 4 提取配合力结果
# 读取数据
dat = read.csv("NCII_dat.csv")
str(dat)

# 转化为因子
for(i in 1:5) dat[,i] = as.factor(as.character(dat[,i]))
str(dat)

# 构建模型
library(sommer)
mod = mmer(yield ~ Block, random = ~ P1 + P2 + Fam, data=dat)
summary(mod)

# 提取配合力结果
gca1 = mod$U$P1
head(gca1)

gca2 = mod$U$P2
gca2

sca = mod$U$Fam
sca

4.其它参数的计算

「配合力的显著性:」 配合力的效应,是效应值(Effect),如果判断显著性,先看方差分析中GCA和SCA是否达到显著性,然后计算每个自交系的GCA的显著性。

因为GCA是总和为0的效应值,如果判断每个品种的GCA是否是显著的,其实是和0进行的比较,其实是进行的是T检验,可以计算LSD0.05和LSD0.01,比较GCA的值,判断其显著水平。

「狭义遗传力与广义遗传力」

不同近交系数,计算方法是不一样的,感兴趣的老师,可以参考我之前写的博客:如何使用NC II 遗传设计估算配合力和遗传力?

5. 数据介绍

「表型数据:」 phenotype

6. 结果怎么看?

上面P11P12P13P14为骨干系(测验种),可以看到P12位1.15,最高,说明P12和所有待测系杂交后代的平均值最高,而且比整体平均值高1.15.

注意,这里的一般配合力都是相对于整体平均值的效应值,他们的整体之和为0,为正表示高于整体平均值,为负表示低于整体平均值。

这里的P20~29为待测系,可以看到,P211的一般配合力为2.08,表现最好。

注意,一般配合力,都是针对于自交系本身的,它的值主要是在本实验内才有意义,重点是排名,为了筛选。

上面是特殊配合力的结果,Fam1_10表现最好,它的亲本分别是:P1和P10,说明他们之间的特殊配合力最高。

注意:特殊配合力为两者之间的值,不能讲某个自交系的特殊配合力是多少,要讲两个自交系间的特殊配合力多高。它也是根植于具体的试验,用于筛选排名。

7. 下载示例数据和代码

**关注公众号:**育种数据分析之放飞自我 回复关键词:「NC」 下载相关数据和代码。

8. 上面数据如何用Excel实现?

确实,Excel是很香的,我被师兄教育很多次了,学什么R语言,Excel能解决的问题,不要增加工具的学习成本,大家都很忙,王者峡谷不香吗???

所以,小伙伴们多多转发点赞吧,点赞数超过40,我就用Excel实现一下!

【手动狗头】:上一次,你这样讲TASSEL做GWAS的教程,说是阅读量达到1000的话,会录制视频教程。教程在哪里???

每一个Flag都在我的心中,从来没有忘却,视频已经安排上了,准备用鬼畜的风格做一期,哈哈……




http://blog.sciencenet.cn/blog-2577109-1250738.html

上一篇:R语音读取txt报错:UTF-16不支持 或者<ff><fe><63>\'多字节字符串有错
下一篇:GWAS分析中:为何有些SNP效应值大却不显著?

0

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

数据加载中...

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

GMT+8, 2020-10-29 05:24

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部