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.

The degrees of freedom for F are r(H ′ o) and (N − r(X)). Note that<br />

y ′ R −1 y − ˆb ′ X ′ R −1 y − û ′ Z ′ R −1 y = y ′ V −1 y − ˆb ′ X ′ V −1 y.<br />

If G and R are not known, then there is no best test because BLUE of b is not possible.<br />

Valid tests exist only under certain circumstances. If estimates of G and R are used to construct<br />

the MME, then the solution for ˆb is not BLUE and the resulting tests are only approximate.<br />

If the estimate of G is considered to be inappropriate, then a test of H ′ ob = c can be<br />

constructed by treating u as a fixed factor, assuming that H ′ ob is estimable in the model with<br />

u as fixed. That is,<br />

( ) ( ) ˆb X<br />

=<br />

′ R −1 X X ′ R −1 − ( )<br />

Z X ′ R −1 y<br />

û Z ′ R −1 X Z ′ R −1 Z Z ′ R −1 ,<br />

y<br />

=<br />

P zx P zz Z ′ R −1 y<br />

(<br />

Pxx P xz<br />

) (<br />

X ′ R −1 y<br />

)<br />

,<br />

and<br />

(<br />

ˆσ e 2 = (y ′ R −1 y − ˆb ′ X ′ R −1 y − û ′ Z ′ R −1 y)/(N − r<br />

X<br />

Z<br />

)<br />

),<br />

s = (H ′ oˆb − c) ′ (H ′ oP xx H o ) −1 (H ′ oˆb − c),<br />

F = (s/r(H ′ o))/ˆσ 2 e.<br />

11 Restrictions on Fixed Effects<br />

There may be functions of b that are known and this knowledge should be incorporated into the<br />

estimation process. For example, in beef cattle, male calves of a particular breed are known to<br />

weigh 25 kg more than female calves of the same breed at 200 days of age. By incorporating a<br />

difference of 25 kg between the sexes in an analysis then all other estimates of fixed and random<br />

effects would be changed accordingly and also their variances.<br />

Let B ′ b = d be the restriction to be placed on b, then the appropriate equations would be<br />

⎛<br />

X ′ R −1 X X ′ R −1 ⎞ ⎛ ⎞ ⎛<br />

Z B ˆb X ′ R −1 ⎞<br />

y<br />

⎜<br />

⎝ Z ′ R −1 X Z ′ R −1 Z + G −1 ⎟ ⎜ ⎟ ⎜<br />

0 ⎠ ⎝ û ⎠ = ⎝ Z ′ R −1 ⎟<br />

y ⎠ .<br />

B ′ 0 0 φ<br />

d<br />

Because B ′ b = d is any general function, then there are three possible effects of this function<br />

on the estimability of K ′ b in the model. The conditions on B ′ are that it<br />

1. must have full row rank, and<br />

2. must not have more than r(X) rows.<br />

13

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

Saved successfully!

Ooh no, something went wrong!