advanced33 min studyLesson 693 in Mathematics

Bernstein, Bennett, and subexponential tails

Use variance and one-sided increment bounds, distinguish exact constants, and specify the local-mgf convention for subexponential variables.

On this page

What you’ll understand

  • State Bennett and Bernstein inequalities with their variance and bound parameters.
  • Explain the quadratic-to-linear tail transition under a specified local-mgf convention.
  • Compare variance-sensitive bounds and derive a valid two-sided inequality.

A range-based bound can waste information when a variable is usually close to its mean and only rarely approaches its maximum. Bennett's inequality retains both the total variance and an upper bound on centered increments. Bernstein replaces Bennett's rate function by a simpler expression with explicit constants. Related subexponential bounds use a quadratic mgf estimate only near zero, so their large-deviation behaviour differs fundamentally from a globally subgaussian bound.

Variance helps only with the correct increment bound

Let independent centered variables X_i have finite variances, satisfy X_i≤b almost surely for one deterministic b>0, and let v be their positive total variance. Bennett bounds the upper tail of their sum by exp(−v h(bt/v)/b²), where h(u)=(1+u)log(1+u)−u. A lower bound on X_i is not needed for this upper-tail theorem. To apply it to the lower tail, bound −X_i from above instead. A two-sided absolute bound follows with a common b only if both directions are controlled.

The proof bounds each centered mgf by variance times the exponential remainder e^(λb)−1−λb. Optimizing over positive λ produces Bennett's h. The elementary lower bound h(u)≥u²/(2(1+u/3)) gives Bernstein's denominator v+bt/3. Thus one third is a specified constant, not an interchangeable convention. The zero-variance case is deterministic and handled separately. If the variables were not centered initially, the almost-sure bound must be checked after centering, rather than applied blindly to their raw maximum.

P ⁣{∑iXi≥t}≤exp⁡ ⁣[−vb2h ⁣(btv)]≤exp⁡ ⁣[−t22(v+bt/3)],h(u)=(1+u)log⁡(1+u)−u\mathbb P\!\left\{\sum_iX_i\ge t\right\}\le\exp\!\left[-\frac v{b^2}h\!\left(\frac{bt}v\right)\right]\le\exp\!\left[-\frac{t^2}{2(v+bt/3)}\right],\quad h(u)=(1+u)\log(1+u)-u
The variables are independent and centered, Xᵢ≤b almost surely, v=ΣE Xᵢ²>0, and t≥0. This is an upper-tail statement.↔ Scroll wide equations horizontally.

Fix a subexponential mgf convention

Here a centered variable is called (ν²,β)-subexponential in the mgf sense if log E exp(λX)≤ν²λ²/2 for every |λ|<1/β, with positive ν and β. This is a local quadratic bound with no singular denominator. Other texts use a subgamma bound or an Orlicz norm and convert parameters with additional constants; they are equivalent classes but not identical numerical definitions. Independence adds ν² parameters for a sum and uses the largest β as a valid common domain scale.

Optimizing within that finite interval gives Gaussian-shaped decay when t≤ν²/β and exponential-shaped decay beyond it. At the boundary, use limits from inside the domain. One convenient two-sided consequence is twice exp of minus one half times the minimum of t²/ν² and t/β. This is a conservative piecewise simplification, not Bennett's exact rate function. The word subexponential in this concentration convention also differs from the heavy-tail theory convention involving convolution tails; it does not describe that class here.

log⁡EeλX≤ν2λ22(∣λ∣<1/β)⟹P{∣X∣≥t}≤2exp⁡ ⁣[−12min⁡(t2ν2,tβ)]\log\mathbb E e^{\lambda X}\le\frac{\nu^2\lambda^2}{2}\quad(|\lambda|<1/\beta)\quad\Longrightarrow\quad\mathbb P\{|X|\ge t\}\le2\exp\!\left[-\frac12\min\left(\frac{t^2}{\nu^2},\frac t\beta\right)\right]
This explicitly defines the centered local-mgf convention used here. Parameters from other subexponential conventions require conversion.↔ Scroll wide equations horizontally.

When the variance-sensitive improvement matters

For a sum of rare Bernoulli successes, centering leaves interval width one, but total variance is np(1−p), much smaller than n/4 when p is small. Bernstein can then improve substantially on Hoeffding without assuming a normal approximation. It still remains an upper bound and can exceed the exact binomial probability by a noticeable factor. Choosing the tight upper centered bound 1−p can improve a directed result further, while b=1 is a simple safe common bound for both centered directions.

The small-deviation regime is controlled mainly by v, while a very large upward deviation feels the bound b. Neither parameter should be estimated from the same data and substituted without a theorem that accounts for its randomness. Empirical Bernstein inequalities exist, but they include additional assumptions and correction terms. Likewise unconditional variance bounds alone do not extend this independent-sum result to an adaptive process; predictable conditional variances are the relevant quantities in a martingale counterpart.

