大家好, 我是邓飞。

今天聊一下广义遗传力的计算方法。

广义遗传力定义

广义遗传力( H 2 H^2 H2)定义为归因于基因型总体遗传变异的表型变异比例。
通常他包括三个解释:(详见我这篇博客的公式推导 回归系数 相关系数 遗传力的关系)

  • (i)它相当于观察到的表型上不可观察的基因型值的线性回归的确定系数;
  • (ii)也是预测表型值和基因型值之间的相关系数的平方;
  • (iii)它表示可以实现为选择响应(RR)的选择差异(SS)的比例(Falconer和Mackay,2005)。

有两个主要原因说明植物育种中对输入平均值的遗传力感兴趣(Schmidt等人,2019a):

  • 1,计算这个参数,可以帮育种家预测选择相应

  • 2,这是一种用于评估品种评估试验结果的有用性和准确性的描述性度量。

2,One stage和Two stage

Two stage,两阶段方法,就是先计算单个地点的育种值,然后再计算一年多点或者多年多点的数据,这里多个年份或者多个地点,都看做是不同的环境。One stage是一步到位的MET的分析方法,可以是一年多点,或者是多年多点,都在一个模型中进行分析。

2.1 Two stage

y i k t = μ + G i + E t + G x E i t + ε i k t y_{ikt} = μ + G_i+E_t+GxE_{it}+ε_{ikt} yikt​=μ+Gi​+Et​+GxEit​+εikt​

这里的表型方差组分为:

σ p 2 = σ g 2 + σ g ∗ e 2 n e + σ e 2 n e ∗ n r \sigma_p^2 = \sigma_g^2 + \frac{\sigma_{g*e}^2}{n_e} + \frac{\sigma_{e}^2}{n_e*n_r} σp2​=σg2​+ne​σg∗e2​​+ne​∗nr​σe2​​

这里的n_e,是环境的个数,n_r是重复的个数。

2.2 One stage

这里,Y是年份,E是地点,公式中包括:

  • 基因型
  • 年份
  • 地点
  • 年份与地点互作
  • 基因型与年份互作
  • 基因型与地点互作
  • 基因型与年份与地点互作
  • 残差

表型方差的计算方法:

每一项,都除以对应的重复数。如果数据有缺失值,用调整均值进行分析,方法参考:

3. 更广义的遗传力计算方法

3.1 标准方法

第一种,标准方法,我们最常用的,Vg/Vp,根据一年多点或者多年多点,计算出表型的方差组分,然后用下面公式计算:

H 2 = σ g 2 σ p 2 H^2 = \frac{\sigma_g^2}{\sigma_p^2} H2=σp2​σg2​​

3.2 Cullis遗传力计算

第二种,Cullis方法

3.3 Piepho遗传力计算

第三种:Piepho方法

4. 代码数据演示

4.1 标准方法演示

4.2 Cullis遗传力演示

4.3 Piepho遗传力演示

