45.02.04 · mathematical-statistics / 02-hypothesis-testing-confidence

Confidence Sets, Pivotal Quantities, and Test-Interval Duality

shipped3 tiersLean: none

Anchor (Master): Lehmann & Romano 2005 Testing Statistical Hypotheses (3rd ed., Springer) Ch. 3 §5 and Ch. 5 (uniformly most accurate confidence sets from UMP tests, unbiased and shortest sets); Casella & Berger 2002 Statistical Inference §9.3 (evaluating interval estimators: size, coverage, optimality)

Intuition Beginner

A point estimate hands you a single number — the average bolt length is millimeters — and then quietly hopes you forget it is almost certainly a little bit wrong. A confidence interval is the honest version of the same report. Instead of one number it hands you a range, like "between and millimeters," together with a promise about how often a range built this way contains the real value.

The promise is the whole point, and it is easy to misread. A confidence interval does not mean "there is a chance the truth is in this particular range." The truth is a fixed number; it is either in your range or it is not. What is random is the range, because it was computed from random data. The promise is about the recipe: if you repeated the whole experiment many times, building a fresh interval each time by the same rule, about out of every of those intervals would capture the true value. The number describes the long-run hit rate of the procedure, not the status of one interval.

There is a beautiful trick for building these ranges, and it ties two ideas together. Suppose for every candidate value you could run a fair test of "is this the truth?" and either reject it or let it stand. Collect all the candidate values your data does not reject. That collection is your confidence interval. A value survives exactly when it is plausible, and the interval is the set of survivors. Testing and interval-building are two views of one machine.

Visual Beginner

Picture a vertical number line of candidate values for the true average. For each candidate you imagine a small hypothesis test: does the data argue against this value? Some candidates get rejected (the data says "no, not that one"); the rest survive. The survivors form a connected band, and that band is the confidence interval. Slide the candidate up or down and you trace out where rejection starts and stops — the two edges of the band are the interval's endpoints.

candidate value does a 5% test reject it? in the interval?
far below the estimate yes no
just below the estimate no yes
at the estimate no yes
just above the estimate no yes
far above the estimate yes no

The takeaway: the confidence interval is the set of candidate values that a level- test refuses to reject. Building an interval and running a family of tests are the same act, seen from two sides. The contiguous band of survivors is what you report as the range.

Worked example Beginner

A lab weighs the same reference mass times. The scale's noise is known to have a standard deviation of milligrams, and the nine readings average to milligrams. You want a confidence interval for the true mass.

Step 1. The spread of the average. Averaging nine independent readings shrinks the noise: the standard deviation of the average is milligram.

Step 2. The magic quantity that does not depend on the unknown. Form , where is the true mass. Whatever the true mass is, this standardized gap behaves like a standard bell curve — its distribution does not depend on at all. That is what makes it useful.

Step 3. Bracket the bell curve. A standard bell curve puts of its area between and . So with probability , the standardized gap lands in that range: .

Step 4. Solve for the unknown. Rearranging, . Plugging in : the interval is to , that is to milligrams.

What this tells us: the data is consistent with any true mass between about and milligrams, at the confidence level. Report a single number and you hide the uncertainty; report this range and you show exactly how much wiggle the nine readings leave. Take more readings and shrinks, the bracket tightens, and the range narrows.

Check your understanding Beginner

Formal definition Intermediate+

Fix a sample space and a family as in 45.02.01. Interval estimation replaces a point estimator with a set-valued rule.

Definition (confidence set). A confidence set is a measurable set-valued map , assigning to each outcome a subset , with for every . When and is an interval, is a confidence interval, written for measurable endpoint statistics .

Definition (coverage probability, confidence coefficient). The coverage probability of at is

The confidence coefficient of is the infimum coverage, . The set is a confidence set if its confidence coefficient is at least . The frequentist content is exactly this: is fixed and non-random; the randomness is in , hence in the set , and is the long-run fraction of realizations whose set contains the fixed .

Definition (acceptance region). For each , let be a level- test of (in the sense of 45.02.01). Its acceptance region is for a non-randomized test, i.e. the outcomes that do not reject . The family has the level property for every .

Definition (pivotal quantity). A pivotal quantity (or pivot) is a function of the data and the parameter whose distribution under is the same for every . For a location family the quantity is a pivot; for a scale family the quantity is a pivot; for the normal model with the sample mean and sample standard deviation, is the pivot with a law, and is the pivot with a law.

Construction (the pivot method). Given a pivot , choose constants with (the same for all , since the law of is parameter-free). Then is a confidence set, because for every .

