Chapman–Robbins bound

From HandWiki

In statistics, the Chapman–Robbins bound or Hammersley–Chapman–Robbins bound is a lower bound on the variance of estimators of a deterministic parameter. It is a generalization of the Cramér–Rao bound; compared to the Cramér–Rao bound, it is both tighter and applicable to a wider range of problems. However, it is usually more difficult to compute.

The bound was independently discovered by John Hammersley in 1950,[1] and by Douglas Chapman and Herbert Robbins in 1951.[2]

Statement

Let Θ be the set of parameters for a family of probability distributions {μθ:θ∈Θ} on Ω.

For any two θ,θ′∈Θ, let χ2(μθ′;μθ) be the χ2-divergence from μθ to μθ′. Then:

Theorem — Given any scalar random variable g^:Ω→ℝ, and any two θ,θ′∈Θ, we have Varθ[g^]≥supθ′≠θ∈Θ(Eθ′[g^]−Eθ[g^])2χ2(μθ′;μθ).

A generalization to the multivariable case is:[3]

Theorem — Given any multivariate random variable g^:Ω→ℝm, and any θ,θ′∈Θ, χ2(μθ′;μθ)≥(Eθ′[g^]−Eθ[g^])TCovθ[g^]−1(Eθ′[g^]−Eθ[g^])

Proof

By the variational representation of chi-squared divergence:[3]χ2(P;Q)=supg(EP[g]−EQ[g])2VarQ[g] Plug in g=g^,P=μθ′,Q=μθ, to obtain: χ2(μθ′;μθ)≥(Eθ′[g^]−Eθ[g^])2Varθ[g^]Switch the denominator and the left side and take supremum over θ′ to obtain the single-variate case. For the multivariate case, we define h=∑i=1mvig^i for any v≠0∈ℝm. Then plug in g=h in the variational representation to obtain: χ2(μθ′;μθ)≥(Eθ′[h]−Eθ[h])2Varθ[h]=⟨v,Eθ′[g^]−Eθ[g^]⟩2vTCovθ[g^]vTake supremum over v≠0∈ℝm, using the linear algebra fact that supv≠0vTwwTvvTMv=wTM−1w, we obtain the multivariate case.

Relation to Cramér–Rao bound

Usually, Ω=𝒳n is the sample space of n independent draws of a 𝒳-valued random variable X with distribution λθ from a by θ∈Θ⊆ℝm parameterized family of probability distributions, μθ=λθ⊗n is its n-fold product measure, and g^:𝒳n→Θ is an estimator of θ. Then, for m=1, the expression inside the supremum in the Chapman–Robbins bound converges to the Cramér–Rao bound of g^ when θ′→θ, assuming the regularity conditions of the Cramér–Rao bound hold. This implies that, when both bounds exist, the Chapman–Robbins version is always at least as tight as the Cramér–Rao bound; in many cases, it is substantially tighter.

The Chapman–Robbins bound also holds under much weaker regularity conditions. For example, no assumption is made regarding differentiability of the probability density function p(x; θ) of λθ. When p(x; θ) is non-differentiable, the Fisher information is not defined, and hence the Cramér–Rao bound does not exist.

See also

References

  1. ↑ "On estimating restricted parameters", Journal of the Royal Statistical Society, Series B 12 (2): 192–240, 1950, doi:10.1111/j.2517-6161.1950.tb00056.x 
  2. ↑ Chapman, D. G. (1951), "Minimum variance estimation without regularity assumptions", Annals of Mathematical Statistics 22 (4): 581–586, doi:10.1214/aoms/1177729548 
  3. ↑ 3.0 3.1 Polyanskiy, Yury (2017). "Lecture notes on information theory, chapter 29, ECE563 (UIUC)". https://people.lids.mit.edu/yp/homepage/data/LN_stats.pdf. 

Further reading

  • Lehmann, E. L.; Casella, G. (1998), Theory of Point Estimation (2nd ed.), Springer, pp. 113–114, ISBN 0-387-98502-6