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.
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