Skip to content
Date: Oct 2, 2026
Builds on Chapter 7.

The oldest theorem of probability says that the average of many independent observations is close to their common mean. Probability Theory states it in its weak form: for independent identically distributed variables with finite variance, the probability that the average of the first \(n\) differs from the mean by more than a given amount tends to zero (math.probability-theory.weak-law-large-numbers). The proof is one application of Chebyshev’s inequality, and it leaves three questions that a first course cannot answer.

Is the variance needed? The statement mentions only the mean, and one suspects that the mean is all that should matter. Is the conclusion true of the sequence of averages actually observed? The weak form speaks of each \(n\) separately; it does not exclude that along the realized sequence the average strays far from the mean again and again, as the typewriter sequence strays from zero. And what happens when there is no mean? If \(X_1,X_2,\dots\) are independent Cauchy variables, the average \((X_1+\dots+X_n)/n\) has exactly the Cauchy law again, for every \(n\) : averaging a million observations is no better than keeping one.

This chapter owns the proofs. It proves the weak law under a first moment only, then the strong law, almost sure convergence of the averages, under the same hypothesis and with independence weakened to pairwise independence, and finally the converse: for independent identically distributed variables the averages converge almost surely to a finite limit if and only if the mean exists. Throughout, \(S_n=X_1+\dots+X_n\) .

8.1 Second moments and a rate

Theorem 8.1 (weak law in mean square). Let \(X_1,X_2,\dots\) be pairwise uncorrelated random variables with common mean \(\mu\) and variances bounded by a constant \(C\) . Then \(S_n/n\to\mu\) in \(L^2\) and in probability, and for every \(\varepsilon>0\) ,

\[ \mathbb P\Bigl(\Bigl\lvert\frac{S_n}{n}-\mu\Bigr\rvert\ge\varepsilon\Bigr)\le\frac{C}{n\varepsilon^2}. \]

Proof. By Corollary 2.7, \(\operatorname{Var}(S_n)=\sum_{k\le n}\operatorname{Var}(X_k)\le nC\) , so \(\mathbb E(S_n/n-\mu)^2=\operatorname{Var}(S_n)/n^2\le C/n\to 0\) . The bound on the probability is Chebyshev’s inequality.

Neither independence nor identical distribution was used, only the additivity of variance. The theorem comes with a rate, and the rate is its limitation: the bound \(C/(n\varepsilon^2)\) is not summable in \(n\) , so the first Borel–Cantelli lemma cannot convert it into almost sure convergence.

Example 8.1: What Chebyshev’s bound is worth

Toss a fair coin \(n=100\) times and let \(S_{100}\) be the number of heads. Each toss has variance \(\tfrac14\) , so Theorem 8.1 gives

\[ \mathbb P\Bigl(\Bigl\lvert\frac{S_{100}}{100}-\frac12\Bigr\rvert\ge 0.1\Bigr)\le\frac{1/4}{100\cdot 0.01}=0.25 . \]

The exact probability is that of at most \(40\) or at least \(60\) heads, a binomial sum equal to \(0.0569\) . For \(n=1000\) and \(\varepsilon=0.05\) the bound is \(0.1\) and the exact probability is \(0.0017\) . The bound has the right qualitative behaviour and is far from sharp, because it uses only the variance, while the true tail of a sum of bounded variables decays exponentially.

The bound is nevertheless enough for planning. To guarantee that the frequency of heads is within \(0.1\) of \(\tfrac12\) with probability at least \(0.95\) , it suffices that \(0.25/(n\cdot 0.01)\le 0.05\) , that is, \(n\ge 500\) .

8.2 Fourth moments and almost sure convergence

A stronger moment assumption gives a summable bound, and with it the strong law, in a few lines.

Theorem 8.2 (strong law under a fourth moment). Let \(X_1,X_2,\dots\) be independent and identically distributed with \(\mathbb E X_1^4<\infty\) and mean \(\mu\) . Then \(S_n/n\to\mu\) almost surely.

Proof. Replacing \(X_k\) by \(X_k-\mu\) we may take \(\mu=0\) . Expand \(\mathbb E S_n^4=\sum\mathbb E[X_iX_jX_kX_l]\) over all quadruples of indices. By independence and \(\mathbb E X_i=0\) , a term vanishes if some index occurs exactly once. The surviving terms are the \(n\) terms \(\mathbb E X_i^4\) and the terms \(\mathbb E X_i^2\,\mathbb E X_j^2\) with \(i\ne j\) , of which there are \(3n(n-1)\) : choose the unordered pair \(\{i,j\}\) in \(\tfrac12n(n-1)\) ways and the two positions of \(i\) among the four in \(6\) ways. With \(\sigma^2=\mathbb E X_1^2\) ,

\[ \mathbb E S_n^4=n\,\mathbb E X_1^4+3n(n-1)\,\sigma^4\le 3n^2\,\mathbb E X_1^4, \]

using \(\sigma^4\le\mathbb E X_1^4\) . By Markov’s inequality, \(\mathbb P(\lvert S_n\rvert>n\varepsilon)\le\mathbb E S_n^4/(n\varepsilon)^4\le 3\,\mathbb E X_1^4/(n^2\varepsilon^4)\) . These bounds are summable in \(n\) , so by the first Borel–Cantelli lemma, for each \(\varepsilon>0\) almost surely \(\lvert S_n/n\rvert\le\varepsilon\) for all large \(n\) ; that is the criterion of Lemma 4.2.

The fourth moment grows like \(n^2\) , not \(n^4\) , because of the cancellation in the cross terms: that is where independence enters, and it is what makes the tail of \(S_n/n\) of order \(n^{-2}\) .

Example 8.2: Random signs and normal numbers

