research29 min studyLesson 649 in Physics

Profile likelihoods and nuisance parameters

Construct a joint likelihood with auxiliary data, profile shared nuisance parameters, and distinguish asymptotic uncertainty from guaranteed coverage.

On this page

What you’ll understand

  • Construct a Poisson likelihood including auxiliary measurements.
  • Distinguish unconditional and conditional maximum-likelihood estimates.
  • State why nuisance modeling, boundaries, and calibration affect interval coverage.

Particle measurements depend on calibration, efficiency, backgrounds, and theory assumptions as well as the parameter of interest. A likelihood represents those dependencies in one joint probability model. Profiling asks how well each proposed signal strength can fit the data after the nuisance parameters are adjusted within their auxiliary constraints. It does not average them away or guarantee that every missing uncertainty has been captured. We will solve a small count model, distinguish profiling from marginalization, and connect likelihood intervals to the repeated-experiment procedure needed for a coverage claim.

Auxiliary measurements are data too

For binned independent Poisson observations, each mean depends on the signal strength mu and nuisance vector theta. Auxiliary measurements such as control counts or calibration estimates supply additional likelihood factors. A Gaussian constraint can represent a measured calibration with an approximately Gaussian sampling distribution, but it should not be confused automatically with a Bayesian prior. The full model must state what is observed, what is fixed, and which variables would fluctuate in repeated experiments.

A nuisance parameter shared between channels creates a physical correlation. Giving every histogram bin an independent shift can remove that relationship, while forcing unrelated uncertainties to share a parameter can overconstrain them. Positive yields and efficiencies may require suitable parameterizations or explicit boundaries. A flexible likelihood can only quantify effects represented by its parameterized model; an omitted process or unsupported functional form is not repaired by the profiling operation itself.

L(μ,θ)=∏iνi(μ,θ)nie−νi(μ,θ)ni! p(a∣θ),λ(μ)=L(μ,θ^^(μ))L(μ^,θ^)L(\mu,\theta)=\prod_i\frac{\nu_i(\mu,\theta)^{n_i}e^{-\nu_i(\mu,\theta)}}{n_i!}\,p(a\mid\theta),\qquad \lambda(\mu)=\frac{L(\mu,\widehat{\widehat\theta}(\mu))}{L(\widehat\mu,\widehat\theta)}
A denotes auxiliary observations. Double-hat theta maximizes the likelihood at fixed mu, while the single-hat pair maximizes it over the declared full parameter space. The profile likelihood ratio is dimensionless and at most one when both optimizations use compatible domains.↔ Scroll wide equations horizontally.

Profiling changes the best nuisance fit at each hypothesis

Consider a signal-region count with mean mu s plus b and an independent control count with mean tau b. Here s and tau are known positive constants, while b is a background nuisance parameter. If the unconstrained fitted signal is positive, the global fit can reproduce both observed counts by setting b to the control count divided by tau and assigning the remaining signal-region excess to mu s. Under the no-signal hypothesis, the two regions instead jointly determine b.

This shift is precisely why a test should not generally fix b to its global best value while varying mu. The conditional maximum can use auxiliary information differently at each proposed hypothesis. Profiling maximizes rather than integrates; Bayesian marginalization instead integrates with a prior and answers a different inferential question. Neither procedure makes parameter degeneracies disappear. Weak auxiliary information can leave a broad profile even when the signal-region count is large.

L(μ,b)=Pois⁡(n∣μs+b)Pois⁡(m∣τb),b^=mτ,μ^=n−m/τs,b^^(0)=n+m1+τL(\mu,b)=\operatorname{Pois}(n\mid\mu s+b)\operatorname{Pois}(m\mid\tau b),\qquad \widehat b=\frac m\tau,\qquad \widehat\mu=\frac{n-m/\tau}{s},\qquad \widehat{\widehat b}(0)=\frac{n+m}{1+\tau}
The global formulas apply to an interior solution with positive fitted mu and b under the physical domain. If the signal-strength maximum is at a boundary, those formulas must be replaced by the constrained fit. The null conditional background uses both counts.↔ Scroll wide equations horizontally.

Intervals need a calibrated sampling distribution

