45.01.01 · mathematical-statistics / 01-decision-theory-estimation

Statistical Decision Theory: Loss, Risk, Admissibility, Minimax, and Bayes Rules

shipped3 tiersLean: none

Anchor (Master): Berger 1985 Statistical Decision Theory and Bayesian Analysis (2nd ed., Springer) Ch. 4-5, Ch. 8 (complete classes); Lehmann & Casella 1998 Theory of Point Estimation (2nd ed., Springer) Ch. 1, Ch. 5; Ferguson 1967 Mathematical Statistics: A Decision Theoretic Approach (Academic Press) Ch. 1-3

Intuition Beginner

Every time you use data to make a call — estimate a quantity, choose between two explanations, decide whether a part is defective — you are picking a rule that turns observations into an action. Two honest people with the same data can use different rules. Statistical decision theory is the framework for asking which rule is better, and what "better" even means once chance is in the picture.

The first ingredient is a way to score mistakes. A loss says how bad it is to take action when the truth is . Guessing a temperature of 19 when it is really 20 should cost less than guessing 5. A common score is the squared miss, but the framework allows any sensible cost. The loss is the value judgment; everything after it is mathematics.

A rule does not have one loss, it has a whole spread of them, because the data are random. So we average the loss over the randomness in the data, holding the truth fixed. That average is the risk of the rule at . A rule with low risk is one whose mistakes are small on average, whatever specific data happen to arrive.

The catch is that risk depends on the unknown truth . One rule might be best when is small, another best when is large. There is rarely a single rule that beats all others everywhere. So the theory offers principles for breaking the tie: avoid rules that are beaten outright (keep only the admissible ones), protect against the worst case (minimax), or weigh the possible truths by a prior belief and minimise the weighted average (Bayes). These three principles, and how they secretly fit together, are the whole subject.

Visual Beginner

Figure: a horizontal axis for the unknown truth and a vertical axis for risk. Two risk curves are drawn. Rule A is a flat line at a moderate height — same risk for every . Rule B dips very low near the centre but rises steeply at the edges. Where B is below A, B is the better rule; where B is above A, A wins. Neither curve lies entirely below the other, so neither rule dominates: both are admissible. A third dashed curve C sits entirely above A everywhere — C is dominated and therefore inadmissible.

 risk
  |                                          .C (dominated: above A everywhere)
  |   .........................................   <- dashed, always worst
  |
  |----A----------------------------------------  rule A: flat risk (minimax-ish)
  |        \                              /
  |         \         rule B            /
  |          \.____________________.../        B: low in the middle, high at edges
  |
  +------------------------------------------------> theta
       A vs B: neither dominates (curves cross)  -> both admissible
       C vs A: A dominates C (A below C for all theta) -> C inadmissible

Worked example Beginner

You observe one draw from a population whose unknown mean is , with the draw scattered around with variance . You want to estimate . Compare two rules: the rule (just report the data), and the shrink-toward-zero rule . Score each by squared error, .

Step 1. Risk of . The risk is the average squared error. For , the error is , whose average square is the variance, . So the risk is for every — a flat line.

Step 2. Risk of . Now the error is . Split it into a typical-offset part and a wobble part. On average sits at , so it is off-target by (this is the bias). The wobble of has variance . The average squared error is bias-squared plus wobble:

Step 3. Compare. When is better? Solve , i.e. , i.e. . So wins for truths near zero and loses for truths far out.

What this tells us. Neither rule beats the other everywhere: their risk curves cross. Both are admissible — you cannot discard either on dominance grounds alone. Choosing between them needs an extra principle: if you fear a large , the flat rule is safer (minimax); if you believe is probably near zero, the shrinker is better on average (Bayes).

Check your understanding Beginner

Formal definition Intermediate+

A statistical decision problem is specified by a parameter space , a family of probability measures on a sample space (the measure-theoretic substrate of 37.01.01), an action space , and a loss function that is measurable in its arguments. A (non-randomized) decision rule is a measurable map ; on observing one takes the action .

Definition (risk function). The risk of a rule at is the expected loss

A rule is identified with the function it induces. For point estimation with and squared-error loss , the risk is the mean squared error, with the bias-variance decomposition [Casella & Berger §7.3].

Definition (dominance and admissibility). A rule dominates if for all with strict inequality for at least one . A rule is inadmissible if some dominates it, and admissible otherwise. Admissibility is a minimal requirement, not a recommendation: it only rules out being beaten outright, and many admissible rules are poor.

Definition (minimax rule). The maximum risk of is . A rule is minimax if it minimises the maximum risk,

the minimax value of the problem. Minimaxity is the conservative, worst-case-protective principle.

Definition (Bayes risk and Bayes rule). Given a prior probability measure on , the Bayes risk of is the prior-averaged risk