For independent signs \(\xi_k=\pm1\) with probability \(\tfrac12\) , \(\mathbb E\xi^4=\sigma^4=1\) and the formula gives \(\mathbb E S_n^4=n+3n(n-1)=3n^2-2n\) ; for \(n=4\) this is \(40\) , which enumeration of the \(16\) sign patterns confirms: \(S_4\) equals \(\pm4\) in \(2\) patterns, \(\pm2\) in \(8\) , and \(0\) in \(6\) , and \((2\cdot 256+8\cdot 16)/16=40\) . Hence \(\mathbb P(\lvert S_n/n\rvert>\varepsilon)\le 3/(n^2\varepsilon^4)\) and \(S_n/n\to 0\) almost surely.

Applied to the binary digits \(D_k\) of a uniform point of \((0,1)\) , which by Proposition 2.8 are independent fair bits, Theorem 8.2 gives

\[ \frac{D_1+\dots+D_n}{n}\longrightarrow\frac12\qquad\text{for Lebesgue-almost every point of }(0,1). \]

This is Borel’s normal number theorem: almost every real number has, in its binary expansion, asymptotically as many ones as zeros. It was the first strong law, and it is a statement of pure analysis that happens to have a probabilistic proof. No explicit example of such a number is produced by the proof.

8.3 First moments: truncation

With only \(\mathbb E\lvert X_1\rvert<\infty\) there is no variance to feed into Chebyshev’s inequality. The remedy is to manufacture one: cut each variable off at a level beyond which it rarely goes, apply second-moment methods to the bounded pieces, and show that the discarded parts do not matter. This device, truncation, turns out to need nothing beyond the mean.

Theorem 8.3 (weak law under a first moment). Let \(X_1,X_2,\dots\) be pairwise independent and identically distributed with \(\mathbb E\lvert X_1\rvert<\infty\) and mean \(\mu\) . Then \(S_n/n\to\mu\) in probability and in \(L^1\) .

