Admissible decision rule

From HandWiki
Short description: Type of "good" decision rule in Bayesian statistics

In statistical decision theory, an admissible decision rule is a rule for making a decision such that there is no other rule that is always "better" than it[1] (or at least sometimes better and never worse), in the precise sense of "better" defined below. This concept is analogous to Pareto efficiency.

Definition

Define sets [math]\displaystyle{ \Theta\, }[/math], [math]\displaystyle{ \mathcal{X} }[/math] and [math]\displaystyle{ \mathcal{A} }[/math], where [math]\displaystyle{ \Theta\, }[/math] are the states of nature, [math]\displaystyle{ \mathcal{X} }[/math] the possible observations, and [math]\displaystyle{ \mathcal{A} }[/math] the actions that may be taken. An observation of [math]\displaystyle{ x \in \mathcal{X}\,\! }[/math] is distributed as [math]\displaystyle{ F(x\mid\theta)\,\! }[/math] and therefore provides evidence about the state of nature [math]\displaystyle{ \theta\in\Theta\,\! }[/math]. A decision rule is a function [math]\displaystyle{ \delta:{\mathcal{X}}\rightarrow {\mathcal{A}} }[/math], where upon observing [math]\displaystyle{ x\in \mathcal{X} }[/math], we choose to take action [math]\displaystyle{ \delta(x)\in \mathcal{A}\,\! }[/math].

Also define a loss function [math]\displaystyle{ L: \Theta \times \mathcal{A} \rightarrow \mathbb{R} }[/math], which specifies the loss we would incur by taking action [math]\displaystyle{ a \in \mathcal{A} }[/math] when the true state of nature is [math]\displaystyle{ \theta \in \Theta }[/math]. Usually we will take this action after observing data [math]\displaystyle{ x \in \mathcal{X} }[/math], so that the loss will be [math]\displaystyle{ L(\theta,\delta(x))\,\! }[/math]. (It is possible though unconventional to recast the following definitions in terms of a utility function, which is the negative of the loss.)

Define the risk function as the expectation

[math]\displaystyle{ R(\theta,\delta)=\operatorname{E}_{F(x\mid\theta)}[{L(\theta,\delta(x))]}.\,\! }[/math]

Whether a decision rule [math]\displaystyle{ \delta\,\! }[/math] has low risk depends on the true state of nature [math]\displaystyle{ \theta\,\! }[/math]. A decision rule [math]\displaystyle{ \delta^*\,\! }[/math] dominates a decision rule [math]\displaystyle{ \delta\,\! }[/math] if and only if [math]\displaystyle{ R(\theta,\delta^*)\le R(\theta,\delta) }[/math] for all [math]\displaystyle{ \theta\,\! }[/math], and the inequality is strict for some [math]\displaystyle{ \theta\,\! }[/math].

A decision rule is admissible (with respect to the loss function) if and only if no other rule dominates it; otherwise it is inadmissible. Thus an admissible decision rule is a maximal element with respect to the above partial order. An inadmissible rule is not preferred (except for reasons of simplicity or computational efficiency), since by definition there is some other rule that will achieve equal or lower risk for all [math]\displaystyle{ \theta\,\! }[/math]. But just because a rule [math]\displaystyle{ \delta\,\! }[/math] is admissible does not mean it is a good rule to use. Being admissible means there is no other single rule that is always as good or better – but other admissible rules might achieve lower risk for most [math]\displaystyle{ \theta\,\! }[/math] that occur in practice. (The Bayes risk discussed below is a way of explicitly considering which [math]\displaystyle{ \theta\,\! }[/math] occur in practice.)

Bayes rules and generalized Bayes rules

Bayes rules

Let [math]\displaystyle{ \pi(\theta)\,\! }[/math] be a probability distribution on the states of nature. From a Bayesian point of view, we would regard it as a prior distribution. That is, it is our believed probability distribution on the states of nature, prior to observing data. For a frequentist, it is merely a function on [math]\displaystyle{ \Theta\,\! }[/math] with no such special interpretation. The Bayes risk of the decision rule [math]\displaystyle{ \delta\,\! }[/math] with respect to [math]\displaystyle{ \pi(\theta)\,\! }[/math] is the expectation

[math]\displaystyle{ r(\pi,\delta)=\operatorname{E}_{\pi(\theta)}[R(\theta,\delta)].\,\! }[/math]

A decision rule [math]\displaystyle{ \delta\,\! }[/math] that minimizes [math]\displaystyle{ r(\pi,\delta)\,\! }[/math] is called a Bayes rule with respect to [math]\displaystyle{ \pi(\theta)\,\! }[/math]. There may be more than one such Bayes rule. If the Bayes risk is infinite for all [math]\displaystyle{ \delta\,\! }[/math], then no Bayes rule is defined.

Generalized Bayes rules

In the Bayesian approach to decision theory, the observed [math]\displaystyle{ x\,\! }[/math] is considered fixed. Whereas the frequentist approach (i.e., risk) averages over possible samples [math]\displaystyle{ x \in \mathcal{X}\,\! }[/math], the Bayesian would fix the observed sample [math]\displaystyle{ x\,\! }[/math] and average over hypotheses [math]\displaystyle{ \theta \in \Theta\,\! }[/math]. Thus, the Bayesian approach is to consider for our observed [math]\displaystyle{ x\,\! }[/math] the expected loss

