⚛ Quantum Self-Study · Technology Engineering
Home / Signal Processing / Modules
Module 21Barkat, chapter 5, pp. 289 to 338

Statistical decision theory

The Bayes criterion and the likelihood ratio, M-ary hypothesis testing and the MAP criterion, the minimax criterion when the priors are unknown, the Neyman–Pearson criterion and ROC curves, composite hypothesis testing and sequential detection (Barkat, chapter 5).

Lecture (slides) Hands-on with the Jupyter Notebook Self-check quiz

Lecture (slides)

MODULE 21 · Barkat, chapter 5, pp. 289 to 338

Statistical decision theory

  • Set up the likelihood ratio test from the Bayes criterion, compute the threshold \eta from costs and priors, and recognise a sufficient statistic.
  • Extend to M-ary hypothesis testing, and derive the maximum a posteriori (MAP) criterion for 0–1 costs.
  • Apply the minimax criterion when the priors are unknown, and interpret the linear Bayes-risk-against-P_1 plot.
  • Derive the Neyman–Pearson test via a Lagrange multiplier, construct ROC curves, and relate the ROC slope to the threshold \eta.
  • Handle composite hypotheses whether the parameter is random or not, recognise a uniformly most powerful (UMP) test, and set up Wald's sequential detection.

Use the buttons, the dots or the ← → keys to move between slides.

PART 1/5
01

The Bayes criterion and the likelihood ratio test

Situation: The binary decision problem: two hypotheses H_0,H_1, costs C_{ij}, and the goal of minimising the average risk \Re=E[C].
Complication: Minimising the risk leads to the likelihood ratio test \Lambda(y)=f_{Y|H_1}/f_{Y|H_0} against a threshold \eta fixed by the costs and priors.
Resolution: A sufficient statistic summarises all the data relevant to the decision; for independent Gaussian noise, the sum of the samples is a sufficient statistic.
  • The four cases and the costs C_{ij}
  • The likelihood ratio and the threshold \eta
  • A Gaussian example and the sufficient statistic
1/5 · The Bayes criterion and the likelihood ratio test

The four outcomes of a binary decision

Each time a decision is made, based on some criterion, for this binary hypothesis testing problem, four possible cases can occur: decide H_0 when H_0 is true; decide H_0 when H_1 is true; decide H_1 when H_0 is true; decide H_1 when H_1 is true. For the first and fourth cases the receiver makes a correct decision, while for the middle two cases it makes an error. From radar nomenclature, the second case is called a miss, the third a false alarm, and the fourth a detection (Barkat, pp. 289 to 291).

1/5 · The Bayes criterion and the likelihood ratio test

Costs, the average risk, and the Bayes goal

The Bayes criterion needs two assumptions: knowledge of the priors P_0,P_1 of the two hypotheses (with P_0+P_1=1), and a cost C_{ij} assigned to each decision D_i when hypothesis H_j is true. The average risk \\Re=E[C] is to be minimised, under the reasonable assumption that a wrong decision costs more than a correct one: C_{01}>C_{11} and C_{10}>C_{00} (Barkat, p. 291). Measured: with C_{00}{=}C_{11}{=}0, C_{01}{=}C_{10}{=}1, P_0{=}P_1{=}0.5 and the error probabilities from the Gaussian example below, the computed risk is 0.1587, exactly the average error probability.

The average Bayes risk
\Re=C_{00}P(D_0,H_0)+C_{01}P(D_0,H_1)+C_{10}P(D_1,H_0)+C_{11}P(D_1,H_1)
1/5 · The Bayes criterion and the likelihood ratio test

The likelihood ratio test and the threshold \eta

Rewriting the risk in terms of the decision region Z_0, the only varying part of the integral is Z_0; minimising the risk means assigning a point y to Z_1 (decide H_1) when P_1(C_{01}-C_{11})f_{Y|H_1}(y|H_1) > P_0(C_{10}-C_{00})f_{Y|H_0}(y|H_0), and otherwise to Z_0. This test is written compactly with the likelihood ratio \\Lambda(y)=f_{Y|H_1}(y|H_1)/f_{Y|H_0}(y|H_0) against the threshold \\eta=P_0(C_{10}-C_{00})/[P_1(C_{01}-C_{11})] (Barkat, pp. 293 to 295). Because the natural logarithm is monotonically increasing, an equivalent test uses \\ln\\Lambda(y) against \\ln\\eta. Measured: with the costs and priors above, \\eta = 1.

The likelihood ratio test
\Lambda(y)=\frac{f_{Y|H_1}(y|H_1)}{f_{Y|H_0}(y|H_0)}\ \underset{H_0}{\overset{H_1}{\gtrless}}\ \eta
1/5 · The Bayes criterion and the likelihood ratio test

Example 5.1: a constant in additive Gaussian noise

