Tutorials › Statistical Theory I - Probability and Distribution Theory › Bennett’s Inequality and Its Bernstein Comparison

Concentration and Exponential Bounds · Tutorial 927 of 994

Bennett's Inequality and Its Bernstein Comparison

Learn how a variance-sensitive exponential bound follows from the moment-generating function, and why Bennett’s deviation function gives a sharper exponent than Bernstein’s.

Advanced 10 min read

What You'll Learn

  • State the assumptions for Bennett’s inequality for independent centered summands
  • Derive the exponential-moment bound that leads to Bennett’s deviation function
  • Optimize the Chernoff bound and interpret the parameter in the deviation function
  • Prove how Bennett’s bound implies Bernstein’s inequality
  • Compare the two bounds in explicit binomial and heterogeneous examples

Why a Variance-Sensitive Bound Has More Than One Form

For a sum of independent centered variables, an upper bound on each summand and the sum’s variance together control its upper tail. Bernstein’s inequality expresses this control through a simple rational function of the deviation. Bennett’s inequality uses a different deviation function: it retains more detail from the exponential-moment calculation and gives a tail bound at least as strong as Bernstein’s under the same assumptions.

The comparison is useful because the two expressions describe different levels of approximation. At moderate deviations their exponents are close, while Bennett’s exponent grows faster for large values of the scaled deviation. The starting point is an upper bound on each centered summand, not necessarily a bound on its absolute value.

Definition: Let \(X_1,\ldots,X_n\) be independent random variables with \(E[X_i]=0\) and finite variances \(\sigma_i^2=\operatorname{Var}(X_i)\). Suppose \(X_i\le b\) almost surely for every \(i\), where \(b>0\). Write \(S=\sum_{i=1}^n X_i\) and \(v=\sum_{i=1}^n\sigma_i^2\). Bennett’s deviation function is \[ h(u)=(1+u)\log(1+u)-u,\qquad u\ge 0. \]

