research33 min studyLesson 730 in Mathematics

Information-theoretic minimax lower bounds

Reduce estimation to testing separated hypotheses and distinguish statistical impossibility from computational limits.

On this page

What you’ll understand

  • Construct a separated finite subset of a parameter class.
  • Combine information bounds with a nearest-hypothesis reduction.
  • Interpret sparse estimation rates without confusing them with support or computational guarantees.

An upper bound describes what a particular estimator can achieve. A minimax lower bound asks what no estimator can uniformly improve over a specified class, regardless of computation. The central method hides the true parameter among many separated alternatives whose induced distributions remain difficult to distinguish. A successful estimator would reveal which alternative generated the data, so a testing impossibility becomes an estimation impossibility.

Choose a loss and a difficult finite experiment

For squared Euclidean loss, minimax risk takes the infimum over estimators of their worst expected error over the parameter class. A lower bound can restrict attention to a finite subset because the worst case over the full class is at least the worst case over that subset. Choose hypotheses separated by at least twice a radius. If an estimator lands within that radius of the truth, nearest-neighbor decoding identifies the correct hypothesis. A testing error therefore implies that the estimation error is at least the radius.

This reduction requires the metric used to separate hypotheses to match the loss being bounded. Prediction loss in a design image is not automatically equivalent to coefficient loss, just as support error is not the same as Euclidean error. The finite subset must also lie entirely in the advertised parameter class: a lower bound using dense alternatives does not establish a sparse-class lower bound. These bookkeeping conditions are part of the argument rather than details to fill in after applying an information inequality.

