๐Ÿงฎ Brain Teaser

The Precision Matrix and Partial Correlations

Let X=(X1,X2,โ€ฆ,Xp)TX = (X_1, X_2, \ldots, X_p)^T be a multivariate Gaussian random vector with covariance matrix ฮฃ\Sigma (assumed invertible). Define the precision matrix ฮ˜=ฮฃโˆ’1\Theta = \Sigma^{-1}.

Claim: XiX_i and XjX_j are conditionally independent given all other variables XrestX_{\text{rest}} if and only if ฮ˜ij=0\Theta_{ij} = 0.

Question: Prove (or give a clean argument for) this claim. In particular, explain why zeros in the precision matrix โ€” not in the covariance matrix โ€” encode conditional independence. What does ฮ˜ij=0\Theta_{ij} = 0 geometrically mean about the conditional distribution?

As a warm-up: why does ฮฃij=0\Sigma_{ij} = 0 (zero covariance) not imply conditional independence in general?

Gaussian graphical modelsconditional independenceprecision matrixpartial correlationSchur complement

Answer: Precision Matrix and Partial Correlations

Key Idea / Intuition

In a Gaussian, all the action of conditional distributions is controlled by partial correlations โ€” correlations after "regressing out" the effect of other variables. The precision matrix ฮ˜=ฮฃโˆ’1\Theta = \Sigma^{-1} encodes exactly these residual relationships: ฮ˜ij\Theta_{ij} measures the direct linear connection between XiX_i and XjX_j after removing all indirect paths through other variables. Zero covariance only says XiX_i and XjX_j are marginally uncorrelated โ€” but conditioning on other variables can create dependence (Berkson's paradox). Zero in ฮ˜\Theta is the right notion because the Gaussian density factors precisely when the off-diagonal precision entry vanishes.


Formal Proof / Solution

Warm-up: Why ฮฃij=0\Sigma_{ij} = 0 is not enough

Consider X1,X2โˆผiidN(0,1)X_1, X_2 \overset{\text{iid}}{\sim} N(0,1) and X3=X1+X2X_3 = X_1 + X_2. Then Cov(X1,X2)=0\text{Cov}(X_1, X_2) = 0, so marginally they are independent. But given X3X_3, knowing X1X_1 tells you exactly X2=X3โˆ’X1X_2 = X_3 - X_1 โ€” they are perfectly conditionally dependent. Conditioning can introduce dependence where none existed marginally (Berkson's paradox / collider effect).


The Gaussian Conditional Distribution

Partition X=(Xi,Xj,Xrest)X = (X_i, X_j, X_{\text{rest}}) and consider the joint Gaussian. For the bivariate case, it suffices to look at (Xi,Xj)โˆฃXrest(X_i, X_j) \mid X_{\text{rest}}.

The conditional distribution of (Xi,Xj)โˆฃXrest(X_i, X_j) \mid X_{\text{rest}} is Gaussian with partial covariance matrix:

ฮฃijโˆฃrest=ฮฃaaโˆ’ฮฃabฮฃbbโˆ’1ฮฃba\Sigma_{ij \mid \text{rest}} = \Sigma_{aa} - \Sigma_{ab}\Sigma_{bb}^{-1}\Sigma_{ba}

where a={i,j}a = \{i,j\} and b=restb = \text{rest}. This is the Schur complement.

Key algebraic fact: By the block matrix inversion formula,

(ฮฃaaโˆ’ฮฃabฮฃbbโˆ’1ฮฃba)โˆ’1=ฮ˜aa(\Sigma_{aa} - \Sigma_{ab}\Sigma_{bb}^{-1}\Sigma_{ba})^{-1} = \Theta_{aa}

where ฮ˜aa\Theta_{aa} is the 2ร—22\times 2 submatrix of ฮ˜=ฮฃโˆ’1\Theta = \Sigma^{-1} corresponding to indices {i,j}\{i,j\}.

So the conditional precision of (Xi,Xj)โˆฃXrest(X_i, X_j) \mid X_{\text{rest}} is exactly ฮ˜aa=(ฮ˜iiฮ˜ijฮ˜jiฮ˜jj)\Theta_{aa} = \begin{pmatrix} \Theta_{ii} & \Theta_{ij} \\ \Theta_{ji} & \Theta_{jj} \end{pmatrix}.


Conditional Independence โ†” ฮ˜ij=0\Theta_{ij} = 0

For jointly Gaussian variables, conditional independence XiโŠฅXjโˆฃXrestX_i \perp X_j \mid X_{\text{rest}} is equivalent to:

Cov(Xi,XjโˆฃXrest)=0\text{Cov}(X_i, X_j \mid X_{\text{rest}}) = 0

i.e., the off-diagonal entry of the conditional covariance ฮฃijโˆฃrest=0\Sigma_{ij\mid\text{rest}} = 0.

The conditional covariance matrix is ฮ˜aaโˆ’1\Theta_{aa}^{-1}. For a 2ร—22\times 2 matrix:

ฮ˜aaโˆ’1=1ฮ˜iiฮ˜jjโˆ’ฮ˜ij2(ฮ˜jjโˆ’ฮ˜ijโˆ’ฮ˜ijฮ˜ii)\Theta_{aa}^{-1} = \frac{1}{\Theta_{ii}\Theta_{jj} - \Theta_{ij}^2} \begin{pmatrix} \Theta_{jj} & -\Theta_{ij} \\ -\Theta_{ij} & \Theta_{ii} \end{pmatrix}

The off-diagonal entry of ฮ˜aaโˆ’1\Theta_{aa}^{-1} is โˆโˆ’ฮ˜ij\propto -\Theta_{ij}.

Therefore:

Cov(Xi,XjโˆฃXrest)=0โ€…โ€ŠโŸบโ€…โ€Šฮ˜ij=0\text{Cov}(X_i, X_j \mid X_{\text{rest}}) = 0 \iff \Theta_{ij} = 0


The Density Factorization Perspective

The multivariate Gaussian log-density is:

logโกp(x)=โˆ’12xTฮ˜โ€‰x+linearย terms+const\log p(x) = -\frac{1}{2} x^T \Theta\, x + \text{linear terms} + \text{const}

The quadratic form expands as:

xTฮ˜โ€‰x=โˆ‘k,lฮ˜klโ€‰xkxlx^T \Theta\, x = \sum_{k,l} \Theta_{kl}\, x_k x_l

The cross-term between xix_i and xjx_j is 2ฮ˜ijxixj2\Theta_{ij} x_i x_j. If ฮ˜ij=0\Theta_{ij} = 0, this cross term vanishes, so the density factors in (xi,xj)(x_i, x_j) given the rest โ€” which is precisely conditional independence.


The Partial Correlation Formula

The partial correlation between XiX_i and XjX_j given the rest is:

ฯijโˆฃrest=โˆ’ฮ˜ijฮ˜iiฮ˜jj\rho_{ij \mid \text{rest}} = -\frac{\Theta_{ij}}{\sqrt{\Theta_{ii}\Theta_{jj}}}

(the minus sign comes from the inverse formula above). So not only does ฮ˜ij=0\Theta_{ij} = 0 encode conditional independence โ€” the magnitude of ฮ˜ij\Theta_{ij} encodes the strength of the direct connection, making ฮ˜\Theta the natural object for Gaussian graphical models (draw an edge iโˆผji \sim j iff ฮ˜ijโ‰ 0\Theta_{ij} \neq 0).


Summary

| Object | What it encodes | |---|---| | ฮฃij=0\Sigma_{ij} = 0 | Marginal independence (only!) | | ฮ˜ij=0\Theta_{ij} = 0 | Conditional independence given all others |

This is why methods like the graphical LASSO penalize entries of ฮ˜\Theta rather than ฮฃ\Sigma โ€” they are recovering the graph of direct dependencies.

Source: The Elements of Statistical Learning, 2nd ed. โ€” Exercise 17.3 and Section 17.3

Type: ML/StatsSource: The Elements of Statistical Learning, 2nd ed. โ€” Exercise 17.3 and Section 17.3Edit on GitHub โ†—