Maximizing the likelihood turns out to be closely related to minimizing the “distance” between Pθ∗ and Pθ. Here's how we formalize that notion of “distance”.
Definition. Given distributions P1 and P2 with PDFs f1(x) and f2(x), their Kullback-Leibler (KL) divergence is:
Here, R is the entire support of P1, which is also the entire support of P2. Also, as usual, the Ef1 notation indicates that the random variable X is to be sampled using the PDF f1.
Note that DKL is not a metric; it is not symmetric, nor does it satisfy the triangle inequality. But Jensen's inequality implies it is at least nonnegative, with DKL(P1∥P2)=0 if and only if P1=P2.
Remark. Jensen's inequality states that for any concave function φ and any random variable Z, we have E[φ(Z)]≤φ(E[Z]). Equality is precisely when Z is almost-surely constant, assuming φ is strictly concave.
Importantly, this implies the following.
Theorem. (KL Divergence Minimization) The map θ↦DKL(Pθ∗∥Pθ) is minimized when θ=θ∗. Also, the population log likelihood ℓ(θ):=Eθ∗[logfθ(X)] is maximized at θ=θ∗.
Proof: The first sentence is true by Jensen's; the point is that log is strictly concave, so:
Equality occurs precisely when fθ∗(X)fθ(X) (with X∼Pθ∗) is almost-surely constant, or when fθ≅fθ∗, as desired.
The second sentence follows because minimizing DKL(Pθ∗∥Pθ) is equivalent to maximizing ∫Rfθ∗(x)logfθ(x)dx; see the far-RHS of the definition of DKL. But ∫Rfθ∗(x)logfθ(x)dx=ℓ(θ) identically. ■
§ MLE Consistency
The upshot of introducing KL Divergence is the second sentence of the “KL Divergence Minimization” theorem. We will leverage this fact to prove that the MLE is a consistent estimator.
Theorem. (MLE Consistency) Under mild regularity conditions (identifiability in particular), the MLE is consistent; that is, if θ^n is the MLE of the parameter θ, then θ^n→Pθ∗.
Proof: Recall the following.
The MLE θ^n is the maximizer of the log likelihood ℓn(θ) by definition.
The true parameter value θ∗ is the maximizer of the population log likelihood ℓ(θ) by the previous theorem.
For any value θ′, by the LLN, the sequence {n1ℓn(θ′)}n=1∞ converges in probability to the constant ℓ(θ′).
These three points together imply θ^n→Pθ∗. Roughly speaking, the third bullet point says n1ℓn→Pℓ, so the maximizer of ℓn should also converge in probability to the maximizer of ℓ. (Yes, this needs to be justified—handwaving a little here.) ■