Tutorials › Statistical Theory I - Probability and Distribution Theory › Poisson Tail Bounds from Its MGF

Probability and Tail Inequalities · Tutorial 908 of 994

Poisson Tail Bounds from Its MGF

Use the Poisson moment-generating function to obtain optimized upper- and lower-tail bounds, and assess their conservatism against exact probabilities.

Intermediate 9 min read

What You'll Learn

  • Derive the moment-generating function of a Poisson random variable.
  • Select the correct exponential Markov inequality for each tail.
  • Optimize the Chernoff parameter and express the result using a rate function.
  • Handle zero, mean, and noninteger thresholds correctly.
  • Compare explicit Poisson Chernoff bounds with exact tail probabilities.

Why the Poisson MGF gives useful tail bounds

A Poisson model is a natural starting point for counts: \(X\sim\operatorname{Pois}(\lambda)\), with \(\lambda>0\), takes values \(0,1,2,\ldots\) and has mean \(\lambda\). Its exact tail probabilities can be computed by summing its probability mass function. For large means or many different thresholds, however, a compact bound can be easier to evaluate and interpret. The moment-generating function (MGF) gives such bounds through the exponential Markov inequality.

The central idea is to replace a count-tail event by an event involving \(e^{tX}\). Choosing \(t>0\) emphasizes large values and bounds an upper tail; choosing \(t<0\) emphasizes small values and bounds a lower tail. Optimizing \(t\) produces an explicit exponential bound. The resulting expression is useful, but it is generally not the exact tail probability.

The Poisson MGF and its exponential form

Definition: If \(X\sim\operatorname{Pois}(\lambda)\), then \(P(X=k)=e^{-\lambda}\lambda^k/k!\) for \(k=0,1,2,\ldots\), where \(\lambda>0\). Its moment-generating function is \(M_X(t)=E[e^{tX}]\), when this expectation is finite.

For a Poisson variable, the MGF is finite for every real \(t\). Directly summing over the probability mass function gives

$$ \begin{aligned} M_X(t) &=\sum_{k=0}^{\infty}e^{tk}e^{-\lambda}\frac{\lambda^k}{k!}\\ &=e^{-\lambda}\sum_{k=0}^{\infty}\frac{(\lambda e^t)^k}{k!}\\ &=e^{-\lambda}e^{\lambda e^t} =\exp\{\lambda(e^t-1)\}. \end{aligned} $$

The series converges for every real \(t\), because it is the exponential series evaluated at the finite number \(\lambda e^t\). Thus any real Chernoff parameter used below is in the MGF domain. The exponential Markov inequality says that if \(Y\geq0\) and \(c>0\), then \(P(Y\geq c)\leq E[Y]/c\). Applying it to \(Y=e^{tX}\) gives the upper-tail bound for \(t>0\). For \(t<0\), the function \(x\mapsto e^{tx}\) decreases, so a lower-tail event instead implies a large value of \(e^{tX}\).

Optimizing the upper and lower tails

For \(u>0\), define the Poisson rate function

$$ h(u)=u\log u-u+1, $$