A rule is a Bayes rule with respect to if it minimises the Bayes risk, , the Bayes value. The posterior obtained by conditioning on (the Bayesian update of 26.07.01) gives the posterior expected loss .

Definition (randomized rule and the risk set). A randomized rule chooses an action by drawing from a distribution on after seeing ; its risk is . When is finite, the risk set is

a convex subset of in the sense of 44.01.01: convexity comes from mixing rules, since the - random choice between and has risk vector the midpoint of theirs.

Counterexamples to common slips

  • Admissible does not mean good. The constant rule (ignore the data, always estimate ) is admissible under squared-error loss on : its risk is at , so no competitor can be uniformly no-worse without also being there, which pins it down. A useless rule can be admissible.
  • Minimax is not unique and need not be Bayes for an obvious prior. Several rules can share the minimax value; the minimax rule is the Bayes rule for a least favourable prior, which is often supported on the boundary of , not the prior a practitioner would write down.
  • A Bayes rule can be randomized only on a null set. When the posterior expected loss has a unique minimiser for almost every , the Bayes rule is essentially non-randomized; randomization buys nothing under strict convexity of the loss. Randomization matters for minimaxity in finite- problems, not for typical Bayes estimation.

Key theorem with proof Intermediate+

The organising theorem is that the Bayes rule is computed pointwise by minimising posterior expected loss, and that Bayes rules with finite Bayes risk are almost the cheapest possible certificate of admissibility.

Theorem (Bayes rule via posterior expected loss; admissibility of Bayes rules). Fix a prior with finite Bayes value .

(i) (Posterior characterization.) Suppose the joint law of admits the disintegration with marginal . If for -almost every the action minimises the posterior expected loss , then is a Bayes rule.

(ii) (Admissibility.) If the Bayes rule is unique up to risk equivalence — any other Bayes rule has the same risk function — then is admissible.

Proof of (i). For any rule , Tonelli's theorem (the integrand is non-negative) lets the Bayes risk be reordered:

The Bayes risk is thus the -integral of the posterior expected loss evaluated at the chosen action. To minimise the integral it suffices to minimise the integrand by , which is what does by hypothesis. Hence for every , so attains and is a Bayes rule.

Proof of (ii). Suppose is inadmissible. Then some dominates it: for all , strictly for some . Integrating against ,

Because is the minimum Bayes risk, the reverse inequality also holds, so and is itself a Bayes rule for . By the uniqueness hypothesis has the same risk function as , contradicting the strict inequality at . Hence no dominating exists and is admissible.

Bridge. This theorem builds toward the complete-class theory and the minimax constructions of the Advanced section, and appears again in the Rao-Blackwell improvement 45.01.02 and the Bayes-optimality of the Neyman-Pearson test 45.02.01. This is exactly the device that converts a global optimisation over rules into a pointwise optimisation over actions, paid for by the prior, and it generalises the elementary posterior-mean estimator of 26.07.01 to an arbitrary loss: under squared error the posterior-expected-loss minimiser is the posterior mean, under absolute error the posterior median, under 0-1 loss the posterior mode. The foundational reason Bayes rules are the workhorse is that uniqueness upgrades them to admissible rules, so the Bayes construction simultaneously delivers optimality-on-average and the no-dominance guarantee. Putting these together with the convexity of the risk set, the minimax rule turns out to be dual to a least-favourable prior — the same prior-averaging integral, now read as a saddle point — which is the connection the Advanced section makes precise via the separating-hyperplane geometry of 44.01.01.

Exercises Intermediate+

Advanced results Master

The geometry of the risk set and admissibility

Take finite and let be the risk set of all randomized rules, a convex set by mixing 44.01.01. Three principles become three geometric operations on . A rule is admissible if and only if its risk vector lies on the lower-left boundary of : the set of points for which no other point of is coordinatewise with a strict inequality, the Pareto frontier of . A rule is Bayes for if and only if its risk vector minimises the linear functional over , that is, has a supporting hyperplane with outer normal at that point; as ranges over the simplex, these supporting points trace out the closed convex lower boundary. A rule is minimax if and only if its risk vector minimises , the lowest corner of the smallest down-translate of the positive orthant meeting . The first identification shows that admissible rules are exactly the exposed and supporting points of the lower boundary; the second that every Bayes rule is a supporting point and hence admissible whenever the support is strict; the third that the minimax point sits where the diagonal meets the lower boundary, which is generically a Bayes point [Ferguson §2].

The minimax theorem and least favourable priors

The minimax-Bayes duality is the statistical face of von Neumann's minimax theorem applied to the loss game between Nature (choosing , or a prior ) and the statistician (choosing ). Under conditions guaranteeing a saddle point — compactness of the (randomized) decision space and risk-set convexity — the value of the game exists:

A prior attaining the right side is least favourable: it is the prior under which the best achievable Bayes risk is largest, and the Bayes rule against is minimax [Berger §5.2]. The standard production recipe inverts this: guess an equalizer rule with constant risk, exhibit a prior for which is Bayes, and conclude is minimax — the constant-risk and Bayes conditions together forcing the saddle. When no proper least favourable prior exists, one uses a least favourable sequence with , the device by which the sample mean is shown minimax for the normal mean via priors with [Lehmann & Casella §5.2].

Stein's phenomenon: admissibility is dimension-sensitive

The sample mean is admissible as an estimator of a normal mean in dimensions one and two under squared-error loss, but Charles Stein proved in 1956 that for with the estimator is inadmissible [Stein 1956]: the James-Stein estimator dominates it everywhere, shrinking toward the origin. The risk improvement is computed via Stein's unbiased risk identity for [Lehmann & Casella §5.5]. The phenomenon shows admissibility is not a local property of a single coordinate: borrowing strength across coordinates strictly reduces total risk despite each coordinate's estimate becoming biased. James-Stein is itself inadmissible (its positive-part truncation dominates it), so the example also illustrates that exhibiting a dominator does not produce an admissible terminus.

Complete class theorems

A class of rules is complete if for every rule there is a dominating it, and essentially complete if the dominating is merely no-worse. The structural payoff is that one may restrict the search for good rules to a complete class without loss. The central results, due to Wald and refined by Stein, Le Cam, and Brown, state that under convexity of the loss and mild regularity the Bayes rules and their limits form a complete class, and conversely every admissible rule is a (generalized or limiting) Bayes rule [Berger §8]. For finite the closed lower boundary of the convex risk set is the minimal complete class, and its supporting hyperplanes are exactly the priors, so "admissible," "supporting point of ," and "Bayes or limit of Bayes" coincide. The complete-class theorems are the converse direction of the admissibility-of-Bayes-rules result: not only is every (unique) Bayes rule admissible, but every admissible rule is, in the limit, Bayes.

Synthesis. The central insight of decision theory is that risk turns the choice of a statistical procedure into a partial order — dominance — whose maximal elements (admissible rules) are exactly the supporting points of the convex risk set, and the three classical principles are three ways of selecting among them. The foundational reason Bayes rules dominate the theory is the posterior-expected-loss characterization: it makes the Bayes rule computable pointwise, makes unique Bayes rules admissible, and — through the complete-class theorems — makes every admissible rule a limit of Bayes rules, so the Bayes construction generalises from one optimality principle into the generating mechanism for the whole admissible frontier. Putting these together, the minimax rule is dual to a least favourable prior in the literal saddle-point sense, the minimax value equalling , which is exactly von Neumann's game value read on the loss matrix and supported by the separating-hyperplane geometry of 44.01.01. The Stein phenomenon shows the order is genuinely global and dimension-sensitive, and the entire apparatus appears again in estimation 45.01.02, testing 45.02.01, and the Bayesian inference of 26.07.01, where these same objects — loss, risk, prior, posterior — recur as the unifying language. The bridge is the convex risk set: feasibility, optimality, admissibility, the Bayes-minimax duality, and the complete-class theorems are one theory about supporting hyperplanes of a single convex object.

Full proof set Master

Proposition 1 (the risk set is convex; Bayes points support it). For finite , the risk set is convex, and is Bayes for in the simplex if and only if .

Proof. Given rules and , the randomized rule that uses with probability and with probability (an independent coin flip before seeing data) has risk at each by linearity of expectation, so ; hence is convex (the mixing operation of 44.01.01). For the Bayes characterization, , so minimising Bayes risk over is minimising the linear functional over . The minimiser is a point at which the hyperplane supports from below.

Proposition 2 (a Bayes rule with constant risk is minimax). If is Bayes for some prior and is constant in , then is minimax and is least favourable.

Proof. For any rule , since is Bayes, . The left side is . The right side is a prior average and so cannot exceed the maximum risk: . Chaining, for every . Taking the infimum over , . But itself has , so . Hence and is minimax. Moreover (the last step is weak duality for every prior), so attains and is least favourable.

Proposition 3 (weak duality: every Bayes value is below the minimax value). For every prior , .

Proof. For any fixed , because an average over is dominated by the supremum. Taking the infimum over on both sides, . (This is the max-min inequality specialized; equality is the minimax theorem, the saddle realized by a least favourable and its Bayes rule.)

Proposition 4 (admissibility of a unique Bayes rule). If is Bayes for with and every Bayes rule for has the same risk function as , then is admissible.