Consider a source whose output under H_1 is a constant voltage m, while under H_0 it is zero, with the received signal corrupted by additive zero-mean Gaussian noise of variance \sigma^2: H_1{:}\,Y=m+N, H_0{:}\,Y=N. The likelihood ratio test reduces to comparing the observation y with the threshold \gamma=(\sigma^2/m)\ln\eta+m/2 (Barkat, pp. 296 to 297). With m=2, \sigma=1, \eta=1 (equal priors), the threshold \gamma = 1. The false-alarm probability P_F=Q(\gamma/\sigma) = 0.1587, the detection probability P_D=Q((\gamma-m)/\sigma) = 0.8413, and the miss probability P_M=1-P_D = 0.1587.

1/5 · The Bayes criterion and the likelihood ratio test

K independent samples: sum then threshold

If K independent samples Y_1,\ldots,Y_K of the same problem are taken, the joint density is the product of the individual densities, and the likelihood ratio test reduces to comparing the sum \sum_{k=1}^K y_k with the threshold (\sigma^2/m)\ln\eta+Km/2 (Barkat, pp. 298 to 299). This means the receiver need only add the K samples and compare against a single threshold — it need not retain the individual samples. Measured: a Monte Carlo simulation with K=5 random samples under each hypothesis finds the decision based on the sum \sum y_k agrees with the decision based on the full likelihood ratio in 1 of all trials.

1/5 · The Bayes criterion and the likelihood ratio test

The sufficient statistic: definition and a Bernoulli example

A statistic is any random variable that can be computed from the observed data. A statistic T is called sufficient if, once T is known, the conditional density of any other statistic T' no longer depends on the parameter being tested or estimated — meaning T has summarised all the data relevant to it (Barkat, p. 300). For K independent Bernoulli samples of parameter p, the joint density f_Y(y,p)=p^{\sum y_k}(1-p)^{K-\sum y_k} depends on the data only through the sum \sum y_k; hence T(Y)=\sum_{k=1}^KY_k is a sufficient statistic for p (Barkat, p. 301). Measured: with K=6 Bernoulli samples and p=0.4, three different binary sequences sharing the same sum (3) numerically give the same relative joint-density value (deviating by 1.25e-16), confirming the density depends only on the sum.

1/5 · The Bayes criterion and the likelihood ratio test

When the variances differ: the sufficient statistic is a sum of squares

If f_{Y|H_0} and f_{Y|H_1} are both zero-mean Gaussian but with different variances, \sigma_0^2<\sigma_1^2, the likelihood ratio test leads to T(Y)=Y^2 against a threshold \gamma; with K independent samples, the sufficient statistic becomes T(Y)=\sum_{k=1}^KY_k^2 (Barkat, pp. 301 to 303). This is an example where the sufficient statistic is not a simple sum but a sum of squares, unlike the Example 5.1/5.2 case. Measured: with \sigma_0^2=1, \sigma_1^2=4, the threshold \gamma=\dfrac{2\sigma_0^2\sigma_1^2}{\sigma_1^2-\sigma_0^2}\ln\eta at \eta=1 gives 0 (zero, since \ln1=0); at \eta=2, \gamma = 1.8484.

1/5 · The Bayes criterion and the likelihood ratio test

Self-check, part 1

  • Why does the Bayes risk \Re depend only on the choice of Z_0, not directly on Z_1?
  • From which quantities is the threshold \eta in the likelihood ratio test determined?
  • What does a sufficient statistic summarise about the observed data?

Hints: because Z_1=Z-Z_0 so only one region varies independently; from the costs C_{ij} and the priors P_0,P_1; all the data information relevant to the decision or the parameter sought.

PART 2/5
02

M-ary hypothesis testing and the MAP criterion

Situation: With M hypotheses there are M^2 possible alternatives and M decision regions Z_0,\ldots,Z_{M-1}; the Bayes risk generalises directly.
Complication: With 0–1 costs (C_{ii}=0, C_{ij}=1), minimising the risk is equivalent to maximising the posterior P(H_i|Y): this is the MAP criterion.
Resolution: With M=3, the decision space is two-dimensional (two likelihood ratios \Lambda_1,\Lambda_2), and the three regions Z_0,Z_1,Z_2 are set by pairwise comparisons of I_i(y).
  • The general Bayes risk for M hypotheses
  • 0–1 costs and the MAP criterion
  • Three hypotheses: a numerical example
2/5 · M-ary hypothesis testing and the MAP criterion

From two to M hypotheses

With M hypotheses H_0,\ldots,H_{M-1}, each trial can have M^2 (decision, true hypothesis) pairs, each carrying a cost C_{ij}. The observation space Z is divided into M disjoint subregions Z_0,\ldots,Z_{M-1}, and the average risk is the sum over all pairs \Re=\sum_{i}\sum_j P_jC_{ij}\int_{Z_i}f_{Y|H_j}(y|H_j)\,dy (Barkat, pp. 303 to 304). Minimising it assigns each point y to the region Z_i for which I_i(y)=\sum_{j\ne i}P_j(C_{ij}-C_{jj})f_{Y|H_j}(y|H_j) is smallest (Barkat, pp. 304 to 305).

