第一节线性模型基础知识第二节BLUP的基本原理第三节BLUP课件.ppt

1、v第一节 线性模型基础知识v第二节 BLUP的基本原理v第三节 BLUP的计算技术v第四节 育种值估计模型 v第五节 多性状BLUP 法的基本原理v第六节 BLUP育种值估计举例v第七节 BLUP育种值估计软件http:/ Best Linear Unbiased PredictionLinear Unbiased Prediction）育种值估计法，育种值估计法，将畜禽遗传育种的理论与实践带入了一个新的发将畜禽遗传育种的理论与实践带入了一个新的发展阶段。展阶段。http:/ 线性模型基础知识线性模型基础知识 v模型（模型（ModelModel）模型是描述观察值与影响观察值变异性的各因子之间的

2、关系的数学方程式 分类 真实模型非常准确地模拟观察值的变异性，模型中不含有未知成分 理想模型根据研究者所掌握的专业知识建立的尽可能接近真实模型的模型 操作模型用于实际统计分析的模型，它通常是理想模型的简化形式 http:/ 离散型通常表现为若干个有限的等级或水平固定因子 有意识地抽取若干个特定的水平，目的是对这些水平的效应进行估计或进行比较，如年效应 随机因子因子的若干水平可看作是来自该因子的所有水平所构成的总体的随机样本，目的是要通过该样本去推断总体，如个体的遗传效应。连续型 它呈现连续性变异，通常是作为影响观察值的协变量（回归变量）http:/ 线性模型（线性模型（Linear model

3、Linear model）线性模型是指在模型中所包含的各个因子是以相加的形式影响观察值，即它们与观察值的关系为线性关系，但对于连续性的协变量也允许出现平方或立方项。一个线性模型应由3个部分组成：1.数学方程式2.方程式中随机变量的期望和方差及协方差3.假设、约束和限制条件 http:/ 数学方程式：数学方程式：v 期望和方差：期望和方差：v 假设和约束条件：假设和约束条件：所有母牛都来自同一品种 所有母牛都在相同的环境下以相同的饲养方式饲养 所有的母牛都来自同一公牛 所有的母牛的母亲对母牛的乳脂量无影响 分组分组 初产年龄（等级）初产年龄（等级）123产犊季节产犊季节 1114 143 145