Proof. This is part (ii) of the Key theorem, reproduced for the proof set. If dominated , integrating the dominance inequality against gives , hence since is minimal, so is also Bayes for . By the uniqueness hypothesis , contradicting the strict inequality that domination requires at some . So no dominator exists.

Connections Master

  • Unbiased estimation and the Rao-Blackwell / Lehmann-Scheffé theorems 45.01.02 are the squared-error specialization of this framework: conditioning an estimator on a sufficient statistic is a convexity-of-loss argument that lowers risk, the same risk and dominance order defined here, and the minimum-variance-unbiased estimator is the admissible representative inside the unbiased class.

  • Hypothesis testing and the Neyman-Pearson lemma 45.02.01 are the two-action decision problem under 0-1 (or weighted 0-1) loss; the most powerful test is the Bayes rule of Exercise 5 for the corresponding prior, and the admissibility and complete-class structure of tests is the same risk-set geometry, now over the power/size plane.

  • Bayesian inference 26.07.01 supplies the prior, posterior, and posterior-expected-loss machinery that the Key theorem turns into the Bayes rule; this unit is the decision-theoretic completion of that survey, replacing "report the posterior" with "choose the action minimising posterior expected loss for a stated loss function."

  • The convex-set and supporting-hyperplane calculus of 44.01.01 is the geometric engine: the risk set is convex by rule-mixing, Bayes rules are its supporting points, admissible rules are its lower boundary, and the minimax-Bayes duality is the separating-hyperplane / saddle-point structure read on this single convex object, exactly as Lagrangian duality reads it on the value set.

Historical & philosophical context Master

The decision-theoretic reframing of statistics is due to Abraham Wald, whose Statistical Decision Functions [Wald 1950] (Wiley, 1950) introduced the loss function, the risk function, randomized decision rules, and the complete-class program, casting estimation and testing as special cases of a single game against Nature; the minimax principle and the proof that Bayes and limiting-Bayes rules form a complete class are Wald's. The game-theoretic substrate traces to John von Neumann's 1928 minimax theorem for zero-sum games, which supplies the saddle-point existence behind the minimax-Bayes duality, and the decision-game reading was sharpened by Wald and by L. J. Savage's subjective-expected-utility foundations in The Foundations of Statistics (1954).

Charles Stein's 1956 Berkeley Symposium paper [Stein 1956] (Proc. Third Berkeley Symposium 1, 197-206) proved the inadmissibility of the sample mean in three or more dimensions, a result so counter to intuition that it reoriented the field's understanding of admissibility; the explicit dominating estimator was given by James and Stein in 1961. The systematic theory of admissibility and complete classes was developed through the 1950s-1970s by Stein, Lucien Le Cam, and Lawrence Brown, whose differential-equation characterization of admissibility for the normal mean closed the dimension-one and -two admissibility questions. Thomas Ferguson's Mathematical Statistics: A Decision Theoretic Approach [Ferguson §2] (1967) gave the geometric risk-set exposition, and James Berger's Statistical Decision Theory and Bayesian Analysis [Berger §4-§5] (1985) and Lehmann and Casella's Theory of Point Estimation [Lehmann & Casella §5] (1998) are the standard modern references, the latter integrating the decision framework with the classical optimality theory of estimation.

Bibliography Master

@book{wald1950statistical,
  author    = {Wald, Abraham},
  title     = {Statistical Decision Functions},
  publisher = {Wiley},
  year      = {1950}
}

@inproceedings{stein1956inadmissibility,
  author    = {Stein, Charles},
  title     = {Inadmissibility of the usual estimator for the mean of a multivariate normal distribution},
  booktitle = {Proceedings of the Third Berkeley Symposium on Mathematical Statistics and Probability},
  volume    = {1},
  pages     = {197--206},
  publisher = {University of California Press},
  year      = {1956}
}

@book{ferguson1967mathematical,
  author    = {Ferguson, Thomas S.},
  title     = {Mathematical Statistics: A Decision Theoretic Approach},
  publisher = {Academic Press},
  year      = {1967}
}

@book{berger1985statistical,
  author    = {Berger, James O.},
  title     = {Statistical Decision Theory and Bayesian Analysis},
  edition   = {2nd},
  series    = {Springer Series in Statistics},
  publisher = {Springer},
  year      = {1985}
}

@book{lehmanncasella1998theory,
  author    = {Lehmann, Erich L. and Casella, George},
  title     = {Theory of Point Estimation},
  edition   = {2nd},
  series    = {Springer Texts in Statistics},
  publisher = {Springer},
  year      = {1998}
}

@book{casellaberger2002statistical,
  author    = {Casella, George and Berger, Roger L.},
  title     = {Statistical Inference},
  edition   = {2nd},
  publisher = {Duxbury},
  year      = {2002}
}