Randomised decision rule

Selection of randomised decision rules

As with nonrandomised decision rules, randomised decision rules may satisfy favourable properties such as admissibility, minimaxity and Bayes. This shall be illustrated in the case of a finite decision problem, i.e. a problem where the parameter space is a finite set of, say, ${\displaystyle k}$ elements. The risk set, henceforth denoted as ${\displaystyle \mathcal{S}}$, is the set of all vectors in which each entry is the value of the risk function associated with a randomised decision rule under a certain parameter: it contains all vectors of the form ${\displaystyle (R({\theta }_1,d^{\ast }),...R({\theta }_k,d^{\ast })),d^{\ast }\in {\mathcal{D}}^{\ast }}$. Note that by the definition of the randomised decision rule, the risk set is the convex hull of the risks ${\displaystyle (R({\theta }_1,d),...R({\theta }_k,d)),d\in \mathcal{D}}$.

In the case where the parameter space has only two elements ${\displaystyle {\theta }_1}$ and ${\displaystyle {\theta }_2}$, this constitutes a subset of ${\displaystyle {\mathbb{R}}^2}$, so it may be drawn with respect to the coordinate axes ${\displaystyle R_1}$ and ${\displaystyle R_2}$ corresponding to the risks under ${\displaystyle {\theta }_1}$ and ${\displaystyle {\theta }_2}$ respectively. An example is shown on the right.

Admissibility

An admissible decision rule is one that is not dominated by any other decision rule, i.e. there is no decision rule that has equal risk as or lower risk than it for all parameters and strictly lower risk than it for some parameter. In a finite decision problem, the risk point of an admissible decision rule has either lower x-coordinates or y-coordinates than all other risk points or, more formally, it is the set of rules with risk points of the form ${\displaystyle (a,b)}$ such that ${\displaystyle \{(R_1,R_2):R_1\le a,R_2\le b\}\cap \mathcal{S}=(a,b)}$. Thus the left side of the lower boundary of the risk set is the set of admissible decision rules.

Minimax

A minimax Bayes rule is one that minimises the supremum risk ${\displaystyle \underset{\theta \in \mathrm{\Theta }}{sup}R(\theta ,d^{\ast })}$ among all decision rules in ${\displaystyle {\mathcal{D}}^{\ast }}$. Sometimes, a randomised decision rule may perform better than all other nonrandomised decision rules in this regard.

In a finite decision problem with two possible parameters, the minimax rule can be found by considering the family of squares ${\displaystyle Q(c)=\{(R_1,R_2):0\le R_1\le c,0\le R_2\le c\}}$. The value of ${\displaystyle c}$ for the smallest of such squares that touches ${\displaystyle \mathcal{S}}$ is the minimax risk and the corresponding point or points on the risk set is the minimax rule.

If the risk set intersects the line ${\displaystyle R_1=R_2}$, then the admissible decision rule lying on the line is minimax. If ${\displaystyle R_2>R_1}$ or ${\displaystyle R_1>R_2}$ holds for every point in the risk set, then the minimax rule can either be an extreme point (i.e. a nonrandomised decision rule) or a line connecting two extreme points (nonrandomised decision rules).

Bayes

A randomised Bayes rule is one that has infimum Bayes risk ${\displaystyle r(\pi ,d^{\ast })}$ among all decision rules. In the special case where the parameter space has two elements, the line ${\displaystyle {\pi }_1{R}_1+(1-{\pi }_1)R_2=c}$, where ${\displaystyle {\pi }_1}$ and ${\displaystyle {\pi }_2}$ denote the prior probabilities of ${\displaystyle {\theta }_1}$ and ${\displaystyle {\theta }_2}$ respectively, is a family of points with Bayes risk ${\displaystyle c}$. The minimum Bayes risk for the decision problem is therefore the smallest ${\displaystyle c}$ such that the line touches the risk set. This line may either touch only one extreme point of the risk set, i.e. correspond to a nonrandomised decision rule, or overlap with an entire side of the risk set, i.e. correspond to two nonrandomised decision rules and randomised decision rules combining the two. This is illustrated by the three situations below:

