1.10 Practice Problems

Problem 1.1. Let \(A\) be partitioned as \(A = \begin {pmatrix} A_{11} & A_{12}\\ A_{21} & A_{22}\end {pmatrix}\) with \(A_{11}\) non-singular. Show that \[\left |A\right | = \left |A_{11}\right |\, \left |A_{22}-A_{21}A_{11}^{-1}A_{12}\right |.\] Where to start: write \(A = BD\) with \(B\) block lower triangular and \(D\) block upper triangular, and use \(|BD| = |B||D|\). The second factor is the Schur complement, which reappears as a conditional covariance in Section 5.1.

Problem 1.2. Prove that \(\operatorname {tr}(AB) = \operatorname {tr}(BA)\) whenever both products are defined, and deduce that \(\operatorname {tr}(P\Lambda P') = \operatorname {tr}(\Lambda )\) for orthogonal \(P\). This is the identity behind Theorem 4.6.

Problem 1.3. The centring matrix is \(H = I - \tfrac {1}{n}\underline {1}\,\underline {1}'\). Show that \(H\) is symmetric and idempotent, that \(H\underline {1}=\underline {0}\), and that \(\operatorname {tr}(H) = n-1\). Explain what each of these three facts says about what \(H\) does to a data matrix.

Show solution

Solution. Symmetry is immediate since \(\underline {1}\,\underline {1}'\) is symmetric. For idempotency, \[H^{2} = I - \tfrac {2}{n}\underline {1}\underline {1}' + \tfrac {1}{n^{2}}\underline {1}\left (\underline {1}'\underline {1}\right )\underline {1}' = I - \tfrac {2}{n}\underline {1}\underline {1}' + \tfrac {1}{n}\underline {1}\underline {1}' = H,\] using \(\underline {1}'\underline {1}=n\). Then \(H\underline {1} = \underline {1}-\tfrac {1}{n}\underline {1}(n) = \underline {0}\), and \(\operatorname {tr}(H) = n - \tfrac {1}{n}\operatorname {tr}(\underline {1}\underline {1}') = n-1\).

Reading them in order: \(H\) subtracts the column means; applying it twice does nothing more, because once centred a column is already centred; a constant column is annihilated entirely; and one degree of freedom is spent on the mean, which is the \(n-1\) that appears in the divisor of \(S\).

Problem 1.4. Show that a symmetric matrix \(A\) is positive definite if and only if every eigenvalue is strictly positive, and that \(\left |A\right | = \prod \lambda _i\).

Problem 1.5. Obtain the spectral decomposition of \(S = \begin {pmatrix} 4 & 1\\ 1 & 4\end {pmatrix}\), and verify both \(\left |S\right | = \lambda _1\lambda _2\) and \(\operatorname {tr}(S) = \lambda _1+\lambda _2\).

Show solution

Solution. \(\left |S-\lambda I\right | = (4-\lambda )^{2}-1 = 0\) gives \(\lambda _1=5\) and \(\lambda _2=3\), with normalised eigenvectors \((1,1)'/\sqrt 2\) and \((1,-1)'/\sqrt 2\). Hence \[S = 5\,\tfrac 12\begin {pmatrix}1&1\\1&1\end {pmatrix} + 3\,\tfrac 12\begin {pmatrix}1&-1\\-1&1\end {pmatrix}.\] Then \(\lambda _1\lambda _2 = 15 = 16-1 = \left |S\right |\) and \(\lambda _1+\lambda _2 = 8 = 4+4 = \operatorname {tr}(S)\).

Problem 1.6. For \(\Sigma \) positive definite, define the symmetric square root \(\Sigma ^{1/2} = P\Lambda ^{1/2}P'\) and verify that \(\Sigma ^{1/2}\Sigma ^{1/2}=\Sigma \) and that \(\Sigma ^{1/2}\) is symmetric. Where is this matrix used later in the course?

Problem 1.7. Show that if \(A\) is \(n\times p\) with \(n<p\) then \(A'A\) is singular. Relate this to the statement in Section 5.3 that a multivariate analysis requires more observations than variables.

Questions on this section

Stuck on something here? Ask below and it stays attached to this topic.