Proof. Fix \(n\) and truncate at level \(n\) : let \(X_{n,k}=X_k\mathbf 1_{\{\lvert X_k\rvert\le n\}}\) for \(k\le n\) and \(S_n'=\sum_{k\le n}X_{n,k}\) . We compare \(S_n/n\) with \(S_n'/n\) and \(S_n'/n\) with its mean.

The truncation is rarely felt. By subadditivity, \(\mathbb P(S_n\ne S_n')\le n\,\mathbb P(\lvert X_1\rvert>n)\le\mathbb E[\lvert X_1\rvert;\lvert X_1\rvert>n]\) , which tends to \(0\) by dominated convergence.

The mean is nearly right. \(\mathbb E S_n'/n=\mathbb E[X_1;\lvert X_1\rvert\le n]\to\mu\) , again by dominated convergence.

The variance is small. The variables \(X_{n,k}\) , \(k\le n\) , are pairwise independent, so \(\operatorname{Var}(S_n'/n)\le\tfrac1n\mathbb E[X_1^2;\lvert X_1\rvert\le n]\) . For any \(K<n\) , bound \(X_1^2\) by \(K^2\) where \(\lvert X_1\rvert\le K\) and by \(n\lvert X_1\rvert\) where \(K<\lvert X_1\rvert\le n\) :

\[ \frac1n\,\mathbb E\bigl[X_1^2;\lvert X_1\rvert\le n\bigr]\le\frac{K^2}{n}+\mathbb E\bigl[\lvert X_1\rvert;\lvert X_1\rvert>K\bigr]. \]

Letting \(n\to\infty\) and then \(K\to\infty\) shows that the left side tends to \(0\) .

Combining, for \(\varepsilon>0\) and \(n\) so large that \(\lvert\mathbb E S_n'/n-\mu\rvert<\varepsilon/2\) ,

\[ \mathbb P\Bigl(\Bigl\lvert\frac{S_n}{n}-\mu\Bigr\rvert>\varepsilon\Bigr)\le\mathbb P(S_n\ne S_n')+\frac{4}{\varepsilon^2}\operatorname{Var}\Bigl(\frac{S_n'}{n}\Bigr)\longrightarrow 0 . \]

For convergence in \(L^1\) : the \(X_k\) are identically distributed and integrable, hence uniformly integrable, so the averages \(S_n/n\) are uniformly integrable by Proposition 5.5, and Theorem 5.6 upgrades convergence in probability to convergence in \(L^1\) .

The three steps, in the order “rarely felt, mean nearly right, variance small”, are the template for every truncation argument. The proof of the strong law has the same three steps, with each strengthened from “tends to zero” to “summable”.

8.4 The strong law

Theorem 8.4 (strong law of large numbers). Let \(X_1,X_2,\dots\) be pairwise independent and identically distributed with \(\mathbb E\lvert X_1\rvert<\infty\) and mean \(\mu\) . Then

\[ \frac{S_n}{n}\longrightarrow\mu\qquad\text{almost surely and in }L^1 . \]

The proof, due to Etemadi, needs two lemmas about the truncated variables. Both are stated for a nonnegative variable; the reduction to that case is the first step of the proof.

Lemma 8.5 (truncated tails and second moments). Let \(X\ge 0\) be integrable. Then

\[ \sum_{k\ge 1}\mathbb P(X>k)\le\mathbb E X\qquad\text{and}\qquad\sum_{k\ge 1}\frac{\mathbb E[X^2;X\le k]}{k^2}\le 2\,\mathbb E X . \]

Proof. The first inequality is contained in Proposition 1.12. For the second, exchange sum and expectation by Tonelli’s theorem:

\[ \sum_{k\ge 1}\frac{\mathbb E[X^2;X\le k]}{k^2}=\mathbb E\Bigl[X^2\sum_{k\ge 1,\ k\ge X}\frac1{k^2}\Bigr]. \]

We claim that \(\sum_{k\ge 1,\,k\ge x}k^{-2}\le 2/x\) for every \(x>0\) . If \(x\le 1\) the sum is \(\sum_{k\ge 1}k^{-2}=\pi^2/6<2\le 2/x\) . If \(x>1\) , let \(m\ge 2\) be the least integer with \(m\ge x\) ; comparing with an integral, \(\sum_{k\ge m}k^{-2}\le m^{-2}+\int_m^\infty t^{-2}\,dt=m^{-2}+m^{-1}\le 2/m\le 2/x\) . Hence the right side of the display is at most \(\mathbb E[X^2\cdot 2/X;X>0]=2\,\mathbb E X\) .

Proof of Theorem 8.4. Step 0: reduction to nonnegative variables. The positive parts \(X_k^{+}\) are pairwise independent, identically distributed, and integrable, and so are the negative parts \(X_k^{-}\) . If the theorem holds for nonnegative variables, then \(\tfrac1n\sum_{k\le n}X_k^{\pm}\to\mathbb E X_1^{\pm}\) almost surely, and subtracting gives the general case. So assume \(X_k\ge 0\) .

Step 1: truncation. Let \(Y_k=X_k\mathbf 1_{\{X_k\le k\}}\) and \(T_n=Y_1+\dots+Y_n\) . By Lemma 8.5 and identical distribution, \(\sum_k\mathbb P(X_k\ne Y_k)=\sum_k\mathbb P(X_1>k)\le\mathbb E X_1<\infty\) . By the first Borel–Cantelli lemma, almost surely \(X_k=Y_k\) for all but finitely many \(k\) , so \(S_n-T_n\) is eventually constant and \((S_n-T_n)/n\to 0\) . It suffices to prove \(T_n/n\to\mu\) almost surely.

Step 2: the means. \(\mathbb E Y_k=\mathbb E[X_1;X_1\le k]\uparrow\mu\) by monotone convergence, and the averages of a convergent sequence of numbers converge to the same limit, so \(\mathbb E T_n/n\to\mu\) .

Step 3: almost sure convergence along a geometric subsequence. Fix \(\alpha>1\) and let \(k_n=\lfloor\alpha^n\rfloor\) for \(n\ge 0\) . The \(Y_k\) are pairwise independent, so \(\operatorname{Var}(T_m)=\sum_{i\le m}\operatorname{Var}(Y_i)\) , and by Chebyshev’s inequality, for \(\varepsilon>0\) ,

\[ \sum_{n\ge 0}\mathbb P\bigl(\lvert T_{k_n}-\mathbb E T_{k_n}\rvert>\varepsilon k_n\bigr)\le\frac1{\varepsilon^2}\sum_{n\ge 0}\frac1{k_n^2}\sum_{i\le k_n}\operatorname{Var}(Y_i)=\frac1{\varepsilon^2}\sum_{i\ge 1}\operatorname{Var}(Y_i)\sum_{n:\,k_n\ge i}\frac1{k_n^2}. \]

To bound the inner sum, note that \(k_n\ge\alpha^n/2\) , because \(\lfloor y\rfloor\ge y/2\) for \(y\ge 1\) . If \(n_0\) is the least \(n\) with \(k_n\ge i\) , then \(\alpha^{n_0}\ge i\) , and

\[ \sum_{n:\,k_n\ge i}\frac1{k_n^2}\le 4\sum_{n\ge n_0}\alpha^{-2n}=\frac{4\,\alpha^{-2n_0}}{1-\alpha^{-2}}\le\frac{4}{(1-\alpha^{-2})\,i^2}. \]

Since \(\operatorname{Var}(Y_i)\le\mathbb E Y_i^2=\mathbb E[X_1^2;X_1\le i]\) , Lemma 8.5 bounds the whole expression by \(8\,\mathbb E X_1/(\varepsilon^2(1-\alpha^{-2}))<\infty\) . By the first Borel–Cantelli lemma, \((T_{k_n}-\mathbb E T_{k_n})/k_n\to 0\) almost surely, and with Step 2, \(T_{k_n}/k_n\to\mu\) almost surely.

Step 4: filling the gaps. Because \(Y_k\ge 0\) , the sums \(T_m\) are nondecreasing in \(m\) . Every \(m\ge 1\) lies between two consecutive terms of the subsequence, and for \(k_n\le m\le k_{n+1}\) ,

\[ \frac{k_n}{k_{n+1}}\cdot\frac{T_{k_n}}{k_n}\le\frac{T_m}{m}\le\frac{k_{n+1}}{k_n}\cdot\frac{T_{k_{n+1}}}{k_{n+1}}. \]

Since \(k_{n+1}/k_n\to\alpha\) , almost surely \(\mu/\alpha\le\liminf_mT_m/m\le\limsup_mT_m/m\le\alpha\mu\) . Applying this with \(\alpha=1+1/j\) for \(j=1,2,\dots\) and intersecting the countably many events of probability one gives \(T_m/m\to\mu\) almost surely.

Convergence in \(L^1\) was proved in Theorem 8.3.

It is worth recording where each hypothesis was used. Identical distribution entered through Lemma 8.5 and Step 2. Independence entered exactly once, as the additivity of variance in Step 3, and for that pairwise independence is enough, just as it was for the second Borel–Cantelli lemma in Theorem 3.5. Nonnegativity made \(T_m\) monotone, which is what allowed a sparse subsequence to control everything in between; without monotonicity one needs a maximal inequality, such as Kolmogorov’s inequality of Corollary 7.9, to control the fluctuations between consecutive terms of the subsequence, and that route requires full independence. The subsequence had to be geometric: the factor \(i^{-2}\) in Step 3 is what Lemma 8.5 can absorb, and a subsequence such as \(k_n=n^2\) would produce \(i^{-3/2}\) instead (Exercise 17).

The first moment is not only sufficient.

Theorem 8.6 (necessity of the first moment). Let \(X_1,X_2,\dots\) be independent and identically distributed with \(\mathbb E\lvert X_1\rvert=\infty\) . Then \(\limsup_n\lvert S_n\rvert/n=\infty\) almost surely. In particular, the averages \(S_n/n\) converge to a finite limit with probability \(0\) .

Proof. Fix \(K>0\) . Since \(\mathbb E\lvert X_1/K\rvert=\infty\) , Exercise 3.5, which is the second Borel–Cantelli lemma applied to the independent events \(\{\lvert X_n\rvert>Kn\}\) , gives \(\lvert X_n\rvert>Kn\) for infinitely many \(n\) , almost surely. Now \(X_n=S_n-S_{n-1}\) , so \(\lvert X_n\rvert>Kn\) forces \(\lvert S_n\rvert>Kn/2\) or \(\lvert S_{n-1}\rvert>Kn/2\) . Hence almost surely \(\lvert S_m\rvert/m>K/2\) for infinitely many \(m\) , and \(\limsup_m\lvert S_m\rvert/m\ge K/2\) . Since \(K\) was arbitrary, the upper limit is infinite.

For independent identically distributed variables there is therefore a clean dichotomy, in keeping with the zero-one law: if \(\mathbb E\lvert X_1\rvert<\infty\) the averages converge almost surely, to the mean; if not, they are almost surely unbounded. Here full independence is used, through the second Borel–Cantelli lemma.

Example 8.3: Cauchy averages do not settle

Let \(X_1,X_2,\dots\) be independent with the Cauchy density \(1/(\pi(1+x^2))\) . The family of Cauchy laws with scale \(a>0\) , with density \(a/(\pi(a^2+x^2))\) , is closed under convolution: the sum of independent Cauchy variables with scales \(a\) and \(b\) is Cauchy with scale \(a+b\) . Chapter 9 will make this a one-line computation with characteristic functions; it can also be verified directly from the convolution integral. Consequently \(S_n\) is Cauchy with scale \(n\) , and \(S_n/n\) is Cauchy with scale \(1\) : the average of \(n\) observations has the law of a single observation. In particular

\[ \mathbb P\Bigl(\Bigl\lvert\frac{S_n}{n}\Bigr\rvert>1\Bigr)=1-\frac2\pi\arctan 1=\frac12\qquad\text{for every }n . \]

No law of large numbers holds, in any mode: \(S_n/n\) does not converge in probability to a constant, and by Theorem 8.6, since \(\mathbb E\lvert X_1\rvert=\infty\) (Example 1.6), \(\limsup_n\lvert S_n\rvert/n=\infty\) almost surely. Even the necessary condition for a weak law with centering constants fails: \(n\,\mathbb P(\lvert X_1\rvert>n)=n(1-\tfrac2\pi\arctan n)\to\tfrac2\pi\) , not \(0\) .

The symmetry of the Cauchy density about \(0\) tempts one to call \(0\) its mean and to expect the averages to approach it. Figure 8.1 shows what happens instead. A single observation of size comparable to \(n\) arrives often enough to move the average by an amount of order \(1\) , at every scale.

Two panels plotting the running average of a simulated sample against the sample size from 1 to 100000 on a logarithmic axis. Left: for exponential observations the average wanders at first and then hugs the dashed line at height 1. Right: for Cauchy observations the average makes sudden jumps of several units at sample sizes near 3, 90, and 2500 and never approaches the dashed line at 0.
Figure 8.1. One simulated path of running averages in each case. With a finite mean the path is eventually trapped in every band around the mean, which is the strong law; without one, a single huge observation keeps resetting the average, as Theorem 8.6 predicts.

Example 8.4: Integration by sampling

To compute \(I=\int_0^1h(x)\,dx\) for an integrable \(h\) , let \(U_1,U_2,\dots\) be independent and uniform on \((0,1)\) . The variables \(h(U_k)\) are independent and identically distributed with mean \(I\) , so by Theorem 8.4

\[ \frac{h(U_1)+\dots+h(U_n)}{n}\longrightarrow I\qquad\text{almost surely.} \]

Take \(h(x)=4/(1+x^2)\) , for which \(I=\pi\) . The strong law guarantees convergence and says nothing about speed. For that, compute the variance: \(\int_0^1(1+x^2)^{-2}\,dx=\tfrac\pi8+\tfrac14\) , so

\[ \operatorname{Var}\,h(U)=16\Bigl(\frac\pi8+\frac14\Bigr)-\pi^2=2\pi+4-\pi^2\approx 0.4136 , \]

and Theorem 8.1 gives \(\mathbb P(\lvert\text{average}-\pi\rvert\ge 0.01)\le 0.4136/(n\cdot 10^{-4})\) , which is at most \(0.05\) once \(n\ge 82{,}717\) . The error of the method is of order \(\sqrt{0.4136/n}\) , about \(0.00064\) for \(n=10^6\) : each further decimal digit costs a hundredfold increase in the sample. The strong law explains why the method works; the variance, and in Chapter 10 the central limit theorem, explain how well.

Example 8.5: The truncation lemma on a law without variance

Let \(X\ge 1\) have the tail \(\mathbb P(X>x)=x^{-3/2}\) for \(x\ge 1\) , so that \(X\) has density \(\tfrac32x^{-5/2}\) , mean \(\mathbb E X=3\) , and infinite variance. Both inequalities of Lemma 8.5 can be evaluated exactly, using the values \(\zeta(\tfrac32)=\sum_kk^{-3/2}\approx 2.6124\) and \(\zeta(2)=\pi^2/6\approx 1.6449\) of the zeta function.

For the first, \(\sum_{k\ge 1}\mathbb P(X>k)=\zeta(\tfrac32)\approx 2.612\le 3\) . For the second, \(\mathbb E[X^2;X\le k]=\int_1^kx^2\cdot\tfrac32x^{-5/2}\,dx=3(\sqrt k-1)\) , which grows without bound, as it must for a law with infinite variance; nevertheless

\[ \sum_{k\ge 1}\frac{3(\sqrt k-1)}{k^2}=3\bigl(\zeta(\tfrac32)-\zeta(2)\bigr)\approx 2.902\le 6=2\,\mathbb E X . \]

The truncated second moments grow like \(\sqrt k\) , and division by \(k^2\) makes them summable with room to spare. This is the quantitative content of the strong law’s proof: the variance of the \(k\) -th truncated variable may tend to infinity, so long as it does so more slowly than \(k\) on average, and a finite first moment guarantees exactly that.

Common mistake. Quoting “the law of large numbers” without its mode or its hypothesis. Two errors follow. The first is to apply it where there is no mean, as to Cauchy observations, on the strength of symmetry or of a finite median; the hypothesis is \(\mathbb E\lvert X_1\rvert<\infty\) , and Theorem 8.6 shows that for independent identically distributed variables nothing less will do. The second is to draw a pathwise conclusion from the weak law: “\(\mathbb P(\lvert S_n/n-\mu\rvert>\varepsilon)\to 0\) , so eventually the running average stays within \(\varepsilon\) of \(\mu\) .” That conclusion is the strong law. It is true under the hypotheses of Theorem 8.4, but it does not follow from convergence in probability, as Chapter 4 showed. When citing a law of large numbers, write the arrow with its mode and state the moment assumption beside it.

8.5 What remains to be said about averages

The strong law is a statement about the limit and is silent about the approach to it. Three kinds of refinement exist, and the rest of the book concerns the first.

Fluctuations. If the variance \(\sigma^2\) is finite, Theorem 8.1 shows that \(S_n/n-\mu\) is of order \(n^{-1/2}\) in mean square. Multiplying by \(\sqrt n\) gives a quantity with constant variance \(\sigma^2\) , and one may ask whether its law settles down. It does, and the limit is the same for all laws with the same variance; that is the central limit theorem, proved in Chapter 10. The convergence there is only in distribution, the weakest mode of Chapter 4, and it requires a separate theory, which is the subject of Chapter 9.

Almost sure rates. The strong law says \(S_n-n\mu=o(n)\) almost surely. The sharp almost sure statement, for finite variance, is the law of the iterated logarithm, \(\limsup_n(S_n-n\mu)/\sqrt{2\sigma^2n\log\log n}=1\) , which we mention without proof.

Dependence. Independence was used through the additivity of variances and through the Borel–Cantelli lemmas. Sequences in which the variables are dependent can obey laws of large numbers for other reasons, as martingale differences or as stationary sequences do, and they can also fail to: if \(X_n=X_1\) for all \(n\) the averages equal \(X_1\) and converge to a random limit. Such results belong to the study of stochastic processes.

The next chapter develops convergence of laws. Its opening question is one of definition: the constants \(1/n\) surely converge in distribution to \(0\) , and yet their distribution functions do not converge to the distribution function of \(0\) at the origin. What is the right notion?

Exercises

  1. Let \(X_1,X_2,\dots\) be pairwise uncorrelated with mean \(\mu\) and variances \(\sigma_k^2\) , not assumed bounded. Show that \(S_n/n\to\mu\) in \(L^2\) if and only if \(n^{-2}\sum_{k\le n}\sigma_k^2\to 0\) , and give an example with \(\sigma_k^2\to\infty\) in which this holds.

    Hint

    \(\mathbb E(S_n/n-\mu)^2=n^{-2}\sum_{k\le n}\sigma_k^2\) exactly, by Corollary 2.7. Example: \(\sigma_k^2=\sqrt k\) , for which \(n^{-2}\sum_{k\le n}\sqrt k\le n^{-1/2}\to 0\) .

  2. Let \(X_1,\dots,X_n\) be independent and identically distributed with mean \(0\) , variance \(\sigma^2\) , and fourth moment \(m_4\) . Prove that \(\mathbb E S_n^4=nm_4+3n(n-1)\sigma^4\) , and verify the formula by enumeration for \(n=3\) when \(X\) takes the values \(-1,0,1\) with probabilities \(\tfrac14,\tfrac12,\tfrac14\) .

    Solution

    Expand \(S_n^4=\sum_{i,j,k,l}X_iX_jX_kX_l\) over all \(n^4\) quadruples and take expectations. If some index value occurs exactly once in the quadruple, the corresponding factor is independent of the product of the others and has mean \(0\) , so the term vanishes by Theorem 2.6. The remaining quadruples are of two kinds. Those with all four indices equal number \(n\) , and each contributes \(m_4\) . Those in which two distinct values each occur twice: the unordered pair of values can be chosen in \(\binom n2\) ways, and the two positions occupied by the smaller value in \(\binom42=6\) ways, giving \(6\binom n2=3n(n-1)\) quadruples, each contributing \(\mathbb E X_i^2\,\mathbb E X_j^2=\sigma^4\) . Hence

    \[ \mathbb E S_n^4=n\,m_4+3n(n-1)\,\sigma^4 . \]

    For the three-point law, \(\sigma^2=\mathbb E X^2=\tfrac12\) and \(m_4=\mathbb E X^4=\tfrac12\) , so the formula gives \(\mathbb E S_3^4=\tfrac32+18\cdot\tfrac14=6\) . Enumerating the \(27\) value patterns with their probabilities and averaging \((X_1+X_2+X_3)^4\) also gives \(6\) . The point of the formula is the order of magnitude: \(n^2\) , not \(n^4\) , because only paired indices survive.

  3. Let \(X_1,X_2,\dots\) be pairwise uncorrelated with mean \(0\) and \(\lvert X_k\rvert\le K\) . Prove that \(S_n/n\to 0\) almost surely by first treating the subsequence \(n=k^2\) with Theorem 8.1 and then bounding \(\lvert S_n-S_{k^2}\rvert\) for \(k^2\le n<(k+1)^2\) .

    Solution

    Along squares. The variances are at most \(K^2\) , so by Theorem 8.1, \(\mathbb P(\lvert S_{k^2}\rvert\ge\varepsilon k^2)\le K^2/(k^2\varepsilon^2)\) . These bounds are summable in \(k\) , and the first Borel–Cantelli lemma gives \(\lvert S_{k^2}\rvert/k^2<\varepsilon\) for all large \(k\) , almost surely. With \(\varepsilon=1/j\) and an intersection over \(j\) , \(S_{k^2}/k^2\to 0\) almost surely.

    Between squares. For \(k^2\le n<(k+1)^2\) the sum \(S_n-S_{k^2}\) has \(n-k^2\le 2k\) terms, each bounded by \(K\) , so \(\lvert S_n-S_{k^2}\rvert\le 2kK\) at every sample point. Therefore

    \[ \frac{\lvert S_n\rvert}{n}\le\frac{\lvert S_{k^2}\rvert}{k^2}+\frac{2kK}{k^2}\longrightarrow 0\qquad\text{almost surely,} \]

    since \(k\to\infty\) as \(n\to\infty\) . The subsequence \(k^2\) is dense enough that bounded terms cannot move the average between consecutive squares, and sparse enough that Chebyshev’s bound becomes summable. For unbounded variables the gaps are the difficulty, which is why the proof of Theorem 8.4 needs truncation and monotonicity.

  4. Let \(f\) be continuous on \([0,1]\) and let \(B_nf(p)=\sum_{k=0}^{n}f(k/n)\binom nkp^k(1-p)^{n-k}\) . Use Theorem 8.1 to prove that \(B_nf\to f\) uniformly on \([0,1]\) , and bound \(\sup_p\lvert B_nf(p)-f(p)\rvert\) for \(f(x)=\lvert x-\tfrac12\rvert\) and \(n=100\) .

    Solution

    Fix \(p\in[0,1]\) and let \(S_n\) be the number of successes in \(n\) independent trials with success probability \(p\) , so that \(B_nf(p)=\mathbb E f(S_n/n)\) . Let \(M=\sup\lvert f\rvert\) and \(\varepsilon>0\) . By uniform continuity of \(f\) on \([0,1]\) there is \(\delta>0\) with \(\lvert f(x)-f(y)\rvert\le\varepsilon\) whenever \(\lvert x-y\rvert\le\delta\) . Splitting according to whether \(\lvert S_n/n-p\rvert\le\delta\) ,

    \[ \lvert B_nf(p)-f(p)\rvert\le\mathbb E\lvert f(S_n/n)-f(p)\rvert\le\varepsilon+2M\,\mathbb P\bigl(\lvert S_n/n-p\rvert>\delta\bigr)\le\varepsilon+\frac{2M\,p(1-p)}{n\delta^2}\le\varepsilon+\frac{M}{2n\delta^2}, \]

    by Theorem 8.1 with variance \(p(1-p)\le\tfrac14\) . The bound does not depend on \(p\) , so \(\limsup_n\sup_p\lvert B_nf(p)-f(p)\rvert\le\varepsilon\) for every \(\varepsilon\) : the convergence is uniform. This is Weierstrass’s approximation theorem, with explicit polynomials.

    For \(f(x)=\lvert x-\tfrac12\rvert\) we have \(M=\tfrac12\) and may take \(\delta=\varepsilon\) , giving the bound \(\varepsilon+1/(4n\varepsilon^2)\) ; with \(n=100\) the best choice \(\varepsilon=200^{-1/3}\approx 0.171\) gives \(0.256\) . The actual maximal error, attained at \(p=\tfrac12\) , is \(\mathbb E\lvert S_{100}/100-\tfrac12\rvert\approx 0.0398\) . As in Example 8.1, the Chebyshev bound proves convergence and overstates the error.

  5. Prove that almost every number in \((0,1)\) is normal in base \(10\) in the strong sense that every finite block of digits occurs with its expected frequency, by applying the strong law to suitable indicator variables that are independent along arithmetic progressions.

    Hint

    The decimal digits of a uniform point are independent and uniform on \(\{0,\dots,9\}\) : the proof of Proposition 2.8 applies word for word with intervals of length \(10^{-n}\) in place of \(2^{-n}\) , as Example 2.5 notes for base three. Fix a block \(w\) of length \(r\) . The indicators that the digits in positions \(jr+s+1,\dots,jr+s+r\) equal \(w\) , for \(j=0,1,2,\dots\) and a fixed offset \(s\in\{0,\dots,r-1\}\) , are functions of disjoint groups of digits, hence independent by Corollary 2.4, with mean \(10^{-r}\) . Apply Theorem 8.2 for each offset, average over the \(r\) offsets, and intersect over the countably many blocks.

  6. A student writes: “The Cauchy density is symmetric about \(0\) , so by the law of large numbers the average of \(n\) independent Cauchy observations is close to \(0\) for large \(n\) .” Compute \(\mathbb P(\lvert S_n/n\rvert>1)\) , identify the hypothesis that fails, and state what Theorem 8.6 says about these averages.

    Solution

    By Example 8.3 the average \(S_n/n\) has the same Cauchy law as a single observation, for every \(n\) . Its distribution function is \(\tfrac12+\tfrac1\pi\arctan x\) , so

    \[ \mathbb P\Bigl(\Bigl\lvert\frac{S_n}{n}\Bigr\rvert>1\Bigr)=1-\frac2\pi\arctan 1=1-\frac2\pi\cdot\frac\pi4=\frac12\qquad\text{for all }n . \]

    The average is not close to \(0\) for large \(n\) ; it is as dispersed as one observation. The hypothesis that fails is integrability: \(\mathbb E\lvert X_1\rvert=\infty\) (Example 1.6), so the Cauchy law has no mean. Symmetry about \(0\) makes \(0\) the median and the centre of symmetry, and neither of these is what a law of large numbers is about. Theorem 8.6 gives the precise statement: since the observations are independent and identically distributed with \(\mathbb E\lvert X_1\rvert=\infty\) , almost surely \(\limsup_n\lvert S_n\rvert/n=\infty\) . The averages are not merely non-convergent but unbounded along almost every sequence of observations, as in the right panel of Figure 8.1.

  7. Let \(X_1,X_2,\dots\) be independent with \(\mathbb P(X_k=2^j)=2^{-j}\) for \(j\ge 1\) . Show that \(\mathbb E X_1=\infty\) and that \(S_n/(n\log_2n)\to 1\) in probability, by truncating at the level \(n\log_2n\) .

    Hint

    \(\mathbb E X_1=\sum_j2^j2^{-j}=\infty\) . With \(b_n=n\log_2n\) and \(X_k'=X_k\mathbf 1_{\{X_k\le b_n\}}\) : \(n\,\mathbb P(X_1>b_n)\le 2n/b_n\to 0\) ; \(\mathbb E X_1'\) is the number of \(j\) with \(2^j\le b_n\) , which is \(\log_2n+\log_2\log_2n+O(1)\) ; and \(\operatorname{Var}(\sum_{k\le n}X_k')\le n\sum_{2^j\le b_n}2^j\le 2nb_n\) , which is small compared with \(b_n^2\) . Conclude as in Theorem 8.3.

  8. Let \(X_1,X_2,\dots\) be independent and identically distributed with \(\mathbb E X_1^{+}=\infty\) and \(\mathbb E X_1^{-}<\infty\) . Prove that \(S_n/n\to\infty\) almost surely.

    Solution

    For \(M>0\) let \(X_k^{M}=X_k\wedge M\) . These variables are independent and identically distributed, and integrable: the positive part is at most \(M\) and the negative part equals \(X_k^{-}\) , which is integrable by hypothesis. By Theorem 8.4,

    \[ \liminf_{n\to\infty}\frac{S_n}{n}\ \ge\ \lim_{n\to\infty}\frac1n\sum_{k\le n}X_k^{M}=\mathbb E[X_1\wedge M]\qquad\text{almost surely,} \]

    because \(X_k\ge X_k^{M}\) . This holds for each \(M\) outside a null event, hence for all integers \(M\) simultaneously outside a null event. As \(M\uparrow\infty\) , the variables \(X_1\wedge M+X_1^{-}\) are nonnegative and increase to \(X_1^{+}\) , so by monotone convergence \(\mathbb E[X_1\wedge M]\uparrow\mathbb E X_1^{+}-\mathbb E X_1^{-}=\infty\) . Therefore \(\liminf_nS_n/n=\infty\) almost surely. Truncation from above produced integrable variables whose averages bound the original averages from below; the strong law for the truncated variables did the rest.

  9. Let \(X_1,X_2,\dots\) be independent and identically distributed with mean \(\mu>0\) . Show that \(S_n\to\infty\) almost surely and that \(\inf_nS_n>-\infty\) almost surely.

    Answer

    By Theorem 8.4, \(S_n/n\to\mu>0\) almost surely, so almost surely \(S_n>n\mu/2\) for all large \(n\) ; hence \(S_n\to\infty\) , and the infimum over \(n\) is the minimum of finitely many finite values.

  10. Let \(\xi_1,\xi_2,\dots\) be independent, identically distributed, and positive with mean \(\mu\in(0,\infty)\) , let \(T_n=\xi_1+\dots+\xi_n\) , and let \(N_t=\max\{n:T_n\le t\}\) . Prove that \(N_t/t\to 1/\mu\) almost surely as \(t\to\infty\) .

    Solution

    Since each \(T_n\) is finite, \(N_t\to\infty\) as \(t\to\infty\) at every sample point, and since \(T_n\to\infty\) almost surely (by the strong law, \(T_n/n\to\mu>0\) ), \(N_t\) is finite for every \(t\) . By the definition of \(N_t\) ,

    \[ T_{N_t}\le t<T_{N_t+1},\qquad\text{so}\qquad\frac{T_{N_t}}{N_t}\le\frac{t}{N_t}<\frac{T_{N_t+1}}{N_t+1}\cdot\frac{N_t+1}{N_t} \]

    once \(N_t\ge 1\) . On the event of probability one where \(T_n/n\to\mu\) , the same limit holds along the integers \(N_t\) and \(N_t+1\) , because they tend to infinity; and \((N_t+1)/N_t\to 1\) . Both outer terms tend to \(\mu\) , so \(t/N_t\to\mu\) , that is, \(N_t/t\to 1/\mu\) almost surely. The argument is purely pathwise once the strong law has been applied: it converts a statement about the time of the \(n\) -th event into a statement about the number of events by time \(t\) . This is the elementary renewal law used in the study of Stochastic Processes.

  11. Let \(X_1,X_2,\dots\) be independent with common distribution function \(F\) and let \(F_n(x)=\tfrac1n\sum_{k\le n}\mathbf 1_{\{X_k\le x\}}\) . Show that \(F_n(x)\to F(x)\) almost surely for each \(x\) , and that when \(F\) is continuous the convergence is almost surely uniform in \(x\) .

    Hint

    For fixed \(x\) the indicators \(\mathbf 1_{\{X_k\le x\}}\) are independent and identically distributed with mean \(F(x)\) . For uniformity with \(F\) continuous, choose points \(x_1<\dots<x_m\) with \(F(x_{i+1})-F(x_i)\le\varepsilon\) and use monotonicity of \(F_n\) and \(F\) between them.

  12. For \(X\) exponential with rate \(1\) evaluate the two sums in Proposition 1.12 that bracket \(\mathbb E X\) , and evaluate \(\sum_{k\ge 1}\mathbb P(X>k)\) as it appears in Lemma 8.5.

    Answer

    \(\sum_{n\ge 1}\mathbb P(X\ge n)=\sum_{n\ge 1}e^{-n}=1/(e-1)\approx 0.582\) and \(\sum_{n\ge 0}\mathbb P(X>n)=e/(e-1)\approx 1.582\) , which bracket \(\mathbb E X=1\) . In Lemma 8.5, \(\sum_{k\ge 1}\mathbb P(X>k)=1/(e-1)\) .

  13. Let \(X_1,X_2,\dots\) be independent and identically distributed and not almost surely \(0\) . Show that the series \(\sum_nX_n\) converges with probability \(0\) .

    Answer

    If the series converged, its terms would tend to \(0\) . But for some \(\varepsilon>0\) , \(\mathbb P(\lvert X_1\rvert>\varepsilon)>0\) , and the independent events \(\{\lvert X_n\rvert>\varepsilon\}\) have this same probability, so by the second Borel–Cantelli lemma they occur infinitely often almost surely.

  14. Let \(X\) be symmetric with \(\mathbb P(\lvert X\rvert>x)=e/(x\log x)\) for \(x\ge e\) . Show that \(\mathbb E\lvert X\rvert=\infty\) while \(n\,\mathbb P(\lvert X\rvert>n)\to 0\) , and adapt the proof of Theorem 8.3 to show that \(S_n/n\to 0\) in probability for independent copies; conclude with Theorem 8.6 that the weak law can hold where the strong law fails.

    Hint

    \(\mathbb E\lvert X\rvert=\infty\) by Exercise 1.16, and \(n\,\mathbb P(\lvert X\rvert>n)=e/\log n\to 0\) . In the proof of Theorem 8.3 the first step needs only \(n\,\mathbb P(\lvert X_1\rvert>n)\to 0\) ; the truncated means vanish by symmetry; and \(\tfrac1n\mathbb E[X_1^2;\lvert X_1\rvert\le n]\le\tfrac1n\int_0^n2y\,\mathbb P(\lvert X_1\rvert>y)\,dy\to 0\) as an average of a function tending to \(0\) .

  15. Prove the \(L^1\) assertion of Theorem 8.4 in full: explain why \((X_k)\) is uniformly integrable, why the averages are, and which theorem then gives \(\mathbb E\lvert S_n/n-\mu\rvert\to 0\) .

    Solution

    The variables. The \(X_k\) all have the law of \(X_1\) , so for every \(k\) , \(\mathbb E[\lvert X_k\rvert;\lvert X_k\rvert>M]=\mathbb E[\lvert X_1\rvert;\lvert X_1\rvert>M]\) , which tends to \(0\) as \(M\to\infty\) because \(X_1\) is integrable. The supremum over \(k\) of these equal quantities tends to \(0\) : the sequence \((X_k)\) is uniformly integrable. No independence is involved.

    The averages. By Proposition 5.5 the averages \(\bar X_n=S_n/n\) of a uniformly integrable sequence are uniformly integrable; the proof used the two faces of Proposition 5.2 and the bound \(\mathbb E[\lvert\bar X_n\rvert;A]\le\sup_k\mathbb E[\lvert X_k\rvert;A]\) .

    The limit. By the almost sure part of Theorem 8.4, \(\bar X_n\to\mu\) almost surely, hence in probability. The Vitali convergence theorem (Theorem 5.6, from (i) to (ii)) now gives \(\mathbb E\lvert\bar X_n-\mu\rvert\to 0\) . Domination would not have been available: \(\sup_n\lvert\bar X_n\rvert\) need not be integrable when \(X_1\) is merely integrable, so the Dominated Convergence Theorem cannot replace Vitali’s theorem here.

  16. An event has unknown probability \(p\) . Use Theorem 8.1 to find a number of independent trials that guarantees, for every \(p\) , that the observed frequency is within \(0.01\) of \(p\) with probability at least \(0.95\) .

    Answer

    The variance of an indicator is \(p(1-p)\le\tfrac14\) , so Theorem 8.1 gives the bound \(1/(4n\cdot 10^{-4})\) , which is at most \(0.05\) when \(n\ge 50{,}000\) .

  17. In the proof of Theorem 8.4, say exactly where identical distribution, pairwise independence, and nonnegativity are used. Then repeat Step 3 with the subsequence \(k_n=n^2\) in place of \(\lfloor\alpha^n\rfloor\) and show which series would have to converge, and why Lemma 8.5 does not provide it.

    Solution

    Identical distribution is used three times: in Step 1, to write \(\sum_k\mathbb P(X_k>k)=\sum_k\mathbb P(X_1>k)\) and apply Lemma 8.5; in Step 2, to get \(\mathbb E Y_k=\mathbb E[X_1;X_1\le k]\uparrow\mu\) ; and in Step 3, to write \(\mathbb E Y_i^2=\mathbb E[X_1^2;X_1\le i]\) and apply Lemma 8.5 again. Pairwise independence is used once, in Step 3, for \(\operatorname{Var}(T_m)=\sum_{i\le m}\operatorname{Var}(Y_i)\) ; the \(Y_i\) are pairwise independent as functions of the \(X_i\) . Nonnegativity is used in Step 2 for monotone convergence and, decisively, in Step 4, where it makes \(T_m\) nondecreasing so that \(T_m/m\) can be sandwiched.

    Now take \(k_n=n^2\) . The inner sum in Step 3 becomes

    \[ \sum_{n:\,n^2\ge i}\frac1{n^4}\ \approx\ \int_{\sqrt i}^\infty\frac{dt}{t^4}=\frac{1}{3\,i^{3/2}}, \]

    so the argument would require \(\sum_i\operatorname{Var}(Y_i)\,i^{-3/2}<\infty\) . Lemma 8.5 controls only \(\sum_i\mathbb E Y_i^2\,i^{-2}\) , and the stronger series can diverge: for the law of Example 8.5, \(\mathbb E Y_i^2=3(\sqrt i-1)\) and \(\sum_i\sqrt i\cdot i^{-3/2}=\sum_i1/i=\infty\) . A polynomial subsequence has too many terms below each index \(i\) ; a geometric one has only about \(\log i\) of them, with reciprocal squares summing to order \(i^{-2}\) . The price of a geometric subsequence is large gaps, and the limits \(\mu/\alpha\) and \(\alpha\mu\) in Step 4, with \(\alpha\downarrow 1\) , are how the proof pays it.

  18. Kronecker’s lemma states that if \(\sum_ka_k/k\) converges then \(\tfrac1n\sum_{k\le n}a_k\to 0\) . Use it with Theorem 7.15 to prove that independent mean-zero variables with \(\sum_k\operatorname{Var}(X_k)/k^2<\infty\) satisfy \(S_n/n\to 0\) almost surely.

    Hint

    The partial sums \(\sum_{k\le n}X_k/k\) form a martingale with \(\mathbb E(\sum_{k\le n}X_k/k)^2=\sum_{k\le n}\operatorname{Var}(X_k)/k^2\) , bounded by hypothesis; Theorem 7.15 gives almost sure convergence of \(\sum_kX_k/k\) , and Kronecker’s lemma applied pathwise with \(a_k=X_k(\omega)\) gives \(S_n/n\to 0\) .

Search knowledge base

Explore

On this page

Note