4、 150 2109 103 163 117 ijkjiijkebayv线性模型举例线性模型举例母牛的乳脂量生产成绩表（Schaeffer L R，1993）jiijkbayE)(0)(ijkeE2)()(iijkijkeVyVhttp:/ 固定效应模型（fixed model）如一个模型中除了随机误差外，其余所有的效应均为固定效应，则称此模型为固定效应模型或固定模型。随机效应模型（random model）若模型中除了总平均数外，其余的所有效应均为随机效应，则称此模型为随机效应模型或随机模型 混合模型（mixed model）若模型中除了总平均数和随机误差之外，既含有固定效应，也含有随机效应，

5、则称之为混合模型 http:/ 传统的选择指数法的基本假设是传统的选择指数法的基本假设是不存在影响观察值的系统环境效应不存在影响观察值的系统环境效应，或者或者在使用前剔除了系统环境效应在使用前剔除了系统环境效应。v 遗憾的是这个基本假设在遗憾的是这个基本假设在几乎所有实际情况下都是不能成立的几乎所有实际情况下都是不能成立的，如乳，如乳用母牛饲养在管理条件不同的牛群中。用母牛饲养在管理条件不同的牛群中。v 为克服以上缺陷，为克服以上缺陷，HendersonHenderson于于19481948年提出了年提出了BLUPBLUP方法，即方法，即最佳线性最佳线性无偏预测无偏预测，这个统计方法可同时估计

6、固定效应（例如系统环境效应），这个统计方法可同时估计固定效应（例如系统环境效应）和育种值。传统的选择指数是具有已知固定效应的和育种值。传统的选择指数是具有已知固定效应的BLUPBLUP方法的一种特方法的一种特殊情形。殊情形。v 随着计算机技术的高速发展，使这一方法的实际应用成为可能，目前随着计算机技术的高速发展，使这一方法的实际应用成为可能，目前BLUPBLUP法已成为世界各国（尤其是发达国家）家畜遗传评定的法已成为世界各国（尤其是发达国家）家畜遗传评定的规范方规范方法法 。http:/ BLUPBLUP的基本原理的基本原理v 一般混合模型可表示为：一般混合模型可表示为：是所有观察值构成的向量

7、 是所有固定效应（包括）构成的向量是固定效应的关联矩阵 是所有随机效应构成的向量 是随机效应的关联矩阵 是随机残差向量 v 随机变量的数学期望：随机变量的数学期望：v 方差方差-协方差矩阵结构：协方差矩阵结构：ybXuZeeZuXbybb)(E，0u)(E，0e)(E，Xby)(ER00GeuVarhttp:/ BLUP 的统计特性的统计特性 可估函数：预测函数：预测误差：BLUP分析的实质是利用观察值的一个线性函数（）对固定效应和随机效应的任意线性可估函数（）进行估计和预测，要求同时满足预测的无偏性和预测误差方差最小（最佳）两个条件，由此得到 的最佳线性无偏估计值最佳线性无偏估计值（BLUE

8、），的最佳线性无最佳线性无偏预测值（偏预测值（BLUPBLUP）。uMbKyLyLuMbKyLuMbKbuhttp:/ BLUP BLUP 估计一般方程估计一般方程 v BLUPBLUP法前提条件法前提条件1.所用的表型信息必须真实可靠，系谱资料必须正确完整2.所用的模型是真实模型；3.模型中的随机效应的方差组分或方差组分的比值已知yVXX)VX(b11)bX(yVZGu1http:/ zzxxVarCG00CubzzzxxzxxVarCCCCuubkdduuCovriiiiiiuueauuuiiuu1/)(),(22222uek 为 中与 个体对应的对角线元素 iudzzCihttp:/ B

9、LUPBLUP的计算技术的计算技术v混合模型方程组的求解混合模型方程组的求解 经典解法 先求出方程组的系数矩阵和等式右边的向量，建立方程组，然后迭代求解 缺点：混合模型方程组往往很大，容易受计算机内存的限制，实际应用范围不广 间接解法 不需建立方程组，直接构建观测数据迭代公式，每次迭代读入原始数据包括性状观测值和系谱记录，并同时计算该次迭代的解 通用性不强，需要构建特定的数据迭代公式http:/ 高斯-赛德尔迭代法（gauss-seidel）雅可比迭代法（jacobi）松弛迭代法（relaxation）v收敛标准收敛标准（convergence criteria)convergence cri

