core26 min studyLesson 242 in Mathematics

Likelihood and sufficient statistics

Use factorization and conditional laws to identify what a statistic preserves about a parameter.

On this page

What you’ll understand

  • Interpret likelihood as a function of a parameter with data fixed.
  • Apply the factorization criterion in a dominated model.
  • Verify sufficiency by a parameter-free conditional distribution.

A binary record can be long while its information about a common success probability is compact. Under independent identical Bernoulli sampling, rearranging the observations changes the time order but not the evidence about that probability. Sufficiency gives a precise version of this claim. It does not say the record is useless for every question; a sequence can reveal an assumption failure such as clustering even when its total is sufficient inside the assumed model. Learn both what a summary preserves and what the model has excluded.

Hold the data fixed

For a model with densities or probability masses fθf_\theta relative to a common dominating measure, the likelihood is L(θ;x)=fθ(x)L(\theta;x)=f_\theta(x) viewed as a function of θ\theta. It need not sum or integrate to one over the parameter space. A density value is not the probability of a continuously distributed observation taking one exact value. Multiplying likelihood by a positive factor depending only on the data preserves all likelihood ratios and maximizers, but multiplying by a parameter-dependent factor generally changes the inference.

If an ordered Bernoulli record has kk successes, its likelihood is pk(1−p)n−kp^k(1-p)^{n-k}. If only the count was recorded, its likelihood also contains the binomial coefficient. This factor counts compatible sequences and is constant as a function of pp. The two experiments consequently yield proportional likelihoods for the common probability. They do not have identical sample spaces, and proportional likelihood alone should not be confused with equality of all possible frequentist procedures under arbitrary sampling plans.

L(p;x)=pT(x)(1−p)n−T(x),T(x)=∑i=1nxiL(p;x)=p^{T(x)}(1-p)^{n-T(x)},\qquad T(x)=\sum_{i=1}^n x_i
Inside the iid Bernoulli model, all parameter dependence in the ordered record enters through the success count.↔ Scroll wide equations horizontally.

Sufficiency as a conditional statement

A statistic TT is sufficient when a version of the conditional distribution of the complete data given TT can be chosen independent of the parameter throughout the model. In a finite Bernoulli sample, conditioning on T=kT=k makes all sequences with that total equally likely. The unknown factors cancel between numerator and denominator. Once the total is known, the remaining ordering uncertainty carries no additional information about the common probability according to this model.

For dominated models, the factorization theorem provides a practical equivalent: the density can be written as a nonnegative parameter-dependent function of T(x)T(x) times a nonnegative data-only function. Supports and zero-density cases must be represented correctly. For a uniform distribution on [0,θ][0,\theta], the sample maximum enters through an indicator constraining θ\theta; discarding that indicator destroys the likelihood. Sufficiency is always relative to the selected family, so adding a time-varying success parameter can make a previously sufficient total insufficient.

fθ(x)=gθ(T(x))h(x),Pp(X=x∣T=k)=(nk)−1(T(x)=k, 0<p<1)f_\theta(x)=g_\theta(T(x))h(x),\qquad P_p(X=x\mid T=k)=\binom nk^{-1}\quad(T(x)=k,\ 0<p<1)
Factorization establishes sufficiency in a dominated family. The displayed conditional probability directly verifies it for interior Bernoulli probabilities; endpoint conditional laws can be extended consistently on impossible totals.↔ Scroll wide equations horizontally.

Compression is not a universal diagnosis

Sufficiency need not mean minimality. The full dataset is trivially sufficient, and appending irrelevant information to a sufficient statistic keeps it sufficient. A minimal sufficient statistic compresses as far as possible in the sense that it is a measurable function of every sufficient statistic, with the customary almost-sure qualifications. In finite positive-density models, comparing which sample likelihood ratios are constant in the parameter often identifies the minimal partition. For Bernoulli data with interior probability, two records lie in the same partition exactly when they have the same total.

Keep raw observations when checking assumptions. A record of ten zeros followed by ten ones and an alternating record have equal likelihoods for a single common Bernoulli probability, yet suggest different temporal mechanisms. Sufficiency says their orders do not help estimate that one parameter if independence and identical distribution hold. It does not certify those assumptions or erase the value of provenance. The practical workflow is to preserve the original data for diagnostics while using sufficient summaries for calculations justified by the declared model.

A conditional distribution without p

Four iid Bernoulli observations have total two. Compute the conditional probability of the ordered record (1,0,1,0), and explain why it does not depend on the success probability.

  1. For 0<p<10<p<1, the record probability is p2(1−p)2p^2(1-p)^2. There are (42)=6\binom42=6 distinct records containing exactly two successes.

  2. The total probability of two successes is 6p2(1−p)26p^2(1-p)^2. Divide the selected record probability by this total to obtain its conditional probability.

  3. All parameter factors cancel, leaving 1/6. The same calculation applies to each of the other five records, and their conditional probabilities sum to one.

The result

The conditional probability is 1/6 for every interior p. At p = 0 or 1, total two has probability zero, so this conditional event cannot be observed under that law.

Explore the model

Interactive lab

Interactive lab · probability

Separate a count law from a sequence law

Reset and set both p=0.5,n=4p=0.5,n=4. Predict the count-two bar: 6/16=0.3756/16=0.375. Reset, set both p=0.25,n=4p=0.25,n=4, and predict 6(0.25)2(0.75)2=0.21093756(0.25)^2(0.75)^2=0.2109375. The bars change with pp, while the conditional probability of each compatible ordered record remains 1/61/6. This instrument displays counts, not sequence order.

Separate a count law from a sequence lawThese are exact binomial probabilities for independent trials with one fixed success probability. Playback tours the distribution; it does not resample outcomes. Current readings: Mean · np: 5; Variance · np(1−p): 2.5; Total probability: 1. All parameters are adjustable with the labeled controls below.03581000.070.150.220.3Successes kP(X = k)k = 0 · P = 9.77e-4
Exact probabilityHighlighted outcomeMean np
Mean · np
5
Variance · np(1−p)
2.5
Total probability
1
t = 0 s

These are exact binomial probabilities for independent trials with one fixed success probability. Playback tours the distribution; it does not resample outcomes. 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 forA sufficient statistic has a parameter-dependent distribution even when the data conditional on it do not.

Check your understanding

Which factor may be removed from a likelihood without changing its parameter maximizers or likelihood ratios?

Try it yourself

For two iid observations from the uniform distribution on [0,θ][0,\theta], θ>0\theta>0, write the joint density and exhibit a sufficient statistic. Explain why the sample sum is not enough by comparing records (1,3)(1,3) and (2,2)(2,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

Include the support indicator before collecting the parameter-dependent terms. Which record permits a parameter value of 2.5?

Compare with the solution

The joint density is θ−21{0≤x1,x2≤θ}\theta^{-2}\mathbf1\{0\le x_1,x_2\le\theta\}. It factors through the maximum, with the nonnegativity restriction in a data-only factor. Thus max⁡(X1,X2)\max(X_1,X_2) is sufficient. Both records sum to four, but at θ=2.5\theta=2.5 the first has zero density and the second positive density, whereas for θ>3\theta>3 their densities agree. Their likelihood ratio therefore depends on the parameter, demonstrating that the sum loses relevant information.

Take it with you

  • Likelihood fixes the observed data and varies the candidate parameter.
  • Sufficiency means a parameter-free conditional law for the remaining data.
  • A sufficient summary for estimation does not replace raw data for model checks.

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.