Pairwise independence

From testwiki
Jump to navigation Jump to search

Template:Short description In probability theory, a pairwise independent collection of random variables is a set of random variables any two of which are independent.[1] Any collection of mutually independent random variables is pairwise independent, but some pairwise independent collections are not mutually independent. Pairwise independent random variables with finite variance are uncorrelated.

A pair of random variables X and Y are independent if and only if the random vector (X, Y) with joint cumulative distribution function (CDF) FX,Y(x,y) satisfies

FX,Y(x,y)=FX(x)FY(y),

or equivalently, their joint density fX,Y(x,y) satisfies

fX,Y(x,y)=fX(x)fY(y).

That is, the joint distribution is equal to the product of the marginal distributions.[2]

Unless it is not clear in context, in practice the modifier "mutual" is usually dropped so that independence means mutual independence. A statement such as " X, Y, Z are independent random variables" means that X, Y, Z are mutually independent.

Example

Pairwise independence does not imply mutual independence, as shown by the following example attributed to S. Bernstein.[3]

Suppose X and Y are two independent tosses of a fair coin, where we designate 1 for heads and 0 for tails. Let the third random variable Z be equal to 1 if exactly one of those coin tosses resulted in "heads", and 0 otherwise (i.e., Z=XY). Then jointly the triple (X, Y, Z) has the following probability distribution:

(X,Y,Z)={(0,0,0)with probability 1/4,(0,1,1)with probability 1/4,(1,0,1)with probability 1/4,(1,1,0)with probability 1/4.

Here the marginal probability distributions are identical: fX(0)=fY(0)=fZ(0)=1/2, and fX(1)=fY(1)=fZ(1)=1/2. The bivariate distributions also agree: fX,Y=fX,Z=fY,Z, where fX,Y(0,0)=fX,Y(0,1)=fX,Y(1,0)=fX,Y(1,1)=1/4.

Since each of the pairwise joint distributions equals the product of their respective marginal distributions, the variables are pairwise independent:

  • X and Y are independent, and
  • X and Z are independent, and
  • Y and Z are independent.

However, X, Y, and Z are not mutually independent, since fX,Y,Z(x,y,z)fX(x)fY(y)fZ(z), the left side equalling for example 1/4 for (x, y, z) = (0, 0, 0) while the right side equals 1/8 for (x, y, z) = (0, 0, 0). In fact, any of {X,Y,Z} is completely determined by the other two (any of X, Y, Z is the sum (modulo 2) of the others). That is as far from independence as random variables can get.

Probability of the union of pairwise independent events

Bounds on the probability that the sum of Bernoulli random variables is at least one, commonly known as the union bound, are provided by the Boole–Fréchet[4][5] inequalities. While these bounds assume only univariate information, several bounds with knowledge of general bivariate probabilities, have been proposed too. Denote by {Ai,i{1,2,...,n}} a set of n Bernoulli events with probability of occurrence (Ai)=pi for each i. Suppose the bivariate probabilities are given by (AiAj)=pij for every pair of indices (i,j). Kounias [6] derived the following upper bound:

(iAi)i=1npimaxj{1,2,..,n}ijpij,

which subtracts the maximum weight of a star spanning tree on a complete graph with n nodes (where the edge weights are given by pij) from the sum of the marginal probabilities ipi.
Hunter-Worsley[7][8] tightened this upper bound by optimizing over τT as follows:

(iAi)i=1npimaxτT(i,j)τpij,

where T is the set of all spanning trees on the graph. These bounds are not the tightest possible with general bivariates pij even when feasibility is guaranteed as shown in Boros et.al.[9] However, when the variables are pairwise independent (pij=pipj), Ramachandra—Natarajan [10] showed that the Kounias-Hunter-Worsley [6][7][8] bound is tight by proving that the maximum probability of the union of events admits a closed-form expression given as:

Template:NumBlk

where the probabilities are sorted in increasing order as 0p1p2pn1. The tight bound in Template:EquationNote depends only on the sum of the smallest n1 probabilities i=1n1pi and the largest probability pn. Thus, while ordering of the probabilities plays a role in the derivation of the bound, the ordering among the smallest n1 probabilities {p1,p2,...,pn1} is inconsequential since only their sum is used.

Comparison with the Boole–Fréchet union bound

It is useful to compare the smallest bounds on the probability of the union with arbitrary dependence and pairwise independence respectively. The tightest Boole–Fréchet upper union bound (assuming only univariate information) is given as:
Template:NumBlk As shown in Ramachandra-Natarajan,[10] it can be easily verified that the ratio of the two tight bounds in Template:EquationNote and Template:EquationNote is upper bounded by 4/3 where the maximum value of 4/3 is attained when

i=1n1pi=1/2, pn=1/2

where the probabilities are sorted in increasing order as 0p1p2pn1. In other words, in the best-case scenario, the pairwise independence bound in Template:EquationNote provides an improvement of 25% over the univariate bound in Template:EquationNote.

Generalization

More generally, we can talk about k-wise independence, for any k ≥ 2. The idea is similar: a set of random variables is k-wise independent if every subset of size k of those variables is independent. k-wise independence has been used in theoretical computer science, where it was used to prove a theorem about the problem MAXEkSAT.

k-wise independence is used in the proof that k-independent hashing functions are secure unforgeable message authentication codes.

See also

References

Template:Reflist

  1. Gut, A. (2005) Probability: a Graduate Course, Springer-Verlag. Template:Isbn. pp. 71–72.
  2. Template:Cite book Definition 2.5.1, page 109.
  3. Template:Cite book Remark 2.6.1, p. 120.
  4. Boole, G. (1854). An Investigation of the Laws of Thought, On Which Are Founded the Mathematical Theories of Logic and Probability. Walton and Maberly, London. See Boole's "major" and "minor" limits of a conjunction on page 299.
  5. Fréchet, M. (1935). Généralisations du théorème des probabilités totales. Fundamenta Mathematicae 25: 379–387.
  6. 6.0 6.1 Template:Cite journal
  7. 7.0 7.1 Template:Cite journal
  8. 8.0 8.1 Template:Cite journal
  9. Template:Cite journal
  10. 10.0 10.1 Template:Cite journal