10、teria)一般标准 改进标准iiijnijkjijkjijikicxcxcrx)(111)1()()()1(1)1()()(kiiinjkjijikixcxcrx)()1()()()(kikikiimprovedixxxx)1()(maxtitixx21)(12)1()()()(nitinititixxxhttp:/ 的父母未知时：个体 的父或母为 时：个体 的父母已知为 或 时：121ti、tttpq1ttaipittiaaa5.0121ti、p1tta0ittiaa pqttaa5.01)(5.0iqipittiaaaa121ti、http:/ 1.构造所有个体的系谱列表，父母亲号先于个

11、体号2.构建三角矩阵 个体 的父母未知时：个体 的父或母为 时：个体的父母已知为 或 ，假设 ，这时：tp1ttl0til121ti、t1210215.0tppipillpiti、ppitittfll25.075.0112L1210215.021)(5.0tqqiqppilpilllqiqipiti、)(25.05.05.01112qppiqitiqjpjttffllllpqqp http:/ 3.令 为 对角线元素组成的对角阵，让4.按以下规则加入已知父母的个体的有关元素构建 DL 2)(11DA1A如果双亲已知为 和 ：如果个体父或母已知 为：pq1A要加入的数值要加入的数值中的位置中的位

12、置p1A要加入的数值要加入的数值中的位置中的位置iia5.0),(),(),(),(qiiqpiipiia25.0),(),(),(),(qqpqqpppiia5.0),(),(),(),(qiiqpiipiia25.0),(),(),(),(qqpqqppphttp:/ 如果双亲已知为 和 ：如果个体父或母已知为 ：p要加入的数值要加入的数值中的位置中的位置2-10.5pq1A),(),(),(),(qiiqpiip),(),(),(),(qqpqqppp),(ii要加入的数值要加入的数值中的位置中的位置3/4-2/31/31A),(ii),(),(piip),(pp1Ahttp:/ 育种值

13、估计模型育种值估计模型 v动物模型动物模型 （animal model animal model）数学方程式：期望和方差：混合模型方程组：eZaXby0a)(E0e)(EXby)(E22eaVarI00AeayZyXabAZZXZZXXX1k22221hhkaehttp:/ （sire model sire model）数学方程式：期望和方差：混合模型方程组：eZsXby0s)(E0e)(EXby)(E2e2ssI00AesVar是公畜间加性遗传相关矩阵 sAyZyXsbAZZXZZXXX1sk22222224hhkssysehttp:/ 假设和约束条件：只可用来估计公畜的育种值 公畜在群体中

14、与母畜的交配是完全随机的 母亲之间没有血源关系 每个母亲只有一个后代，即一个公畜的所有后代都是父系半同胞。http:/ （sire-dam model sire-dam model）数学方程式：期望和方差：混合模型方程组：edZsZXbyds0s)(E0d)(E0e)(EXby)(E222edsVarI000A000Aedsds是公畜间加性遗传相关矩阵 sA是母畜间加性遗传相关矩阵 dAyZyZyXdsbAZZZZXZZZAZZXZZXZXXXds1dddsddds1ssssds21kk222214hhkse2222224hhkdehttp:/ 假设和约束条件：只适用于后裔测定的父、母亲育种值

15、预测，而且主要适用于猪、鸡等母畜繁殖力高的畜禽 动物只有一个记录 有记录的动物不是其它动物的双亲 双亲无记录 http:/ （maternal grandsire model maternal grandsire model）数学方程式：期望和方差：混合模型方程组：egZsZXbygs0s)(E0g)(E0e)(EXby)(E222edgsVarI000A000Aegss是公畜间加性遗传相关矩阵 sA是外祖父间加性遗传相关矩阵 gAyZyZyXgsbAZZZZXZZZAZZXZZXZXXXgs1gggsgggs1ssssgs21kk222214hhkse22222)516(hhkgehttp:

16、/ 假设和约束条件：主要适用于种公牛评定 动物只有一个记录 有记录的动物不是其它动物的双亲 双亲无记录 每个母畜只有一个后代，且外祖母只有一个女儿 母畜在外祖父所有女儿中随机抽样，http:/ 多性状多性状BLUP BLUP 法的基本原理法的基本原理v BLUPBLUP原理同样可使用于对多个性状进行育种值估计。原理同样可使用于对多个性状进行育种值估计。v 当我们要对个体在多个性状上的育种值进行估计时，一种当我们要对个体在多个性状上的育种值进行估计时，一种方法可以方法可以分别对每一性状单独进行估计分别对每一性状单独进行估计，然后根据性状之，然后根据性状之间的经济重要性间的经济重要性进行综合进行综

17、合。v 另一种方法可以另一种方法可以利用一个多性状模型对多个性状同时进行利用一个多性状模型对多个性状同时进行估计估计。由于同时进行估计时考虑了性状间的相关，利用了。由于同时进行估计时考虑了性状间的相关，利用了更多的信息，同时可校正由于对某些性状进行了选择而产更多的信息，同时可校正由于对某些性状进行了选择而产生的偏差，因而可提高估计的准确度。生的偏差，因而可提高估计的准确度。http:/ 两性状线性模型两性状线性模型v 合并的矩阵形式：合并的矩阵形式：222222111111euZbXyeuZbXy212121212121eeeuuaZ00ZZbbbX00XXyyy,eZuXby0u)(E0e)

