17.01.2015 Views

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

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!