The Precision Matrix and Partial Correlations
Let be a multivariate Gaussian random vector with covariance matrix (assumed invertible). Define the precision matrix .
Claim: and are conditionally independent given all other variables if and only if .
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 geometrically mean about the conditional distribution?
As a warm-up: why does (zero covariance) not imply conditional independence in general?
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 encodes exactly these residual relationships: measures the direct linear connection between and after removing all indirect paths through other variables. Zero covariance only says and are marginally uncorrelated โ but conditioning on other variables can create dependence (Berkson's paradox). Zero in is the right notion because the Gaussian density factors precisely when the off-diagonal precision entry vanishes.
Formal Proof / Solution
Warm-up: Why is not enough
Consider and . Then , so marginally they are independent. But given , knowing tells you exactly โ they are perfectly conditionally dependent. Conditioning can introduce dependence where none existed marginally (Berkson's paradox / collider effect).
The Gaussian Conditional Distribution
Partition and consider the joint Gaussian. For the bivariate case, it suffices to look at .
The conditional distribution of is Gaussian with partial covariance matrix:
where and . This is the Schur complement.
Key algebraic fact: By the block matrix inversion formula,
where is the submatrix of corresponding to indices .
So the conditional precision of is exactly .
Conditional Independence โ
For jointly Gaussian variables, conditional independence is equivalent to:
i.e., the off-diagonal entry of the conditional covariance .
The conditional covariance matrix is . For a matrix:
The off-diagonal entry of is .
Therefore:
The Density Factorization Perspective
The multivariate Gaussian log-density is:
The quadratic form expands as:
The cross-term between and is . If , this cross term vanishes, so the density factors in given the rest โ which is precisely conditional independence.
The Partial Correlation Formula
The partial correlation between and given the rest is:
(the minus sign comes from the inverse formula above). So not only does encode conditional independence โ the magnitude of encodes the strength of the direct connection, making the natural object for Gaussian graphical models (draw an edge iff ).
Summary
| Object | What it encodes | |---|---| | | Marginal independence (only!) | | | Conditional independence given all others |
This is why methods like the graphical LASSO penalize entries of rather than โ they are recovering the graph of direct dependencies.
Source: The Elements of Statistical Learning, 2nd ed. โ Exercise 17.3 and Section 17.3