Your American History Reference Guide!
- Schur complement

HistoryMania Information Site on Schur complement American History American History Search        American History Browse welcome to our free resource site for all enthusiasts!

Schur complement

In linear algebra and the theory of matrices, the Schur complement (named after Issai Schur) of a block of a matrix within the larger matrix is defined as follows. Suppose A, B, C, D are respectively p×p, p×q, q×p and q×q matrices, and D is invertible. Let

M=\left[\begin{matrix} A & B \\ C & D \end{matrix}\right]

so that M is a (p+q)×(p+q) matrix.

Then the Schur complement of the block D of the matrix M is the p×p matrix

A - BD - 1C.

The Schur complement arises as the result of performing a "partial" Gaussian elimination by multiplying the matrix M from the right with the "lower triangular" block matrix

LT=\left[\begin{matrix} E_p & 0 \\ -D^{-1}C & D^{-1} \end{matrix}\right].

Here Ep denotes a p×p unit matrix. After multiplication with the matrix LT the Schur complement appears in the upper p×p block. The product matrix is

M\cdot LT=\left[\begin{matrix} A-BD^{-1}C & BD^{-1} \\ 0 & E_q \end{matrix}\right].

If M is a positive definite symmetric matrix, then so is the Schur complement of D in M.

Applications to probability theory and statistics

Suppose the random column vectors X, Y live in Rn and Rm respectively, and the vector (X, Y) in 'Rn+m has a multivariate normal distribution whose variance is the symmetric positive-definite matrix

V=\left[\begin{matrix} A & B \\ B^T & C \end{matrix}\right].

Then the conditional variance of X given Y is the Schur complement of C in V:

\operatorname{var}(X\mid Y)=A-BC^{-1}B^T.

If we take the matrix V above to be, not a variance of a random vector, but a sample variance, then it may have a Wishart distribution. In that case, the Schur complement of C in V also has a Wishart distribution.

The contents of this article are licensed from Wikipedia.org under the
GNU Free Documentation License. How to see transparent copy
Search | Browse | Contact | Legal info