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
The question for today
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.
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.
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.
Work it through
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.
The control region gives the interior global estimate . The signal-region excess relative to that value is 2.
The global signal-strength estimate is . Both fitted parameters are positive, so the interior solution lies within the declared domain and reproduces the two observed means.
Under mu = 0, maximizing the combined likelihood gives . It differs from 10 because the signal-region count now also pulls the background fit.
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.
The living notebook
Explore the model
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.
- y at chosen x
- 0.608
- Tangent slope · 2ax
- 1.35
- Signed area · ax³/3
- 0.182
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.
Pause & reason
Check your understanding
Your turn
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 .
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 , giving and b = 6.8. The conditional cost is , 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.
Keep the idea
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.