and set \(h(0)=1\), its continuous limiting value. For \(u\geq0\), \(h(u)\geq0\): for \(u>0\), \(h'(u)=\log u\), so \(h\) decreases up to \(u=1\), increases after \(u=1\), and attains its minimum \(h(1)=0\).

Theorem: Let \(X\sim\operatorname{Pois}(\lambda)\), with \(\lambda>0\). For \(a>\lambda\), \(P(X\geq a)\leq\exp\{-\lambda h(a/\lambda)\}\). For \(0<a<\lambda\), \(P(X\leq a)\leq\exp\{-\lambda h(a/\lambda)\}\). At \(a=\lambda\), the same expression equals \(1\), giving only the trivial bound. At \(a=0\), \(P(X\leq0)=e^{-\lambda}\) exactly.

Thresholds \(a\) need not be integers. The event notation has its usual meaning for the integer-valued variable \(X\); for instance, \(X\geq a\) is \(X\geq\lceil a\rceil\). The theorem remains valid for real \(a\), because the exponential Markov argument uses the threshold \(a\) itself.

Proof. First suppose \(a>\lambda\). For every \(t>0\), if \(X\geq a\), then \(e^{tX}\geq e^{ta}\). Exponential Markov therefore gives

$$ P(X\geq a)\leq e^{-ta}E[e^{tX}] =\exp\{\lambda(e^t-1)-ta\}. $$

The exponent \(g(t)=\lambda(e^t-1)-ta\) has derivative \(g'(t)=\lambda e^t-a\) and second derivative \(g''(t)=\lambda e^t>0\). Its unique minimum over \(t>0\) is at \(t^\ast=\log(a/\lambda)>0\). Substitution yields

$$ g(t^\ast)=a-\lambda-a\log(a/\lambda) =-\lambda h(a/\lambda), $$

which proves the upper-tail claim.

Now let \(0<a<\lambda\). For \(t<0\), the event \(X\leq a\) implies \(e^{tX}\geq e^{ta}\). Thus

$$ P(X\leq a)\leq e^{-ta}E[e^{tX}] =\exp\{\lambda(e^t-1)-ta\}. $$

The same strictly convex function \(g(t)\) has its unique minimum at \(t^\ast=\log(a/\lambda)<0\), which is now an admissible lower-tail parameter. Substituting this value again gives \(g(t^\ast)=-\lambda h(a/\lambda)\). Finally, because \(X\) is nonnegative and integer-valued, \(X\leq0\) means \(X=0\), whose probability is \(e^{-\lambda}\). This establishes the endpoint directly. \(\square\)

For a lower threshold \(a=0\), one can also see the limiting expression in the rate function: \(\exp\{-\lambda h(0)\}=e^{-\lambda}\), which matches the exact probability. At \(a=\lambda\), the optimizing parameter is \(t=0\), and the exponential bound is \(1\); there is no nontrivial tail bound at the mean from this optimization.

Worked comparisons with exact probabilities

Worked Example: An Upper Tail with Mean 4

Let \(X\sim\operatorname{Pois}(4)\). Bound \(P(X\geq8)\), then compare it with the exact tail.

Since \(8>4\), the optimized upper-tail bound applies:

$$ \begin{aligned} P(X\geq8) &\leq \exp\{-4h(2)\}\\ &=\exp\{-4(2\log2-2+1)\}\\ &=\exp(4-8\log2) \approx 0.2131. \end{aligned} $$

For an exact comparison, use \(p_k=P(X=k)=e^{-4}4^k/k!\). The first omitted mass in the tail is \(p_8\), and later masses can be obtained from \(p_{k+1}=p_k\,4/(k+1)\). In particular, \(p_8\approx0.02977\), \(p_9\approx0.01323\), \(p_{10}\approx0.00529\), and subsequent terms decrease by the same recurrence. Summing \(p_8+p_9+\cdots\) gives \(P(X\geq8)\approx0.0511\). The bound is above the exact probability, as required, and is about 4.2 times as large.

Worked Example: A Lower Tail with Mean 5

Let \(X\sim\operatorname{Pois}(5)\). Bound \(P(X\leq2)\) and compare the bound with the exact sum.

Here \(0<2<5\), so the lower-tail version gives

$$ \begin{aligned} P(X\leq2) &\leq \exp\{-5h(2/5)\}\\ &=\exp\{-5[(2/5)\log(2/5)-2/5+1]\}\\ &=\exp\{-5(0.233483\ldots)\} \approx0.3112. \end{aligned} $$

The exact probability includes the three masses at 0, 1, and 2:

$$ \begin{aligned} P(X\leq2) &=e^{-5}\left(1+5+\frac{5^2}{2!}\right)\\ &=18.5e^{-5} \approx0.12465. \end{aligned} $$

Thus the optimized exponential bound is valid but not exact: in this case it is about 2.5 times the exact lower-tail probability.

Worked Example: A Noninteger Upper Threshold

Let \(X\sim\operatorname{Pois}(10)\) and consider \(P(X\geq15.5)\). Since \(X\) is integer-valued, this event is exactly \(X\geq16\), but the theorem can be applied directly with \(a=15.5\).

The bound using the stated threshold is

$$ P(X\geq15.5) \leq \exp\{-10h(1.55)\} =\exp\{-10[1.55\log(1.55)-1.55+1]\} \approx0.2745. $$

For an exact comparison, sum from \(k=16\) upward using \(p_{10}=e^{-10}10^{10}/10!\approx0.12511\) and the recurrence \(p_{k+1}=10p_k/(k+1)\). It gives \(p_{16}\approx0.02170\), \(p_{17}\approx0.01276\), \(p_{18}\approx0.00709\), and the remaining terms decrease successively; their full sum from 16 onward is approximately \(0.0487\). The bound is again larger than the exact probability.

One may instead use the integer threshold \(16\) in the bound, which yields \(\exp\{-10h(1.6)\}\approx0.2187\). This is also valid and is tighter here. Applying the theorem with the real threshold \(15.5\) is valid, but replacing it by the least integer threshold can improve the bound for an integer-valued variable.

How to read the bound and avoid common errors

The same rate-function expression describes both tails, but the parameter constraints differ: an upper-tail deviation uses \(t>0\), while a lower-tail deviation uses \(t<0\). The optimized exponent is nonpositive because \(h(u)\geq0\), and it is strictly negative whenever \(u\ne1\). Consequently, the bound decreases exponentially as the rate-function cost of the threshold grows.

A common mistake is to use \(t>0\) for a lower-tail event and claim that \(X\leq a\) implies \(e^{tX}\geq e^{ta}\). For positive \(t\), the exponential is increasing, so the implication goes in the opposite direction: \(X\leq a\) gives \(e^{tX}\leq e^{ta}\). The correction is to choose \(t<0\), making the exponential decreasing; then \(X\leq a\) does imply \(e^{tX}\geq e^{ta}\), and Markov's inequality applies.

Another important distinction is between a bound and an exact probability. Optimizing the exponential Markov expression finds the best bound within this one-parameter family of exponential Markov bounds. It does not make the resulting number equal to the Poisson tail in general. Exact probabilities still come from summing the mass function, and the examples show that the bound can be conservative by different amounts in different settings.

Key takeaway: For a Poisson count, the MGF \(M_X(t)=\exp\{\lambda(e^t-1)\}\) gives the optimized tail bound \(\exp\{-\lambda h(a/\lambda)\}\). Use it above the mean for an upper tail and below the mean for a lower tail, with the sign of \(t\) chosen accordingly.

Check Your Understanding

Work through these questions using the Poisson MGF and the rate function.

  1. For \(X\sim\operatorname{Pois}(6)\), which sign of \(t\) should be used to bound \(P(X\geq9)\), and what is the optimizing value of \(t\)?
  2. Write the optimized bound for \(P(X\leq3)\) when \(X\sim\operatorname{Pois}(6)\), leaving the answer in terms of \(h\).
  3. Why does the optimized expression equal 1 when the threshold is \(a=\lambda\)?
  4. For integer-valued \(X\), express the event \(X\geq7.2\) as an event with an integer threshold. Which threshold may give a tighter Chernoff bound?
  5. Explain why an optimized Chernoff bound need not equal the exact Poisson tail probability.