Counterexamples to common slips

  • Coverage is the worst case over , not the value at one point. A set can have coverage at the true parameter you happen to consider yet drop below elsewhere; the confidence coefficient is the infimum, and a procedure advertised as must clear at every .
  • The confidence level is not the probability the parameter lies in the realized interval. Once is observed, is a fixed set and is a fixed number; the statement "" is a category error in the frequentist framework. The probability lives in the procedure, before the data are seen.
  • A pivot is not the same as a statistic. A pivot depends on ; a statistic does not. The point of the pivot is precisely that it mixes data and parameter so that the dependence on cancels out of its distribution, which is what lets you invert the probability bracket into a parameter set.

Key theorem with proof Intermediate+

The organizing result is the exact correspondence between tests and confidence sets: inverting a family of level- acceptance regions produces a confidence set, and slicing a confidence set produces a family of level- tests. The two objects carry the same information.

Theorem (test-inversion duality). Let be a family of acceptance regions of level- tests of , so for every . Define

Then is a confidence set: for every . Conversely, if is a confidence set, then the tests with acceptance regions have level , and the two constructions are mutually inverse.

Proof. Membership unfolds symmetrically: for any and any ,

directly from the definition of . Fix an arbitrary and take in this equivalence; intersecting with the event over the random ,

Apply to both sides:

the inequality being the level property of the test at . Since was arbitrary, , so is a confidence set.

For the converse, suppose has confidence coefficient and set . Then the type I error of the corresponding test at is

so the test has level . Applying the forward map to this acceptance family returns , and applying the backward map to an acceptance family returns ; the two maps are inverse.

Bridge. This duality builds toward the optimal-interval theory of the Master section, where the optimality of a test transfers to the dual set, and it appears again in the large-sample Wald, score, and likelihood-ratio intervals of 45.04.09, each the inversion of an asymptotic level- test. This is exactly the same membership bookkeeping that the pivot method exploits: a pivot supplies, for every , the acceptance region , and inverting that family is the pivot construction read through the duality, so the pivot method is dual to a family of pivot-based tests. The central insight is that a confidence set and a family of tests are one object indexed two ways — by outcomes giving sets of parameters, or by parameters giving sets of outcomes — and the foundational reason the coverage equals one minus the size is the pointwise identity . Putting these together, every method of finding tests (Neyman-Pearson thresholds, likelihood ratios, score statistics) becomes a method of finding intervals, and every optimality criterion for tests generalises to an accuracy criterion for sets; the bridge is the inversion map .

Exercises Intermediate+

Advanced results Master

Uniformly most accurate sets from uniformly most powerful tests

False coverage is the interval analogue of type II error. For a confidence set and a wrong value , the false-coverage probability is , the chance the set captures a value that is not the truth. A set is uniformly most accurate (UMA) if it minimizes over all sets, simultaneously for every pair in the relevant range. The duality theorem turns UMA-ness into a statement about tests: by the identity , minimizing false coverage at is minimizing the acceptance probability , that is, maximizing the power of the dual test of at the alternative . Inverting a family of uniformly most powerful tests therefore yields a UMA confidence set [Lehmann & Romano §3.5]. In a one-parameter model with monotone likelihood ratio, the one-sided UMP tests of 45.02.02 invert to the UMA one-sided confidence bounds; for two-sided problems no UMP test exists, and the construction is repaired by restricting to unbiased tests.

Unbiased sets and the shortest-length problem

A confidence set is unbiased if its false-coverage probability never exceeds its coverage: for all , so the set is at least as likely to cover the truth as any falsehood. Unbiased sets are exactly the inversions of unbiased tests, and inverting the UMP-unbiased tests produces UMA-unbiased sets, the canonical optimum for two-sided normal problems where no unrestricted UMP test exists [Lehmann & Romano Ch. 5]. Shortest length is a different criterion and need not coincide with UMA-unbiasedness. By the Pratt identity (Exercise 8), expected length equals integrated false coverage, so a UMA-unbiased set minimizes expected length within the unbiased class whenever the dominated family makes the two criteria agree; the equal-tailed interval for the normal variance is not the shortest, and the minimum-length variance interval solves the stationarity at the pivot density rather than splitting the tail probability evenly.

The relation of the Cramér-Rao bound to interval width

