Odwrotna granica Chernoffa


31

Czy istnieje odwrotna granica Chernoffa, która ogranicza, że ​​prawdopodobieństwo ogona jest co najmniej tak duże.

tj. jeśli X 1 , X 2 , , X nX1,X2,,Xn są niezależnymi dwumianowymi zmiennymi losowymi, a μ = E [ n i = 1 X i ]μ=E[ni=1Xi] . Czy możemy zatem udowodnić P r [ n i = 1 X i( 1 + δ ) μ ] f ( μ , δ , n )Pr[ni=1Xi(1+δ)μ]f(μ,δ,n) dla niektórych funkcji ff .


1
Twój przykład o zbyt wiele: z p = n - 2 / 3p=n2/3 , standardową nierówność chernoffa pokazuje, że Pr [ | T S 1 | 1,1 n 1 / 3 ]Pr[|TS1|1.1n1/3]iPr[| TS2| 1,1n 1 / 3 ]Pr[|TS2|1.1n1/3]jest najwyżejexp(-cn 1 / 3 ),exp(cn1/3)przez pewienCc.
Colin McQuillan

You are right, I got confused about which term in chernoff bound has the square. I have changed the question to reflect a weaker bound. I don't think it will help me in my current application, but it might be interesting for other reasons.
Ashwinkumar B V

Odpowiedzi:


28

Here is an explicit proof that a standard Chernoff bound is tight up to constant factors in the exponent for a particular range of the parameters. (In particular, whenever the variables are 0 or 1, and 1 with probability 1/2 or less, and ϵ(0,1/2)ϵ(0,1/2), and the Chernoff upper bound is less than a constant.)

If you find a mistake, please let me know.

Lemma 1. (tightness of Chernoff bound) Let XX be the average of kk independent, 0/1 random variables (r.v.). For any ϵ(0,1/2]ϵ(0,1/2] and p(0,1/2]p(0,1/2], assuming ϵ2pk3ϵ2pk3,

(i) If each r.v. is 1 with probability at most pp, then Pr[X(1ϵ)p]  exp(9ϵ2pk).

Pr[X(1ϵ)p]  exp(9ϵ2pk).

(ii) If each r.v. is 1 with probability at least pp, then Pr[X(1+ϵ)p]  exp(9ϵ2pk).

Pr[X(1+ϵ)p]  exp(9ϵ2pk).

Proof. We use the following observation:

Claim 1. If 1k11k1, then (k)  1e2π(k)(kk)k(k)  1e2π(k)(kk)k

Proof of Claim 1. By Stirling's approximation, i!=2πi(i/e)ieλi!=2πi(i/e)ieλ where λ[1/(12i+1),1/12i].λ[1/(12i+1),1/12i].

Thus, (k)(k), which is k!!(k)!k!!(k)!, is at least 2πk(ke)k2π(e)  2π(k)(ke)kexp(112k+1112112(k))

2πk(ke)k2π(e)  2π(k)(ke)kexp(112k+1112112(k))
  12π(k)(kk)ke1.
  12π(k)(kk)ke1.
QED

Proof of Lemma 1 Part (i). Without loss of generality assume each 0/1 random variable in the sum XX is 1 with probability exactly pp. Note Pr[X(1ϵ)p]Pr[X(1ϵ)p] equals the sum (1ϵ)pki=0Pr[X=i/k](1ϵ)pki=0Pr[X=i/k], and Pr[X=i/k]=(ki)pi(1p)kiPr[X=i/k]=(ki)pi(1p)ki.

Fix =(12ϵ)pk+1=(12ϵ)pk+1. The terms in the sum are increasing, so the terms with index ii each have value at least Pr[X=/k]Pr[X=/k], so their sum has total value at least (ϵpk2)Pr[X=/k](ϵpk2)Pr[X=/k]. To complete the proof, we show that (ϵpk2)Pr[X=/k]  exp(9ϵ2pk).

(ϵpk2)Pr[X=/k]  exp(9ϵ2pk).

The assumptions ϵ2pk3ϵ2pk3 and ϵ1/2ϵ1/2 give ϵpk6ϵpk6, so the left-hand side above is at least 23ϵpk(k)p(1p)k23ϵpk(k)p(1p)k. Using Claim 1, to bound (k)(k), this is in turn at least ABAB where A=23eϵpk/2πA=23eϵpk/2π and B=(k)(kk)kp(1p)k.B=(k)(kk)kp(1p)k.