Rn(Θ)=inf⁡θ^sup⁡θ∈ΘEθ∥θ^−θ∥22,KL⁡(Pθ∥Pθ′)=n∥θ−θ′∥222σ2 for Pθ=Np(θ,σ2Ip/n).\mathcal R_n(\Theta)=\inf_{\widehat\theta}\sup_{\theta\in\Theta}\mathbb E_\theta\|\widehat\theta-\theta\|_2^2,\qquad \operatorname{KL}(P_\theta\|P_{\theta'})=\frac{n\|\theta-\theta'\|_2^2}{2\sigma^2}\ \text{for }P_\theta=N_p(\theta,\sigma^2I_p/n).
The Gaussian experiment is one noisy vector with covariance sigma squared divided by n; it is equivalent to the sample mean of n independent Gaussian observations.↔ Scroll wide equations horizontally.

Information limits decoding

Let the hypothesis index be uniform on M≥2M\geq2 alternatives and let the data follow the corresponding distribution. Fano’s inequality bounds average decoding error below by one minus mutual information plus a binary entropy constant divided by the logarithm of the number of alternatives. All logarithms here are natural. Mutual information can be bounded above by the average divergence from any fixed reference distribution when those divergences are finite. Choosing the reference well avoids unnecessary calculations with the mixture distribution.

Combining Fano with the separation argument gives the displayed expected-risk lower bound, with a nonpositive bracket understood as the trivial zero bound. There is a balance: increasing hypothesis separation improves the loss scale but also increases divergence and makes testing easier. Increasing the packing size supplies more uncertainty about the index, but the parameter class limits how many well-separated points fit. A valid proof chooses amplitude and packing together rather than inserting arbitrary large signals into a convenient entropy formula.

Rn(Θ)≥r2[1−I(J;Y)+log⁡2log⁡M]+,I(J;Y)≤1M∑j=1MKL⁡(Pj∥Q),∥θj−θk∥2≥2r.\mathcal R_n(\Theta)\geq r^2\left[1-\frac{I(J;Y)+\log2}{\log M}\right]_+,\qquad I(J;Y)\leq\frac1M\sum_{j=1}^M\operatorname{KL}(P_j\|Q),\quad \|\theta_j-\theta_k\|_2\geq2r.
The packing points belong to Theta, Q is a common finite-divergence reference, and the index is uniform. The positive part only removes a vacuous negative bound.↔ Scroll wide equations horizontally.

What a sparse minimax rate actually says

In the Gaussian sequence model over vectors with at most ss nonzero coordinates, for 1≤s≤p/21\leq s\leq p/2, the squared Euclidean minimax risk is of order σ2slog⁡(ep/s)/n\sigma^2s\log(ep/s)/n, up to universal constants. A sparse packing captures both unknown values and unknown support. Corresponding upper bounds show that this is a statistical benchmark for that experiment. Transporting the rate to a regression design requires conditions connecting its measurement geometry to the coefficient metric; a singular unrestricted design can invalidate such a transfer.

An information-theoretic lower bound includes computationally unbounded estimators and therefore cannot by itself establish a computational hardness gap. Nor does a Euclidean minimax rate imply exact support recovery: arbitrarily small nonzero coefficients remain compatible with that class. A fixed-dimensional two-point construction can also miss a logarithmic dimension penalty, so many-hypothesis arguments are often needed. The practical value of a lower bound is its precise scope: it tells whether improving an estimator’s rate is statistically possible under the stated experiment, class, and loss.

A sixteen-way sparse testing problem

In Y∼N16(θ,I16/100)Y\sim N_{16}(\theta,I_{16}/100), take sixteen one-sparse alternatives θj=aej\theta_j=ae_j with a2=log⁡16/200a^2=\log16/200. Use the zero-mean experiment as reference to obtain a numerical squared-risk lower bound.

  1. Two distinct alternatives have squared separation 2a22a^2, so choose radius squared r2=a2/2r^2=a^2/2. Every alternative is one-sparse, and all are included in the one-sparse parameter class.

  2. For the zero-mean reference, each divergence is 100a2/2=log⁡16/4=log⁡2100a^2/2=\log16/4=\log2. Thus mutual information is at most log⁡2\log2, and the Fano bracket is at least 1−(log⁡2+log⁡2)/log⁡16=1/21-(\log2+\log2)/\log16=1/2.

  3. The worst expected squared error of any estimator over this class is at least r2/2=a2/4=log⁡16/800≈0.00346574r^2/2=a^2/4=\log16/800\approx0.00346574. This is a lower bound from a particular packing, not an assertion that its constant is the exact minimax constant.

The result

Even an estimator with unlimited computation cannot uniformly push the squared risk below this valid finite-packing lower bound.

Explore the model

Interactive lab

Interactive lab · probability

Overlapping experiments are hard to distinguish

Pause/Reset, set Bernoulli probability 0.4 and trials 10, and predict count mean 4 and variance 2.4. For a fresh comparison Pause/Reset, set both probability 0.6 and trials 10; the mean becomes 6 with the same variance. In both runs the probability bar at count 5 is 252 times 0.24 to the fifth, about 0.20066.

Overlapping experiments are hard to distinguishThese are exact binomial probabilities for independent trials with one fixed success probability. Playback tours the distribution; it does not resample outcomes. Current readings: Mean · np: 5; Variance · np(1−p): 2.5; Total probability: 1. All parameters are adjustable with the labeled controls below.03581000.070.150.220.3Successes kP(X = k)k = 0 · P = 9.77e-4
Exact probabilityHighlighted outcomeMean np
Mean · np
5
Variance · np(1−p)
2.5
Total probability
1
t = 0 s

These are exact binomial probabilities for independent trials with one fixed success probability. Playback tours the distribution; it does not resample outcomes. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forDifferent parameters can produce substantially overlapping data laws. These exact binomial distributions illustrate testing ambiguity, not the Gaussian packing or an automatic calculation of mutual information.

Check your understanding

What does a minimax lower bound under squared coefficient loss establish?

Try it yourself

Consider the one-dimensional alternatives θ=±a\theta=\pm a with a=1/(2n)a=1/(2\sqrt n) and observation variance 1/n1/n. Use the optimal equal-prior sign test to obtain a squared-risk lower bound, then evaluate it at n=100n=100.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

The likelihood-ratio threshold is zero. Under the positive mean, a wrong sign has probability Phi of minus one half; decoding error forces absolute estimation error at least a.

Compare with the solution

The optimal average testing error is Φ(−an)=Φ(−1/2)≈0.3085375\Phi(-a\sqrt n)=\Phi(-1/2)\approx0.3085375. Nearest-sign decoding from any estimator gives worst squared risk at least a2Φ(−1/2)=Φ(−1/2)/(4n)a^2\Phi(-1/2)=\Phi(-1/2)/(4n). At sample size one hundred this is approximately 0.0007713440.000771344. The argument is valid but one-dimensional: unlike a large sparse packing, it cannot by itself explain a logarithmic dependence on an ambient dimension that never appeared among its alternatives.

Take it with you

  • Lower bounds must match an experiment, a parameter class, and a loss.
  • Packing separation and information must be balanced to make a Fano bound nonvacuous.
  • Statistical impossibility, computational hardness, and exact support recovery are different conclusions.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.