Section 3.4

Normal Linear Models

[21/0%]
Show
LaTeX
Problem 3.4.1
?
(a)

Suppose Xi:N(ξi,σ2)X_{i}: N\left(\xi_{i}, \sigma^{2}\right) with ξi=α+βti\xi_{i}=\alpha +\beta t_{i}. If the first column of the matrix CC leading to the canonical form (4.7) is (1/n,…,1/n)′(1 / \sqrt{n}, \ldots , 1 / \sqrt{n})^{\prime }, find the second column of CC.

(b)

If Xi:N(ξi,σ2)X_{i}: N\left(\xi_{i}, \sigma^{2}\right) with ξi,=α+βti+γti2\xi_{i},=\alpha +\beta t_{i}+\gamma t_{i}^{2}, and the first two columns of CC are those of (a), find the third column under the simplifying assumptions Σti=0,Σti2=1\Sigma t_{i}=0, \Sigma t_{i}^{2}=1.

Remark.
?
(b)

The orthogonal polynomials that are progressively built up in this way are frequently used to simplify regression analysis.

Problem 3.4.2

Write out explicit expressions for the transformations (4.10) when ΠΩ\Pi_{\Omega } is given by

(a) ξi=α+βti\xi_{i}=\alpha +\beta t_{i} and (b) ξi=α+βti+γti2\xi_{i}=\alpha +\beta t_{i}+\gamma t_{i}^{2}.

?
Problem 3.4.3

Use Problem 3.10 to prove (iii) of Theorem 4.3.

?
Problem 3.4.4
?
(a)

In Example 4.7, determine α^,β^\hat{\alpha }, \hat{\beta }, and hence ξ^i\hat{\xi }_{i} by minimizing Σ(Xi−α−βti)2\Sigma \left(X_{i}-\alpha -\beta t_{i}\right)^{2}.

(b)

Verify the expressions (4.12) for α\alpha and β\beta, and the corresponding expressions for α^\hat{\alpha } and β^\hat{\beta }.

Problem 3.4.5

In Example 4.2, find the UMVU estimators of α,β,γ\alpha , \beta , \gamma, and σ2\sigma^{2} when Σti=0\Sigma t_{i}=0 and Σti2=1\Sigma t_{i}^{2}=1.

?
Problem 3.4.6

Let XijX_{i j} be independent N(ξij,σ2)N\left(\xi_{i j}, \sigma^{2}\right) with ξij=αi+βtij\xi_{i j}=\alpha_{i}+\beta t_{i j}. Find the UMVU estimators of the αi\alpha_{i} and β\beta.

?
Problem 3.4.7
?
(a)

In Example 4.9, show that the vectors of the coefficients in the α^i\hat{\alpha }_{i} are not orthogonal to the vector of the coefficients of μ^\hat{\mu }.

(b)

Show that the conclusion of (a) is reversed if α^i\hat{\alpha }_{i} and μ^\hat{\mu } are replaced by α^^i\hat{\hat{\alpha }}_{i} and μ^^\hat{\hat{\mu }}.

Problem 3.4.8

In Example 4.9, find the UMVU estimator of μ\mu when the αi\alpha_{i} are known to be zero and compare it with μ^\hat{\mu }.

?
Problem 3.4.9

The coefficient vectors of the XijkX_{i j k} given by (4.32) for μ^,α^i\hat{\mu }, \hat{\alpha }_{i}, and β^j\hat{\beta }_{j} are orthogonal to the coefficient vectors for the γ^ij\hat{\gamma }_{i j} given by (4.33).

?
Problem 3.4.10

In the model defined by (4.26) and (4.27), determine the UMVU estimators of αi\alpha_{i}, βj\beta_{j}, and σ2\sigma^{2} under the assumption that the γij\gamma_{i j} are known to zero.

?
Problem 3.4.11
?
(a)

In Example 4.11, show that