As different priors result in different slopes, the set of all rules that are Bayes with respect to some prior are the same as the set of admissible rules.

Note that no situation is possible where a nonrandomised Bayes rule does not exist but a randomised Bayes rule does. The existence of a randomised Bayes rule implies the existence of a nonrandomised Bayes rule. This is also true in the general case, even with infinite parameter space, infinite Bayes risk, and regardless of whether the infimum Bayes risk can be attained. This supports the intuitive notion that the statistician need not utilise randomisation to arrive at statistical decisions.

In practice

As randomised Bayes rules always have nonrandomised alternatives, they are unnecessary in Bayesian statistics. However, in frequentist statistics, randomised rules are theoretically necessary under certain situations, and were thought to be useful in practice when they were first invented: Egon Pearson forecast that they 'will not meet with strong objection'. However, few statisticians actually implement them nowadays.

Randomised test

Randomized tests should not be confused with permutation tests.

In the usual formulation of the likelihood ratio test, the null hypothesis is rejected whenever the likelihood ratio ${\displaystyle \mathrm{\Lambda }}$ is smaller than some constant ${\displaystyle K}$, and accepted otherwise. However, this is sometimes problematic when ${\displaystyle \mathrm{\Lambda }}$ is discrete under the null hypothesis, when ${\displaystyle \mathrm{\Lambda }=K}$ is possible.

A solution is to define a test function ${\displaystyle \varphi (x)}$, whose value is the probability at which the null hypothesis is accepted:

${\displaystyle \varphi (x)=\{\begin{array}{ll}1& \text{ if }\mathrm{\Lambda }>K\\ p(x)& \text{ if }\mathrm{\Lambda }=K\\ 0& \text{ if }\mathrm{\Lambda }<K\end{array}}$

This can be interpreted as flipping a biased coin with a probability ${\displaystyle p(x)}$ of returning heads whenever ${\displaystyle \mathrm{\Lambda }=k}$ and rejecting the null hypothesis if a heads turns up.

A generalised form of the Neyman–Pearson lemma states that this test has maximum power among all tests at the same significance level ${\displaystyle \alpha }$, that such a test must exist for any significance level ${\displaystyle \alpha }$, and that the test is unique under normal situations.

As an example, consider the case where the underlying distribution is Bernoulli with probability ${\displaystyle p}$, and we would like to test the null hypothesis ${\displaystyle p\le \lambda }$ against the alternative hypothesis ${\displaystyle p>\lambda }$. It is natural to choose some ${\displaystyle k}$ such that ${\displaystyle P(\hat{p}>k|H_0)=\alpha }$, and reject the null whenever ${\displaystyle \hat{p}>k}$, where ${\displaystyle \hat{p}}$ is the test statistic. However, to take into account cases where ${\displaystyle \hat{p}=k}$ , we define the test function:

${\displaystyle \varphi (x)=\{\begin{array}{ll}1& \text{ if }\hat{p}>k\\ \gamma & \text{ if }\hat{p}=k\\ 0& \text{ if }\hat{p}<k\end{array}}$

where ${\displaystyle \gamma }$ is chosen such that ${\displaystyle P(\hat{p}>k|H_0)+\gamma P(\hat{p}=k|H_0)=\alpha }$.

Randomised confidence intervals

An analogous problem arises in the construction of confidence intervals. For instance, the Clopper-Pearson interval is always conservative because of the discrete nature of the binomial distribution. An alternative is to find the upper and lower confidence limits ${\displaystyle U}$ and ${\displaystyle L}$ by solving the following equations:

${\displaystyle \{\begin{array}{ll}Pr(\hat{p}<k|p=U)+\gamma P(\hat{p}=k|p=U)& =\alpha /2\\ Pr(\hat{p}>k|p=L)+\gamma P(\hat{p}=k|p=L)& =\alpha /2\end{array}}$

where ${\displaystyle \gamma }$ is a uniform random variable on (0, 1).