Rare successes improve the exponent

Let K be binomial with n=100 and p=0.05. Bound P(K≥10) using Bennett and Bernstein with b=1, and compare with Hoeffding and the exact tail.

  1. Center each Bernoulli observation: the sum is K−5, with total variance v=100·0.05·0.95=4.75. The centered increments lie between −0.05 and 0.95, so b=1 is valid. The required upper deviation is t=5.

  2. Bennett gives exp(−4.75 h(5/4.75))≈0.133795. Bernstein gives exp(−25/[2(4.75+5/3)])≈0.142551. Bennett is the stronger of these two because Bernstein relaxed its rate function.

  3. Hoeffding uses width total 100 and gives exp(−2·25/100)=exp(−0.5)≈0.606531. Direct summation of the binomial tail gives about 0.0281883, below all three bounds. Variance sensitivity improves the guarantee without making it exact.

The result

Using the safe bound b=1, Bennett≈0.1338 and Bernstein≈0.1426 both improve on Hoeffding≈0.6065; the true tail is about 0.02819.

Explore the model

Interactive lab

Interactive lab · sampling

Small variance behind the same unit width

Pause and Reset, then set p = 0.05 and n = 100. Predict sample-mean variance 0.000475 and add bars at K/n≥0.1 to obtain the worked tail. Reset and set both p = 0.5 and n = 100: the same unit observation range now gives sample-mean variance 0.0025. Compare the variance readouts while retaining the same range-based Hoeffding inputs. The diagram computes exact binomial masses and a separate normal approximation; it does not plot Bennett or Bernstein bounds.

Small variance behind the same unit widthExact binomial masses for the mean of n independent Bernoulli trials, compared with normal mass between adjacent half-count boundaries; normal tails outside the finite support remain outside the bars. Playback highlights outcomes, not random samples. The Berry–Esseen bound applies to standardized CDFs, not individual bars or continuity-corrected errors. At p = 0 or 1, normal comparison and that bound are disabled. Normal-CDF evaluation uses a numerical approximation with absolute error below 7.5×10⁻⁸ (at most 1.5×10⁻⁷ per interval, before rounding). Current readings: Mean sample proportion: 0.5; Sample-mean variance: 0.0125; Uniform CDF error bound: ≤ 0.1062. All parameters are adjustable with the labeled controls below.Exact distribution of the mean of 20 independent Bernoulli variables with success probability 0.5; playback highlights outcome 0.00.250.50.75100.0580.1150.1730.23Sample mean · k / nProbability massNormal bin masses are an approximationk = 0 · exact probability = 9.54e-7
Exact binomial massNormal interval massHighlighted outcome / mean
Mean sample proportion
0.5
Sample-mean variance
0.0125
Uniform CDF error bound
≤ 0.1062
t = 0 s

Exact binomial masses for the mean of n independent Bernoulli trials, compared with normal mass between adjacent half-count boundaries; normal tails outside the finite support remain outside the bars. Playback highlights outcomes, not random samples. The Berry–Esseen bound applies to standardized CDFs, not individual bars or continuity-corrected errors. At p = 0 or 1, normal comparison and that bound are disabled. Normal-CDF evaluation uses a numerical approximation with absolute error below 7.5×10⁻⁸ (at most 1.5×10⁻⁷ per interval, before rounding). 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 forIdentical observation widths can hide very different variances and therefore very different variance-sensitive guarantees.

Check your understanding

Independent centered variables satisfy only Xᵢ≤b. What does the displayed Bennett bound directly control?

Try it yourself

A centered X satisfies the stated (ν²,β) local-mgf condition with ν²=4 and β=1. Give the two-sided simplified bounds at t=1 and t=8, identify their regimes, and cap either answer at one if needed.

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

Give me a hint

Compare t²/4 with t before multiplying the smaller value by minus one half. The transition occurs at t=4.

Compare with the solution

At t=1, the minimum is 1/4, giving 2e−1/8>12e^{-1/8}>1; cap this uninformative answer at one. This is the quadratic regime. At t=8, the minimum of sixteen and eight is eight, giving 2e−4≈0.03663132e^{-4}\approx0.0366313, in the linear regime. The transition t=ν²/β=4 follows from equating the two expressions. These bounds use only the local-mgf convention, not an assertion that X is bounded or has variance exactly four. The parameter ν² is an exponential-moment proxy and need not equal the actual variance.

Take it with you

  • Bennett uses total variance and a directed centered-increment bound; Bernstein retains explicit one-third constants.
  • A local quadratic mgf yields a quadratic-to-linear tail transition.
  • Different subexponential definitions require parameter conversion rather than silent reuse of constants.

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.