Prediction Theory 1 Introduction 2 General Linear Mixed Model
Prediction Theory 1 Introduction 2 General Linear Mixed Model
Prediction Theory 1 Introduction 2 General Linear Mixed Model
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
17.2.3 G and R<br />
The covariance matrix of the means of residuals is R. The variance of a mean of random<br />
variables is the variance of individual variables divided by the number of variables in the mean.<br />
Let n ij equal the number of progeny in a sire by contemporary group subclass, then the variance<br />
of the subclass mean is σe/n 2 ij . Thus,<br />
R =<br />
⎛<br />
⎜<br />
⎝<br />
σ 2 e/3 0 0 0 0<br />
0 σ 2 e/4 0 0 0<br />
0 0 σ 2 e/5 0 0<br />
0 0 0 σ 2 e/6 0<br />
0 0 0 0 σ 2 e/3<br />
⎞<br />
.<br />
⎟<br />
⎠<br />
The matrix G is similarly partitioned into two submatrices, one for contemporary groups<br />
and one for sires.<br />
( )<br />
Gc 0<br />
G =<br />
,<br />
0 G s<br />
where<br />
and<br />
G c =<br />
G s =<br />
⎛<br />
⎜<br />
⎝<br />
(<br />
σ<br />
2<br />
c 0<br />
0 σ 2 c<br />
)<br />
σ 2 s 0 0<br />
0 σ 2 s 0<br />
0 0 σ 2 s<br />
⎞<br />
= Iσ 2 c = I σ2 e<br />
6.0 ,<br />
⎟<br />
⎠ = Iσ 2 s<br />
= I σ2 e<br />
11.5 .<br />
and<br />
The inverses of G and R are needed for the MME.<br />
⎛<br />
⎞<br />
3 0 0 0 0<br />
0 4 0 0 0<br />
R −1 1<br />
=<br />
0 0 5 0 0<br />
⎜<br />
⎟ σ<br />
⎝ 0 0 0 6 0 ⎠<br />
2 ,<br />
e<br />
0 0 0 0 3<br />
G −1 =<br />
⎛<br />
⎜<br />
⎝<br />
6 0 0 0 0<br />
0 6 0 0 0<br />
0 0 11.5 0 0<br />
0 0 0 11.5 0<br />
0 0 0 0 11.5<br />
⎞<br />
1<br />
⎟ σ<br />
⎠ e<br />
2 .<br />
Because both are expressed in terms of the inverse of σ 2 e, then that constant can be ignored.<br />
The relative values between G and R are sufficient to get solutions to the MME.<br />
21