The information floor on variance is an information floor on length. For a regular model the score-based and likelihood-based intervals of 45.04.09 have asymptotic half-width governed by the inverse Fisher information: the Wald interval is , whose squared half-width scales as . Because the Cramér-Rao lower bound of 45.01.05 caps the precision of any unbiased estimator at , no asymptotically-efficient interval centered at an unbiased estimator can be narrower in expectation than the scale permits; the inverse information is simultaneously the variance floor and, through the Pratt identity, the expected-length floor. Exact finite-sample attainment of both floors is the exponential-family rigidity of 45.01.05: only there does an affine-in-the-score statistic give both a minimum-variance estimator and the corresponding minimal-length pivot interval.

Bayesian credible sets contrasted with frequentist confidence sets

A Bayesian credible set at level satisfies , where is the posterior; the probability statement is over , conditional on the observed . A frequentist confidence set's statement is over , conditional on a fixed . These are distinct probabilities on distinct sample spaces: the credible set answers "given this data, where is the parameter likely to be?" while the confidence set answers "does this procedure, run repeatedly, capture the fixed parameter at the advertised rate?". The highest-posterior-density credible set is the Bayesian shortest-set analogue, minimizing Lebesgue measure at fixed posterior mass exactly as the condition minimizes length at fixed coverage. The two coincide numerically in special conjugate cases — a flat prior on a normal mean returns the credible set , identical to the frequentist interval — but the agreement is a coincidence of that model, and the interpretations remain incompatible: the frequentist quantity is calibrated over hypothetical repetitions, the Bayesian quantity over the prior-to-posterior update, and a credible set need not have frequentist coverage nor conversely.

Synthesis. The central insight is that a confidence set is a family of tests indexed the other way, so the entire optimality theory of testing maps onto interval estimation through one membership identity, . This is exactly why coverage is one minus size and why false coverage is one minus power: the dual test's two error rates are the set's coverage and accuracy, and the foundational reason uniformly most accurate sets exist precisely where uniformly most powerful tests exist is that inversion is a bijection preserving these probabilities pointwise. Putting these together, the one-sided MLR theory of 45.02.02 generalises to UMA one-sided bounds, the two-sided unbiased theory to UMA-unbiased intervals, and the asymptotic tests of 45.04.09 to the Wald, score, and likelihood-ratio intervals, every optimality crossing the duality intact. The Pratt identity is dual to the power calculation — expected length is integrated false coverage — so short intervals and powerful tests are one demand, and the Cramér-Rao bound of 45.01.05 caps both the estimator variance and, through that identity, the achievable expected length, with exact attainment confined to the exponential family. The Bayesian credible set sits beside this structure as a different probability over a different space, numerically tangent in conjugate-flat-prior models yet interpretively distinct; the bridge across the whole chapter is that inversion turns optimal tests into optimal sets, so testing and estimation are not two subjects but one viewed through two indices.

Full proof set Master

Proposition 1 (the inversion bijection). The map , , is a bijection between families of acceptance regions and set-valued maps , with inverse , . It carries level- test families to confidence sets and back.

Proof. Both maps are defined by the same biconditional , read in the two directions, so each is determined by the relation . A family is the -slices of , and a map is the -slices of the same ; the correspondence is a bijection because slicing a fixed subset of a product in either coordinate is invertible. The level/coverage transfer is the Key theorem: holds iff for every .

Proposition 2 (the pivot construction has exact coverage). If is a pivot with law (independent of ) and is a Borel set with , then has coverage for every .

Proof. By the definition of , for the fixed under consideration. Hence , the middle equality because the law of under is for every (the pivot property). Coverage is exact and constant, not merely bounded below.

Proposition 3 (UMA from UMP). Suppose for each the test is uniformly most powerful at level for against the relevant alternatives, with acceptance region . Then is uniformly most accurate among sets: for any other set and any in range, .

Proof. Let be any set, with dual acceptance regions , which by Proposition 1 are the regions of level- tests . Fix . By the membership identity,

where is the power of the test of evaluated at the true ; likewise . Because is UMP at level and is some level- test, , so , giving . Uniformity over is the uniformity of the UMP family.

Proposition 4 (Pratt expected-length identity). For an interval-valued set with ,

Proof. For each , . The integrand is nonnegative and jointly measurable, so Tonelli's theorem permits exchanging (integration in against ) with the Lebesgue integral in :

Proposition 5 (shortest interval from a unimodal symmetric pivot). Let be a pivot with a continuous density that is symmetric about and strictly decreasing in -argument, and let the confidence set be with . Among all such the length-minimizing choice in the parameter (when is affine and decreasing in with constant slope) is the symmetric .