A likelihood-ratio curve orders parameter values by fit quality, but converting it into a confidence interval requires a threshold with a justified sampling distribution. Wilks-type chi-square approximations need regularity, identifiability, sufficiently informative data, and appropriate treatment of boundaries. Small counts, nearly degenerate nuisance effects, or a parameter appearing only under an alternative can violate familiar asymptotic assumptions. Simulated repeated experiments can test or calibrate the chosen construction.

Frequentist coverage means that a specified procedure includes the fixed true parameter in the stated fraction of repeated datasets. It is not a posterior probability that the fixed parameter lies inside this already observed interval. Auxiliary fluctuations and the nuisance treatment belong to that procedure. Conditional or unconditional simulation choices can produce different questions, so a coverage study should state which parameters and observations vary. A successful fit and small optimizer error do not by themselves establish coverage or physical model adequacy.

A background fit changes under the null

In a synthetic two-region model observe n = 12 and m = 20, with s = 4 and tau = 2. All means are nonnegative and mu is constrained nonnegative. Find the global fitted b and mu, then the conditional fitted background under mu = 0.

  1. The control region gives the interior global estimate b^=20/2=10\widehat b=20/2=10. The signal-region excess relative to that value is 2.

  2. The global signal-strength estimate is μ^=(12−10)/4=0.5\widehat\mu=(12-10)/4=0.5. Both fitted parameters are positive, so the interior solution lies within the declared domain and reproduces the two observed means.

  3. Under mu = 0, maximizing the combined likelihood gives b^^(0)=(12+20)/(1+2)=32/3≃10.667\widehat{\widehat b}(0)=(12+20)/(1+2)=32/3\simeq10.667. It differs from 10 because the signal-region count now also pulls the background fit.

The result

The global fit is (mu, b) = (0.5, 10); the no-signal conditional fit uses b = 32/3. A profile test must allow that nuisance adjustment.

Explore the model

Interactive lab

Interactive lab · function

A local quadratic likelihood analogy

Interpret the parabola only as a possible local quadratic cost. Reset, set Coefficient a = 0.5 and Position x = 1, and predict y = 0.5. For a fresh run Reset, set a = 2 and x = 1, and predict y = 2 with a steeper tangent.

A local quadratic likelihood analogyThis lab uses only y = ax². Playback traces the curve within the fixed integration interval from 0 to x; it does not estimate numerical error or change the chosen point. Current readings: y at chosen x: 0.608; Tangent slope · 2ax: 1.35; Signed area · ax³/3: 0.182. All parameters are adjustable with the labeled controls below.−2.4−1.201.22.4−9−4.504.59xy(x, ax²)
y = ax²TangentSigned integral
y at chosen x
0.608
Tangent slope · 2ax
1.35
Signed area · ax³/3
0.182
t = 0 s

This lab uses only y = ax². Playback traces the curve within the fixed integration interval from 0 to x; it does not estimate numerical error or change the chosen point. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forThis is one chosen parabola, not a profile optimizer or confidence-interval calculator. It has no Poisson data, auxiliary measurement, nuisance parameter, or boundary; a real likelihood need not be locally quadratic over the relevant range.

Check your understanding

What defines profiling over nuisance parameters at a fixed signal strength?

Try it yourself

Use a separate synthetic Gaussian model: independent measurements y = 10 with standard deviation 2 and a = 6 with standard deviation 1 have means mu + b and b. Mu and b are unrestricted real parameters. Find the global fit and the profiled b at mu = 0 by minimizing (10−μ−b)2/4+(6−b)2(10-\mu-b)^2/4+(6-b)^2.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

The global fit can make both residuals zero. At fixed mu = 0, differentiate the weighted sum with respect to b.

Compare with the solution

The global fit is b = 6 and mu = 4. At mu = 0 the derivative condition is (b−10)/2+2(b−6)=0(b-10)/2+2(b-6)=0, giving 2.5b=172.5b=17 and b = 6.8. The conditional cost is 3.22/4+0.82=3.23.2^2/4+0.8^2=3.2, compared with zero globally. This continuous Gaussian exercise illustrates profiling; it is not a substitute for a low-count Poisson likelihood or a generic exact confidence threshold.

Take it with you

  • Likelihoods combine primary observations with explicitly modeled auxiliary data.
  • Profiling maximizes nuisance parameters conditionally and differs from marginalization.
  • Model completeness, boundaries, and sampling calibration determine what uncertainty and coverage claims are justified.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.