Chernoff Bound

Interact

If the form of a distribution is intractable in that it is difficult to find exact probabilities by integration, then good estimates and bounds become important. Bounds on the tails of the distribution of a random variable help us quantify roughly how close to the mean the random variable is likely to be.

We already know two such bounds. Let X be a random variable with expectation μ and SD σ.

Markov’s Bound on the Right Hand Tail

If X is non-negative,

P(Xc)  μc

Chebychev’s Bound on Two Tails

P(|Xμ|c)  σ2c2

Moment generating functions can help us improve upon these bounds in many cases. In what follows, we will assume that the moment generating function of X is finite over the whole real line. If it is finite only over a smaller interval around 0, the calculations of the mgf below should be confined to that interval.

Chernoff Bound on the Right Tail

Observe that if g is an increasing function, then the event Xc is the same as the event g(X)g(c).

For any fixed t>0, the function defined by g(x)=etx is increasing as well as non-negative. So for each t>0,

P(Xc) =P(etXetc) E(etX)etc    (Markov's bound)= MX(t)etc

This is the first step in developing a Chernoff bound on the right hand tail.

For the next step, notice that you can choose t to be any positive number. Some choices of t will give sharper bounds than others. Because these are upper bounds, the sharpest among all of the bounds will correspond to the value of t that minimizes the right hand side. So the Chernoff bound has an optimized form:

P(Xc)  mint>0MX(t)etc

Application to the Normal Distribution

Suppose X has the normal (μ,σ2) distribution and we want to get a sense of how far X can be above the mean. Fix c>0. The exact chance that the value of X is at least c above the mean is

P(Xμc) = 1Φ(c/σ)

because the distribution of Xμ is normal (0,σ2). This exact answer looks neat and tidy, but the standard normal cdf Φ is not easy to work with analytically. Sometimes we can gain more insight from a good bound.

The optimized Chernoff bound is

P(Xμc)  mint>0MXμ(t)etc= mint>0eσ2t2/2etc= mint>0ect+σ2t2/2

The curve below is the graph of exp(ct+σ2t2/2) as a function of t, in the case σ=2 and c=5. The flat line is the exact probability of P(Xμc). The curve is always above the flat line: no matter what t is, the bound is an upper bound. The sharpest of all the upper bounds corresponds to the minimizing value t which is somewhere in the 1.2 to 1.3 range.

png

To find the minimizing value of t analytically, we will use the standard calculus method of minimization. But first we will simplify our calculations by observing that finding the point at which a positive function is minimized is the same as finding the point at which the log of the function is minimized. This is because log is an increasing function.

So the problem reduces to finding the value of t that minimizes the function h(t)=ct+σ2t2/2. By differentiation, the minimizing value of t solves

c = σ2t

and hence

t = cσ2

So the Chernoff bound is

P(Xμc)  ect+σ2t2/2 = ec2/2σ2

Compare this with the bounds we already have. Markov’s bound can’t be applied directly as Xμ can have negative values. Because the distribution of Xμ is symmetric about 0, Chebychev’s bound becomes

P(Xμc)  σ22c2

When c is large, the optimized Chernoff bound is quite a bit sharper than Chebychev’s. In the case σ=2, the graph below shows the exact value of P(Xμc) as a function of c, along with the Chernoff and Chebychev bounds.

png

Chernoff Bound on the Left Tail

By an analogous argument we can derive a Chernoff bound on the left tail of a distribution. For a fixed t>0, the function g(x)=etx is decreasing and non-negative. So for t>0 and any fixed c,

P(Xc) = P(etXetc)  E(etX)etc = MX(t)etc

and therefore

P(Xc)  mint>0MX(t)etc

Sums of Independent Random Variables

The Chernoff bound is often applied to sums of independent random variables. Let X1,X2,,Xn be independent and let Sn=X1+X2++Xn. Fix any number c. For every t>0,

P(Snc)  MSn(t)etc = ni=1MXi(t)etc

This result is useful for finding bounds on binomial tails because the moment generating function of a Bernoulli random variable has a straightforward form. It is also used for bounding tails of sums of independent indicators with possibly different success probabilities. We will leave all this for a subsequent course.