ΣΣΣ(Xijk−μ−αi−βj−γij)2=S2+Sμ2+Sα2+Sβ2+Sγ2 \Sigma \Sigma \Sigma \left(X_{i j k}-\mu -\alpha _{i}-\beta _{j}-\gamma _{i j}\right)^{2}=S^{2}+S_{\mu }^{2}+S_{\alpha }^{2}+S_{\beta }^{2}+S_{\gamma }^{2}

where S2=ΣΣΣ(Xijk−Xij.)2,Sμ2=IJm⁡(X…−μ)2,Sα2=Jm⁡Σ(X1..−X…−αi)2S^{2}=\Sigma \Sigma \Sigma \left(X_{i j k}-X_{i j.}\right)^{2}, S_{\mu }^{2}=\operatorname {IJm}\left(X_{\ldots }-\mu \right)^{2}, S_{\alpha }^{2}=\operatorname {Jm} \Sigma \left(X_{1. .}-X_{\ldots }-\alpha_{i}\right)^{2}, and Sβ2,S2γS_{\beta }^{2}, S^{2} \gamma are defined analogously.

(b)

Use the decomposition of (a) to show that the least squares estimators of μ,αi,…\mu , \alpha_{i}, \ldots are given by (4.32) and (4.33).

(c)

Show that the error sum of squares S2S^{2} is equal to ΣΣΣ(Xijk−ξ^ij)2\Sigma \Sigma \Sigma \left(X_{i j k}-\hat{\xi }_{i j}\right)^{2} and hence in the canonical form to Σj=s+1nYj2\Sigma_{j=s+1}^{n} Y_{j}^{2}.

Problem 3.4.12
?
(a)

Show how the decomposition in Problem 4.11(a) must be modified when it is known that the γij\gamma_{i j} are zero.

(b)

Use the decomposition of (a) to solve Problem 4.10.

Problem 3.4.13

Let Xijk(i=1,…,I,j=1,…,J,k=1,…,K)X_{i j k}(i=1, \ldots , I, j=1, \ldots , J, k=1, \ldots , K) be N(ξijk,σ2)N\left(\xi_{i j k}, \sigma^{2}\right) with

ξijk=μ+αi+βj+γk \xi _{i j k}=\mu +\alpha _{i}+\beta _{j}+\gamma _{k}

where Σαi=Σβj=Σγk=0\Sigma \alpha_{i}=\Sigma \beta_{j}=\Sigma \gamma_{k}=0. Express μ,αi,βj\mu , \alpha_{i}, \beta_{j}, and γk\gamma_{k} in terms of the ξ\xi 's and find their UMVU estimators. Viewed as a special case of (4.4), what is the value of ss?

?
Problem 3.4.14

Extend the results of the preceding problem to the model

ξijk=μ+αi+βj+γk+δij+εik+λjk \xi _{i j k}=\mu +\alpha _{i}+\beta _{j}+\gamma _{k}+\delta _{i j}+\varepsilon _{i k}+\lambda _{j k}

where

∑iδij=∑jδij=∑iεik=∑kεik=∑jλjk=∑kλjk=0. \sum _{i} \delta _{i j}=\sum _{j} \delta _{i j}=\sum _{i} \varepsilon _{i k}=\sum _{k} \varepsilon _{i k}=\sum _{j} \lambda _{j k}=\sum _{k} \lambda _{j k}=0.
?
Problem 3.4.15

In the preceding problem, if it is known that the λ\lambda 's are zero, determine whether the UMVU estimators of the remaining parameters remain unchanged.

?
Problem 3.4.16
?
(a)

Show that under assumptions (4.35), if ξ=θA\xi =\theta A, then the least squares estimate of θ\theta is xA(AA′)−1\mathbf{x} A\left(A A^{\prime }\right)^{-1}.

(b)

If (X,A)(X, A) is multivariate normal with all parameters unknown, show that the least squares estimator of part (a) is a function of the complete sufficient statistic and, hence, prove part (a) of Theorem 4.14.

Problem 3.4.17

A generalization of the order statistics, to vectors, is given by the following definition.