Proof. When is affine in with constant scale , the inverted set is of length , so minimizing parameter length is minimizing subject to . Introduce a multiplier: stationarity of gives and , so . Strict unimodality and symmetry of force , and the constraint then determines as the upper quantile. The second-order condition holds because decreasing away from makes the constrained length strictly convex at the symmetric point. For a skewed pivot density no longer forces symmetry, and the equal-tail set () differs from the shortest set.

Connections Master

  • The most-powerful and uniformly-most-powerful test theory 45.02.01 is the source of every optimal interval here: a level- test of inverts to a confidence set, and the most-accurate sets are exactly the inversions of the most-powerful tests, so the Neyman-Pearson and Karlin-Rubin optimality proved upstream becomes minimal false-coverage probability downstream by the membership identity .

  • The Cramér-Rao lower bound 45.01.05 caps interval width: the inverse Fisher information that floors the variance of any unbiased estimator also floors, through the Pratt expected-length identity, the achievable expected length of any asymptotically-efficient interval centered at that estimator, since the Wald interval's squared half-width scales as and exact attainment of both floors is the same exponential-family rigidity.

  • The large-sample Wald, score, and likelihood-ratio intervals 45.04.09 are the asymptotic instances of the test-inversion construction: each inverts the corresponding asymptotic level- test of , the Wald interval from the asymptotically-normal estimator, the score interval from the mean-zero score, and the likelihood-ratio interval from the chi-squared limit of the ratio, and the duality theorem is what licenses reading each test off as an interval.

  • The one-sided composite testing via monotone likelihood ratio 45.02.02 supplies the one-sided uniformly-most-powerful tests whose inversions are the uniformly-most-accurate one-sided confidence bounds; the rejection region with monotone is exactly the structure that makes the inverted set a single bound rather than a disconnected region.

Historical & philosophical context Master

The confidence-set framework is due to Jerzy Neyman, who introduced the confidence interval and the frequentist confidence coefficient in a 1934 discussion and developed the theory fully in his 1937 paper Outline of a theory of statistical estimation based on the classical theory of probability in the Philosophical Transactions of the Royal Society [Neyman 1937] (Phil. Trans. R. Soc. Lond. A 236, 333-380). Neyman framed the interval as a procedure with a guaranteed long-run coverage rate, deliberately separating his frequentist account from R. A. Fisher's contemporaneous fiducial argument, which sought to assign a probability distribution to the parameter directly from the data without a prior. The fiducial program was never made consistent and is now largely abandoned, but the Fisher-Neyman dispute over whether an interval statement is about the parameter or about the procedure shaped the foundations of estimation for decades.

The duality between tests and confidence sets was made explicit in the development running from Neyman and Pearson through Erich Lehmann, whose Testing Statistical Hypotheses organized uniformly most accurate and unbiased confidence sets as the inversions of uniformly most powerful and unbiased tests [Lehmann & Romano Ch. 5]. John W. Pratt's 1961 identity relating the expected length of an interval to the integrated probability of false coverage [Pratt 1961] (J. Amer. Statist. Assoc. 56, 549-567) tied the geometric notion of a short interval to the operational notion of a powerful dual test. The contrast with Bayesian credible sets traces to the Bayesian revival of the mid-twentieth century, where the posterior-probability interval was advanced as the answer to the question Neyman's coverage statement does not address, and the conditions under which the two numerically agree — flat priors in location models, matching priors more generally — remain an active subject in the theory of probability-matching priors.

Bibliography Master

@article{neyman1937outline,
  author    = {Neyman, Jerzy},
  title     = {Outline of a theory of statistical estimation based on the classical theory of probability},
  journal   = {Philosophical Transactions of the Royal Society of London. Series A},
  volume    = {236},
  pages     = {333--380},
  year      = {1937}
}

@book{lehmannromano2005testing,
  author    = {Lehmann, Erich L. and Romano, Joseph P.},
  title     = {Testing Statistical Hypotheses},
  edition   = {3rd},
  series    = {Springer Texts in Statistics},
  publisher = {Springer},
  year      = {2005}
}

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

@article{pratt1961length,
  author    = {Pratt, John W.},
  title     = {Length of confidence intervals},
  journal   = {Journal of the American Statistical Association},
  volume    = {56},
  number    = {295},
  pages     = {549--567},
  year      = {1961}
}

@book{schervish1995theory,
  author    = {Schervish, Mark J.},
  title     = {Theory of Statistics},
  series    = {Springer Series in Statistics},
  publisher = {Springer},
  year      = {1995}
}