植物育种中广义遗传力的定义相关推荐

  1. GS在植物育种中的成本和遗传增益分析

    这是2010年的文章, 文章的亮点是基因组选择在植物育种中的优势分析, 从时间成本, 金钱成本, 效益成本进行分析, 今天看来还是有很强的参考性. 1. 参考文献 Heffner E L , Lore ...

  2. ISME:全基因组关联研究揭示了控制根际微生物组遗传力的植物基因位点

    编译:微科盟阿Z,编辑:微科盟木木夕.江舜尧. 微科盟原创微文,欢迎转发转载,转载须注明来源<微生态>公众号. 导读 最近已证明宿主遗传学是植物微生物组组成的驱动因素之一.然而,确定控制微 ...

  3. 育种中一般配合力和特殊配合力的计算方法

    研究生期间,学习数量遗传学时,对配合力的概念感到很神奇,曾经和同学讨论过,按道理说曹操的一般配合力很高,几个孩子都很出色. 那么问题来了,什么是一般配合力?一般配合力就是可以稳定遗传的因素,即曹操和不 ...

  4. EJSS: 南土所梁玉婷组-生境决定了土壤植物系统中细菌、抗性基因和可移动元件间的关系...

    编译:微科盟黄瑞林,编辑:微科盟木木夕.江舜尧. 导读 土壤抗生素抗性主要是由细菌群落组成决定的.然而,植物微生物群的抗生素耐药性及其与土壤-植物系统中土壤微生物组的关系在很大程度上仍是未知的.在这里 ...

  5. 不同试验设计遗传力的计算方法

    前言 农业和林业, 经常涉及到要计算遗传力的问题, 这是一个数量遗传学问题. 和动物计算遗传力不同, 植物和林木计算遗传力时, 一般是使用家系遗传力, 动物计算遗传力一般是计算个体遗传力. 他们的主要 ...

  6. 遗传进展系列 | 1. 回归系数 相关系数 遗传力的关系

    大家好,我是飞哥. 前几天群里有老师问如何计算遗传进展​? 怎么根据公式计算遗传进展? 遗传进展不同的元素是怎么来的? 如何手动计算遗传进展? 如何使用软件计算遗传进展? 我计划写几篇博客总结学习一下 ...

  7. 细胞工程在植物组织培养中的应用

    细胞工程在植物组织培养中的应用 [文章摘要] 目前,生物技术正在世界突飞猛进地发展,而且在医学.农业.食品工业.能源工业.环境保护各个领域显示出极大的生产潜力.作为生物技术有力手段的组织培养,也日益受 ...

  8. PJ|韩方普实验室在小麦远缘杂交及抗赤霉病育种中取得新进展

    小麦-长穗偃麦草易位系的系统发展及其在小麦赤霉病抗性育种中的应用 Systemic development of wheat-Thinopyrum elongatum translocation li ...

  9. stk中天体坐标系的定义

    stk中天体坐标系的定义 中心天体坐标系统 每个中心天体都有相应的坐标系统.坐标系原点均位于中心天体的质心,区别在于坐标轴的定义. 每个中心天体均支持固定坐标系(Fixed coordinate sy ...

最新文章

  1. python 功能 代码_挑战“不可能”的代码:你不知道的Python功能
  2. 权威赛事来了!千言-文本生成评测启动,聚焦NLG技术痛点
  3. 消息队列(2):常见的消息队列协议
  4. mysql如何创建用户代码_MySQL创建用户与受权方法_mysql
  5. Activiti工作流的简单介绍
  6. .html好 还是.asp好,各位说说在ASP.net里 用静态函数的好 还是实例函数出处HTML好???...
  7. yolov3目标检测android,目标检测 | YOLOv3训练自己的数据全流程
  8. Atitit java ftp client sumup apache common net jar 1.1. 协议解析 1 1.2. 读取文件 1 2. } 3 1. 下载文件 ftpCli
  9. javaparser - java源码分析修改框架
  10. 2021数学建模学习笔记
  11. 修改win10更新服务器,修改win10更新服务器地址
  12. HEIC图片如何批量转换成jpg格式
  13. 学习整理软测(八)-----数据库命令与操作
  14. SQL总结 学期前8周学习内容
  15. leetcode 19. 删除链表的倒数第 N 个结点(c++)
  16. Taro 小程序开发大型实战(六):尝鲜微信小程序云(上篇)
  17. Java实现O(nlogn)最长上升子序列
  18. vue 使用百度地图(全景图)
  19. 在Windows 10下搭建Android开发环境
  20. 网站域名被拦截,网站被流量劫持,该怎么做?

热门文章

  1. 软件分享:最好用的免费图片转漫画风图片软件,真的不要吗?
  2. pari的用法(stl)
  3. 今天在csdn看到一博主今年大三,初中就喜欢上编程,学习了各种知识,初中…挺震惊。不禁细数我的编程史,挺有感。
  4. 面向 Java 人员的 Python
  5. MoveIt! 学习笔记14 - Kinematics Configuration/Fake Controller Manager
  6. UI设计专业术语大汇总啦!
  7. odoo 如何解决 由于附件过多 备份超时的现象
  8. 统计蚂蚁(ants)
  9. uniapp 页面登录控制
  10. 沥青都市7723java_“横看成岭侧成峰,远近高低各不同”——线激光车辙测试仪(LIR)及其在沥青路面车辙检测中的应用...