Definition 8.1 The cjc_{j}-order statistics of a sample of vectors are the vectors arranged in increasing order according to their jj th components.

Let Xi,i=1,…,n\mathbf{X}_{i}, i=1, \ldots , n, be an iid sample of p×1p \times 1 vectors, and let X=(X1,…,Xn)X=\left(X_{1}, \ldots , X_{n}\right) be a p×np \times n matrix.

Theorem 4.14 states: under assumptions (4.35), with ξ=θA\xi =\theta A, the following hold.

?
(a)

If (X,A)(X, A) are jointly multivariate normal with all parameters unknown, then the least squares estimator γiξ^i\gamma_{i}\hat{\xi }_{i} is the UMVU estimator of γiξi\gamma_{i}\xi_{i}.

(b)

If the distribution of AA is unknown, then the least squares estimator γiξ^i\gamma_{i}\hat{\xi }_{i} is UMVU among all linear estimators of γiξi\gamma_{i}\xi_{i}.

(c)

If E[AA′]\mathbb {E}\left[AA'\right] is known, no best linear unbiased estimator of γiξi\gamma_{i}\xi_{i} exists.

  1. If the distribution of Xi\mathbf{X}_{i} is completely unknown, show that, for any j,j=1,…,pj, j=1, \ldots , p, the cjc_{j}-order statistics of (X1,…,Xn)(X_{1}, \ldots , X_{n}) are complete sufficient. (That is, the vectors X1,…,XnX_{1}, \ldots , X_{n} are ordered according to their jj th coordinate.)

  2. Let Y1×nY_{1 \times n} be a random variable with unknown distribution (possibly different from Xi\mathbf{X}_{i}). Form the (p−1)×n(p-1) \times n matrix (xy)\binom {\mathbf{x}}{\mathbf{y}}, and for any j=1,…,pj=1, \ldots , p, calculate the cjc_{j}-order statistics based on the columns of (xy)\binom {\mathbf{x}}{\mathbf{y}}. Show that these cjc_{j}-order statistics are sufficient.

  3. Use parts (a) and (b) to prove Theorem 4.14(b).

Problem 3.4.18

The proof of Theorem 4.14(c) is based on two results. Establish that:

?
(a)

For large values of θ\theta, the unconditional variance of a linear unbiased estimator will be greater than that of the least squares estimator.

(b)

For θ=0\theta =0, the variance of XA(AA′)−1X A\left(A A^{\prime }\right)^{-1} is greater than that of XA[E[AA′]]−1X A\left[\mathbb {E}\left[A A^{\prime }\right]\right]^{-1}.

(c)

Parts (a) and (b) imply that no best linear unbiased estimator of Σγiξi\Sigma \gamma_{i} \xi_{i} exists if E[AA′]\mathbb {E}\left[A A^{\prime }\right] is known.

Problem 3.4.19
?
(a)

Under the assumptions of Example 4.15, find the variance of ΣλiSi2\Sigma \lambda_{i} S_{i}^{2}.

(b)

Show that the variance of (a) is minimized by the values stated in the example.

Problem 3.4.20

In the linear model (4.4), a function Σciξi\Sigma c_{i} \xi_{i} with Σci=0\Sigma c_{i}=0 is called a contrast. Show that a linear function Σdiξi\Sigma d_{i} \xi_{i} is a contrast if and only if it is translation invariant, that is, satisfies Σdi(ξi+a)=Σdiξi\Sigma d_{i}\left(\xi_{i}+a\right)=\Sigma d_{i} \xi_{i} for all aa, and hence if and only if it is a function of the differences ξi−ξj\xi_{i}-\xi_{j}.

?
Problem 3.4.21

Determine which of the following are contrasts:

?
(a)

The regression coefficients α,β\alpha , \beta, or γ\gamma of (4.2).

(b)

The parameters μ,αi,βj\mu , \alpha_{i}, \beta_{j}, or γij\gamma_{i j} of (4.27).

(c)

The parameters μ\mu or αi\alpha_{i} of (4.23) and (4.24).