18、(EXby)(E IR00AGeu00Var22211211gggg0G22211211rrrr0Rhttp:/ 令令v 两性状混合模型方程组的简化形式两性状混合模型方程组的简化形式:v 获得综合育种值获得综合育种值 得到各个个体两个性状的估计育种值后，可用性状经济重要性进行加权计算综合育种值，或者将估计育种值转化为标准化的估计育种值，然后再加权计算综合育种值。即：2212121122121211rrrrgggg1010RG，2221221221112221221221112222121222121212111112112212221212111211ryryryryryryryrygrgrr

19、rgrgrrrrrrrrrrr221122112121122112221212111121112212221221112111ZZZZXXXXaabbAZZAZZXZXZAZZAZZXZXZZXZXXXXXZXZXXXXX22112211iiiiiiVEBwVEBwIEBVwEBVwIAEBVVEBhttp:/ BLUPBLUP育种值估计举例育种值估计举例v单性状动物模型单性状动物模型BLUPBLUP育种值估计育种值估计 某种猪场有如下种猪性能测定资料，测定性状为达100 kg日龄，已知该性状的遗传力为0.33，试对该性状资料进行个体育种值估计。种猪达100kg日龄记录猪场猪场个体个体父亲父亲

20、母亲母亲达达100kg日龄日龄111401215213113524121432532160http:/ 个体间加性遗传相关矩阵的计算个体间加性遗传相关矩阵的计算1375.05.05.025.0375.0125.05.05.05.025.0105.05.05.001025.05.05.001A111a122a02112 aa133a5.05.0113113aaa15.013255aa25.0)05.0(5.0)(5.012131551aaaav 个体间加性遗传相关矩阵逆矩阵的计算个体间加性遗传相关矩阵逆矩阵的计算 构建 ：http:/ 个体间加性遗传相关矩阵逆矩阵的计算个体间加性遗传相关矩阵逆矩

21、阵的计算 构建对角矩阵：令：http:/ 构建线性模型构建线性模型 根据资料性质，可对种猪达100kg日龄写出如下动物模型：用矩阵形式表示，则对于该资料有：ijjiijeahy2524131211543212110000010000010000010000011010010101160143135152140eeeeeaaaaahhyhttp:/ 构建混合模型方程组构建混合模型方程组 因此有：2003XX1100000111ZX)ZX(XZ1000001000001000001000001ZZ303427yX160143135152140yZ0003.233.033.01122hhk16014

22、313515214030342740220100402210206667.413333.1012215101023333.116667.401110002000111035432121aaaaahhv求解混合模型方程组有：求解混合模型方程组有：http:/ kg日龄和达100kg背膘，试以两个性状资料进行个体育种值估计。种猪达100kg日龄和达100kg背膘厚测定记录猪场猪场个体个体父亲父亲母亲母亲达达100kg日龄日龄(d)达达100kg背膘厚背膘厚（mm）1114013121521413113512241214313253216016http:/ 根据资料性质，可对种猪达100kg日龄和达

23、100kg背膘厚写出如下动物模型：是第 性状，第 猪场，第 个体的观测值 是第 性状，第 猪场的效应 是第 性状，第 个体的育种值 是随机残差 猪两个性状的表型、遗传参数和经济加权值（表中右边2项的右上角为表型相关，左下角为遗传相关）性状性状单位单位达达100kg日龄（日龄（）d-0.60.332250.55达达100kg背膘厚（背膘厚（）mm-0.80.501.440.451X2X2Xw2h2P1Xijkikijijkeahyijkyijkijhijikaikijkev 因为所有个体两个性状都有记录，因此有因为所有个体两个性状都有记录，因此有http:/ 加性遗传相关矩阵的逆矩阵（单性状例子

24、获得）：加性遗传相关矩阵的逆矩阵（单性状例子获得）：201100201110611213211212210132216111Ahttp:/ 由遗传参数表可计算出性状间的遗传和误差方差及协方差由遗传参数表可计算出性状间的遗传和误差方差及协方差为：为：v 逆矩阵为：逆矩阵为：7200.0002500.1467200.03885.33885.37500.7800RG，1.3889000.00681.74160.0749-0.0749-0.01591100RG，http:/ 混合模型方程组：混合模型方程组：18.055718.055716.666819.444618.05571.08800.97240

25、.91801.03360.952036.111454.16712.06042.90364.87210000.01.7416-1.7416-0000.0 0.1498-0000.00.07490.07490000.01.38890000.00000.00000.00.00004.87210000.01.7416-1.7416-0000.0 0.1498-0000.00.07490.07491.38890000.00000.00000.01.7416-0000.04.58180.87081.1611-0.07490000.00.1373-0.0374-0.04990000.01.38890000.

26、00000.01.7416-1.7416-0.87084.87210.87080.07490.07490.0374-0.1498-0.0374-0000.01.38890000.00000.00000.01.7416-1.1611-0.87084.58180000.00.07490.04990.0374-0.1373-0000.01.38890000.00000.00.1498-0000.00.07490.07490000.00.03860000.00.0159-0.0159-0000.00000.00000.00.00680000.00.0000 0.1498-0000.00.07490.0

27、7490000.00.03860000.00.0159-0.0159-0000.00000.00.00680000.00.07490.00000.1373-0.0374-0.04990.0159-0000.00.03590.00800.0106-0000.00000.00000.00.00680.07490.07490.0374-0.1498-0.0374-0.0159-0.0159-0.00800.03860.00800000.00000.00000.00.00680000.00.07490.04990.0374-0.1373-0000.00.0159-0.0106-0.00800.0359

28、0000.00000.00000.00.00681.38891.38890000.00000.00000.00000.00000.00000.00000.00000.02.77780000.00000.00000.00000.00000.01.38891.38891.38890000.00000.00000.00000.00000.00000.04.16670000.00000.00000.00000.00000.00000.00000.00.0068 0.0068 0000.00000.00000.00000.00000.00.01360000.00000.00000.00000.00000

29、.00000.00000.00000.00.0068 0.0068 0.0068 0000.00000.00000.00.02042524232221151413121122211211aaaaaaaaaahhhhhttp:/ 方程组的解为：方程组的解为：v 综合育种值为：综合育种值为：12.936913.0859150.9676143.080722211211hhhh2.15441.0895-3.7480-4.37172.8657-1514131211aaaaa0.06390.06230.5402-0.51020.2278-2524232221aaaaa1.3438-0.60392.6810

30、 3.0312-1.901754321aaaaahttp:/ 如果说模型是BLUP法的关键，那么，计算问题则是BLUP法的难点。从前面例子我们可以看出，仅仅5个个体两个性状就产生很大的方程组，而对于猪、鸡等畜禽在BLUP法中所涉及的线性方程组是非常大的，对一些跨群（场）的遗传评定，方程组个数可达几万至几十万，如此大数量的方程组用手工计算是根本不可能的。v 近年来，世界各国育种学家在BLUP法的计算问题上做了大量的工作，已开发出相应的电脑软件，如国外的PEST 和国内的NETPIG。下面对一些常用的遗传评估软件进行介绍。第七节第七节 BLUP育种值估计软件育种值估计软件http:/ PEST P

31、EST是由美国Illinois大学的Groeneveld、Kovac和Wang（1990）开发研制的多性状遗传评估软件，其英文全文名为（Multivariate Prediction and ESTimation，目前已在世界各国广泛应用。根据性能测定和生产数据，PEST提供了基于30多种数学模型的单性状或多性状BLUP育种值的计算，包括固定模型、个体动物模型、公畜模型、公畜母畜模型和外祖父模型等 。为了满足实际育种的需要，系统还提供了可自行定义性状、修改模型和设定参数的余地。PEST可以在不同的操作系统下运行。http:/ PIGBLUPPIGBLUP PIGBLUP软件是由澳大利亚New

32、England大学编写的，在1989发行了第一版，从那时起版本不断更新，到目前为止，最新的版本5.10已发展成为WINDOWS界面操作 PIGBLUP是一种专为育种猪场设计使用的现代遗传评估系统，PIGBLUP主要包含种猪评估、遗传进展分析、选配和遗传审计四个模块 目前，PIGBLUP已国际化，正在多个国家使用 v GBSGBS GBS是“猪场生产管理与育种数据分析系统”的英文缩写，是中国农业大学动物科学技术学院GBS软件创作小组开发研制的系统软件。GBS是中文WINDOWS95/WINDOWNT下的管理信息系统。它集种猪、商品猪生产和育种数据的采集与分析于一体，十分适合大型种猪生产集团使用，

33、并支持联合育种方案 http:/ NETPIGNETPIG NETPIG（种猪场网络管理系统）是四川农业大学动物科技学院和重庆市养猪科学研究院（系统指导：李学伟、王金勇；程序设计：徐顺来）联合研制开发WINGDOWS界面的种猪场网络管理系统 该系统借鉴了加拿大、丹麦等国的成功经验，应用先进的数学模型进行育种值估计，非常易于实现“联合育种”系统主要包含“生产管理”和“育种管理”、两大模块，模块之间相辅相成，数据共享，完全无缝街接 该软件有两个版本，一个用于遗传评估中心，遗传评估中心通过该模块对整个地区的所有种猪进行统一遗传评估并将结果在网上发布。另一个在猪场，各种猪场用于进行场内遗传评估，以解决场内选种与遗传评估中心发布遗传评估结果的时间差问题，为选种作一定参考 Add your company slogan