[math]\displaystyle{ \rho(\pi,\delta \mid x)=\operatorname{E}_{\pi(\theta \mid x)} [ L(\theta,\delta(x)) ]. \,\! }[/math]

where the expectation is over the posterior of [math]\displaystyle{ \theta\,\! }[/math] given [math]\displaystyle{ x\,\! }[/math] (obtained from [math]\displaystyle{ \pi(\theta)\,\! }[/math] and [math]\displaystyle{ F(x\mid\theta)\,\! }[/math] using Bayes' theorem).

Having made explicit the expected loss for each given [math]\displaystyle{ x\,\! }[/math] separately, we can define a decision rule [math]\displaystyle{ \delta\,\! }[/math] by specifying for each [math]\displaystyle{ x\,\! }[/math] an action [math]\displaystyle{ \delta(x)\,\! }[/math] that minimizes the expected loss. This is known as a generalized Bayes rule with respect to [math]\displaystyle{ \pi(\theta)\,\! }[/math]. There may be more than one generalized Bayes rule, since there may be multiple choices of [math]\displaystyle{ \delta(x)\,\! }[/math] that achieve the same expected loss.

At first, this may appear rather different from the Bayes rule approach of the previous section, not a generalization. However, notice that the Bayes risk already averages over [math]\displaystyle{ \Theta\,\! }[/math] in Bayesian fashion, and the Bayes risk may be recovered as the expectation over [math]\displaystyle{ \mathcal{X} }[/math] of the expected loss (where [math]\displaystyle{ x\sim\theta\,\! }[/math] and [math]\displaystyle{ \theta\sim\pi\,\! }[/math]). Roughly speaking, [math]\displaystyle{ \delta\,\! }[/math] minimizes this expectation of expected loss (i.e., is a Bayes rule) if and only if it minimizes the expected loss for each [math]\displaystyle{ x \in \mathcal{X} }[/math] separately (i.e., is a generalized Bayes rule).

Then why is the notion of generalized Bayes rule an improvement? It is indeed equivalent to the notion of Bayes rule when a Bayes rule exists and all [math]\displaystyle{ x\,\! }[/math] have positive probability. However, no Bayes rule exists if the Bayes risk is infinite (for all [math]\displaystyle{ \delta\,\! }[/math]). In this case it is still useful to define a generalized Bayes rule [math]\displaystyle{ \delta\,\! }[/math], which at least chooses a minimum-expected-loss action [math]\displaystyle{ \delta(x)\!\, }[/math] for those [math]\displaystyle{ x\,\! }[/math] for which a finite-expected-loss action does exist. In addition, a generalized Bayes rule may be desirable because it must choose a minimum-expected-loss action [math]\displaystyle{ \delta(x)\,\! }[/math] for every [math]\displaystyle{ x\,\! }[/math], whereas a Bayes rule would be allowed to deviate from this policy on a set [math]\displaystyle{ X \subseteq \mathcal{X} }[/math] of measure 0 without affecting the Bayes risk.

More important, it is sometimes convenient to use an improper prior [math]\displaystyle{ \pi(\theta)\,\! }[/math]. In this case, the Bayes risk is not even well-defined, nor is there any well-defined distribution over [math]\displaystyle{ x\,\! }[/math]. However, the posterior [math]\displaystyle{ \pi(\theta\mid x)\,\! }[/math]—and hence the expected loss—may be well-defined for each [math]\displaystyle{ x\,\! }[/math], so that it is still possible to define a generalized Bayes rule.

Admissibility of (generalized) Bayes rules

According to the complete class theorems, under mild conditions every admissible rule is a (generalized) Bayes rule (with respect to some prior [math]\displaystyle{ \pi(\theta)\,\! }[/math]—possibly an improper one—that favors distributions [math]\displaystyle{ \theta\,\! }[/math] where that rule achieves low risk). Thus, in frequentist decision theory it is sufficient to consider only (generalized) Bayes rules.

Conversely, while Bayes rules with respect to proper priors are virtually always admissible, generalized Bayes rules corresponding to improper priors need not yield admissible procedures. Stein's example is one such famous situation.

Examples

The James–Stein estimator is a nonlinear estimator of the mean of Gaussian random vectors and can be shown to dominate the ordinary least squares technique with respect to a mean-squared-error loss function.[2] Thus least squares estimation is not an admissible estimation procedure in this context. Some others of the standard estimates associated with the normal distribution are also inadmissible: for example, the sample estimate of the variance when the population mean and variance are unknown.[3]

Notes

  1. Dodge, Y. (2003) The Oxford Dictionary of Statistical Terms. OUP. ISBN:0-19-920613-9 (entry for admissible decision function)
  2. Cox & Hinkley 1974, Section 11.8
  3. Cox & Hinkley 1974, Exercise 11.7

References

  • Cox, D. R.; Hinkley, D. V. (1974). Theoretical Statistics. Wiley. ISBN 0-412-12420-3. 
  • Berger, James O. (1980). Statistical Decision Theory and Bayesian Analysis (2nd ed.). Springer-Verlag. ISBN 0-387-96098-8. 
  • DeGroot, Morris (2004). Optimal Statistical Decisions. Wiley Classics Library. ISBN 0-471-68029-X. 
  • Robert, Christian P. (1994). The Bayesian Choice. Springer-Verlag. ISBN 3-540-94296-3.