To finish we show Aexp(ϵ2pk)Aexp(ϵ2pk) and Bexp(8ϵ2pk)Bexp(8ϵ2pk).

Claim 2. Aexp(ϵ2pk)Aexp(ϵ2pk)

Proof of Claim 2. The assumptions ϵ2pk3ϵ2pk3 and ϵ1/2ϵ1/2 imply (i) pk12pk12.

By definition, pk+1pk+1. By (i), pk12pk12. Thus, (ii) 1.1pk1.1pk.

Substituting the right-hand side of (ii) for in AA gives (iii) A23eϵpk/2.2πA23eϵpk/2.2π.

The assumption, ϵ2pk3ϵ2pk3, implies ϵpk3ϵpk3, which with (iii) gives (iv) A23e3/2.2π0.1A23e3/2.2π0.1.

From ϵ2pk3ϵ2pk3 it follows that (v) exp(ϵ2pk)exp(3)0.04exp(ϵ2pk)exp(3)0.04.

(iv) and (v) together give the claim. QED

Claim 3. Bexp(8ϵ2pk)Bexp(8ϵ2pk).

Proof of Claim 3. Fix δδ such that =(1δ)pk=(1δ)pk.
The choice of implies δ2ϵδ2ϵ, so the claim will hold as long as Bexp(2δ2pk)Bexp(2δ2pk). Taking each side of this latter inequality to the power 1/1/ and simplifying, it is equivalent to pk(k(1p)k)k/1  exp(2δ2pk).

pk(k(1p)k)k/1  exp(2δ2pk).
Substituting =(1δ)pk=(1δ)pk and simplifying, it is equivalent to (1δ)(1+δp1p)1(1δ)p1  exp(2δ21δ).
(1δ)(1+δp1p)1(1δ)p1  exp(2δ21δ).
Taking the logarithm of both sides and using ln(1+z)zln(1+z)z twice, it will hold as long as δ+δp1p(1(1δ)p1)  2δ21δ.
δ+δp1p(1(1δ)p1)  2δ21δ.
The left-hand side above simplifies to δ2/(1p)(1δ)δ2/(1p)(1δ), which is less than 2δ2/(1δ)2δ2/(1δ) because p1/2p1/2. QED

Claims 2 and 3 imply ABexp(ϵ2pk)exp(8ϵ2pk)ABexp(ϵ2pk)exp(8ϵ2pk). This implies part (i) of the lemma.

Proof of Lemma 1 Part (ii). Without loss of generality assume each random variable is 11 with probability exactly pp.

Note Pr[X(1+ϵ)p]=ni=(1ϵ)pkPr[X=i/k]Pr[X(1+ϵ)p]=ni=(1ϵ)pkPr[X=i/k]. Fix ˆ=(1+2ϵ)pk1^=(1+2ϵ)pk1.

The last ϵpkϵpk terms in the sum total at least (ϵpk2)Pr[X=ˆ/k](ϵpk2)Pr[X=^/k], which is at least exp(9ϵ2pk)exp(9ϵ2pk). (The proof of that is the same as for (i), except with replaced by ˆ^ and δδ replaced by ˆδδ^ such that ˆ=(1+ˆδ)pk^=(1+δ^)pk.) QED


Several [math processing error]s -- any chance of fixing them?
Aryeh

Those math expressions used to display just fine. For some reason the \choose command is not working in mathjax. Neither is \binom. E.g. $a \choose b$ gives (ab). Presumably this is a bug in the mathjax configuration. Hopefully it will be fixed soon. Meanwhile see Lemma 5.2 in the appendix of arxiv.org/pdf/cs/0205046v2.pdf or cs.ucr.edu/~neal/Klein15Number.
Neal Young

22

The Berry-Esseen theorem can give tail probability lower bounds, as long as they are higher than n1/2.

Another tool you can use is the Paley-Zygmund inequality. It implies that for any even integer k, and any real-valued random variable X,

Pr[|X|>=12(E[Xk])1/k]E[Xk]24E[X2k]

Together with the multinomial theorem, for X a sum of n rademacher random variables Paley-Zygmund can get you pretty strong lower bounds. Also it works with bounded-independence random variables. For example you easily get that the sum of n 4-wise independent ±1 random variables is Ω(n) with constant probability.


14

If you are indeed okay with bounding sums of Bernoulli trials (and not, say, bounded random variables), the following is pretty tight.

