首页
学习
活动
专区
圈层
工具
发布
首页
学习
活动
专区
圈层
工具
MCP广场
社区首页 >专栏 >配合力为负值咋整啊???

配合力为负值咋整啊???

作者头像
邓飞
发布于 2024-03-25 01:38:51
发布于 2024-03-25 01:38:51
17110
代码可运行
举报
运行总次数:0
代码可运行

大家好,我是邓飞,作为育种数据分析的段子手,数据分析水平不敢保证,但是段子确实香。经常有小伙伴抛出各种育种数据分析相关的问题,当问题比较多了,我觉得应该写一篇博客了,毕竟,别人问你一个问题,再没有抛过去一篇博客,什么话都不说更爽的了,暗爽。

作为一个资深的水电工,身上带着锤子扳手之类的很正常,作为一个资深的数据分析时,手里面有各种博客用于答疑也很正常。

常见问题:

1,gca和sca都为0,是什么原因?

答:这是因为亲本之间没有显著性差异,配合力分析的前提是方差分析中亲本之间达到显著性差异,要不然gca和sca都没有意义。

2,模型迭代失败,是什么原因?

答:NCII试验,需要有重复数据,否则不能计算gca和sca

3,方差分析和混合线性模型LMM模型有什么区别

答:LMM模型结果更靠谱,而且可以直接计算方差组分,计算广义遗传力和狭义遗传力。当有缺失值存在时,推荐使用LMM模型,相关的包有:lme4,sommer,asreml等。

4,配合力为何有正有负?

答:配合力之和为0,所以有正有负,正的就是高于均值,负的就是低于均值,0就是均值。所以,选择配合力高的,就选择最大的就行。

常见的问题都有了,举个栗子介绍一下如何进行配合力分析吧,毕竟说了一大堆有钱真幸福的话,没有告诉大家如何搞钱显然不是一个合格的安利文。

…………

…………

1. 玉米测交试验

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

试验设计如下:

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

「计算目的:」

  • 1 计算待测系的一般配合力(选择最优品系)
  • 2 计算待测系和骨干系的特殊配合力(选出最优组合)
  • 3 了解材料的遗传特性

2. 概念解释

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

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

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

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

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

3. 计算流程

  • 1 读取数据
  • 2 查看数据结构,转化为因子
  • 3 构建模型
  • 4 提取配合力结果
代码语言:javascript
代码运行次数:0
运行
AI代码解释
复制
# 读取数据
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的显著性。e

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

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

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

5. 数据介绍

「表型数据:」

6. 结果怎么看?

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

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

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

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

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

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

7. 下载示例数据和代码

