06.06.2013 Views

Theory of Statistics - George Mason University

Theory of Statistics - George Mason University

Theory of Statistics - George Mason University

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

We now write the original model as<br />

Now, for the transformations. Let<br />

Y = δ1 T n + ɛ.<br />

Z = HmXH T n ,<br />

5.5 Applications 433<br />

δ = Hmδ,<br />

and<br />

ɛ = HmɛH T n .<br />

We first <strong>of</strong> all note that the transformations are all nonsingular and<br />

Z = H1 T n + ɛ.<br />

Next, we see because <strong>of</strong> the orthonormality <strong>of</strong> the Helmert matrices that<br />

the distributions <strong>of</strong> δ and ɛ are the same as those <strong>of</strong> δ and ɛ and they are still<br />

independent. Furthermore, the Zij are independent, and we have<br />

and<br />

Zi1<br />

iid ∼ N(0, σ 2 a + σ 2 ), for i = 1, . . ., m<br />

Zij iid ∼ N(0, σ 2 ), for i = 1, . . ., m; j = 2, . . ., n0.<br />

To continue with the analysis, we follow the same steps as in Example 5.28,<br />

and get the same decomposition <strong>of</strong> the “adjusted total sum <strong>of</strong> squares” as in<br />

equation (5.94):<br />

m n<br />

i=1 j=1<br />

(Zij − Z) 2 = SSA + SSE. (5.99)<br />

Again, we get chi-squared distributions, but the distribution involving SSA<br />

is not the same as in expression (5.95) for the fixed-effects model.<br />

Forming<br />

MSA = SSA/(m − 1)<br />

and<br />

we see that<br />

and<br />

MSE = SSE/(m(n − 1)),<br />

E(MSA) = nσ 2 δ + σ 2 ɛ<br />

E(MSE) = σ 2 ɛ.<br />

Unbiased estimators <strong>of</strong> σ 2 δ and σ2 ɛ are therefore<br />

<strong>Theory</strong> <strong>of</strong> <strong>Statistics</strong> c○2000–2013 James E. Gentle<br />

s 2 δ = (MSA − MSE)/n (5.100)

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

Saved successfully!

Ooh no, something went wrong!