Best Linear Unbiased Prediction of Maize Single-Cross Performance

文章目录

数据分析

对每个杂种优势模式分别进行数据分析,以便很好地确定估计遗传方差和确定杂种之间协方差所必需的概念基础群体(Stuber和Cockerham,1966)。
假设单交由1组的1个自交系和2组的2个自交系组成。在不同的产量试验中,对单交种和nc检验杂种进行了评价,得出p个总数据点。在单交叉预测的混合模型方程中,假设可忽略上位性的线性模型是:
y = X β + Z 0 c + Z 1 g 1 + Z 2 g 2 + Z d + e y=X\beta+Z_0c+Z_1g_1+Z_2g_2+Zd+e y=Xβ+Z0c+Z1g1+Z2g2+Zd+e
其中 y y y是给定性状的观测值向量, β \beta β是产量试验效应, c c c是对照杂交种效应, g 1 g_1 g1是组合1自交系的一般配合力(GCA)效应, g 2 g_2 g2组合2自交系的一般配合力效应, d d d是特殊配合力效应, e e e是残差; X , Z 0 , Z 1 , Z 2 和 Z X,Z_0,Z_1,Z_2和Z XZ0Z1Z2Z分别是与y相关的1和0的关联矩阵。多地点产量试验效应被认为是固定的,而模型中的所有其他效应被认为是随机的。因为数据集包括单个杂交的平均表现和不同地点的杂交检验,虽然个别地点的影响被认为是随机的,但它们不在模型中。产量试验效应被认为是固定的,因此在方差分量估计和单交叉预测之前,每个数据点都可以调整其产量试验效应的细节。
假设 A A A A ′ A' A是组1中的两个近交系, B B B B ′ B' B是组合2中的两个近交系。在 A × B A \times B A×B A ′ × B ′ A' \times B' A×B的协变量矩阵是:
C o v [ ( A × B ) , ( A ′ × B ′ ) ] = f A A ′ V G C A ( 1 ) + f B B ′ V G C A ( 2 ) + f A A ′ f B B ′ V S C A ( 1 ) Cov[(A \times B),(A' \times B')]=f_{AA'}V_{GCA_(1)}+f_{BB'}V_{GCA_(2)}+f_{AA'}f_{BB'}V_{SCA} \qquad(1) Cov[(A×B),(A×B)]=fAAVGCA(1)+fBBVGCA(2)+fAAfBBVSCA(1)

其中 V G C A ( 1 ) V_{GCA_(1)} VGCA(1)是组合1的等位基因的一般配合力方差和, V G C A ( 2 ) V_{GCA_(2)} VGCA(2)是组合2的等位基因的一般配合力方差和, V S C A V_{SCA} VSCA是组合1和组合2的等位基因的特殊配合力方差和; f A A ′ f_{AA'} fAA是近交系 A A A A ′ A' A之间的同源系数, f B B ′ f_{BB'} fBB是近交系 B B B B ′ B' B之间的同源系数.随机效应 c , g 1 , g 2 , d , e c,g_1,g_2,d,e c,g1,g2,d,e的均值为0,遵循以下协方差:
V a r [ c g 1 g 2 d e ] = [ I V c 0 0 0 0 0 G 1 V G C A ( 1 ) 0 0 0 0 0 G 2 V G C A ( 2 ) 0 0 0 0 0 D V S C A 0 0 0 0 0 R V R ] V_{ar} \begin{bmatrix} c \\ g_1 \\ g_2 \\ d \\ e \\ \end{bmatrix}=\begin{bmatrix} IV_c &0 &0 &0 &0 \\ 0& G_1V_{GCA(1)}&0 &0 &0\\ 0 & 0 & G_2V_{GCA(2)} &0 &0 \\ 0 & 0 & 0 & DV_{SCA} &0\\ 0 & 0& 0& 0 & RV_R \\ \end{bmatrix} Varcg1g2de=IVc00000G1VGCA(1)00000G2VGCA(2)00000DVSCA00000RVR
其中 V c V_c Vc是对照杂交种之间的方差, V R V_R VR是残差。

矩阵元素分别等于 G 1 G_1 G1中的 f A A ′ f_{AA'} fAA G 2 G_2 G2中的 f B B ′ f_{BB'} fBB和D中的 f A A ′ f B B ′ f_{AA'}f_{BB'} fAAfBB。假定对照杂种不相关, c c c的协方差为 I V C IV_C IVC,其中 i i i为单位矩阵。如果一些对照杂交种的关系数据可用,则可以使用适当的遗传关系矩阵代替I进行检验。在R中,假定非对角线元素为零,并且第i个对角元素是y中第i个数据点的数目或f个观测(位置)的倒数。例如,如果y中的第一个元素是在四个位置的试验中单个杂交的平均表现,则R中的第一个对角元素是1/4。
β 、 c 、 g 1 、 g 2 和 d \beta、c、g_1、g_2和d βcg1g2d的解是通过求解以下MME得到的:
[ β c g 1 g 2 d ] = [ x ′ R 1 X x ′ R 1 Z 0 x ′ R 1 Z 1 x ′ R 1 Z 2 x ′ R 1 Z Z 4 ′ R − 1 X x − 1 R 1 X x − 1 R 1 X x − 1 R 1 X ] \begin{bmatrix} \beta \\ c \\ g_1 \\ g_2 \\ d \\ \end{bmatrix} =\begin{bmatrix} x' R^1X & x'R^1Z_0 &x'R^1Z_1 & x'R^1Z_2 &x'R^1Z \\ Z_4'R^{-1}X & x^{-1} R^1X & x^{-1} R^1X & x^{-1} R^1X \\ \end{bmatrix} βcg1g2d=[xR1XZ4R1XxR1Z0x1R1XxR1Z1x1R1XxR1Z2x1R1XxR1Z]

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值