本文参与 腾讯云自媒体同步曝光计划,分享自微信公众号。
原始发表:2024-03-23,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 育种数据分析之放飞自我 微信公众号,前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体同步曝光计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
1 条评论
热度
最新
现在根本找不到这个镜像!
现在根本找不到这个镜像!
回复回复点赞举报
推荐阅读
编辑精选文章
换一批
如何通过Genstat软件来计算遗传力和配合力?
对于不经常编程的小伙伴,Genstat窗口的选择模型,方差分析、混合线性模型、配合力分析、基因与环境互作等功能,特别适合农业方向的数据分析。这个软件功能很强大,通过鼠标点击的形式就可以进行复杂模型的分析,非常666。
邓飞
2024/03/25
6830
如何通过Genstat软件来计算遗传力和配合力?
育种中一般配合力和特殊配合力的计算方法
研究生期间,学习数量遗传学时,对配合力的概念感到很神奇,曾经和同学讨论过,按道理说曹操的一般配合力很高,几个孩子都很出色。
邓飞
2019/12/05
1.7K0
NCII试验 计算配合力和遗传力的方法
育种中,有几个必须要掌握的概念,配合力是其中之一。配合力包括一般配合力和特殊配合力。这个概念很抽象,下面用曹操的例子,解释一下。
邓飞
2022/12/13
4270
NCII试验 计算配合力和遗传力的方法
GS在植物育种中的成本和遗传增益分析
文章的亮点是基因组选择在植物育种中的优势分析, 从时间成本, 金钱成本, 效益成本进行分析, 今天看来还是有很强的参考性.
邓飞
2019/06/13
1.3K0
如何使用NC II 遗传设计估算配合力和遗传力?
非常令人疑惑,于是我就多看了一些教材,发现是有一个万能公式的,即是考虑近交系数,整体结论:
邓飞
2020/08/24
1.6K0
如何使用NC II 遗传设计估算配合力和遗传力?
如何计算配合力+方差组分+遗传力
育种中,有几个必须要掌握的概念,配合力是其中之一。配合力包括一般配合力和特殊配合力。这个概念很抽象,下面用曹操的例子,解释一下。
邓飞
2022/12/13
6920
如何计算配合力+方差组分+遗传力
我的农学转数据分析之花花草草
数据分析界育种知识最好、育种界编程最扎实、段子讲得最好的数据分析师,所以:编程+数据分析+育种,就是我的日常工作了。
邓飞
2022/12/12
5650
我的农学转数据分析之花花草草
通过基因组选择预测杂交水稻的表现(数据挖掘)
基因组选择,预测杂种优势,在水稻可以用,在玉米,高粱中也可以用,在动物选择配套系时也可以用,根据加性效应和非加性效应进行预测,前景广阔。
邓飞
2019/06/13
1K0
GLMM:广义线性混合模型(遗传参数评估)
这篇文章,主要是介绍了抗性数据,如何利用GLMM模型进行的分析,文中,他将9级分类性状变为了二分类性状,进行分析。
邓飞
2020/09/30
2.2K0
GLMM:广义线性混合模型(遗传参数评估)
全同胞家系如何计算遗传力及育种值
大家好,我是飞哥,今天是农历七月初七,中国的情人节,祝大家节日快乐,有情人终成眷属,单身的尽快脱单。
邓飞
2022/12/13
8490
全同胞家系如何计算遗传力及育种值
当我们在说方差分析时,我们在说些什么?
方差分析或变方分析(Analysis of variance,简称ANOVA)为数据分析中常见的统计模型,主要为探讨连续型(Continuous)资料型态之因变量(Dependent variable)与类别型资料型态之自变量(Independent variable)的关系,当自变项的因子中包含等于或超过三个类别情况下,检定其各类别间平均数是否相等的统计模式,广义上可将T检定中方差相等(Equality of variance)的合并T检定(Pooled T-test)视为是方差分析的一种,基于T检定为分析两组平均数是否相等,并且采用相同的计算概念,而实际上当方差分析套用在合并T检定的分析上时,产生的F值则会等于T检定的平方项。
邓飞
2019/09/25
1.6K0
当我们在说方差分析时,我们在说些什么?
文献阅读: 林木中遗传参数评估
这篇文章, 非常具有代表性, 可以为林木和作物的数据分析提供思路, 不一定非要有系谱才可以计算育种值和遗传相关, 混合线性模型代替一般线性模型进行育种值的筛选, 是大势所趋, 必须推而广之.
邓飞
2019/06/13
6620
统计遗传学:第二章,统计分析概念
前几天推荐了这本书,可以领取pdf和配套数据代码。这里,我将各个章节介绍一下,总结也是学习的过程。
邓飞
2022/07/27
8000
统计遗传学:第二章,统计分析概念
环境遗传相关 | 育种中的基因与环境互作
基因与环境互作,植物中同一个品种多年多点种植,评价基因与环境互作,找到品种最适合推广的区域。
邓飞
2022/12/13
1.1K0
环境遗传相关 | 育种中的基因与环境互作
完结篇 | GWAS计算BLUE值4--联合方差分析演示
本篇,用书籍中的数据和结论,用R语言的一般线性模型和混合线性模型,做一下一年多点的联合方差分析的演示。
邓飞
2021/12/20
8110
完结篇 | GWAS计算BLUE值4--联合方差分析演示
多年多点数据如何计算遗传力以及BLUP值
之前写过一篇博客, 介绍领导安利我哔哩哔哩的故事, 介绍到我将我从YouTube上收集的关于混合线性模型, 关于GWAS, 关于GS, 关于农业数据分析相关的视频, 上传到了哔哩哔哩上面. 今天我们看一下介绍多年多点遗传力及BLUP值计算的视频内容. 阅读原文可以查看视频, 这里我用文字和代码进行重演.
邓飞
2019/06/13
4.9K1
GWAS计算BLUE值3--LMM考虑残差异质计算BLUE值
本节,介绍如何使用R语言的asreml包拟合混合线性模型,定义残差异质,计算最佳线性无偏估计(blue)
邓飞
2021/12/20
8960
GWAS计算BLUE值3--LMM考虑残差异质计算BLUE值
动物育种统计发展的百年--翻译版
Daniel Gianola1–5 and Guilherme J.M. Rosa1,2
邓飞
2023/09/06
4040
动物育种统计发展的百年--翻译版
农学的为何要学习神经网络???
在我刚开始学习GS的时候,我是从混合线性模型(LMM)的基础上理解的,因为动物模型BLUP,所以基因组选择GBLUP,再所以一步法ssGBLUP。
邓飞
2022/12/13
3430
农学的为何要学习神经网络???
MP长篇综述 | 植物泛基因组及其应用
2022年12月15日,中山大学史俊鹏副教授、中国科学院遗传与发育生物学研究所田志喜研究员、中国农业大学赖锦盛教授和上海师范大学黄学辉教授共同撰文,在Molecular Plant杂志发表了题为“Plant pan-genomics and its applications”的长篇综述。该论文对植物中泛基因组产生的源动力、表征、分析方法以及在植物遗传学和育种中的应用进行了系统总结和展望。
生信宝典
2023/08/30
6920
MP长篇综述 | 植物泛基因组及其应用
推荐阅读
相关推荐
如何通过Genstat软件来计算遗传力和配合力?
更多 >
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档
本文部分代码块支持一键运行,欢迎体验
本文部分代码块支持一键运行,欢迎体验