# Law of large numbers

(Redirected from Borel's law of large numbers)
An illustration of the law of large numbers using a particular run of rolls of a single die. As the number of rolls in this run increases, the average of the values of all the results approaches 3.5. While different runs would show a different shape over a small number of throws (at the left), over a large number of rolls (to the right) they would be extremely similar.

In probability theory, the law of large numbers (LLN) is a theorem that describes the result of performing the same experiment a large number of times. According to the law, the average of the results obtained from a large number of trials should be close to the expected value, and will tend to become closer as more trials are performed.

The LLN is important because it "guarantees" stable long-term results for the averages of some random events. For example, while a casino may lose money in a single spin of the roulette wheel, its earnings will tend towards a predictable percentage over a large number of spins. Any winning streak by a player will eventually be overcome by the parameters of the game. It is important to remember that the LLN only applies (as the name indicates) when a large number of observations is considered. There is no principle that a small number of observations will coincide with the expected value or that a streak of one value will immediately be "balanced" by the others (see the gambler's fallacy).

## Examples

For example, a single roll of a fair, six-sided die produces one of the numbers 1, 2, 3, 4, 5, or 6, each with equal probability. Therefore, the expected value of a single die roll is

${\displaystyle {\frac {1+2+3+4+5+6}{6}}=3.5}$

According to the law of large numbers, if a large number of six-sided dice are rolled, the average of their values (sometimes called the sample mean) is likely to be close to 3.5, with the precision increasing as more dice are rolled.

It follows from the law of large numbers that the empirical probability of success in a series of Bernoulli trials will converge to the theoretical probability. For a Bernoulli random variable, the expected value is the theoretical probability of success, and the average of n such variables (assuming they are independent and identically distributed (i.i.d.)) is precisely the relative frequency.

For example, a fair coin toss is a Bernoulli trial. When a fair coin is flipped once, the theoretical probability that the outcome will be heads is equal to 1/2. Therefore, according to the law of large numbers, the proportion of heads in a "large" number of coin flips "should be" roughly 1/2. In particular, the proportion of heads after n flips will almost surely converge to 1/2 as n approaches infinity.

Though the proportion of heads (and tails) approaches 1/2, almost surely the absolute difference in the number of heads and tails will become large as the number of flips becomes large. That is, the probability that the absolute difference is a small number, approaches zero as the number of flips becomes large. Also, almost surely the ratio of the absolute difference to the number of flips will approach zero. Intuitively, expected absolute difference grows, but at a slower rate than the number of flips, as the number of flips grows.

## History

Diffusion is an example of the law of large numbers, applied to chemistry. Initially, there are solute molecules on the left side of a barrier (purple line) and none on the right. The barrier is removed, and the solute diffuses to fill the whole container.
Top: With a single molecule, the motion appears to be quite random.
Middle: With more molecules, there is clearly a trend where the solute fills the container more and more uniformly, but there are also random fluctuations.
Bottom: With an enormous number of solute molecules (too many to see), the randomness is essentially gone: The solute appears to move smoothly and systematically from high-concentration areas to low-concentration areas. In realistic situations, chemists can describe diffusion as a deterministic macroscopic phenomenon (see Fick's laws), despite its underlying random nature.

The Italian mathematician Gerolamo Cardano (1501–1576) stated without proof that the accuracies of empirical statistics tend to improve with the number of trials.[1] This was then formalized as a law of large numbers. A special form of the LLN (for a binary random variable) was first proved by Jacob Bernoulli.[2] It took him over 20 years to develop a sufficiently rigorous mathematical proof which was published in his Ars Conjectandi (The Art of Conjecturing) in 1713. He named this his "Golden Theorem" but it became generally known as "Bernoulli's Theorem". This should not be confused with Bernoulli's principle, named after Jacob Bernoulli's nephew Daniel Bernoulli. In 1837, S.D. Poisson further described it under the name "la loi des grands nombres" ("The law of large numbers").[3][4] Thereafter, it was known under both names, but the "Law of large numbers" is most frequently used.

After Bernoulli and Poisson published their efforts, other mathematicians also contributed to refinement of the law, including Chebyshev,[5] Markov, Borel, Cantelli and Kolmogorov and Khinchin. Markov showed that the law can apply to a random variable that does not have a finite variance under some other weaker assumption, and Khinchin showed in 1929 that if the series consists of independent identically distributed random variables, it suffices that the expected value exists for the weak law of large numbers to be true.[6][7] These further studies have given rise to two prominent forms of the LLN. One is called the "weak" law and the other the "strong" law, in reference to two different modes of convergence of the cumulative sample means to the expected value; in particular, as explained below, the strong form implies the weak.[6]

## Forms

Two different versions of the law of large numbers are described below; they are called the strong law of large numbers, and the weak law of large numbers. Stated for the case where X1, X2, ... is an infinite sequence of i.i.d. Lebesgue integrable random variables with expected value E(X1) = E(X2) = ...= µ, both versions of the law state that – with virtual certainty – the sample average

${\displaystyle {\overline {X}}_{n}={\frac {1}{n}}(X_{1}+\cdots +X_{n})}$

converges to the expected value

${\displaystyle {\begin{matrix}{}\\{\overline {X}}_{n}\,\to \,\mu \qquad {\textrm {for}}\qquad n\to \infty ,\\{}\end{matrix}}}$

(law. 1)

(Lebesgue integrability of Xj means that the expected value E(Xj) exists according to Lebesgue integration and is finite.)

An assumption of finite variance Var(X1) = Var(X2) = ... = σ2 < ∞ is not necessary. Large or infinite variance will make the convergence slower, but the LLN holds anyway. This assumption is often used because it makes the proofs easier and shorter.

The difference between the strong and the weak version is concerned with the mode of convergence being asserted. For interpretation of these modes, see Convergence of random variables.

### Weak law

Simulation illustrating the law of large numbers. Each frame, a coin that is red on one side and blue on the other is flipped, and a dot is added in the corresponding column. A pie chart shows the proportion of red and blue so far. Notice that while the proportion varies significantly at first, it approaches 50% as the number of trials increases.

The weak law of large numbers (also called Khintchine's law) states that the sample average converges in probability towards the expected value[8]

${\displaystyle {\begin{matrix}{}\\{\overline {X}}_{n}\ {\xrightarrow {P}}\ \mu \qquad {\textrm {when}}\ n\to \infty .\\{}\end{matrix}}}$

(law. 2)

That is to say that for any positive number ε,

${\displaystyle \lim _{n\to \infty }\Pr \!\left(\,|{\overline {X}}_{n}-\mu |>\varepsilon \,\right)=0.}$

Interpreting this result, the weak law essentially states that for any nonzero margin specified, no matter how small, with a sufficiently large sample there will be a very high probability that the average of the observations will be close to the expected value; that is, within the margin.

Convergence in probability is also called weak convergence of random variables. This version is called the weak law because random variables may converge weakly (in probability) as above without converging strongly (almost surely) as below.

As mentioned earlier, the weak law applies in the case of independent identically distributed random variables having an expected value. But it also applies in some other cases. For example, the variance may be different for each random variable in the series, keeping the expected value constant. If the variances are bounded, then the law applies, as shown by Chebyshev as early as 1867. (If the expected values change during the series, then we can simply apply the law to the average deviation from the respective expected values. The law then states that this converges in probability to zero.) In fact, Chebyshev's proof works so long as the variance of the average of the first n values goes to zero as n goes to infinity.[7] As an example, assume that each random variable in the series follows a Gaussian distribution with mean zero, but with variance equal to ${\displaystyle 2n/\log(n+1).}$ At each stage, the average will be normally distributed (since it is the average of a set of normally distributed variables). The variance of the sum is equal to the sum of the variances, which is asymptotic to ${\displaystyle n^{2}/\log n}$. The variance of the average is therefore asymptotic to ${\displaystyle 1/\log n}$ and goes to zero.

An example where the law of large numbers does not apply is the Cauchy distribution. Let the random numbers equal the tangent of an angle uniformly distributed between −90° and +90°. The median is zero, but the expected value does not exist, and indeed the average of n such variables has the same distribution as one such variable. It does not tend toward zero as n goes to infinity.

There are also examples of the weak law applying even though the expected value does not exist. See #Differences between the weak law and the strong law.

### Strong law

The strong law of large numbers states that the sample average converges almost surely to the expected value[9]

${\displaystyle {\begin{matrix}{}\\{\bar {X}}_{n}\ {\xrightarrow {\text{a.s.}}}\ \mu \qquad {\textrm {when}}\ n\to \infty .\\{}\end{matrix}}}$

(law. 3)

That is,

${\displaystyle \Pr \!\left(\lim _{n\to \infty }{\bar {X}}_{n}=\mu \right)=1.}$

What this means is that as the number of trials n goes to infinity, the probability that the average of the observations is equal to the expected value will be equal to one.

The proof is more complex than that of the weak law.[10] This law justifies the intuitive interpretation of the expected value (for Lebesgue integration only) of a random variable when sampled repeatedly as the "long-term average".

Almost sure convergence is also called strong convergence of random variables. This version is called the strong law because random variables which converge strongly (almost surely) are guaranteed to converge weakly (in probability). The strong law implies the weak law but not vice versa, when the strong law conditions hold the variable converges both strongly (almost surely) and weakly (in probability). However the weak law may hold in conditions where the strong law does not hold and then the convergence is only weak (in probability).[clarification needed]

To date it has not been possible to prove that the strong law conditions are the same as those of the weak law.[citation needed]

The strong law of large numbers can itself be seen as a special case of the pointwise ergodic theorem.

The strong law applies to independent identically distributed random variables having an expected value (like the weak law). This was proved by Kolmogorov in 1930. It can also apply in other cases. Kolmogorov also showed, in 1933, that if the variables are independent and identically distributed, then for the average to converge almost surely on something (this can be considered another statement of the strong law), it is necessary that they have an expected value (and then of course the average will converge almost surely on that).[11]

If the summands are independent but not identically distributed, then

${\displaystyle {\bar {X}}_{n}-\operatorname {E} {\big [}{\bar {X}}_{n}{\big ]}\ {\xrightarrow {\text{a.s.}}}\ 0,}$

provided that each Xk has a finite second moment and

${\displaystyle \sum _{k=1}^{\infty }{\frac {1}{k^{2}}}\operatorname {Var} [X_{k}]<\infty .}$

This statement is known as Kolmogorov's strong law, see e.g. Sen & Singer (1993, Theorem 2.3.10).

An example of a series where the weak law applies but not the strong law is when Xk is plus or minus ${\displaystyle {\sqrt {k/\log \log \log k}}}$ (starting at sufficiently large k so that the denominator is positive) with probability 1/2 for each.[11] The variance of Xk is then ${\displaystyle k/\log \log \log k.}$ Kolmogorov's strong law does not apply because the partial sum in his criterion up to k=n is asymptotic to ${\displaystyle \log n/\log \log \log n}$ and this is unbounded.

If we replace the random variables with Gaussian variables having the same variances, namely ${\displaystyle {\sqrt {k/\log \log \log k}},}$ then the average at any point will also be normally distributed. The width of the distribution of the average will tend toward zero (standard deviation asymptotic to ${\displaystyle 1/{\sqrt {2\log \log \log n}}}$), but for a given ε, there is probability which does not go to zero with n that the average sometime after the nth trial will come back up to ε. Since this probability does not go to zero, it must have a positive lower bound p(ε), which means there is a probability of at least p(ε) that the average will attain ε after n trials. It will happen with probability p(ε)/2 before some m which depends on n. But even after m, there is still a probability of at least p(ε) that it will happen. (This seems to indicate that p(ε)=1 and the average will attain ε an infinite number of times.)

### Differences between the weak law and the strong law

The weak law states that for a specified large n, the average ${\displaystyle {\overline {X}}_{n}}$ is likely to be near μ. Thus, it leaves open the possibility that ${\displaystyle |{\overline {X}}_{n}-\mu |>\varepsilon }$ happens an infinite number of times, although at infrequent intervals. (Not necessarily ${\displaystyle |{\overline {X}}_{n}-\mu |\neq 0}$ for all n ).

The strong law shows that this almost surely will not occur. In particular, it implies that with probability 1, we have that for any ε > 0 the inequality ${\displaystyle |{\overline {X}}_{n}-\mu |<\varepsilon }$ holds for all large enough n.[12]

The strong law does not hold in the following cases, but the weak law does.[13][14][15]

1. Let x be an exponentially distributed random variable with parameter 1. The random variable ${\displaystyle {\frac {\sin(x)e^{x}}{x}}}$ has no expected value according to Lebesgue integration, but using conditional convergence and interpreting the integral as a Dirichlet integral, which is an improper Riemann integral, we can say:

${\displaystyle E\left({\frac {\sin(x)e^{x}}{x}}\right)=\ \int _{0}^{\infty }{\frac {\sin(x)e^{x}}{x}}e^{-x}dx={\frac {\pi }{2}}}$

2. Let x be geometric distribution with probability 0.5. The random variable ${\displaystyle {\frac {2^{x}(-1)^{x}}{x}}}$ does not have an expected value in the conventional sense because the infinite series is not absolutely convergent, but using conditional convergence, we can say:

${\displaystyle E\left({\frac {2^{x}(-1)^{x}}{x}}\right)=\ \sum _{1}^{\infty }{\frac {2^{x}(-1)^{x}}{x}}2^{-x}=-\ln(2)}$

3. If the cumulative distribution function of a random variable is

${\displaystyle 1-F(x)={\frac {e}{2x\ln(x)}},x\geq e}$
${\displaystyle F(x)={\frac {e}{-2x\ln(-x)}},x\leq -e}$
then it has no expected value, but the weak law is true.[16][17]

### Uniform law of large numbers

Suppose f(x,θ) is some function defined for θ ∈ Θ, and continuous in θ. Then for any fixed θ, the sequence {f(X1,θ), f(X2,θ), …} will be a sequence of independent and identically distributed random variables, such that the sample mean of this sequence converges in probability to E[f(X,θ)]. This is the pointwise (in θ) convergence.

The uniform law of large numbers states the conditions under which the convergence happens uniformly in θ. If[18][19]

1. Θ is compact,
2. f(x,θ) is continuous at each θ ∈ Θ for almost all xs, and measurable function of x at each θ.
3. there exists a dominating function d(x) such that E[d(X)] < ∞, and
${\displaystyle \left\|f(x,\theta )\right\|\leq d(x)\quad {\text{for all}}\ \theta \in \Theta .}$

Then E[f(X,θ)] is continuous in θ, and

${\displaystyle \sup _{\theta \in \Theta }\left\|{\frac {1}{n}}\sum _{i=1}^{n}f(X_{i},\theta )-\operatorname {E} [f(X,\theta )]\right\|{\xrightarrow {\mathrm {a.s.} }}\ 0.}$

This result is useful to derive consistency of a large class of estimators (see Extremum estimator).

### Borel's law of large numbers

Borel's law of large numbers, named after Émile Borel, states that if an experiment is repeated a large number of times, independently under identical conditions, then the proportion of times that any specified event occurs approximately equals the probability of the event's occurrence on any particular trial; the larger the number of repetitions, the better the approximation tends to be. More precisely, if E denotes the event in question, p its probability of occurrence, and Nn(E) the number of times E occurs in the first n trials, then with probability one,[20]

${\displaystyle {\frac {N_{n}(E)}{n}}\to p{\text{ as }}n\to \infty .\,}$

Chebyshev's inequality. Let X be a random variable with finite expected value μ and finite non-zero variance σ2. Then for any real number k > 0,

${\displaystyle \Pr(|X-\mu |\geq k\sigma )\leq {\frac {1}{k^{2}}}.}$

This theorem makes rigorous the intuitive notion of probability as the long-run relative frequency of an event's occurrence. It is a special case of any of several more general laws of large numbers in probability theory.

## Proof of the weak law

Given X1, X2, ... an infinite sequence of i.i.d. random variables with finite expected value E(X1) = E(X2) = ... = µ < ∞, we are interested in the convergence of the sample average

${\displaystyle {\overline {X}}_{n}={\tfrac {1}{n}}(X_{1}+\cdots +X_{n}).}$

The weak law of large numbers states:

Theorem: ${\displaystyle {\begin{matrix}{}\\{\overline {X}}_{n}\ {\xrightarrow {P}}\ \mu \qquad {\textrm {when}}\ n\to \infty .\\{}\end{matrix}}}$

(law. 2)

### Proof using Chebyshev's inequality assuming finite variance

This proof uses the assumption of finite variance ${\displaystyle \operatorname {Var} (X_{i})=\sigma ^{2}}$ (for all ${\displaystyle i}$). The independence of the random variables implies no correlation between them, and we have that

${\displaystyle \operatorname {Var} ({\overline {X}}_{n})=\operatorname {Var} ({\tfrac {1}{n}}(X_{1}+\cdots +X_{n}))={\frac {1}{n^{2}}}\operatorname {Var} (X_{1}+\cdots +X_{n})={\frac {n\sigma ^{2}}{n^{2}}}={\frac {\sigma ^{2}}{n}}.}$

The common mean μ of the sequence is the mean of the sample average:

${\displaystyle E({\overline {X}}_{n})=\mu .}$

Using Chebyshev's inequality on ${\displaystyle {\overline {X}}_{n}}$ results in

${\displaystyle \operatorname {P} (\left|{\overline {X}}_{n}-\mu \right|\geq \varepsilon )\leq {\frac {\sigma ^{2}}{n\varepsilon ^{2}}}.}$

This may be used to obtain the following:

${\displaystyle \operatorname {P} (\left|{\overline {X}}_{n}-\mu \right|<\varepsilon )=1-\operatorname {P} (\left|{\overline {X}}_{n}-\mu \right|\geq \varepsilon )\geq 1-{\frac {\sigma ^{2}}{n\varepsilon ^{2}}}.}$

As n approaches infinity, the expression approaches 1. And by definition of convergence in probability, we have obtained

${\displaystyle {\begin{matrix}{}\\{\overline {X}}_{n}\ {\xrightarrow {P}}\ \mu \qquad {\textrm {when}}\ n\to \infty .\\{}\end{matrix}}}$

(law. 2)

### Proof using convergence of characteristic functions

By Taylor's theorem for complex functions, the characteristic function of any random variable, X, with finite mean μ, can be written as

${\displaystyle \varphi _{X}(t)=1+it\mu +o(t),\quad t\rightarrow 0.}$

All X1, X2, ... have the same characteristic function, so we will simply denote this φX.

Among the basic properties of characteristic functions there are

${\displaystyle \varphi _{{\frac {1}{n}}X}(t)=\varphi _{X}({\tfrac {t}{n}})\quad {\text{and}}\quad \varphi _{X+Y}(t)=\varphi _{X}(t)\varphi _{Y}(t)\quad }$ if X and Y are independent.

These rules can be used to calculate the characteristic function of ${\displaystyle \scriptstyle {\overline {X}}_{n}}$ in terms of φX:

${\displaystyle \varphi _{{\overline {X}}_{n}}(t)=\left[\varphi _{X}\left({t \over n}\right)\right]^{n}=\left[1+i\mu {t \over n}+o\left({t \over n}\right)\right]^{n}\,\rightarrow \,e^{it\mu },\quad {\text{as}}\quad n\rightarrow \infty .}$

The limit  eitμ  is the characteristic function of the constant random variable μ, and hence by the Lévy continuity theorem, ${\displaystyle \scriptstyle {\overline {X}}_{n}}$ converges in distribution to μ:

${\displaystyle {\overline {X}}_{n}\,{\xrightarrow {\mathcal {D}}}\,\mu \qquad {\text{for}}\qquad n\to \infty .}$

μ is a constant, which implies that convergence in distribution to μ and convergence in probability to μ are equivalent (see Convergence of random variables.) Therefore,

${\displaystyle {\begin{matrix}{}\\{\overline {X}}_{n}\ {\xrightarrow {P}}\ \mu \qquad {\textrm {when}}\ n\to \infty .\\{}\end{matrix}}}$

(law. 2)

This shows that the sample mean converges in probability to the derivative of the characteristic function at the origin, as long as the latter exists.

## Notes

1. ^ Mlodinow, L. The Drunkard's Walk. New York: Random House, 2008. p. 50.
2. ^ Jakob Bernoulli, Ars Conjectandi: Usum & Applicationem Praecedentis Doctrinae in Civilibus, Moralibus & Oeconomicis, 1713, Chapter 4, (Translated into English by Oscar Sheynin)
3. ^ Poisson names the "law of large numbers" (la loi des grands nombres) in: S.D. Poisson, Probabilité des jugements en matière criminelle et en matière civile, précédées des règles générales du calcul des probabilitiés (Paris, France: Bachelier, 1837), p. 7. He attempts a two-part proof of the law on pp. 139–143 and pp. 277 ff.
4. ^ Hacking, Ian. (1983) "19th-century Cracks in the Concept of Determinism", Journal of the History of Ideas, 44 (3), 455-475 JSTOR 2709176
5. ^ Tchebichef, P. (1846). "Démonstration élémentaire d'une proposition générale de la théorie des probabilités". Journal für die reine und angewandte Mathematik (Crelles Journal). 1846 (33): 259–267. doi:10.1515/crll.1846.33.259.
6. ^ a b
7. ^ a b Yuri Prohorov. "Law of large numbers". Encyclopedia of Mathematics.
8. ^ Loève 1977, Chapter 1.4, p. 14
9. ^ Loève 1977, Chapter 17.3, p. 251
10. ^ "The strong law of large numbers – What's new". Terrytao.wordpress.com. Retrieved 2012-06-09.
11. ^ a b Yuri Prokhorov. "Strong law of large numbers". Encyclopedia of Mathematics.
12. ^ Ross (2009)
13. ^
14. ^ "A NOTE ON THE WEAK LAW OF LARGE NUMBERS FOR EXCHANGEABLE RANDOM VARIABLES" (PDF). Dguvl Hun Hong and Sung Ho Lee.
15. ^
16. ^ Mukherjee, Sayan. "Law of large numbers" (PDF).
17. ^ J. Geyer, Charles. "Law of large numbers" (PDF).
18. ^ Newey & McFadden 1994, Lemma 2.4
19. ^ Jennrich, Robert I. (1969). "Asymptotic Properties of Non-Linear Least Squares Estimators". The Annals of Mathematical Statistics. 40 (2): 633–643. doi:10.1214/aoms/1177697731.
20. ^ An Analytic Technique to Prove Borel's Strong Law of Large Numbers Wen, L. Am Math Month 1991

## References

• Grimmett, G. R.; Stirzaker, D. R. (1992). Probability and Random Processes, 2nd Edition. Clarendon Press, Oxford. ISBN 0-19-853665-8.
• Richard Durrett (1995). Probability: Theory and Examples, 2nd Edition. Duxbury Press.
• Martin Jacobsen (1992). Videregående Sandsynlighedsregning (Advanced Probability Theory) 3rd Edition. HCØ-tryk, Copenhagen. ISBN 87-91180-71-6.
• Loève, Michel (1977). Probability theory 1 (4th ed.). Springer Verlag.
• Newey, Whitney K.; McFadden, Daniel (1994). Large sample estimation and hypothesis testing. Handbook of econometrics, vol. IV, Ch. 36. Elsevier Science. pp. 2111–2245.
• Ross, Sheldon (2009). A first course in probability (8th ed.). Prentice Hall press. ISBN 978-0-13-603313-4.
• Sen, P. K; Singer, J. M. (1993). Large sample methods in statistics. Chapman & Hall, Inc.
• Seneta, Eugene (2013), "A Tricentenary history of the Law of Large Numbers", Bernoulli, 19 (4): 1088–1121, doi:10.3150/12-BEJSP12