Slud's Inequality*. Let {Xi}ni=1 be i.i.d. draws from a Bernoulli r.v. with E(X1)=p, and let integer kn be given. If either (a) p1/4 and npk, or (b) npkn(1p), then Pr[iXik]1Φ(knpnp(1p)),

where Φ is the cdf of a standard normal.

(Treating the argument to Φ as transforming the standard normal, this agrees exactly with what the CLT tells you; in fact, it tells us that Binomials satisfying the conditions of the theorem will dominate their corresponding Gaussians on upper tails.)

From here, you can use bounds on Φ to get something nicer. For instance, in Feller's first book, in the section on Gaussians, it is shown for every z>0 that z1+z2φ(z)<1Φ(z)<1zφ(z),

where φ is the density of a standard normal. There are similar bounds in the Wikipedia article for "Q-function" as well.

Other than that, and what other people have said, you can also try using the Binomial directly, perhaps with some Stirling.

(*) Some newer statements of Slud's inequality leave out some of these conditions; I've reproduced the one in Slud's paper.


7

The de Moivre-Laplace Theorem shows that variables like |TS1|, after being suitably normalised and under certain conditions, will converge in distribution to a normal distribution. That's enough if you want constant lower bounds.

For lower bounds like nc, you need a slightly finer tool. Here's one reference I know of (but only by accident - I've never had the opportunity to use such an inequality myself). Some explicit lower bounds on tail probabilities of binomial distributions are given as Theorem 1.5 the book Random graphs by Béla Bollobás, Cambridge, 2nd edition, where further references are given to An introduction to probability and its applications by Feller and Foundations of Probability by Rényi.


4

The Generalized Littlewood-Offord Theorem isn't exactly what you want, but it gives what I think of as a "reverse Chernoff" bound by showing that the sum of random variables is unlikely to fall within a small range around any particular value (including the expectation). Perhaps it will be useful.

Formally, the theorem is as follows.

Generalized Littlewood-Offord Theorem: Let a1,,an, and s>0 be real numbers such that |ai|s for 1in and let X1,,Xn be independent random variables that have values zero and one. For 0<p12, suppose that pPr[Xi=0]1p for all 1in. Then, for any rR, Pr[rni=1aiXi<r+s]cpn

Where cp is a constant depending only on p.

3
It may be helpful to others to know that this type of result is also known as a "small ball inequality" and Nguyen and Vu have a terrific survey people.math.osu.edu/nguyen.1261/cikk/LO-survey.pdf. My perspective here slightly differs from yours. I think of a "reverse Chernoff" bound as giving a lower estimate of the probability mass of the small ball around 0. I think of a small ball inequality as qualitatively saying that the small ball probability is maximized by the ball at 0. In this sense reverse Chernoff bounds are usually easier to prove than small ball inequalities.
Sasho Nikolov

3

The exponent in the standard Chernoff bound as it is stated on Wikipedia is tight for 0/1-valued random variables. Let 0<p<1 and let X1,X2, be a sequence of independent random variables such that for each i, Pr[Xi=1]=p and Pr[Xi=0]=1p. Then for every ε>0, 2D(p+εp)nn+1Pr[ni=1Xi(p+ε)n]2D(p+εp)n.

Here, D(xy)=xlog2(x/y)+(1x)log2((1x)/(1y)), which is the Kullback-Leibler divergence between Bernoulli random variables with parameters x and y.

As mentioned, the upper bound in the inequality above is proved on Wikipedia (https://en.wikipedia.org/wiki/Chernoff_bound) under the name "Chernoff-Hoeffding Theorem, additive form". The lower bound can be proved using e.g. the "method of types". See Lemma II.2 in [1]. Also, this is covered in the classic textbook on information theory by Cover and Thomas.

[1] Imre Csiszár: The Method of Types. IEEE Transactions on Information Theory (1998). http://dx.doi.org/10.1109/18.720546


It is also worth noting that D(p+δpp)=p22pδ2+O(δ3), and for common case of p=1/2 it is 12δ2+O(δ4). This shows that when δ=O(n1/3) the typical eCδ2 bound is sharp. (And when δ=O(n1/4) for p=1/2).
Thomas Ahle
Korzystając z naszej strony potwierdzasz, że przeczytałeś(-aś) i rozumiesz nasze zasady używania plików cookie i zasady ochrony prywatności.
Licensed under cc by-sa 3.0 with attribution required.