2/5 · M-ary hypothesis testing and the MAP criterion

0–1 costs: minimising risk becomes maximising the posterior

In communication problems it is common to set C_{ii}=0 (a correct decision costs nothing) and C_{ij}=1 for i\ne j (every wrong decision costs the same). Then minimising the risk is equivalent to minimising I_i(y)=[1-P(H_i|Y)]f_Y(y), which is equivalent to maximising the posterior P(H_i|Y) (Barkat, p. 305). The receiver computes the posteriors P(H_i|Y), i=0,\ldots,M-1, and picks the hypothesis with the largest value — this is called the maximum a posteriori (MAP) receiver. Measured: for three Gaussian hypotheses of means 0,1,2 (variance 1, equal priors), at observation y=1.5 the three posteriors are 0.1554, 0.4223, 0.4223, summing to 1.

The MAP criterion
\text{chọn } H_i \text{ sao cho } P(H_i|Y)=\max_j P(H_j|Y)
2/5 · M-ary hypothesis testing and the MAP criterion

Three hypotheses: verifying MAP against the direct risk

With M=3 and 0–1 costs, instead of computing the posteriors we can directly compare the densities f_{Y|H_j}(y|H_j)P_j for the three hypotheses and pick the largest — the two routes must give the same decision because P(H_i|Y)\propto f_{Y|H_i}(y|H_i)P_i (Bayes' rule). Measured: scanning y continuously and comparing the decision from the \"largest posterior\" rule against the decision from the \"largest prior-weighted density\" rule over 2001 grid points, the agreement rate is 1.

2/5 · M-ary hypothesis testing and the MAP criterion

The decision boundary sits midway between two peaks

With three equal-variance, equal-prior Gaussians of means 0,1,2, the two decision boundaries are the points where two neighbouring densities are equal — by symmetry, exactly the midpoints between the means, i.e. y=0.5 and y=1.5 (Barkat, pp. 307 to 308, the Fig. 5.6 illustration). Measured: numerically solving f_{Y|H_0}(y|H_0)=f_{Y|H_1}(y|H_1) gives y = 0.5, and f_{Y|H_1}(y|H_1)=f_{Y|H_2}(y|H_2) gives y = 1.5, exactly matching the symmetric prediction.

2/5 · M-ary hypothesis testing and the MAP criterion

Two likelihood ratios suffice for a three-hypothesis test

With M=3, there are only two independent likelihood ratios, \Lambda_1(y)=f_{Y|H_1}/f_{Y|H_0} and \Lambda_2(y)=f_{Y|H_2}/f_{Y|H_0}, and the decision space is the two-dimensional (\Lambda_1,\Lambda_2)-plane rather than three-dimensional (Barkat, p. 308). The regions Z_0, Z_1, Z_2 are set by pairwise comparisons of I_i(y) as in (5.51)–(5.52). Measured: at y=1.5 (the theoretical boundary between H_1 and H_2), \Lambda_1(1.5) = 2.7183 and \Lambda_2(1.5) = 2.7183; since neighbouring hypotheses have equal f\cdot P at the boundary, I_1(y)=I_2(y) there, differing by only 1.00e-16.

2/5 · M-ary hypothesis testing and the MAP criterion

MAP-test risk against an off-centre test

Since the MAP test minimises the risk (the average error probability under 0–1 costs), any other choice of decision boundary — say shifting both boundaries right by 0.3 — must give an equal or larger risk. Measured: the average error probability of the correctly centred MAP test (y=0.5,1.5) is 0.4114; the error probability of the shifted test (0.3) is 0.4217, exceeding the MAP test by 0.0103, confirming MAP is optimal.

2/5 · M-ary hypothesis testing and the MAP criterion

Connection to the parameter-estimation module

The MAP criterion here chooses discretely among M known hypotheses; in the next chapter on parameter estimation (module 22), a closely related idea — MAP estimation — will choose a continuous value of a parameter \theta to maximise the posterior density f_{\Theta|Y}(\theta|y), instead of choosing among a finite number of hypotheses. Both rest on the same Bayes' rule inverting f_{Y|\Theta} into f_{\Theta|Y}, differing only in whether \theta ranges over a discrete or a continuous set (Barkat, p. 305, connecting to chapter 6).

2/5 · M-ary hypothesis testing and the MAP criterion

Self-check, part 2

  • With 0–1 costs, minimising the risk is equivalent to maximising what quantity?
  • With M=3 hypotheses, how many dimensions does the decision space have?
  • Where does the decision boundary between two equal-variance, equal-prior Gaussians lie?

Hints: the posterior probability P(H_i|Y); two dimensions (two independent likelihood ratios); exactly midway between the two means.

PART 3/5
03

The minimax criterion

Situation: When the priors are unknown, minimax picks the P_1 that maximises the risk, then designs the Bayes test minimising that risk — a "minimum of the maximum."
Complication: The Bayes risk \Re(P_1) is linear in P_1; the minimax equation sets the derivative to zero, giving (C_{11}{-}C_{00})+(C_{01}{-}C_{11})P_M-(C_{10}{-}C_{00})P_F=0.
Resolution: With 0–1 costs, the equation reduces to P_F=P_M: the minimax threshold is the one equalising the two error probabilities.
  • Why minimax is needed
  • The risk linear in P_1
  • Example: the same result whether or not P_1 is known
3/5 · The minimax criterion

When there is no reliable prior

The Bayes criterion needs knowledge of the priors; in many situations this information is not available. One approach is to choose the value of P_1 that maximises the risk, then minimise that risk — the principle of minimising the largest average cost over the chosen P_1 is called the minimax criterion (Barkat, p. 313). Substituting P_0=1-P_1 into the risk formula (5.15), the risk becomes a linear function of P_1: \Re(P_1)=C_{00}(1-P_F)+C_{10}P_F+P_1[(C_{11}-C_{00})+(C_{01}-C_{11})P_M-(C_{10}-C_{00})P_F] (Barkat, pp. 313 to 314).

3/5 · The minimax criterion

The two extremes: P_1=0 and P_1=1

If P_1=0, the threshold becomes infinite and H_0 is always decided; then P_F=0, P_M=1, and the risk reduces to \Re=C_{00}. If P_1=1, the threshold is zero and H_1 is always decided; then P_F=1, P_M=0, and the risk reduces to \Re=C_{11} (Barkat, pp. 314 to 315). The plot of risk against P_1 is therefore a straight line connecting (0,C_{00}) and (1,C_{11}) geometrically at the two extremes, but the true Bayes curve (optimal P_1 at each point) must be concave below that line. Measured: with C_{00}=0, C_{11}=0 (0–1 costs), the two extreme values of \Re(P_1) are exactly 0 and 0.

3/5 · The minimax criterion

The minimax equation and the 0–1-cost case

Differentiating the risk with respect to P_1 and setting it to zero at the optimal P_1^* gives the minimax equation (C_{11}-C_{00})+(C_{01}-C_{11})P_M-(C_{10}-C_{00})P_F=0 (Barkat, p. 315). If the cost of a correct decision is zero (C_{00}=C_{11}=0), the equation reduces to C_{01}P_M=C_{10}P_F; and if in addition the cost of a wrong decision is one (C_{01}=C_{10}=1), then P_F=P_M — the minimax threshold is exactly the one that equalises the two error probabilities (Barkat, p. 316). Measured: for the Gaussian example m=2, \sigma=1, the minimax threshold \gamma^*=m/2 = 1; there P_F = 0.1587 and P_M = 0.1587, agreeing to within 0.00e+00.

The minimax equation
(C_{11}-C_{00})+(C_{01}-C_{11})P_M-(C_{10}-C_{00})P_F=0
3/5 · The minimax criterion

Example 5.6: knowing or not knowing P_0,P_1 gives the same result

In the constant-source-in-Gaussian-noise example, part (a) assumes P_0=P_1=1/2 is known, giving the threshold \gamma=m/2 and average error probability P(\varepsilon)=Q(m/2\sigma). Part (b) assumes P_0,P_1 are completely unknown, using minimax with the condition P_F=P_M: this also gives the threshold \gamma^*=m/2 and the same error probability Q(m/2\sigma) (Barkat, pp. 315 to 317). This is a notable result: for this symmetric problem, not knowing the priors costs nothing in performance. Measured: with m=2, \sigma=1, both routes give P(\varepsilon) = 0.1587, differing by only 0.00e+00.

3/5 · The minimax criterion

Minimax is not always "free"

The result of Example 5.6 is special because of the problem's symmetry (two equal-variance Gaussian densities, symmetric about the midpoint). In general, when a test designed for a fixed P_1 is applied under a different actual P_1, the actual risk exceeds the optimal Bayes risk for that actual P_1 — this is exactly why the Bayes risk curve must be concave below the line joining the two extreme points. Measured: using the threshold designed for P_1=0.5 (i.e. \gamma=m/2=1) but applied when the actual P_1 is 0.8: the actual risk is 0.1587, larger than the optimal Bayes risk for P_1=0.8, which is 0.1121.

3/5 · The minimax criterion

Reading the risk-against-P_1 plot

On the plot of \Re against P_1 (the book's Fig. 5.10), the tangent at P_1^* is horizontal because that is the maximum of the concave Bayes curve — meaning the derivative is exactly zero at the minimax point. The minimax risk \Re^*(P_1^*) is exactly the largest value of the optimal Bayes risk over the whole range P_1\in[0,1] (Barkat, p. 315). Measured: scanning P_1 from 0 to 1, computing the optimal Bayes risk (the threshold \gamma(P_1) optimal for each P_1) and finding the P_1 that maximises it; the scanned result is P_1^* = 0.5, matching the symmetry-predicted 0.5 to within the grid error 0.

3/5 · The minimax criterion

Self-check, part 3

  • Why is minimax needed instead of Bayes in some practical situations?
  • With 0–1 costs, what does the minimax equation reduce to?
  • In Example 5.6, why does not knowing P_0,P_1 not degrade performance?

Hints: when the priors are not reliably known; P_F=P_M; because the problem's symmetry makes the optimal Bayes threshold coincide with the minimax threshold.

PART 4/5
04

The Neyman–Pearson criterion and ROC curves

Situation: When realistic costs or priors cannot be assigned, Neyman–Pearson fixes P_F=\alpha and maximises P_D via a Lagrange multiplier — still leading to a likelihood ratio test.
Complication: The ROC curve plots P_D against P_F with the threshold as a parameter; it depends only on the conditional densities, not on costs or priors.
Resolution: The slope of the ROC at a point is exactly the threshold \eta needed to reach that (P_F,P_D) point: dP_D/dP_F=\eta.
  • Maximising P_D with P_F fixed
  • A numerical example: exponential and uniform densities
  • The ROC curve and the meaning of its slope
4/5 · The Neyman–Pearson criterion and ROC curves

When costs and priors are not realistic

In many physical situations, such as radar detection, it is very difficult to assign realistic costs and priors. To overcome this, the Neyman–Pearson criterion uses the conditional probabilities P_F and P_D directly: fix P_F=\alpha while maximising P_D (Barkat, p. 317). Since P_M=1-P_D, maximising P_D is equivalent to minimising P_M subject to P_F=\alpha; using constrained-extremum calculus, we form the objective function J=P_M+\lambda(P_F-\alpha) with \lambda\ge0 the Lagrange multiplier (Barkat, pp. 317 to 318).

4/5 · The Neyman–Pearson criterion and ROC curves

Still leads to a likelihood ratio test

Rewriting J in terms of the decision region Z_0, we get J=\lambda(1-\alpha)+\int_{Z_0}[f_{Y|H_1}(y|H_1)-\lambda f_{Y|H_0}(y|H_0)]\,dy; J is minimised when the points with f_{Y|H_1}>f_{Y|H_0} (a negative difference) are assigned to Z_1. This gives exactly the likelihood ratio test \Lambda(y)\gtrless\lambda, with \lambda playing the role of the equivalent Bayes threshold \eta (Barkat, pp. 318 to 319). So the Neyman–Pearson and Bayes criteria use the same form of test — differing only in how the threshold is fixed.

4/5 · The Neyman–Pearson criterion and ROC curves

Example 5.7: two differently-shaped densities

With f_{Y|H_0}(y|H_0)=\dfrac{1}{2(1-e^{-1})}e^{-|y|} and f_{Y|H_1}(y|H_1)=\tfrac12 (both on [-1,1]), equally likely hypotheses. The minimum-error-probability test has a boundary at \pm0.459, deciding H_1 when 0.459\le|y|\le1 (Barkat, pp. 319 to 320). Measured: the exact boundary solved numerically gives 0.4587, matching the book's 0.459; the false-alarm probability at this boundary is P_F = 0.418 and the miss probability P_M = 0.4587, giving average error probability P(\varepsilon) = 0.4383.

4/5 · The Neyman–Pearson criterion and ROC curves

The same example, a Neyman–Pearson threshold with P_F=0.5

If instead of minimising the error probability we fix P_F=0.5 (the Neyman–Pearson criterion), the threshold changes to a \gamma solved from P_F=\int_{-1}^{-\gamma}f_{Y|H_0}\,dy+\int_{\gamma}^{1}f_{Y|H_0}\,dy=0.5, giving \gamma\approx0.38; the corresponding detection probability is P_D=2[(1-\gamma)(1/2)]\approx0.62 (Barkat, pp. 320 to 321). Measured: solving numerically gives \gamma = 0.3799, and P_D = 0.6201, matching the book's values.

4/5 · The Neyman–Pearson criterion and ROC curves

The ROC curve: P_D against P_F with the threshold as a parameter

The ROC curve plots P_D against P_F as the threshold varies; it depends only on the conditional densities f_{Y|H_j}(y|H_j), not on the costs or priors (Barkat, p. 321). For the example of a sum of K Gaussian samples, the normalised sufficient statistic has detection parameter d=\sqrt{K}\,m/\sigma (a signal-to-noise ratio), and P_F=Q(\ln\eta/d+d/2), P_D=Q(\ln\eta/d-d/2) (Barkat, pp. 321 to 322). Measured: with K=4, m=2, \sigma=1, the parameter d = 4; at \ln\eta=0, P_F = 0.0228 and P_D = 0.9772 — P_D far exceeds P_F, showing a large d gives good discrimination.

Figure 1. Gaussian ROC curves with detection parameter d equal to 0.5, 1 and 2: a larger d pulls the curve closer to the top-left corner.
Figure 1. Gaussian ROC curves with detection parameter d equal to 0.5, 1 and 2: a larger d pulls the curve closer to the top-left corner.
4/5 · The Neyman–Pearson criterion and ROC curves

The slope of the ROC is the threshold \eta

Differentiating P_D and P_F with respect to the threshold \eta, one can show dP_D/dP_F=\eta (Barkat, pp. 323 to 325): the tangent slope at a point on the ROC is exactly the Neyman–Pearson threshold needed to reach that (P_F,P_D) point. In the Bayes criterion, \eta is set by the costs and priors, so (P_F,P_D) sits at the point on the ROC whose tangent has that slope \eta. Measured: the numerical derivative dP_D/dP_F at \ln\eta=0 (i.e. \eta=1) gives 1, matching \eta=1; at \ln\eta=1 (i.e. \eta=e\approx2.718), the numerical slope is 2.7183, matching 2.7183.

The ROC slope
\frac{dP_D}{dP_F}=\eta
4/5 · The Neyman–Pearson criterion and ROC curves

Example 5.8: the ROC of two exponential densities

With f_{Y|H_0}(y|H_0)=e^{-y} and f_{Y|H_1}(y|H_1)=\alpha e^{-\alpha y} (y\ge0, \alpha>1), the likelihood ratio gives the threshold \gamma=\tfrac{1}{1-\alpha}\ln(\eta/\alpha), and P_D=1-e^{-\alpha\gamma}, P_F=1-e^{-\gamma} (Barkat, pp. 325 to 326). Measured: with \alpha=2 and \gamma=1, P_D = 0.8647, P_F = 0.6321; the theoretical slope at this point, \eta=\alpha e^{-(\alpha-1)\gamma} (from directly differentiating dP_D/dP_F), gives 0.7358, matching the numerically computed slope dP_D/dP_F of 0.7358.

4/5 · The Neyman–Pearson criterion and ROC curves

Self-check, part 4

  • What quantity does the Neyman–Pearson criterion fix, and what does it maximise?
  • What does the ROC curve depend on, and what does it not depend on?
  • What does the slope of the ROC curve at a point equal?

Hints: fixes P_F=\alpha, maximises P_D; depends on the conditional densities, not on costs or priors; the threshold \eta needed to reach that point.

PART 5/5
05

Composite hypotheses and sequential detection

Situation: A composite hypothesis: one or more unknown parameters. If the parameter is random with a known density, average over it, then apply Bayes as usual.
Complication: If the parameter is nonrandom but unknown, use Neyman–Pearson; a test independent of the true parameter value is called uniformly most powerful (UMP).
Resolution: Wald's sequential detection: samples taken sequentially, the likelihood ratio compared with two thresholds \eta_0,\eta_1; the sample count K is random, with expectation smaller than a fixed-sample-size test of the same performance.
  • A random parameter: average then apply Bayes
  • A nonrandom parameter: the uniformly most powerful test
  • Wald's sequential detection and the expected sample count
5/5 · Composite hypotheses and sequential detection

A composite hypothesis: when the parameter is not fully known

In the problems considered so far, the parameter characterising each hypothesis was known exactly. In many situations this parameter is not known; the hypothesis is then called composite. For example, if H_1{:}\,Y=m+N with m unknown, then H_1 is composite while H_0{:}\,Y=N remains simple (Barkat, pp. 326 to 327). There are two cases: the parameter \Theta is a random variable with a known density, or \theta is an unknown but nonrandom constant.

5/5 · Composite hypotheses and sequential detection

A random parameter: average over its density

If \Theta is random with a known density f_{\Theta|H_j}(\theta|H_j), the Bayes criterion still applies, one need only average the conditional density over the parameter: f_{Y|H_j}(y|H_j)=\int f_{Y|\Theta,H_j}(y|\theta,H_j)f_{\Theta|H_j}(\theta|H_j)\,d\theta (Barkat, p. 327). With M (capitalised, to distinguish it from the constant value m) zero-mean Gaussian of variance \sigma_m^2, the marginal density f_{Y|H_1}(y|H_1) turns out still Gaussian but with the combined variance \sigma_m^2+\sigma^2 (Barkat, pp. 327 to 329). Measured: the numerical integral of \int f_{Y|M,H_1}(y|m,H_1)f_M(m)\,dm at y=1.3 gives 0.1507, matching the closed form \dfrac{1}{\sqrt{2\pi(\sigma_m^2+\sigma^2)}}\exp\!\left(-\dfrac{y^2}{2(\sigma_m^2+\sigma^2)}\right) = 0.1507.

5/5 · Composite hypotheses and sequential detection

A nonrandom parameter: the uniformly most powerful (UMP) test

If \theta is not a random variable but unknown, the Bayes criterion does not apply since there is no density to average over; instead Neyman–Pearson is used, maximising P_D with P_F fixed for each assumed value of \theta. The plot of P_D against \theta is called the power function. A test that maximises P_D for every possible value of \theta, i.e. independent of the true value of \theta, is called uniformly most powerful (UMP) (Barkat, p. 329).

5/5 · Composite hypotheses and sequential detection

Example 5.11: when a UMP test exists, and when it does not

With H_1{:}\,Y=m+N, m unknown but known to be positive, the likelihood ratio test reduces to y\gtrless\gamma_1, with \gamma_1 solved from P_F=\int_{\gamma_1}^\infty f_{Y|H_0}(y|H_0)\,dy — entirely independent of the specific value of m. So a best test exists independent of m, and a UMP test exists (Barkat, pp. 330 to 331). But if m could be positive or negative, the test designed for m>0 differs from that for m<0 (the inequality reverses), so no single test is optimal for every m — a UMP test does not exist in that case (Barkat, pp. 331 to 332). Measured: with \sigma=1 and P_F=0.1, the threshold \gamma_1 solved numerically gives 1.2816, independent of the trial value of m (1 or 3 both give the same \gamma_1, differing by 0).

5/5 · Composite hypotheses and sequential detection

Sequential detection: deciding after each observation

In earlier sections, the number of observations K was always fixed in advance. In sequential detection, after each observation, one of three decisions is made: choose H_1, choose H_0, or not enough information yet (take another observation). The test compares the cumulative likelihood ratio \Lambda(y_K)=\prod_{k=1}^K f_{Y_k|H_1}/f_{Y_k|H_0} with two thresholds \eta_0<\eta_1: if \Lambda\ge\eta_1 choose H_1; if \Lambda\le\eta_0 choose H_0; if in between, take another observation (Barkat, pp. 332 to 333). The number of observations K becomes a random variable.

The two Wald sequential-test thresholds
\eta_1\le\frac{1-\beta}{\alpha},\qquad \eta_0\ge\frac{\beta}{1-\alpha}
5/5 · Composite hypotheses and sequential detection

The expected sample count: Example 5.12

With the sequential-sampling receiver of Example 5.2, \sigma^2=1, m=1, and the target P_F=\alpha=0.1=P_M=\beta, the two log-thresholds are \ln\eta_1=\ln[(1-\beta)/\alpha] and \ln\eta_0=\ln[\beta/(1-\alpha)]; and the expected number of observations under each hypothesis is E[K|H_j]=\dfrac{(1-\beta)\ln\eta_1+\beta\ln\eta_0}{E[L(y)|H_j]} or a similar formula (Barkat, pp. 334 to 337). Measured: \ln\eta_1 = 2.1972, \ln\eta_0 = -2.1972; E[K|H_1] = 3.5156 and E[K|H_0] = 3.5156 — meaning on average about four samples suffice for performance P_F=P_M=0.1, fewer than a larger fixed sample size would need for the same non-sequential performance.

5/5 · Composite hypotheses and sequential detection

Looking back: one thread running through decision theory

From Bayes to MAP, minimax, Neyman–Pearson, and sequential detection, every criterion in this chapter reduces to the same form of likelihood ratio test \Lambda(y)\gtrless a threshold — differing only in how that threshold is set: from costs and priors (Bayes), from the condition P_F=P_M (minimax with 0–1 costs), from a fixed constraint P_F=\alpha (Neyman–Pearson), or from two thresholds changing at each observation step (sequential). This is the chapter's core lesson: the likelihood ratio test is the common structure, and the criterion only decides the threshold.

5/5 · Composite hypotheses and sequential detection

Self-check, part 5

  • How does a composite hypothesis differ from a simple hypothesis?
  • What special property does a uniformly most powerful (UMP) test have?
  • In sequential detection, which random variable is not fixed in advance, unlike a fixed-sample-size test?

Hints: the parameter characterising the hypothesis is not fully known; it maximises P_D for every possible true value of the parameter; the observation count K.

SUMMARY

Takeaways

  • The Bayes criterion minimising the risk \Re=E[C] leads to the likelihood ratio test \Lambda(y)=f_{Y|H_1}/f_{Y|H_0}\gtrless\eta, with \eta set by the costs and priors.
  • With 0–1 costs, minimising the risk is equivalent to maximising the posterior P(H_i|Y): the MAP criterion, extending naturally to M hypotheses.
  • Minimax picks the threshold making P_F=P_M when the priors are unknown; Neyman–Pearson fixes P_F=\alpha and maximises P_D, giving an ROC curve with slope dP_D/dP_F=\eta.
  • Composite hypotheses: a random parameter uses Bayes after averaging; a nonrandom parameter uses Neyman–Pearson, and a test independent of the true parameter is called UMP.
  • Wald's sequential detection compares the cumulative likelihood ratio with two thresholds \eta_0,\eta_1; the sample count K is random, with expectation smaller than a fixed-sample-size test of the same performance.

📜 Theory & historical background

Statistical decision theory originates in the work of Neyman and Pearson in the 1930s on statistical hypothesis testing, and of Wald in the 1940s on sequential analysis (Wald, 1947). Its application to radar signal detection and digital communication in the 1950s–1960s (Van Trees; Helstrom) shaped the theoretical framework presented in this chapter (Barkat, pp. 289, 317, 332 to 338, the end-of-chapter references).

🔎 Real-world case study

One source, five views of the same decision. With H_1{:}Y=m+N, H_0{:}Y=N, m=2, \sigma=1: the Bayes criterion with equal priors gives threshold \gamma = 1, P_F = 0.1587, P_D = 0.8413. The minimax criterion, needing no knowledge of the priors, gives exactly the same threshold \gamma^* = 1 and the same error probability P(\varepsilon) = 0.1587 (differing by only 0.00e+00) — thanks to the problem's symmetry. With K=4 independent samples, the detection parameter d = 4 and the ROC curve at \ln\eta=0 gives P_F = 0.0228, P_D = 0.9772, with tangent slope exactly 1 matching \eta=1. If m is replaced by a random Gaussian parameter of variance \sigma_m^2, the marginal density becomes Gaussian with the combined variance — measured at y=1.3 it matches the closed form to 0.1507 against 0.1507. And with sequential sampling instead of a fixed sample size, only 3.5156 samples on average are needed to reach P_F=P_M=0.1.

🧮 Formula sheet for this module

The average Bayes risk
\Re=C_{00}P(D_0,H_0)+C_{01}P(D_0,H_1)+C_{10}P(D_1,H_0)+C_{11}P(D_1,H_1)
The likelihood ratio test
\Lambda(y)=\frac{f_{Y|H_1}(y|H_1)}{f_{Y|H_0}(y|H_0)}\ \underset{H_0}{\overset{H_1}{\gtrless}}\ \eta
The MAP criterion
\text{chọn } H_i \text{ sao cho } P(H_i|Y)=\max_j P(H_j|Y)
The minimax equation
(C_{11}-C_{00})+(C_{01}-C_{11})P_M-(C_{10}-C_{00})P_F=0
The ROC slope
\frac{dP_D}{dP_F}=\eta
The two Wald sequential-test thresholds
\eta_1\le\frac{1-\beta}{\alpha},\qquad \eta_0\ge\frac{\beta}{1-\alpha}

🧪 Hands-on with the Jupyter Notebook

Steps
  1. Open the notebook and run the setup cell.
  2. Task 1: repeat Example 5.1 with m=3, \sigma=1.5 and recompute \gamma, P_F, P_D.
  3. Task 2: build a MAP test for M=4 Gaussian hypotheses of means 0,1,2,3 and find the three decision boundaries.
  4. Task 3: for the exponential example (Example 5.8), plot the ROC for \alpha=3 and \alpha=5, comparing with \alpha=2.
  5. Task 4: recompute the expected sample count of sequential detection with \alpha=\beta=0.01 instead of 0.1, and compare.

Open in Google Colab Download notebook (.ipynb)

The notebook generates its own data in code, so no external files are needed. Every number on this page is printed by the notebook and cross-checked with two independent methods.

⚠️ Common misreadings

"The minimax criterion always gives the same performance as Bayes with known priors."

True only in symmetric problems like Example 5.6; in general, a test designed for an assumed P_1 then applied under a different actual P_1 gives a risk larger than the optimal Bayes risk — a numerical illustration shows a mismatched risk of 0.1587 against the optimal 0.1121.

"The ROC curve depends on the chosen priors."

No: the ROC depends only on the conditional densities f_{Y|H_j}(y|H_j); the priors and costs only decide which point on the ROC is selected (via the threshold \eta), never the shape of the curve itself.

"A uniformly most powerful (UMP) test always exists for any composite hypothesis."

No: Example 5.11 shows a UMP test exists when m can only be positive (or only negative), but does not exist when m can take either sign, since the optimal tests for m>0 and m<0 have opposite inequalities.

✅ Self-check quiz

Click an answer to grade it immediately and read the explanation. Results are not saved.

References

← Previous module