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.

18 EXERCISES<br />

1. Below are data on progeny of 6 rams used in 5 sheep flocks (for some trait). The rams<br />

were unrelated to each other and to any of the ewes to which they were mated. The first<br />

number is the number of progeny in the herd, and the second (within parentheses) is the<br />

sum of the observations.<br />

Let the model equation be<br />

Ram<br />

Flocks<br />

ID 1 2 3 4 5<br />

1 6(638) 8(611) 6(546) 5(472) 0(0)<br />

2 5(497) 5(405) 5(510) 0(0) 4(378)<br />

3 15(1641) 6(598) 5(614) 6(639) 5(443)<br />

4 6(871) 11(1355) 0(0) 3(412) 3(367)<br />

5 2(235) 4(414) 8(874) 4(454) 6(830)<br />

6 0(0) 0(0) 4(460) 12(1312) 5(558)<br />

y ijk = µ + F i + R j + e ijk<br />

where F i is a flock effect, R j is a ram effect, and e ijk is a residual effect. There are a total<br />

of 149 observations and the total sum of squares was equal to 1,793,791. Assume that<br />

σ 2 e = 7σ 2 f = 1.5σ2 r when doing the problems below.<br />

(a) Set up the mixed model equations and solve. Calculate the SEPs and reliabilities of<br />

the ram solutions.<br />

(b) Repeat the above analysis, but assume that flocks are a fixed factor (i.e. do not add<br />

any variance ratio to the diagonals of the flock equations). How do the evaluations,<br />

SEP, and reliabilities change from the previous model<br />

(c) Assume that rams are a fixed factor, and flocks are random.<br />

similarly to the previous two models<br />

Do the rams rank<br />

2. When a model has one random factor and its covariance matrix is an identity matrix times<br />

a scalar constant, then prove the the solutions for that factor from the MME will sum to<br />

zero. Try to make the proof as general as possible.<br />

26

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

Saved successfully!

Ooh no, something went wrong!