文章目录
数据分析
对每个杂种优势模式分别进行数据分析,以便很好地确定估计遗传方差和确定杂种之间协方差所必需的概念基础群体(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
X,Z0,Z1,Z2和Z分别是与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′)]=fAA′VGCA(1)+fBB′VGCA(2)+fAA′fBB′VSCA(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}
Var⎣⎢⎢⎢⎢⎡cg1g2de⎦⎥⎥⎥⎥⎤=⎣⎢⎢⎢⎢⎡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'}
fAA′fBB′。假定对照杂种不相关,
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
β、c、g1、g2和d的解是通过求解以下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⎦⎥⎥⎥⎥⎤=[x′R1XZ4′R−1Xx′R1Z0x−1R1Xx′R1Z1x−1R1Xx′R1Z2x−1R1Xx′R1Z]