The function satisfies \(h(0)=0\) and \(h'(u)=\log(1+u)>0\) for \(u>0\), so it is increasing away from zero. Its local behavior is quadratic: \(h(u)=u^2/2+o(u^2)\) as \(u\) tends to zero. For large \(u\), it grows on the order of \(u\log u\). These features help explain how the resulting tail bound responds to different deviation sizes.

Deriving Bennett’s Bound

The key step is to control the moment-generating function of a single centered summand. For \(\lambda\ge 0\) and real \(x\), Taylor’s formula with integral remainder gives \[ e^{\lambda x}-1-\lambda x =\lambda^2x^2\int_0^1(1-s)e^{s\lambda x}\,ds. \] For fixed \(\lambda\ge0\), the integral on the right is nondecreasing in \(x\), since its integrand is nondecreasing in \(x\) for every \(s\in[0,1]\). Thus, whenever \(x\le b\), \[ e^{\lambda x}\le 1+\lambda x+ \frac{x^2}{b^2}\bigl(e^{\lambda b}-1-\lambda b\bigr). \] This includes negative \(x\); no lower bound on \(X_i\) is needed.

Theorem (Bennett’s inequality): Under the conditions in the definition, if \(v>0\), then for every \(t>0\), \[ P(S\ge t)\le \exp\left\{-\frac{v}{b^2}h\left(\frac{bt}{v}\right)\right\}. \] If \(v=0\), then \(S=0\) almost surely and \(P(S\ge t)=0\) for \(t>0\).

To prove the result, take expectations in the pointwise inequality with \(x=X_i\). Since \(E[X_i]=0\) and \(E[X_i^2]=\sigma_i^2\), \[ E[e^{\lambda X_i}] \le 1+\frac{\sigma_i^2}{b^2}\bigl(e^{\lambda b}-1-\lambda b\bigr) \le \exp\left\{\frac{\sigma_i^2}{b^2}\bigl(e^{\lambda b}-1-\lambda b\bigr)\right\}. \] The last step uses \(1+y\le e^y\) for \(y\ge0\). Independence gives \[ E[e^{\lambda S}] =\prod_{i=1}^n E[e^{\lambda X_i}] \le \exp\left\{\frac{v}{b^2}\bigl(e^{\lambda b}-1-\lambda b\bigr)\right\}. \] For \(\lambda>0\), Markov’s inequality applied to \(e^{\lambda S}\) now yields \[ P(S\ge t)\le \exp\left\{-\lambda t+\frac{v}{b^2} \bigl(e^{\lambda b}-1-\lambda b\bigr)\right\}. \] The exponent is minimized over \(\lambda\ge0\) at the solution of \[ -t+\frac{v}{b}(e^{\lambda b}-1)=0, \qquad \lambda_*=\frac{1}{b}\log\left(1+\frac{bt}{v}\right). \] Set \(u=bt/v\). Substituting \(\lambda_*\) gives \[ -\lambda_*t+\frac{v}{b^2} \bigl(e^{\lambda_*b}-1-\lambda_*b\bigr) =-\frac{v}{b^2}\bigl((1+u)\log(1+u)-u\bigr). \] This is the stated Bennett bound. If \(v=0\), every \(X_i\) has zero variance and mean zero, so each equals zero almost surely; the asserted tail probability follows. \(\square\)

The optimized parameter \(u=bt/v\) is dimensionless. It compares the requested deviation \(t\), scaled by the per-summand upper bound \(b\), with the total variance \(v\). Small \(u\) corresponds to the range where a quadratic approximation is accurate; larger \(u\) exposes the stronger growth of \(h(u)\).

Worked Example: A Binomial Tail and the Correct Numerical Comparison

Let \(K\sim\operatorname{Binomial}(8,0.25)\). Write \(K-2=\sum_{i=1}^8 X_i\), where the \(X_i\) are independent Bernoulli\((0.25)\) variables centered by subtracting \(0.25\). Each \(X_i\le 0.75\), so take \(b=0.75\). Each variance is \(0.25(0.75)=0.1875\), giving \(v=8(0.1875)=1.5\). The event \(K\ge4\) is \(K-2\ge2\), so \(t=2\) and \(bt/v=0.75(2)/1.5=1\).

Since \(h(1)=2\log 2-1\approx0.3862943611\) and \(v/b^2=1.5/(0.75)^2=8/3\), Bennett’s bound is \[ P(K\ge4)\le \exp\left\{-\frac{8}{3}(0.3862943611)\right\} =\exp(-1.0301182963)\approx 0.357. \]

Bernstein’s bound, derived below, gives exponent \(t^2/[2(v+bt/3)]\). Here that exponent is \(4/[2(1.5+0.75(2)/3)]=4/4=1\), so the Bernstein bound is \(e^{-1}\approx0.368\). For comparison, the exact probability is \[ P(K\ge4)=1-\sum_{k=0}^3 {8\choose k}(0.25)^k(0.75)^{8-k} =1-(0.1001129150+0.2669677734+0.3114624023+0.2076416016) \approx0.113815. \] Both bounds are valid upper bounds; Bennett is slightly smaller in this case.

Why Bennett Implies Bernstein

Bernstein’s inequality has a more compact denominator. The comparison follows from a lower bound on \(h\), which we prove rather than assume.

Theorem (Bernstein comparison): For \(u\ge0\), \[ h(u)\ge \frac{u^2}{2(1+u/3)}. \] Consequently, under the assumptions of Bennett’s inequality, for \(t>0\), \[ P(S\ge t)\le \exp\left\{-\frac{t^2}{2(v+bt/3)}\right\}. \]

For the function inequality, define \[ g(u)=h(u)-\frac{3u^2}{2(u+3)}. \] We have \(g(0)=0\). Differentiating gives \[ g'(u)=\log(1+u)-\frac{3u(u+6)}{2(u+3)^2}. \] Call the right-hand side \(q(u)\). Then \(q(0)=0\), and direct differentiation gives \[ q'(u)=\frac{1}{1+u}-\frac{27}{(u+3)^3} =\frac{u^2(u+9)}{(1+u)(u+3)^3}\ge0. \] Therefore \(q(u)\ge0\) and \(g'(u)\ge0\) for \(u\ge0\). Since \(g(0)=0\), it follows that \(g(u)\ge0\), proving the claimed inequality. Substituting \(u=bt/v\) into it gives \[ \frac{v}{b^2}h\left(\frac{bt}{v}\right) \ge \frac{v}{b^2}\frac{(bt/v)^2}{2(1+bt/(3v))} =\frac{t^2}{2(v+bt/3)}. \] Applying Bennett’s inequality proves the Bernstein comparison. \(\square\)

For \(u>0\), the comparison is strict, so Bennett’s exponent is strictly larger and its upper bound strictly smaller. Near zero, \(h(u)\) is approximately \(u^2/2\). For large \(u\), Bennett’s deviation function grows roughly as \(u\log u\), whereas the Bernstein exponent’s corresponding function \(u^2/[2(1+u/3)]\) grows roughly as \(3u/2\). Thus the two bounds can differ more substantially at large scaled deviations.

Worked Example: Unequal Bernoulli Summands

Let \(X_1=B_1-0.2\), \(X_2=B_2-0.6\), and \(X_3=B_3-0.4\), where the \(B_i\) are independent Bernoulli variables with success probabilities \(0.2,0.6,0.4\), respectively. The upper bounds on the centered variables are \(0.8,0.4,0.6\), so a common bound is \(b=0.8\). Their variances sum to \[ v=0.2(0.8)+0.6(0.4)+0.4(0.6)=0.16+0.24+0.24=0.64. \] For \(t=1\), the scaled deviation is \(u=bt/v=0.8/0.64=1.25\), and \[ h(1.25)=2.25\log(2.25)-1.25\approx0.574592986. \] Since \(v/b^2=0.64/0.64=1\), Bennett gives \(P(S\ge1)\le e^{-0.574592986}\approx0.563\). Bernstein’s exponent is \[ \frac{1^2}{2(0.64+0.8/3)}=\frac{1}{1.813333\ldots}\approx0.551471, \] giving \(P(S\ge1)\le e^{-0.551471}\approx0.576\). The exact event requires all three Bernoulli variables to equal one: \(S=B_1+B_2+B_3-1.2\ge1\) means \(B_1+B_2+B_3\ge2.2\), hence the integer-valued sum must be at least \(3\). Its probability is \(0.2(0.6)(0.4)=0.048\).

Worked Example: A Deviation for a Sample Proportion

Let \(K\sim\operatorname{Binomial}(100,0.5)\). Its centered sum \(S=K-50\) has \(v=100(0.5)(0.5)=25\), and each centered Bernoulli summand is at most \(b=0.5\). Consider the upper deviation \(K\ge60\), or \(S\ge10\). Then \(u=bt/v=0.5(10)/25=0.2\).

Here \(h(0.2)=1.2\log(1.2)-0.2\approx0.018785868\), so the Bennett exponent is \(25/(0.5)^2\) times this value, or \(100(0.018785868)\approx1.878587\). Thus Bennett gives \[ P(K\ge60)\le e^{-1.878587}\approx0.1528. \] For Bernstein, the exponent is \[ \frac{10^2}{2(25+0.5(10)/3)} =\frac{100}{53.333333\ldots}=1.875, \] and the bound is \(e^{-1.875}\approx0.1534\). The bounds are close because the scaled deviation \(u=0.2\) is small, but Bennett remains the sharper one.

Using the Bound Carefully

Bennett’s inequality is a one-sided statement: its assumption is \(X_i\le b\), and its conclusion bounds \(P(S\ge t)\). If an absolute two-sided bound \(|X_i|\le b\) is available, apply the same result to \(X_i\) and to \(-X_i\), then use the union bound: \[ P(|S|\ge t)\le 2\exp\left\{-\frac{v}{b^2}h\left(\frac{bt}{v}\right)\right\}. \] The factor of two comes from combining the upper and lower tails; it is not part of the one-sided Bennett bound.

A common mistake is to replace Bennett’s deviation function by its quadratic approximation \(h(u)\approx u^2/2\) for every \(u\). That approximation describes the behavior near zero, not the full function. For larger \(u\), using it as an equality can substantially misstate the exponent. Another error is to use \(b\) as though it were an upper bound on \(|X_i|\): for the upper-tail theorem, only \(X_i\le b\) is required, but the same \(b\) does not automatically bound the lower tail. Check the direction and size of the bound before applying the result.

Key takeaway: Bennett’s inequality follows by bounding each centered moment-generating function using the summand’s upper bound and variance, then optimizing the Chernoff parameter. Its deviation function \(h(u)\) yields a sharper upper-tail exponent than Bernstein’s rational approximation under the same assumptions.

Check Your Understanding

Use the assumptions and results above to answer the following questions.

  1. Why does the Bennett derivation require an upper bound on each centered summand for the upper tail, but not a lower bound?
  2. For independent centered summands with total variance \(v\), upper bound \(b\), and threshold \(t\), write the scaled deviation \(u\) and the Bennett exponent.
  3. What inequality relating \(h(u)\) to a rational function of \(u\) leads from Bennett’s inequality to Bernstein’s inequality?
  4. If \(|X_i|\le b\) for every summand, how does the one-sided Bennett result yield a two-sided bound?
  5. Why is replacing \(h(u)\) by \(u^2/2\) unreliable when \(u\) is not small?