Mathematics · Book 5 · Bachelor Year 3

University Mathematics — Year 3

University Mathematics — Year 3 · Bachelor Year 3

22Probability: Foundations and the Law of Large Numbers

Year 2 built probability on countable spaces; measure theory now removes every restriction. A probability space is a measure space of total mass 11, random variables are measurable maps, expectation is the Lebesgue integral — and at once the whole analytic arsenal (Chapters 9, 10 and 11) applies to chance. This chapter installs the dictionary, constructs infinite sequences of independent random variables (on [0,1]\intcc01, from binary digits: randomness is hiding inside Lebesgue measure), proves the Borel–Cantelli lemmas and Kolmogorov’s zero–one law, sorts out the modes of convergence, and proves the law of large numbers — the theorem that makes frequencies converge to probabilities and statistics possible. The weekend problem gives Etemadi’s proof of the strong law in its definitive L1L^1 form.

22.1 The dictionary

Definition 22.1

A probability space is a measure space (Ω,A,P)(\Omega, \mathcal A, \P) with P(Ω)=1\P(\Omega) = 1; elements of A\mathcal A are events, and a property holds almost surely (a.s.) if its event has probability 11. A random variable is a measurable map X ⁣:ΩRX \colon \Omega \to \R (or Rd\R^d: a random vector); its law is the pushforward probability measure PX=XP\P_X = X_*\P on R\R (Exercise 11.9), determined by the distribution function FX(t)=P(Xt)F_X(t) = \P(X \leq t) (Exercise 9.3). XX has density ff if PX=f ⁣dλ\P_X = f\,\dd\lambda; it is discrete if PX\P_X is a countable combination of Dirac masses. The expectation is

E[X]=ΩX ⁣dP(X0 or XL1(P)),\E[X] = \int_\Omega X\,\dd\P \qquad (X \geq 0 \text{ or } X \in L^1(\P)),

and the transfer theorem (Exercise 11.9) computes it in the law: E[g(X)]=Rg ⁣dPX\E[g(X)] = \int_\R g\,\dd\P_X=g(xk)pk= \sum g(x_k)p_k in the discrete case, =g(x)f(x) ⁣dx= \int g(x)f(x)\dd x in the density case: Year 2’s formulas, now theorems of one theory. The variance is V(X)=E[(XEX)2]=E[X2](EX)2\V(X) = \E[(X - \E X)^2] = \E[X^2] - (\E X)^2 for XL2X \in L^2.

Example 22.2

The standard laws and their transforms of note: Bernoulli B(p)\mathcal B(p), binomial B(n,p)\mathcal B(n, p), geometric, Poisson P(λ)\mathcal P(\lambda) (discrete: Year 2’s tables remain valid); uniform on [0,1]\intcc01 (Lebesgue measure itself); exponential E(λ)\mathcal E(\lambda) (density λeλx1x>0\lambda\eu^{-\lambda x}\mathbf 1_{x>0}); the Gaussian N(m,σ2)\mathcal N(m, \sigma^2) with density 1σ2πexp((xm)22σ2)\frac1{\sigma\sqrt{2\pi}}\exp\bigl(-\frac{(x - m)^2}{2\sigma^2}\bigr) — a probability density by Problem 10.1, with mean mm and variance σ2\sigma^2 (Gaussian moments, Exercise 11.10).

Proposition 22.3 (Markov and Chebyshev)

For X0X \geq 0 and a>0a > 0: P(Xa)EXa\P(X \geq a) \leq \frac{\E X}{a}; for XL2X \in L^2: P(XEXa)V(X)a2\P\bigl(\abs{X - \E X} \geq a\bigr) \leq \frac{\V(X)}{a^2}.

Proof. Exercise 10.5(a); Chebyshev is Markov applied to (XEX)2(X - \E X)^2.

22.2 Independence

Definition 22.4

Sub-σ\sigma-algebras A1,,AnA\mathcal A_1, \dots, \mathcal A_n \subseteq \mathcal A are independent if P(A1An)=P(Ai)\P(A_1\cap\dots\cap A_n) = \prod\P(A_i) for all AiAiA_i \in \mathcal A_i; events are independent if the σ\sigma-algebras {,Ai,Aic,Ω}\{\varnothing, A_i, A_i^c, \Omega\} are; random variables X1,,XnX_1, \dots, X_n if the σ\sigma-algebras σ(Xi)=Xi1(B(R))\sigma(X_i) = X_i^{-1}(\mathcal B(\R)) are. An infinite family is independent if every finite subfamily is.

Theorem 22.5

X1,,XnX_1, \dots, X_n are independent iff the law of the vector (X1,,Xn)(X_1, \dots, X_n) is the product measure PX1PXn\P_{X_1}\otimes\cdots\otimes\P_{X_n}. In that case, for gi0g_i \geq 0 (or such that the products are integrable):

E[igi(Xi)]=iE[gi(Xi)],\E\Bigl[\prod_ig_i(X_i)\Bigr] = \prod_i\E[g_i(X_i)],

in particular E[XY]=EXEY\E[XY] = \E X\,\E Y and V(X1++Xn)=V(Xi)\V(X_1 + \dots + X_n) = \sum\V(X_i) for independent L2L^2 variables.

Proof. If the XiX_i are independent, the two probability measures P(X1,,Xn)\P_{(X_1,\dots,X_n)} and PXi\bigotimes\P_{X_i} agree on all products B1××BnB_1\times\dots\times B_n of Borel sets — a π\pi-system generating B(Rn)\mathcal B(\R^n) (Proposition 11.2(b)) — hence everywhere (Theorem 9.7). Conversely, a product law factorizes all events iXi1(Bi)\bigcap_iX_i^{-1}(B_i): independence. The expectation formula is then Tonelli/Fubini (Theorem 11.5) through the transfer theorem; E[XY]=EXEY\E[XY] = \E X\E Y is the case gi=idg_i = \mathrm{id}, and expanding the square gives the additivity of variances (cross terms E[(XiEXi)(XjEXj)]=0\E[(X_i - \E X_i)(X_j - \E X_j)] = 0).

Theorem 22.6 (Existence of independent sequences)

On ([0,1],L,λ)\bigl(\intcc01, \mathcal L, \lambda\bigr) there exists a sequence (Un)n1(U_n)_{n\geq1} of independent random variables, each uniform on [0,1]\intcc01. Consequently, for any prescribed laws (μn)(\mu_n) on R\R there exist independent (Xn)(X_n) with PXn=μn\P_{X_n} = \mu_n.

Proof. Digits. For ω[0,1]\omega \in \intcc01, let (bk(ω))(b_k(\omega)) be its binary digits (ω=bk2k\omega = \sum b_k2^{-k}; choose the expansion not ending in all 11’s — ambiguity concerns only a countable, hence null, set). Each bkb_k is a random variable ({bk=1}\{b_k = 1\} is a finite union of dyadic intervals) and the vector (b1,,bm)(b_1, \dots, b_m) takes each value in {0,1}m\{0,1\}^m on a dyadic interval of length 2m2^{-m}: the bkb_k are independent Bernoulli(12)(\frac12).

Regrouping. Split N\N^* into infinitely many disjoint infinite sets (In)(I_n) (e.g. by prime powers, or diagonals); let (kjn)j(k^n_j)_j enumerate InI_n and set

Un=j1bkjn2j.U_n = \sum_{j\geq1} b_{k^n_j}\,2^{-j} .

Each UnU_n is uniform: its binary digits are independent fair bits, so P(Un[l2m,(l+1)2m))=2m\P(U_n \in [l2^{-m}, (l+1)2^{-m})) = 2^{-m} for every dyadic interval, and dyadic intervals determine the law (Theorem 9.7). The UnU_n are independent: they are functions of disjoint blocks of the independent family (bk)(b_k) — formally, events {UnDn}\{U_n \in D_n\} for dyadic DnD_n depend on finitely many digits from disjoint sets, and factorize; the π\pi-system argument upgrades to all Borel sets.

Arbitrary laws. Let Gn(u)=inf{t:Fμn(t)u}G_n(u) = \inf\{t : F_{\mu_n}(t) \geq u\} (the quantile function of the distribution function FμnF_{\mu_n}); the key equivalence Gn(u)t    uFμn(t)G_n(u) \leq t \iff u \leq F_{\mu_n}(t) (right-continuity of FF, monotonicity) shows Xn=Gn(Un)X_n = G_n(U_n) is measurable with P(Xnt)=P(UnFμn(t))=Fμn(t)\P(X_n \leq t) = \P(U_n \leq F_{\mu_n}(t)) = F_{\mu_n}(t): law μn\mu_n; independence is inherited (functions of independent variables, Exercise 22.3).

Example 22.7 (The birthday problem, honestly)

Among nn people with independent, uniform birthdays over N=365N = 365 days, the probability that all birthdays differ is

pn=k=1n1(1kN),p_n = \prod_{k=1}^{n-1}\Bigl(1 - \frac kN\Bigr),

by iterated conditioning (or directly: the favorable N(N1)(Nn+1)N(N-1)\cdots(N - n + 1) over the total NnN^n, a counting argument the product formula of independence makes rigorous). Taking logarithms and using ln(1x)=x+O(x2)-\ln(1 - x) = x + O(x^2):

lnpn=n(n1)2N+O(n3N2),sopnen2/2N.\ln p_n = -\frac{n(n-1)}{2N} + O\Bigl(\frac{n^3}{N^2}\Bigr), \qquad\text{so}\qquad p_n \approx \eu^{-n^2/2N} .

The tipping point pn=12p_n = \frac12 sits at n2Nln21.18Nn \approx \sqrt{2N\ln2} \approx 1.18\sqrt N: for N=365N = 365, n=23n = 23 (p23=0.4927p_{23} = 0.4927). Two morals. First, collisions among nn items in NN boxes appear at the scale nNn \sim \sqrt N, not nNn \sim N — the birthday scaling that governs hash collisions and the N\sqrt N cost of birthday attacks in cryptography. Second, the computation is a template: the (n2)\binom n2 pair-collision events are not independent, yet the answer behaves as if they were (e(n2)/N\eu^{-\binom n2/N} is exactly the independent-pairs heuristic) — a first instance of the Poisson approximation made rigorous in Chapter 23’s weekend problem (Le Cam’s inequality).

22.3 Borel–Cantelli and the zero–one law

Theorem 22.8 (Borel–Cantelli)

Let (An)(A_n) be events and lim supAn=NnNAn\limsup A_n = \bigcap_N \bigcup_{n\geq N}A_n (“AnA_n occurs infinitely often”).

  1. If P(An)<\sum\P(A_n) < \infty, then P(lim supAn)=0\P(\limsup A_n) = 0.
  2. If P(An)=\sum\P(A_n) = \infty and the AnA_n are independent, then P(lim supAn)=1\P(\limsup A_n) = 1.

Proof. (1) is Exercise 9.4. (2): for NMN \leq M, independence of complements (Exercise 22.3) gives

P(n=NMAnc)=n=NM(1P(An))exp(n=NMP(An))M0\P\Bigl(\bigcap_{n=N}^{M}A_n^c\Bigr) = \prod_{n=N}^M\bigl(1 - \P(A_n)\bigr) \leq \exp\Bigl(-\sum_{n=N}^M\P(A_n)\Bigr) \xrightarrow[M \to \infty]{} 0

(1xex1 - x \leq \eu^{-x}; the series diverges). So P(nNAn)=1\P\bigl(\bigcup_{n\geq N}A_n\bigr) = 1 for every NN, and the decreasing intersection over NN still has probability 11 (continuity from above, Proposition 9.6).

Theorem 22.9 (Kolmogorov’s zero–one law)

Let (Xn)(X_n) be independent and T=Nσ(XN,XN+1,)\mathcal T = \bigcap_N\sigma(X_N, X_{N+1}, \dots) the tail σ\sigma-algebra (events insensitive to any finite number of the XnX_n: convergence of Xn\sum X_n, of Snn\frac{S_n}n, values of lim sup\limsup’s, …). Then every TTT \in \mathcal T has P(T){0,1}\P(T) \in \{0, 1\}.

Proof. Fix NN. The σ\sigma-algebras σ(X1,,XN)\sigma(X_1, \dots, X_N) and σ(XN+1,)\sigma(X_{N+1}, \dots) are independent: events depending on disjoint blocks factorize on the generating π\pi-systems (cylinders iN{XiBi}\bigcap_{i\leq N}\{X_i \in B_i\}, resp. finite conditions on later variables), and Dynkin (Theorem 9.4, applied twice, one side at a time) extends the factorization. A tail event TT lies in σ(XN+1,)\sigma(X_{N+1}, \dots) for every NN: TT is independent of every σ(X1,,XN)\sigma(X_1, \dots, X_N), hence of the σ\sigma-algebra they generate, σ(X1,X2,)\sigma(X_1, X_2, \dots) (Dynkin once more: the union of the σ(X1,,XN)\sigma(X_1,\dots,X_N) is a π\pi-system generating it). But Tσ(X1,X2,)T \in \sigma(X_1, X_2, \dots) too: TT is independent of itself, P(T)=P(TT)=P(T)2\P(T) = \P(T\cap T) = \P(T)^2: P(T){0,1}\P(T) \in \{0, 1\}.

22.4 Modes of convergence

Definition 22.10

XnXX_n \to X almost surely if P(XnX)=1\P(X_n \to X) = 1; in probability if P(XnXε)0\P(\abs{X_n - X} \geq \varepsilon) \to 0 for every ε>0\varepsilon > 0; in LpL^p if EXnXp0\E\abs{X_n - X}^p \to 0.

Proposition 22.11

(a) a.s. convergence implies convergence in probability; (b) LpL^p convergence implies convergence in probability; (c) convergence in probability implies a.s. convergence along a subsequence; (d) no other implication holds in general.

Proof. (a) P(XnXε)P(supmnXmXε)P(lim sup{XmXε})=0\P(\abs{X_n - X} \geq \varepsilon) \leq \P\bigl(\sup_{m\geq n}\abs{X_m - X} \geq \varepsilon\bigr) \downarrow \P\bigl(\limsup\{\abs{X_m - X} \geq \varepsilon\}\bigr) = 0 under a.s. convergence (continuity from above; the limsup event excludes convergence). (b) Markov: P(XnXε)εpEXnXp\P(\abs{X_n - X} \geq \varepsilon) \leq \varepsilon^{-p}\,\E\abs{X_n - X}^p. (c) Pick nkn_k with P(XnkX2k)2k\P(\abs{X_{n_k} - X} \geq 2^{-k}) \leq 2^{-k}; Borel–Cantelli (1) makes XnkX<2k\abs{X_{n_k} - X} < 2^{-k} eventually, a.s. (d) The typewriter (Exercise 12.3) on ([0,1],λ)(\intcc01, \lambda) converges in L1L^1 and in probability but nowhere pointwise; n1(0,1/n)0n\mathbf 1_{\intoo0{1/n}} \to 0 a.s. but not in L1L^1; details and the remaining counterexamples in Exercise 22.6.

22.5 The law of large numbers

Throughout, (Xn)(X_n) are independent with the same law (i.i.d.), Sn=X1++XnS_n = X_1 + \dots + X_n.

Theorem 22.12 (Weak law of large numbers)

If X1L2X_1 \in L^2, with m=EX1m = \E X_1:

P(Snnmε)V(X1)nε2n0:\P\Bigl(\Bigl|\frac{S_n}{n} - m\Bigr| \geq \varepsilon\Bigr) \leq \frac{\V(X_1)}{n\,\varepsilon^2} \xrightarrow[n\to\infty]{} 0 :

Snnm\frac{S_n}n \to m in probability (and in L2L^2).

Proof. ESnn=m\E\frac{S_n}n = m and V(Snn)=nV(X1)n2\V\bigl(\frac{S_n}n\bigr) = \frac{n\V(X_1)}{n^2} (Theorem 22.5); Chebyshev.

Theorem 22.13 (Strong law of large numbers)

If X1L1X_1 \in L^1, then

Snnna.s.E[X1].\frac{S_n}{n} \xrightarrow[n\to\infty]{\text{a.s.}} \E[X_1].

We prove it here under the stronger hypothesis X1L4X_1 \in L^4; the general case (L1L^1: Etemadi’s proof) is the weekend problem.

Proof under EX14<\E X_1^4 < \infty. Centering (XiXimX_i \mapsto X_i - m), assume m=0m = 0. Expand:

E[Sn4]=i,j,k,lE[XiXjXkXl]=nE[X14]+3n(n1)(E[X12])2Cn2,\E[S_n^4] = \sum_{i,j,k,l}\E[X_iX_jX_kX_l] = n\,\E[X_1^4] + 3n(n-1)\,\bigl(\E[X_1^2]\bigr)^2 \leq C\,n^2 ,

since independence and centering kill every term containing an isolated factor (E[XiXjXkXl]=E[Xi]E[]=0\E[X_iX_jX_kX_l] = \E[X_i]\E[\cdots] = 0 unless the indices pair up: the only survivors are the nn terms i=j=k=li=j=k=l and the 3n(n1)3n(n-1) terms with two distinct pairs). Markov:

P(Snnε)=P(Sn4n4ε4)Cn2n4ε4=Cε4n2,\P\Bigl(\Bigl|\frac{S_n}n\Bigr| \geq \varepsilon\Bigr) = \P\bigl(S_n^4 \geq n^4\varepsilon^4\bigr) \leq \frac{Cn^2}{n^4\varepsilon^4} = \frac{C}{\varepsilon^4n^2},

summable: Borel–Cantelli (1) gives, for each rational ε\varepsilon, that Sn/n<ε\abs{S_n/n} < \varepsilon eventually, a.s.; intersecting over εQ+\varepsilon \in \Q_+^* (countably many probability-11 events): Sn/n0S_n/n \to 0 a.s.

Example 22.14 (What the strong law buys)

(a) Frequencies: for i.i.d. coin flips, the observed frequency of heads converges a.s. to pp — the empirical justification of probability itself. (b) Monte Carlo: for gL1([0,1])g \in L^1(\intcc01) and (Un)(U_n) i.i.d. uniform (Theorem 22.6), 1nkng(Uk)01g\frac1n\sum_{k\leq n}g(U_k) \to \int_0^1g a.s.: integrals by sampling, in any dimension, at the dimension-independent rate n1/2\sim n^{-1/2} made precise in Chapter 23. (c) Normal numbers: almost every real number has, in its binary expansion, asymptotic frequency 12\frac12 of ones (apply the strong law to the digit variables of Theorem 22.6) — Borel’s theorem, a statement about everyday numbers proved by measure: Problem 22.1 completes it in all bases.

Method 22.15

The working order for asymptotic statements about random sequences: (1) Is the event a tail event? Then its probability is 00 or 11 (Theorem 22.9) and one only has to decide which. (2) To prove a.s. statements: Borel–Cantelli — summable probabilities for the “bad” events, via Markov/Chebyshev-type bounds on whatever moments exist; independence only needed for the converse direction. (3) Subsequence + sandwich: prove convergence along a manageable subsequence, control the oscillation in between by monotonicity or maximal inequalities — the skeleton of Etemadi’s proof. (4) For distributional limits, wait for Chapter 23.

22.6 Exercises

Exercise 22.1

(a) Let XX have continuous strictly increasing distribution function FF. Show that F(X)F(X) is uniform on [0,1]\intcc01, and that G(U)FG(U) \sim F for UU uniform, G=F1G = F^{-1}: simulation by inversion. (b) Compute the distribution function and density of X2X^2 for XX uniform on [1,1]\intcc{-1}1, and of 1λlnU-\frac1\lambda\ln U for UU uniform on (0,1)\intoo01.

Solution

Solution of Exercise 22.1.

(a) For u(0,1)u \in \intoo01: P(F(X)u)=P(XF1(u))=F(F1(u))=u\P(F(X) \leq u) = \P(X \leq F^{-1}(u)) = F(F^{-1}(u)) = u (continuity and strict monotonicity make FF a bijection onto (0,1)\intoo01 with {F(X)u}={XF1(u)}\{F(X) \leq u\} = \{X \leq F^{-1}(u)\}): F(X)F(X) is uniform. Conversely P(G(U)t)=P(UF(t))=F(t)\P(G(U) \leq t) = \P(U \leq F(t)) = F(t): to simulate a law, apply the inverse distribution function to a uniform sample.

(b) Y=X2Y = X^2, XX uniform on [1,1]\intcc{-1}1: for t[0,1]t \in \intcc01, FY(t)=P(tXt)=tF_Y(t) = \P(-\sqrt t \leq X \leq \sqrt t) = \sqrt t: density 12t1(0,1)\frac1{2\sqrt t}\mathbf 1_{\intoo01}. And P(1λlnUt)=P(Ueλt)=1eλt\P\bigl(-\frac1\lambda\ln U \leq t\bigr) = \P(U \geq \eu^{-\lambda t}) = 1 - \eu^{-\lambda t}: the exponential E(λ)\mathcal E(\lambda) — inversion in action.

Exercise 22.2

(a) Compute mean and variance of the Poisson P(λ)\mathcal P(\lambda) and geometric laws via the transfer theorem. (b) Show that a positive random variable TT with P(T>t)>0\P(T > t) > 0 for all tt satisfies the memoryless property P(T>t+sT>t)=P(T>s)\P(T > t + s \mid T > t) = \P(T > s) for all s,t0s, t \geq 0 iff TT is exponential. (The survival function satisfies Cauchy’s functional equation; monotonicity replaces continuity.)

Solution

Solution of Exercise 22.2.

(a) Poisson: EX=k0keλλkk!=λ\E X = \sum_{k\geq0}k\,\eu^{-\lambda} \frac{\lambda^k}{k!} = \lambda, E[X(X1)]=λ2\E[X(X-1)] = \lambda^2, so V=λ2+λλ2=λ\V = \lambda^2 + \lambda - \lambda^2 = \lambda. Geometric (P(X=k)=p(1p)k1\P(X = k) = p(1-p)^{k-1}): EX=1p\E X = \frac1p, V=1pp2\V = \frac{1-p}{p^2} (differentiate the geometric series twice).

(b) G(t)=P(T>t)G(t) = \P(T > t) is nonincreasing with G(0+)G(0^+)\dots G ⁣:[0,)(0,1]G \colon \intco0\infty \to \intoc01; memorylessness reads G(t+s)=G(t)G(s)G(t + s) = G(t)G(s). Then G(nt)=G(t)nG(n t) = G(t)^n and G(t/n)=G(t)1/nG(t/n) = G(t)^{1/n}: G(q)=G(1)qG(q) = G(1)^q for rational q0q \geq 0; writing G(1)=eλG(1) = \eu^{-\lambda} ((0,1)\in \intoo01: G(1)=1G(1) = 1 would force G1G \equiv 1, impossible for a finite random variable; G(1)=0G(1) = 0 is excluded by hypothesis) and squeezing an arbitrary tt between rationals (monotonicity): G(t)=eλtG(t) = \eu^{-\lambda t} — the exponential law. The converse is a computation.

Exercise 22.3 ★★

(a) Show that if X1,,XnX_1, \dots, X_n are independent and fif_i are Borel functions, the fi(Xi)f_i(X_i) are independent. (b) Show that events A1,,AnA_1, \dots, A_n are independent iff their complements are, iff the indicators 1Ai\mathbf 1_{A_i} are independent random variables. (c) (Pairwise is weaker) Two fair coins: A=A = first is heads, B=B = second is heads, C=C = the two agree. Show A,B,CA, B, C are pairwise independent but not independent.

Solution

Solution of Exercise 22.3.

(a) σ(fi(Xi))=fi(Xi)1(B)Xi1(B)=σ(Xi)\sigma(f_i(X_i)) = f_i(X_i)^{-1}(\mathcal B) \subseteq X_i^{-1}(\mathcal B) = \sigma(X_i) (fif_i Borel), and sub-σ\sigma-algebras of independent σ\sigma-algebras are independent (the defining identity holds a fortiori).

(b) σ(Ai)={,Ai,Aic,Ω}=σ(Aic)=σ(1Ai)\sigma(A_i) = \{\varnothing, A_i, A_i^c, \Omega\} = \sigma(A_i^c) = \sigma(\mathbf 1_{A_i}): all three statements assert independence of the same σ\sigma-algebras. (That factorization over the AiA_i propagates to complements is the λ\lambda-system argument inside Definition 22.4’s equivalence — or direct inclusion-exclusion.)

(c) P(A)=P(B)=P(C)=12\P(A) = \P(B) = \P(C) = \frac12; AB=AC=BCA\cap B = A\cap C = B\cap C on pairs: each intersection is “both heads” or analogous, of probability 14\frac14: pairwise independent. But P(ABC)=P(HH)=1418\P(A\cap B\cap C) = \P(\text{HH}) = \frac14 \neq \frac18: not independentCC is determined by AA and BB.

Exercise 22.4 ★★

(a) (Infinite monkey) An i.i.d. sequence of uniform keystrokes on a finite alphabet a.s. contains every finite text infinitely often: prove it with Borel–Cantelli (2) on disjoint blocks. (b) (Runs) For i.i.d. fair bits, let RnR_n be the length of the run of ones starting at position nn. Show that a.s. Rn(1+ε)log2nR_n \geq (1+\varepsilon)\log_2n finitely often, and Rnlog2nR_n \geq \log_2 n infinitely often (both halves of Borel–Cantelli; for the second, pass to disjoint blocks to gain independence): the longest run in the first nn digits grows like log2n\log_2n.

Solution

Solution of Exercise 22.4.

(a) Let the text TT have length LL and q=aLq = a^{-L} (aa the alphabet size). The events Ek={E_k = \{positions kL+1,,(k+1)LkL+1, \dots, (k+1)L spell T}T\} are independent (disjoint blocks of i.i.d. letters), each of probability q>0q > 0: P(Ek)=\sum\P(E_k) = \infty, and Borel–Cantelli (2) gives infinitely many occurrences a.s.

(b) Upper: P(Rn(1+ε)log2n)2(1+ε)log2n=n(1+ε)\P\bigl(R_n \geq (1+\varepsilon)\log_2n\bigr) \leq 2^{-(1+\varepsilon)\log_2n} = n^{-(1+\varepsilon)}, summable: by Borel–Cantelli (1), a.s. only finitely many such nn. Lower: pack disjoint blocks — the jj-th of length j=log2sj\ell_j = \lceil\log_2s_j\rceil starting at sj=i<jis_j = \sum_{i<j}\ell_i; the events “block jj is all ones” are independent with probability 2j1sj1jlog2j2^{-\ell_j} \asymp \frac1{s_j} \asymp \frac1{j\log_2 j}, whose sum diverges: Borel–Cantelli (2) gives infinitely many all-ones blocks, i.e. Rsjlog2sjR_{s_j} \geq \log_2 s_j infinitely often. Together: the maximal run length in the first nn digits is (1+o(1))log2n(1 + o(1))\log_2n a.s.

Exercise 22.5 ★★

Let (Xn)(X_n) be independent. (a) Show that the radius of convergence of Xnzn\sum X_n z^n is an a.s. constant (possibly 00 or \infty). (b) Show that P(Xn converges){0,1}\P(\sum X_n \text{ converges}) \in \{0, 1\} and P(Sn/nm){0,1}\P(S_n/n \to m) \in \{0,1\}. (c) Give an event about (Xn)(X_n) that is not a tail event, and check that the zero–one law can fail for it.

Solution

Solution of Exercise 22.5.

(a) R=(lim supXn1/n)1R = \bigl(\limsup\abs{X_n}^{1/n}\bigr)^{-1} is unchanged if finitely many XnX_n are modified: for every NN, RR is σ(XN,XN+1,)\sigma(X_N, X_{N+1}, \dots)-measurable, i.e. tail-measurable. Then each event {Rc}\{R \leq c\} has probability 00 or 11 (Theorem 22.9), so the distribution function of RR takes only the values 0,10, 1: it jumps at a single point c0[0,+]c_0 \in \intcc0{+\infty}, and R=c0R = c_0 a.s.

(b) Convergence of Xn\sum X_n and of Snn\frac{S_n}n are insensitive to changing finitely many terms (for the second: the modified terms contribute O(1/n)0O(1/n) \to 0): tail events; zero–one law.

(c) {X1>0}\{X_1 > 0\} depends on X1X_1: for i.i.d. signs (P(X1=±1)=12\P(X_1 = \pm1) = \frac12), its probability is 12{0,1}\frac12 \notin \{0,1\} — no contradiction, it is not a tail event.

Exercise 22.6 ★★

On ([0,1],λ)(\intcc01, \lambda), exhibit — with proofs — random variables such that: (a) Xn0X_n \to 0 in probability and in every LpL^p, but nowhere a.s.; (b) Xn0X_n \to 0 a.s. but in no LpL^p; (c) Xn0X_n \to 0 in L1L^1 but not in L2L^2; (d) and show: if XnXX_n \to X in probability and XnYL1\abs{X_n} \leq Y \in L^1, then XnXX_n \to X in L1L^1 (subsequences + dominated convergence + the subsubsequence trick).

Solution

Solution of Exercise 22.6.

Work on ([0,1],λ)(\intcc01, \lambda). (a) The typewriter 1In\mathbf 1_{I_n} (Exercise 12.3): Xnpp=λ(In)0\norm{X_n}_p^p = \lambda(I_n) \to 0 (all p<p < \infty), hence also in probability; at every ω\omega the values 00 and 11 both recur: no pointwise convergence anywhere. (b) Xn=n1(0,1/n)0X_n = n\mathbf 1_{\intoo0{1/n}} \to 0 off 00, but Xnpn11/p1\norm{X_n}_p \geq n^{1 - 1/p} \geq 1. (c) Xn=n1(0,1/n)X_n = \sqrt n\,\mathbf 1_{\intoo0{1/n}}: EXn=n1/20\E\abs{X_n} = n^{-1/2} \to 0, EXn2=1\E X_n^2 = 1. (d) From any subsequence extract (convergence in probability) a further subsequence converging a.s. (Proposition 22.11(c)); dominated convergence gives L1L^1 convergence along it, with the same limit XX. Thus every subsequence of the numerical sequence EXnX\E\abs{X_n - X} has a subsubsequence tending to 00: the whole sequence tends to 00.

Exercise 22.7 ★★

An opinion poll estimates an unknown proportion pp by the empirical frequency p^n\hat p_n of nn independent draws. (a) Chebyshev: show P(p^npε)14nε2\P(\abs{\hat p_n - p} \geq \varepsilon) \leq \frac1{4n\varepsilon^2} (use p(1p)14p(1-p) \leq \frac14). (b) How many draws guarantee an error 3%\leq 3\% with probability 95%\geq 95\% by this bound? (The true answer, via Chapter 23, is about 10701070: Chebyshev is honest but crude.)

Solution

Solution of Exercise 22.7.

(a) p^n=Snn\hat p_n = \frac{S_n}n with SnS_n binomial: V(p^n)=p(1p)n14n\V(\hat p_n) = \frac{p(1-p)}n \leq \frac1{4n}, and Chebyshev (Proposition 22.3) gives the bound. (b) Solve 14n(0.03)20.05\frac1{4n(0.03)^2} \leq 0.05: n140.00090.055556n \geq \frac{1}{4\cdot0.0009\cdot0.05} \approx 5556. The central limit theorem will justify n1070n \approx 1070 for the same guarantee: Chebyshev pays for its generality with a factor 5\approx 5.

Exercise 22.8 ★★★

(Bernstein) For fC([0,1])f \in \mathcal C(\intcc01) define the Bernstein polynomial Bnf(x)=k=0n(nk)xk(1x)nkf(kn)B_nf(x) = \sum_{k=0}^n\binom nkx^k(1-x)^{n-k}f\bigl(\frac kn\bigr). (a) Recognize Bnf(x)=E[f(Snn)]B_nf(x) = \E\bigl[f\bigl(\frac {S_n}n\bigr)\bigr] for SnS_n binomial B(n,x)\mathcal B(n, x). (b) Prove BnffB_nf \to f uniformly on [0,1]\intcc01: split on {Snnxδ}\{\abs{\frac{S_n}n - x} \leq \delta\} and its complement, using uniform continuity and Chebyshev with the uniform bound V(Snn)14n\V(\frac{S_n}n) \leq \frac1{4n}. (c) Conclude: a second, probabilistic proof of the Weierstrass approximation theorem (Corollary 7.16), with the explicit rate Bnff32ωf(n1/2)\norm{B_nf - f}_\infty \leq \frac32\,\omega_f(n^{-1/2}) for the modulus of continuity ωf\omega_f — prove at least the O(ωf(n1/2))O(\omega_f(n^{-1/2})) form.

Solution

Solution of Exercise 22.8.

(a) If SnB(n,x)S_n \sim \mathcal B(n, x), the transfer theorem gives E[f(Snn)]=k(nk)xk(1x)nkf(kn)=Bnf(x)\E\bigl[f(\frac{S_n}n)\bigr] = \sum_k\binom nkx^k(1-x)^{n-k}f(\frac kn) = B_nf(x).

(b)–(c) Let ω=ωf\omega = \omega_f be the modulus of continuity (f(u)f(v)ω(uv)\abs{f(u) - f(v)} \leq \omega(\abs{u - v}), and ω(cδ)(1+c)ω(δ)\omega(c \delta) \leq (1 + c)\,\omega(\delta) by chaining steps). Then, for any δ>0\delta > 0,

f(u)f(x)(1+(ux)2δ2)ω(δ)\abs{f(u) - f(x)} \leq \Bigl(1 + \frac{(u - x)^2}{\delta^2}\Bigr)\omega(\delta)

(if uxδ\abs{u - x} \leq \delta, clear; otherwise ω(ux)(1+uxδ)ω(δ)(1+(ux)2δ2)ω(δ)\omega(\abs{u-x}) \leq (1 + \frac{\abs{u-x}}\delta) \omega(\delta) \leq (1 + \frac{(u-x)^2}{\delta^2}) \omega(\delta)). Take expectations at u=Snnu = \frac{S_n}n:

Bnf(x)f(x)(1+V(Sn/n)δ2)ω(δ)(1+14nδ2)ω(δ);\abs{B_nf(x) - f(x)} \leq \Bigl(1 + \frac{\V(S_n/n)}{\delta^2}\Bigr)\omega(\delta) \leq \Bigl(1 + \frac{1}{4n\delta^2}\Bigr)\omega(\delta) ;

with δ=n1/2\delta = n^{-1/2}: Bnff54ω(n1/2)32ω(n1/2)0\norm{B_nf - f}_\infty \leq \frac54\,\omega\bigl(n^{-1/2}\bigr) \leq \frac32\,\omega\bigl(n^{-1/2}\bigr) \to 0 (uniform continuity on the compact): a probabilistic Weierstrass theorem, with an explicit and uniform rate.

Exercise 22.9 ★★★

(Coupon collector) Cards of nn types are drawn uniformly with replacement; let TnT_n be the number of draws until all types are seen. (a) Write Tn=k=1nτkT_n = \sum_{k=1}^{n}\tau_k with τk\tau_k geometric of parameter nk+1n\frac{n - k + 1}n, the τk\tau_k independent, and deduce ETn=nHnnlnn\E T_n = n\,H_n \sim n\ln n (HnH_n the harmonic number) and V(Tn)π26n2\V(T_n) \leq \frac{\pi^2}6n^2. (b) Chebyshev: Tnnlnn1\frac{T_n}{n\ln n} \to 1 in probability. (c) Refine with Borel–Cantelli: show directly P(Tn>βnlnn)n1β\P(T_n > \beta n\ln n) \leq n^{1 - \beta} for β>1\beta > 1 (union bound on the event that some type is missed after βnlnn\beta n\ln n draws, using 1xex1 - x \leq \eu^{-x}), and deduce that along n=2mn = 2^m, a.s. TnβnlnnT_n \leq \beta n\ln n eventually, for every β>2\beta > 2.

Solution

Solution of Exercise 22.9.

(a) After k1k - 1 types are collected, each draw is new with probability pk=nk+1np_k = \frac{n-k+1}n: τk\tau_k is geometric (pk)(p_k), and the τk\tau_k are independent (the draws are). Sums: ETn=knnk+1=nHnnlnn\E T_n = \sum_k\frac n{n-k+1} = nH_n \sim n\ln n; V(Tn)=1pkpk2n2j=1n1j2π26n2\V(T_n) = \sum\frac{1 - p_k}{p_k^2} \leq n^2\sum_{j=1}^n\frac1{j^2} \leq \frac{\pi^2}6n^2.

(b) Chebyshev: P(TnnHnεnlnn)π2n2/6ε2n2ln2n0\P\bigl(\abs{T_n - nH_n} \geq \varepsilon n\ln n\bigr) \leq \frac{\pi^2n^2/6}{\varepsilon^2n^2\ln^2n} \to 0, and nHnnlnn1\frac{nH_n}{n\ln n} \to 1: Tnnlnn1\frac{T_n}{n\ln n} \to 1 in probability.

(c) Union bound: Tn>tT_n > t means some type is unseen after t\lceil t\rceil draws, so P(Tn>t)n(11n)tnet/n\P(T_n > t) \leq n(1 - \frac1n)^{t} \leq n\,\eu^{-t/n}; at t=βnlnnt = \beta n\ln n: n1β\leq n^{1 - \beta}. For β>1\beta > 1, m2m(1β)<\sum_m 2^{m(1-\beta)} < \infty: Borel–Cantelli gives, along n=2mn = 2^m, a.s. TnβnlnnT_n \leq \beta n\ln n eventually — in particular for every β>2\beta > 2 as stated (any β>1\beta > 1 works along the subsequence).

Exercise 22.10 ★★

Using the digit construction (Theorem 22.6): (a) verify by direct computation that U=b2k2kU = \sum b_{2k}2^{-k} (even-indexed digits of a uniform ω\omega) is uniform and independent of V=b2k12kV = \sum b_{2k-1}2^{-k}; (b) deduce a measurable bijection-up-to-null-sets between [0,1]\intcc01 and [0,1]2\intcc01^2 preserving measure, and comment: one uniform random number contains two (and countably many) independent ones — compare with the Peano curve (Problem 6.1), which achieved surjectivity but not measure-preservation or injectivity.

Solution

Solution of Exercise 22.10.

(a) The even-indexed digits (b2k)k(b_{2k})_k are i.i.d. fair bits (a subfamily of the independent digit family), so U=kb2k2kU = \sum_kb_{2k}2^{-k} gives every dyadic interval its correct probability (as in Theorem 22.6): uniform; likewise VV; and (U,V)(U, V) depend on disjoint digit blocks: independent (factorization on dyadic rectangles, then Dynkin).

(b) Φ(ω)=(U(ω),V(ω))\Phi(\omega) = (U(\omega), V(\omega)) is measurable with Φλ=λλ=λ2\Phi_*\lambda = \lambda\otimes\lambda = \lambda_2 (agreement on dyadic rectangles + uniqueness). Interleaving digits defines an inverse defined off the (null) set of dyadic rationals in either factor: a measure-preserving bijection between full-measure subsets of [0,1]\intcc01 and [0,1]2\intcc01^2. Contrast with Peano (Problem 6.1): continuity forced surjectivity without injectivity; dropping continuity for mere measurability buys a measure-isomorphism — dimension is invisible to measure theory, visible to topology.

Exercise 22.11 ★★

(Records) Let (Xn)n1(X_n)_{n\geq1} be i.i.d. with continuous distribution function, and say a record occurs at time nn if Xn>max(X1,,Xn1)X_n > \max(X_1, \dots, X_{n-1}) (time 11 is a record). Let RnR_n be the record indicator. (a) Show P(Rn=1)=1n\P(R_n = 1) = \frac1n (by symmetry, each of the n!n! orderings of X1,,XnX_1, \dots, X_n is equally likely and ties have probability 00). (b) Show that the RnR_n are independent (count orderings compatible with prescribed record positions, or argue that the relative order of X1,,Xn1X_1, \dots, X_{n-1} is independent of the rank of XnX_n among them). (c) Deduce from Borel–Cantelli (Theorem 22.8, both halves) that infinitely many records occur a.s., but records at consecutive times n,n+1n, n+1 occur infinitely often with probability — decide which! — and compute nP(Rn=1,Rn+1=1)\sum_n\P(R_n = 1, R_{n+1} = 1).

Solution

Solution of Exercise 22.11.

(a) Continuity of the distribution makes ties null events (as in the chapter’s order-statistics arguments), and the n!n! relative orderings of (X1,,Xn)(X_1, \dots, X_n) are exchangeable, hence equally likely. Rn=1R_n = 1 means the maximum sits in the last position: probability (n1)!n!=1n\frac{(n-1)!}{n!} = \frac1n.

(b) Fix nn and condition on the relative ordering of X1,,Xn1X_1, \dots, X_{n-1}: inserting XnX_n into the nn possible rank slots is uniform and independent of that ordering (exchangeability of the nn-tuple). Hence RnR_n (the event “XnX_n takes the top slot”) is independent of the whole record history (R1,,Rn1)(R_1, \dots, R_{n-1}), which is a function of the relative ordering of the first n1n - 1 variables. Induction gives full independence with P(Rn=1)=1n\P(R_n = 1) = \frac1n.

(c) P(Rn=1)=1n=\sum\P(R_n = 1) = \sum\frac1n = \infty with independence: the second Borel–Cantelli half gives records infinitely often a.s. (records never stop — but they thin out logarithmically: E[#recordsn]=Hnlnn\E[\#\text{records} \leq n] = H_n \approx \ln n). Consecutive records: P(Rn=Rn+1=1)=1n(n+1)\P(R_n = R_{n+1} = 1) = \frac1{n(n+1)} (independence), and

n1n(n+1)=n(1n1n+1)=1<:\sum_n\frac1{n(n+1)} = \sum_n\Bigl(\frac1n - \frac1{n+1}\Bigr) = 1 < \infty :

the first Borel–Cantelli half applies — only finitely many consecutive-record pairs occur, a.s.

Exercise 22.12 ★★

(Longest head run) Flip a fair coin infinitely often, and let LnL_n be the length of the longest run of consecutive heads within the first nn flips. (a) Show that for every ε>0\varepsilon > 0, a.s. Ln(1+ε)log2nL_n \leq (1 + \varepsilon)\log_2n eventually (the probability that some run of length \ell starts among the first nn flips is at most n2n2^{-\ell}; Borel–Cantelli along n=2kn = 2^k). (b) Show that a.s. Ln(1ε)log2nL_n \geq (1 - \varepsilon)\log_2n eventually (chop the first nn flips into n/\lfloor n/\ell\rfloor disjoint blocks of length =(1ε)log2n\ell = \lceil(1 - \varepsilon)\log_2n\rceil; the blocks are independent, each all-heads with probability 22^{-\ell}, and the probability that none is all-heads is at most exp(n2/)\exp(-n2^{-\ell}/\ell); sum along n=2kn = 2^k again). (c) Conclude Lnlog2n1\frac{L_n}{\log_2n} \to 1 a.s.: in a million fair flips one should expect a run of about 2020 heads — and a dataset without one is probably fabricated.

Solution

Solution of Exercise 22.12.

(a) A run of length \ell starting at position ini \leq n has probability 22^{-\ell}; union bound: P(Ln)n2\P(L_n \geq \ell) \leq n2^{-\ell}. With n=(1+ε)log2n\ell_n = (1 + \varepsilon)\log_2n: P(Lnn)nε\P(L_n \geq \ell_n) \leq n^{-\varepsilon}. Along n=2kn = 2^k: k2kε<\sum_k2^{-k\varepsilon} < \infty, so a.s. L2k<(1+ε)kL_{2^k} < (1+\varepsilon)k eventually (Borel–Cantelli); for general nn pick 2k1<n2k2^{k-1} < n \leq 2^k and use monotonicity of LnL_n plus log22k1log2n\log_22^{k-1} \leq \log_2n: LnL2k<(1+ε)k(1+ε)kk1log2nL_n \leq L_{2^k} < (1 + \varepsilon)k \leq (1 + \varepsilon)\frac{k}{k-1} \log_2n, and the extra factor is absorbed by enlarging ε\varepsilon slightly.

(b) With =(1ε)log2n\ell = \lceil(1 - \varepsilon)\log_2n\rceil and m=n/m = \lfloor n/\ell\rfloor disjoint blocks: the blocks are independent, each all-heads with probability 2n(1ε)/22^{-\ell} \geq n^{-(1-\varepsilon)}/2, so

P(Ln<)(12)mexp(m2)exp(cnεlog2n)\P(L_n < \ell) \leq \bigl(1 - 2^{-\ell}\bigr)^{m} \leq \exp\bigl(-m2^{-\ell}\bigr) \leq \exp\Bigl(-c\,\frac{n^{\varepsilon}}{\log_2n}\Bigr)

for a constant c>0c > 0 and nn large. These probabilities are summable along n=2kn = 2^k (indeed along all nn): Borel–Cantelli gives a.s. Ln(1ε)log2nL_n \geq (1 - \varepsilon)\log_2n eventually (monotonicity fills in between the 2k2^k as in (a), harmlessly).

(c) Both bounds along a sequence ε=1j\varepsilon = \frac1j, intersecting countably many full-measure events: Lnlog2n1\frac{L_n}{\log_2n} \to 1 a.s. For n=106n = 10^6: log2n19.9\log_2n \approx 19.9 — a run of 20\approx 20 heads is not a suspicious anomaly but a mathematical certainty, and its absence is evidence of a human faking “randomness” (humans rarely dare write more than 55 or 66 heads in a row).

22.7 Problem: Etemadi’s proof of the strong law

Problem 22.1

Weekend problem — the strong law of large numbers for i.i.d. integrable variables

Kolmogorov’s strong lawSnnEX1\frac{S_n}n \to \E X_1 a.s. for i.i.d. XnL1X_n \in L^1 — long had only intricate proofs; in 1981 N. Etemadi found one of striking economy, using nothing beyond this chapter (and even weakening independence to pairwise independence). We follow it. Let (Xn)(X_n) be pairwise independent, identically distributed, integrable; m=EX1m = \E X_1, Sn=X1++XnS_n = X_1 + \dots + X_n.

Part I — Reductions.

  1. Show that it suffices to treat Xn0X_n \geq 0 (split Xn=Xn+XnX_n = X_n^+ - X_n^-: check the two halves are again pairwise independent i.i.d. integrable). Assume henceforth Xn0X_n \geq 0.
  2. (Truncation) Let Yn=Xn1XnnY_n = X_n\,\mathbf 1_{X_n \leq n} and Sn=Y1++YnS_n^* = Y_1 + \dots + Y_n. Show

    n1P(XnYn)=n1P(X1>n)E[X1]<\sum_{n\geq1}\P(X_n \neq Y_n) = \sum_{n\geq1}\P(X_1 > n) \leq \E[X_1] < \infty

    (Exercise 11.3), and deduce via Borel–Cantelli that SnSnn0\frac{S_n - S_n^*}{n} \to 0 a.s.: it suffices to prove Snnm\frac{S^*_n}n \to m a.s.

  3. Show EYn=E[X11X1n]m\E Y_n = \E\bigl[X_1\mathbf 1_{X_1\leq n}\bigr] \to m (monotone convergence), hence 1nknEYkm\frac1n\sum_{k\leq n}\E Y_k \to m (Cesàro): it suffices to prove SnESnn0\frac{S_n^* - \E S_n^*}{n} \to 0 a.s.

Part II — The variance estimate.

  1. Show

    V(Yn)E[Yn2]=E[X121X1n]\V(Y_n) \leq \E[Y_n^2] = \E\bigl[X_1^2\,\mathbf 1_{X_1 \leq n}\bigr]

    and, using the layer cake (Proposition 11.8), the key bound

    n1V(Yn)n2n11n2E[X121X1n]CE[X1]<\sum_{n\geq1}\frac{\V(Y_n)}{n^2} \leq \sum_{n\geq1}\frac1{n^2}\, \E\bigl[X_1^2\mathbf 1_{X_1\leq n}\bigr] \leq C\,\E[X_1] < \infty

    (exchange the sum and the expectation — Tonelli for series — and bound nx1n22max(x,1)\sum_{n \geq x}\frac1{n^2} \leq \frac2{\max(x,1)} for the inner estimate x2nxn22xx^2\sum_{n\geq x}n^{-2} \leq 2x).

Part III — Convergence along geometric subsequences. Fix α>1\alpha > 1 and let kj=αjk_j = \lfloor\alpha^j\rfloor.

  1. Using pairwise independence (variances add, Theorem 22.5 — check that additivity of variances needs only pairwise independence) and Chebyshev, show for every ε>0\varepsilon > 0:

    j1P(SkjESkjkjε)1ε2j11kj2nkjV(Yn)=1ε2n1V(Yn)j:kjn1kj2.\sum_{j\geq1}\P\Bigl(\Bigl| \frac{S^*_{k_j} - \E S^*_{k_j}}{k_j}\Bigr| \geq \varepsilon\Bigr) \leq \frac1{\varepsilon^2}\sum_{j\geq1}\frac1{k_j^2} \sum_{n\leq k_j}\V(Y_n) = \frac1{\varepsilon^2}\sum_{n\geq1}\V(Y_n) \sum_{j\,:\,k_j\geq n}\frac1{k_j^2} .
  2. Show j:kjnkj2Cαn2\sum_{j : k_j \geq n}k_j^{-2} \leq \frac{C_\alpha}{n^2} (geometric series; beware the floor: kjαj2k_j \geq \frac{\alpha^j}2 for αj2\alpha^j \geq 2-type care), and conclude with question 4 and Borel–Cantelli:

    SkjESkjkjja.s.0,henceSkjkjm a.s.\frac{S^*_{k_j} - \E S^*_{k_j}}{k_j} \xrightarrow[j\to\infty]{\text{a.s.}} 0, \qquad\text{hence}\qquad \frac{S^*_{k_j}}{k_j} \to m \ \text{a.s.}

Part IV — Sandwich and conclusion.

  1. For kjnkj+1k_j \leq n \leq k_{j+1}, use the monotonicity of SnS^*_n (nonnegative summands!) to show

    kjkj+1Skjkj    Snn    kj+1kjSkj+1kj+1,\frac{k_j}{k_{j+1}}\,\frac{S^*_{k_j}}{k_j} \;\leq\; \frac{S^*_n}{n} \;\leq\; \frac{k_{j+1}}{k_j}\,\frac{S^*_{k_{j+1}}}{k_{j+1}},

    and deduce, a.s.:

    mαlim infSnnlim supSnnαm.\frac m\alpha \leq \liminf\frac{S^*_n}n \leq \limsup\frac{S^*_n}n \leq \alpha\,m .
  2. Let α1\alpha \downarrow 1 along a sequence and conclude Snnm\frac{S_n^*}n \to m a.s., hence (Part I) the strong law of large numbers:

     Snnna.s.E[X1]. \boxed{\ \frac{S_n}{n} \xrightarrow[n\to\infty]{\text{a.s.}} \E[X_1].\ }
  3. Where exactly was pairwise independence (rather than full independence) sufficient? List the three places where independence-type hypotheses were invoked.

Part V — Dividends.

  1. (Borel’s normal numbers) Show that λ\lambda-almost every x[0,1]x \in \intcc01 is normal in every base b2b \geq 2: each digit 0,,b10, \dots, b-1 appears with asymptotic frequency 1b\frac1b (fix bb and a digit, apply the strong law to the indicator variables — justify that base-bb digits of a uniform variable are i.i.d. uniform on {0,,b1}\{0,\dots,b-1\} as in Theorem 22.6 — then intersect the countably many probability-one events). Exhibit one explicit non-normal number, and reflect: the theorem asserts normality of almost all numbers, yet proving normality of 2\sqrt2 or π\pi remains open.
  2. (Monte Carlo, guaranteed) Justify completely the method of Example 22.14(b) for gL1([0,1]d)g \in L^1(\intcc01^d): construct the i.i.d. uniform sample on [0,1]d\intcc01^d from Theorem 22.6 and Exercise 22.10, and state what the strong law delivers.

Part VI — What full independence buys: maximal inequalities and random series. Etemadi spends only pairwise independence; the remaining parts exploit the full (mutual) version. Let (Zn)(Z_n) be independent centered variables of L2L^2 and Sk=Z1++ZkS_k = Z_1 + \dots + Z_k (a fresh notation, unrelated to the XnX_n above).

  1. (Kolmogorov’s maximal inequality) For ε>0\varepsilon > 0 prove

    P(max1knSkε)    1ε2k=1nV(Zk):\P\Bigl(\max_{1\leq k\leq n}\abs{S_k} \geq \varepsilon\Bigr) \;\leq\; \frac1{\varepsilon^2}\sum_{k=1}^n\V(Z_k) :

    Chebyshev’s price buys the maximum (partition the event according to the first index kk with Skε\abs{S_k} \geq \varepsilon; on that piece write Sn2Sk2+2Sk(SnSk)S_n^2 \geq S_k^2 + 2S_k(S_n - S_k) and use the independence of the coalitions (Z1,,Zk)(Z_1, \dots, Z_k) and (Zk+1,,Zn)(Z_{k+1}, \dots, Z_n), Theorem 22.5). Point out the step where pairwise independence would no longer suffice.

  2. (Khinchin–Kolmogorov one-series theorem) Deduce: if nV(Zn)<\sum_n\V(Z_n) < \infty, then nZn\sum_nZ_n converges almost surely (show that a.s. the partial sums form a Cauchy sequence: let mm \to \infty in the maximal inequality applied to ZN+1,,ZN+mZ_{N+1}, \dots, Z_{N+m}, then let NN \to \infty).
  3. (Rademacher series) Let (εn)(\varepsilon_n) be i.i.d. signs, P(εn=±1)=12\P(\varepsilon_n = \pm1) = \frac12 (Theorem 22.6), and let (xn)(x_n) be real numbers. Show that nxnεn\sum_nx_n\varepsilon_n converges a.s. as soon as nxn2<\sum_nx_n^2 < \infty; show also that, whatever (xn)(x_n), the probability that nxnεn\sum_nx_n\varepsilon_n converges is 00 or 11 (Theorem 22.9).
  4. The converse, elementarily. Set Tn=knxkεkT_n = \sum_{k\leq n}x_k\varepsilon_k and sn2=knxk2s_n^2 = \sum_{k\leq n}x_k^2, and suppose sns_n \to \infty. (a) Prove the Paley–Zygmund inequality: for Z0Z \geq 0 with EZ2<\E Z^2 < \infty and 0<θ<10 < \theta < 1,

    P(Z>θEZ)    (1θ)2(EZ)2EZ2\P\bigl(Z > \theta\,\E Z\bigr) \;\geq\; (1 - \theta)^2\,\frac{(\E Z)^2}{\E Z^2}

    (split EZ\E Z at the level θEZ\theta\E Z and apply Cauchy–Schwarz to the upper piece). (b) Show ETn43sn4\E T_n^4 \leq 3s_n^4. (c) Deduce P(Tn>sn2)316\P\bigl(\abs{T_n} > \frac{s_n}2\bigr) \geq \frac3{16} and conclude that nxnεn\sum_nx_n\varepsilon_n diverges a.s.; hence the dichotomy

    nxnεn converges a.s.    nxn2<.\sum_nx_n\varepsilon_n\ \text{converges a.s.} \iff \sum_nx_n^2 < \infty .
  5. (Random harmonic series) Conclude that nεnns\sum_n\frac{\varepsilon_n}{n^s} converges a.s. if and only if s>12s > \frac12. For 12<s1\frac12 < s \leq 1 the series converges a.s. while nns=\sum_nn^{-s} = \infty: random signs produce square-root-strength cancellation — compare with the alternating series n(1)nns\sum_n\frac{(-1)^n}{n^s}, which converges for every s>0s > 0.

Part VII — Concentration: Hoeffding’s inequality. The strong law says Snnm\frac{S_n}n \to m; concentration inequalities say how unlikely a deviation is at each fixed nn.

  1. (Hoeffding’s lemma) (a) Show coshλeλ2/2\cosh\lambda \leq \eu^{\lambda^2/2} for all λR\lambda \in \R, by comparing the two series termwise. (b) Let ZZ be centered with aZba \leq Z \leq b, a<ba < b. Show

    EeλZexp(λ2(ba)28)\E\,\eu^{\lambda Z} \leq \exp\Bigl(\frac{\lambda^2(b - a)^2}8\Bigr)

    (bound eλz\eu^{\lambda z} on [a,b]\intcc ab by its chord, take expectations, and study φ(t)=pt+log(1p+pet)\varphi(t) = -pt + \log(1 - p + p\eu^t) with p=abap = \frac{-a}{b-a} and t=λ(ba)t = \lambda(b - a): show φ(0)=φ(0)=0\varphi(0) = \varphi'(0) = 0 and φ14\varphi'' \leq \frac14).

  2. (Hoeffding’s inequality) Let X1,,XnX_1, \dots, X_n be independent with aiXibia_i \leq X_i \leq b_i and Sn=X1++XnS_n = X_1 + \dots + X_n. Prove, for t>0t > 0,

    P(SnESnt)exp(2t2i=1n(biai)2),\P\bigl(S_n - \E S_n \geq t\bigr) \leq \exp\Bigl(\frac{-2t^2}{\sum_{i=1}^n(b_i - a_i)^2}\Bigr),

    and the same bound for the lower tail (exponential Chebyshev: bound Eeλ(SnESn)\E\,\eu^{\lambda(S_n - \E S_n)} using independence and question 17, then optimize over λ>0\lambda > 0).

  3. (The strong law, bounded case, with a rate) Let the XiX_i be i.i.d. with values in [a,b]\intcc ab and m=EX1m = \E X_1. Show

    P(Snnmε)2exp(2nε2(ba)2)\P\Bigl(\Bigl|\frac{S_n}n - m\Bigr| \geq \varepsilon\Bigr) \leq 2\exp\Bigl(\frac{-2n\varepsilon^2}{(b - a)^2}\Bigr)

    and recover Snnm\frac{S_n}n \to m a.s. by Borel–Cantelli: a second proof of the strong law for bounded variables — no truncation, an exponential rate at every finite nn, but bounded summands and full independence. Compare the hypotheses with Etemadi’s.

  4. (Monte Carlo, guaranteed at fixed nn) Let g ⁣:[0,1]d[0,1]g \colon \intcc01^d \to \intcc01 be measurable and (Uk)(U_k) the i.i.d. uniform sample of question 11. Given ε,δ>0\varepsilon, \delta > 0, show

    nlog(2/δ)2ε2    P(1nk=1ng(Uk)g ⁣dλdε)δ,n \geq \frac{\log(2/\delta)}{2\varepsilon^2} \implies \P\Bigl(\Bigl|\frac1n\sum_{k=1}^ng(U_k) - \int g\,\dd\lambda_d\Bigr| \geq \varepsilon\Bigr) \leq \delta,

    and evaluate the threshold for ε=δ=102\varepsilon = \delta = 10^{-2}. The bound does not involve dd: compare with question 11 and with deterministic grids.

Part VIII — How big is a random walk? Toward the iterated logarithm. Let Sn=ε1++εnS_n = \varepsilon_1 + \dots + \varepsilon_n be the simple random walk built from i.i.d. fair signs.

  1. (Sub-Gaussian tails) Show EeλSn=(coshλ)nenλ2/2\E\,\eu^{\lambda S_n} = (\cosh\lambda)^n \leq \eu^{n\lambda^2/2} and deduce, for x>0x > 0,

    P(Snx)ex2/(2n),P(Snx)2ex2/(2n).\P(S_n \geq x) \leq \eu^{-x^2/(2n)}, \qquad \P(\abs{S_n} \geq x) \leq 2\,\eu^{-x^2/(2n)} .
  2. Deduce, via Borel–Cantelli,

    lim supnSn2nlogn1a.s.\limsup_{n\to\infty}\frac{\abs{S_n}} {\sqrt{2n\log n}} \leq 1 \quad\text{a.s.}

    (for η>0\eta > 0, sum the tail bounds at x=(1+η)2nlognx = (1 + \eta)\sqrt{2n\log n}, then intersect over η=1p\eta = \frac1p). In particular the walk lives on the CLT scale n\sqrt n up to a logarithmic factor — far below the crude bound Snn\abs{S_n} \leq n.

  3. Along the doubling subsequence nj=2jn_j = 2^j, show

    lim supjSnj2njloglognj1a.s.,\limsup_{j\to\infty}\frac{S_{n_j}} {\sqrt{2n_j\log\log n_j}} \leq 1 \quad\text{a.s.},

    and reflect: the law of the iterated logarithm (Khinchin; Hartman–Wintner for general centered L2L^2 summands) states that

    lim supnSn2nloglogn=1a.s.\limsup_{n\to\infty}\frac{S_n} {\sqrt{2n\log\log n}} = 1 \quad\text{a.s.}

    Explain precisely what separates the subsequence estimate just proved from the upper half of this statement (one must control maxnjnnj+1Sn\max_{n_j \leq n \leq n_{j+1}}S_n inside each block, which requires a maximal inequality at the exponential scale) and check quantitatively that question 12’s inequality is too weak for that purpose. The lower half rests on the second Borel–Cantelli lemma applied to independent blocks; both halves are honest Year 3 material for a dedicated probability course.

  4. (Uniform deviation over a finite class) Let A1,,ANA_1, \dots, A_N be events in a repeatable experiment, and estimate each probability by its empirical frequency p^i\hat p_i over nn i.i.d. repetitions. Combining Hoeffding’s inequality with a union bound, show

    P(maxiNp^iP(Ai)>ε)    2Ne2nε2,\P\Bigl(\max_{i\leq N}\,\abs{\hat p_i - \P(A_i)} > \varepsilon\Bigr) \;\leq\; 2N\,\eu^{-2n\varepsilon^2},

    and deduce the sample-size rule: nln(2N/δ)2ε2n \geq \frac{\ln(2N/\delta)}{2\varepsilon^2} guarantees all NN estimates simultaneously ε\varepsilon-accurate with probability 1δ\geq 1 - \delta. Compute nn for N=106N = 10^6, ε=0.01\varepsilon = 0.01, δ=0.05\delta = 0.05: the logarithmic price of uniformity.

  5. (The random harmonic window) Combining the two halves of the random series theory, show that for i.i.d. signs (εn)(\varepsilon_n) the series nεnnα\sum_n\frac{\varepsilon_n}{n^\alpha} converges a.s. if α>12\alpha > \frac12 and diverges a.s. if α12\alpha \leq \frac12; contrast with absolute convergence (which requires α>1\alpha > 1): on the window α(12,1]\alpha \in \intoc{\frac12}1, convergence is a genuinely probabilistic phenomenon — cancellation, not size.
Solution

Solution of Problem 22.1.

1. Xn±X_n^{\pm} are Borel functions of XnX_n: they remain pairwise independent (Exercise 22.3(a)) and identically distributed, integrable, with EX1=EX1+EX1\E X_1 = \E X_1^+ - \E X_1^-. If the theorem holds for nonnegative variables, apply it to both halves and subtract: Snn=Sn+nSnnEX1+EX1=m\frac{S_n}n = \frac{S_n^+}n - \frac{S_n^-}n \to \E X_1^+ - \E X_1^- = m a.s.

2. P(XnYn)=P(Xn>n)=P(X1>n)\P(X_n \neq Y_n) = \P(X_n > n) = \P(X_1 > n) (identical laws), and nP(X1>n)nP(X1n)EX1<\sum_n\P(X_1 > n) \leq \sum_n\P(X_1 \geq n) \leq \E X_1 < \infty (Exercise 11.3(a)). Borel–Cantelli (1): a.s. Xn=YnX_n = Y_n for all large nn, so SnSnS_n - S_n^* is eventually constant in nn: SnSnn0\frac{S_n - S_n^*}n \to 0 a.s., and the two normalized sums share their asymptotic behavior.

3. X11X1nX1X_1\mathbf 1_{X_1 \leq n} \nearrow X_1: MCT gives EYnm\E Y_n \to m; Cesàro means of a convergent sequence converge to the same limit: ESnn=1nknEYkm\frac{\E S_n^*}n = \frac1n\sum_{k\leq n}\E Y_k \to m. Hence it suffices to prove SnESnn0\frac{S^*_n - \E S^*_n}{n} \to 0 a.s.

4. V(Yn)EYn2=E[X121X1n]\V(Y_n) \leq \E Y_n^2 = \E[X_1^2\mathbf 1_{X_1\leq n}]. By Tonelli for series,

nE[X121X1n]n2=E[X12 ⁣ ⁣nmax(X1,1) ⁣1n2]E[X124max(X1,1)]4E[X1]<,\sum_n\frac{\E[X_1^2\mathbf 1_{X_1\leq n}]}{n^2} = \E\Bigl[X_1^2\!\!\sum_{n \geq \max(X_1, 1)}\!\frac1{n^2} \Bigr] \leq \E\Bigl[X_1^2\cdot\frac{4}{\max(X_1,1)}\Bigr] \leq 4\,\E[X_1] < \infty,

using nxn24x\sum_{n\geq x}n^{-2} \leq \frac4x for x1x \geq 1 (for x2x \geq 2: 1x12x\leq \frac1{x-1} \leq \frac2x; for 1x<21 \leq x < 2: π264x\leq \frac{\pi^2}6 \leq \frac4x since 4x>2\frac4x > 2), and X12/max(X1,1)X1X_1^2/\max(X_1, 1) \leq X_1 in both cases X11X_1 \gtrless 1.

5. Pairwise independence gives E[(YiEYi)(YjEYj)]=0\E[(Y_i - \E Y_i)(Y_j - \E Y_j)] = 0 for iji \neq j (the product formula for two variables), so variances add: V(Sk)=nkV(Yn)\V(S^*_k) = \sum_{n\leq k}\V(Y_n). Chebyshev on each kjk_j and summing:

jP(SkjESkjεkj)1ε2j1kj2nkjV(Yn)=1ε2nV(Yn) ⁣ ⁣j:kjn ⁣1kj2\sum_j\P\Bigl(\abs{S^*_{k_j} - \E S^*_{k_j}} \geq \varepsilon k_j\Bigr) \leq \frac1{\varepsilon^2}\sum_j\frac1{k_j^2}\sum_{n\leq k_j}\V(Y_n) = \frac1{\varepsilon^2}\sum_n\V(Y_n)\!\!\sum_{j : k_j\geq n}\!\frac1{k_j^2}

(Tonelli for the nonnegative double series).

6. kj=αjαj2k_j = \lfloor\alpha^j\rfloor \geq \frac{\alpha^j}2 (valid once αj1\alpha^j \geq 1, i.e. all j0j \geq 0: xx2\lfloor x\rfloor \geq \frac x2 for x1x \geq 1). Hence

j:kjn1kj24j:αjnα2j41α21n2=Cαn2,\sum_{j : k_j \geq n}\frac1{k_j^2} \leq 4\sum_{j : \alpha^j \geq n}\alpha^{-2j} \leq \frac{4}{1 - \alpha^{-2}}\cdot\frac1{n^2} = \frac{C_\alpha}{n^2},

(geometric series from the first jj with αjn\alpha^j \geq n). Combining with questions 4–5, the double sum is finite; Borel–Cantelli (1), applied for each rational ε\varepsilon and intersected, gives SkjESkjkj0\frac{S^*_{k_j} - \E S^*_{k_j}}{k_j} \to 0 a.s., and with question 3: Skjkjm\frac{S^*_{k_j}}{k_j} \to m a.s.

7. Yn0Y_n \geq 0 makes nSnn \mapsto S^*_n nondecreasing: for kjnkj+1k_j \leq n \leq k_{j+1},

Skjkj+1SnnSkj+1kj,\frac{S^*_{k_j}}{k_{j+1}} \leq \frac{S^*_n}{n} \leq \frac{S^*_{k_{j+1}}}{k_j},

which is the displayed sandwich after inserting kjkj+1\frac{k_j}{k_{j+1}} and kj+1kj\frac{k_{j+1}}{k_j}. Since kj+1kjα\frac{k_{j+1}}{k_j} \to \alpha, question 6 gives a.s.

mαlim infnSnnlim supnSnnαm.\frac m\alpha \leq \liminf_n\frac{S^*_n}n \leq \limsup_n\frac{S^*_n}n \leq \alpha m .

8. Apply question 7 for α=1+1p\alpha = 1 + \frac1p, pNp \in \N^*: countably many a.s. events; on their intersection, letting pp \to \infty: limSnn=m\lim\frac{S^*_n}n = m a.s. With questions 1–3, SnnEX1\frac{S_n}n \to \E X_1 a.s.: the strong law of large numbers, under pairwise independence.

9. Independence-type hypotheses appeared three times: (i) additivity of variances (question 5) — pairwise suffices; (ii) identical distribution, in the truncation sums (question 2) and the mean computation (question 3) — no independence at all; (iii) Borel–Cantelli (1) (questions 2 and 6) — valid without any independence. Full mutual independence was never invoked: Etemadi’s observation.

10. Fix a base bb and a digit rr. The base-bb digits (dk)(d_k) of a uniform ω\omega are i.i.d. uniform on {0,,b1}\{0, \dots, b-1\} (each digit vector value occupies an interval of length bmb^{-m}: the argument of Theorem 22.6 verbatim). The strong law applied to the i.i.d. bounded variables 1dk=r\mathbf 1_{d_k = r} gives: a.s., the frequency of the digit rr tends to 1b\frac1b. Intersecting over the countably many pairs (b,r)(b, r): almost every number is simply normal in every base. An explicit non-normal number: x=0.1001001002x = 0.100100100\ldots_2 (frequency of ones 1312\frac13 \neq \frac12). The contrast is humbling: almost all numbers are normal, yet for 2\sqrt2, e\eu, or π\pi normality remains unproved — measure theory counts without exhibiting.

11. By Exercise 22.10 iterated, a single uniform variable yields a sequence of i.i.d. uniform vectors UkU_k on [0,1]d\intcc01^d (split the digit set of each UnU_n of Theorem 22.6 into dd subfamilies). For gL1([0,1]d)g \in L^1(\intcc01^d), the variables g(Uk)g(U_k) are i.i.d. integrable with mean g ⁣dλd\int g\,\dd\lambda_d (transfer): the strong law gives

1nk=1ng(Uk)na.s.[0,1]dg ⁣dλd:\frac1n\sum_{k=1}^ng(U_k) \xrightarrow[n\to\infty]{\text{a.s.}} \int_{\intcc01^d}g\,\dd\lambda_d :

Monte Carlo integration converges almost surely, in every dimension — the error size is the business of the central limit theorem (Chapter 23).

12. Let Ak={Skε}j<k{Sj<ε}A_k = \{\abs{S_k} \geq \varepsilon\} \cap \bigcap_{j<k}\{\abs{S_j} < \varepsilon\}: the AkA_k are disjoint with union A={maxknSkε}A = \{\max_{k\leq n}\abs{S_k} \geq \varepsilon\}. Then

ESn2k=1nE[Sn21Ak]=k=1nE[(Sk2+2Sk(SnSk)+(SnSk)2)1Ak]k=1nE[Sk21Ak],\E S_n^2 \geq \sum_{k=1}^n\E\bigl[S_n^2\mathbf 1_{A_k}\bigr] = \sum_{k=1}^n\E\Bigl[\bigl(S_k^2 + 2S_k(S_n - S_k) + (S_n - S_k)^2\bigr)\mathbf 1_{A_k}\Bigr] \geq \sum_{k=1}^n\E\bigl[S_k^2\mathbf 1_{A_k}\bigr],

because the cross term vanishes: Sk1AkS_k\mathbf 1_{A_k} is a Borel function of the coalition (Z1,,Zk)(Z_1, \dots, Z_k), which is independent of SnSkS_n - S_k, a function of (Zk+1,,Zn)(Z_{k+1}, \dots, Z_n) (Theorem 22.5), so E[Sk1Ak(SnSk)]=E[Sk1Ak]E[SnSk]=0\E[S_k\mathbf 1_{A_k}(S_n - S_k)] = \E[S_k\mathbf 1_{A_k}]\,\E[S_n - S_k] = 0. On AkA_k, Sk2ε2S_k^2 \geq \varepsilon^2, whence ESn2ε2kP(Ak)=ε2P(A)\E S_n^2 \geq \varepsilon^2\sum_k\P(A_k) = \varepsilon^2\P(A); and ESn2=knV(Zk)\E S_n^2 = \sum_{k\leq n}\V(Z_k) (variances add). The decisive step is the factorization: Sk1AkS_k\mathbf 1_{A_k} is a nonlinear function of the whole first block, and its independence from the second block is coalition independence — pairwise independence of the ZiZ_i only decorrelates pairs and would not justify it.

13. Fix NN and apply question 12 to ZN+1,,ZN+mZ_{N+1}, \dots, Z_{N+m}:

P(maxN<kN+mSkSN>ε)1ε2j=N+1N+mV(Zj)rNε2,rN=j>NV(Zj).\P\Bigl(\max_{N < k \leq N+m}\abs{S_k - S_N} > \varepsilon\Bigr) \leq \frac1{\varepsilon^2}\sum_{j=N+1}^{N+m}\V(Z_j) \leq \frac{r_N}{\varepsilon^2}, \qquad r_N = \sum_{j>N}\V(Z_j) .

The events increase with mm; continuity from below gives P(supk>NSkSN>ε)rN/ε2\P(\sup_{k>N}\abs{S_k - S_N} > \varepsilon) \leq r_N/\varepsilon^2, and rN0r_N \to 0 by hypothesis. Hence for each pNp \in \N^*, P(N{supk>NSkSN>1p})infNp2rN=0\P\bigl(\bigcap_N\{\sup_{k>N} \abs{S_k - S_N} > \frac1p\}\bigr) \leq \inf_Np^2r_N = 0: almost surely, for every pp there is NN with supk>NSkSN1p\sup_{k>N}\abs{S_k - S_N} \leq \frac1p (intersect the countably many a.s. events over pp), so that SkSl2p\abs{S_k - S_l} \leq \frac2p for all k,l>Nk, l > N: the partial sums are a.s. Cauchy, hence a.s. convergent.

14. The variables Zn=xnεnZ_n = x_n\varepsilon_n are independent (Borel functions of independent variables, Exercise 22.3(a)), centered, with V(Zn)=xn2\V(Z_n) = x_n^2: question 13 applies when nxn2<\sum_nx_n^2 < \infty and gives a.s. convergence. In general, for each NN the convergence of nxnεn\sum_nx_n\varepsilon_n is unaffected by the values of ε1,,εN\varepsilon_1, \dots, \varepsilon_N: the convergence event lies in the tail σ\sigma-algebra of the independent sequence (εn)(\varepsilon_n), so Kolmogorov’s zero–one law (Theorem 22.9) forces its probability to be 00 or 11.

15. (a) Splitting at the level θEZ\theta\E Z and using Cauchy–Schwarz on the upper piece,

EZ=E[Z1ZθEZ]+E[Z1Z>θEZ]θEZ+EZ2P(Z>θEZ),\E Z = \E\bigl[Z\mathbf 1_{Z \leq \theta\E Z}\bigr] + \E\bigl[Z\mathbf 1_{Z > \theta\E Z}\bigr] \leq \theta\,\E Z + \sqrt{\E Z^2}\, \sqrt{\P(Z > \theta\E Z)} ,

so (1θ)EZEZ2P(Z>θEZ)(1 - \theta)\E Z \leq \sqrt{\E Z^2\,\P(Z > \theta\E Z)}; square. (b) Expand Tn4=i,j,k,lxixjxkxlE[εiεjεkεl]T_n^4 = \sum_{i,j,k,l}x_ix_jx_kx_l\, \E[\varepsilon_i\varepsilon_j\varepsilon_k\varepsilon_l]: the expectation is 11 when the indices pair off (all four equal, or two distinct pairs, the latter in 33 arrangements) and 00 otherwise (an unpaired sign has zero mean and factors out by independence). Hence

ETn4=kxk4+3ijxi2xj2=3sn42kxk43sn4.\E T_n^4 = \sum_kx_k^4 + 3\sum_{i\neq j}x_i^2x_j^2 = 3s_n^4 - 2\sum_kx_k^4 \leq 3s_n^4 .

(c) Paley–Zygmund with Z=Tn2Z = T_n^2, EZ=sn2\E Z = s_n^2, θ=14\theta = \frac14:

P(Tn>sn2)=P(Tn2>sn24)(34)2sn43sn4=316.\P\Bigl(\abs{T_n} > \frac{s_n}2\Bigr) = \P\Bigl(T_n^2 > \frac{s_n^2}4\Bigr) \geq \Bigl(\frac34\Bigr)^2 \frac{s_n^4}{3s_n^4} = \frac3{16} .

If the series converged with positive probability, it would converge a.s. (question 14), so supnTn<\sup_n\abs{T_n} < \infty a.s., and some MM would satisfy P(supnTn>M)<316\P(\sup_n\abs{T_n} > M) < \frac3{16}; but as soon as sn>2Ms_n > 2M, P(Tn>M)P(Tn>sn2)316\P(\abs{T_n} > M) \geq \P(\abs{T_n} > \frac{s_n}2) \geq \frac3{16}: contradiction. So divergence is almost sure, and with question 14 the dichotomy is complete.

16. Here xn=nsx_n = n^{-s} and nn2s<\sum_nn^{-2s} < \infty exactly when s>12s > \frac12: by questions 14–15, nεnns\sum_n\frac{\varepsilon_n}{n^s} converges a.s. if and only if s>12s > \frac12 (for s12s \leq \frac12, a.s. divergence). For 12<s1\frac12 < s \leq 1 the convergence is never absolute. The comparison is instructive: perfectly alternating signs cancel at strength nsn^{-s} for every s>0s > 0, while typical random signs cancel only at square-root strength — the random walk of question 21 grows like n\sqrt n, and Abel summation converts exactly that growth into convergence of εnns\sum\varepsilon_nn^{-s} for s>12s > \frac12.

17. (a) coshλ=kλ2k(2k)!\cosh\lambda = \sum_k\frac{\lambda^{2k}}{(2k)!} and eλ2/2=kλ2k2kk!\eu^{\lambda^2/2} = \sum_k\frac{\lambda^{2k}}{2^kk!}; and (2k)!2kk!(2k)! \geq 2^kk! holds termwise, because (2k)!k!=i=1k(k+i)i=1k(2i)=2kk!\frac{(2k)!}{k!} = \prod_{i=1}^k(k + i) \geq \prod_{i=1}^k(2i) = 2^kk! (each factor satisfies k+i2ik + i \geq 2i for iki \leq k), so that in fact (2k)!2k(k!)22kk!(2k)! \geq 2^k(k!)^2 \geq 2^kk!. (b) Note a0ba \leq 0 \leq b (ZZ is centered), and by convexity of zeλzz \mapsto \eu^{\lambda z}, for z[a,b]z \in \intcc ab:

eλzbzbaeλa+zabaeλb,soEeλZbeλaaeλbba=(1p)ept+pe(1p)t=eφ(t)\eu^{\lambda z} \leq \frac{b - z}{b - a}\,\eu^{\lambda a} + \frac{z - a}{b - a}\,\eu^{\lambda b}, \qquad\text{so}\qquad \E\,\eu^{\lambda Z} \leq \frac{b\,\eu^{\lambda a} - a\,\eu^{\lambda b}}{b - a} = (1 - p)\eu^{-pt} + p\,\eu^{(1-p)t} = \eu^{\varphi(t)}

with p=aba[0,1]p = \frac{-a}{b-a} \in \intcc01, t=λ(ba)t = \lambda(b - a), φ(t)=pt+log(1p+pet)\varphi(t) = -pt + \log(1 - p + p\eu^t). Then φ(0)=0\varphi(0) = 0, φ(t)=p+pet1p+pet\varphi'(t) = -p + \frac{p\eu^t}{1 - p + p\eu^t} vanishes at 00, and φ(t)=ρ(1ρ)14\varphi''(t) = \rho(1 - \rho) \leq \frac14 for ρ=pet1p+pet[0,1]\rho = \frac{p\eu^t}{1 - p + p\eu^t} \in \intcc01: Taylor at order 22 gives φ(t)t28=λ2(ba)28\varphi(t) \leq \frac{t^2}8 = \frac{\lambda^2(b-a)^2}8.

18. For λ>0\lambda > 0, Markov applied to the positive variable eλ(SnESn)\eu^{\lambda(S_n - \E S_n)} (Proposition 22.3) and the product formula for independent variables give

P(SnESnt)eλti=1nEeλ(XiEXi)exp(λt+λ28i(biai)2),\P(S_n - \E S_n \geq t) \leq \eu^{-\lambda t}\prod_{i=1}^n\E\,\eu^{\lambda(X_i - \E X_i)} \leq \exp\Bigl(-\lambda t + \frac{\lambda^2}8\sum_i(b_i - a_i)^2\Bigr),

by question 17(b) applied to each centered XiEXi[aiEXi,biEXi]X_i - \E X_i \in \intcc{a_i - \E X_i}{b_i - \E X_i} (same width). Minimizing the exponent at λ=4tD\lambda = \frac{4t}{D}, D=i(biai)2D = \sum_i(b_i - a_i)^2, yields 2t2D-\frac{2t^2}D. The lower tail follows by applying the result to (Xi)(-X_i).

19. Take t=nεt = n\varepsilon and D=n(ba)2D = n(b - a)^2:

P(Snnmε)2exp(2n2ε2n(ba)2)=2exp(2nε2(ba)2),\P\Bigl(\Bigl|\frac{S_n}n - m\Bigr| \geq \varepsilon\Bigr) \leq 2\exp\Bigl(\frac{-2n^2\varepsilon^2}{n(b-a)^2}\Bigr) = 2\exp\Bigl(\frac{-2n\varepsilon^2}{(b-a)^2}\Bigr),

which is summable in nn (a geometric-type series): Borel–Cantelli (Theorem 22.8) gives that a.s. Snnm<ε\abs{\frac{S_n}n - m} < \varepsilon eventually; intersecting over ε=1p\varepsilon = \frac1p yields Snnm\frac{S_n}n \to m a.s. Comparison: Etemadi asks only X1L1X_1 \in L^1 and pairwise independence, and delivers no rate; Hoeffding asks boundedness and full independence, and delivers an explicit exponential guarantee at every finite nn — the two theorems answer different questions about the same limit.

20. The g(Uk)g(U_k) are i.i.d. with values in [0,1]\intcc01 and mean g ⁣dλd\int g\,\dd\lambda_d (transfer), so question 18 with biai=1b_i - a_i = 1, t=nεt = n\varepsilon gives the two-sided bound 2e2nε2δ2\eu^{-2n\varepsilon^2} \leq \delta as soon as e2nε22δ\eu^{2n\varepsilon^2} \geq \frac2\delta, i.e. nlog(2/δ)2ε2n \geq \frac{\log(2/\delta)}{2\varepsilon^2}. For ε=δ=102\varepsilon = \delta = 10^{-2}:

nlog2002104=5.29830.000226492:n \geq \frac{\log 200}{2\cdot10^{-4}} = \frac{5.2983\ldots}{0.0002} \approx 26\,492 :

about 2650026\,500 samples guarantee a 1%1\% accuracy with 99%99\% confidence — in every dimension dd, for every measurable integrand with values in [0,1]\intcc01. Question 11’s strong law promised convergence with no finite-nn guarantee; a deterministic grid with kk points per axis costs kdk^d evaluations, exponential in dd. Concentration is what makes Monte Carlo a method rather than a hope.

21. Independence and the product formula: EeλSn=(Eeλε1)n=(coshλ)nenλ2/2\E\,\eu^{\lambda S_n} = (\E\,\eu^{\lambda\varepsilon_1})^n = (\cosh\lambda)^n \leq \eu^{n\lambda^2/2} by question 17(a). Markov on eλSn\eu^{\lambda S_n}:

P(Snx)eλx+nλ2/2=ex2/(2n)at the optimum λ=xn,\P(S_n \geq x) \leq \eu^{-\lambda x + n\lambda^2/2} = \eu^{-x^2/(2n)} \qquad\text{at the optimum } \lambda = \frac xn,

and the symmetric bound for Sn-S_n (same law) doubles the constant for Sn\abs{S_n}.

22. Fix η>0\eta > 0 and set xn=(1+η)2nlognx_n = (1 + \eta)\sqrt{2n\log n} for n2n \geq 2:

P(Snxn)2exp((1+η)2logn)=2n(1+η)2,\P(\abs{S_n} \geq x_n) \leq 2\exp\bigl(-(1 + \eta)^2\log n\bigr) = \frac{2}{n^{(1+\eta)^2}},

summable since (1+η)2>1(1 + \eta)^2 > 1. Borel–Cantelli: a.s. Sn<(1+η)2nlogn\abs{S_n} < (1 + \eta)\sqrt{2n\log n} for all large nn, so lim supnSn2nlogn1+η\limsup_n\frac{\abs{S_n}}{\sqrt{2n\log n}} \leq 1 + \eta a.s.; intersecting the a.s. events for η=1p\eta = \frac1p, pNp \in \N^*, gives the claim. The walk of size nn has typical amplitude n\sqrt n (its variance), and even its worst excursions exceed that scale by at most 2logn\sqrt{2\log n}.

23. With nj=2jn_j = 2^j and x=(1+η)2njloglognjx = (1 + \eta)\sqrt{2n_j\log\log n_j} (defined for j2j \geq 2), question 21 gives

P(Snjx)exp((1+η)2loglognj)=(jlog2)(1+η)2,\P\bigl(S_{n_j} \geq x\bigr) \leq \exp\bigl(-(1 + \eta)^2\log\log n_j\bigr) = (j\log 2)^{-(1+\eta)^2},

summable in jj since (1+η)2>1(1 + \eta)^2 > 1: Borel–Cantelli and η=1p\eta = \frac1p give lim supjSnj/2njloglognj1\limsup_jS_{n_j}/\sqrt{2n_j \log\log n_j} \leq 1 a.s. What is missing for the full upper half is the bridge between checkpoints: one must show maxnjnnj+1Sn\max_{n_j\leq n\leq n_{j+1}}S_n exceeds (1+η)2njloglognj(1+\eta)\sqrt{2n_j\log\log n_j} only finitely often, which demands a maximal inequality with Gaussian tails (Lévy’s reflection inequality or Ottaviani’s inequality, not proved here). Question 12 is quantitatively too weak: it bounds the probability by

nj(1+η)22njloglognj=12(1+η)2log(jlog2),\frac{n_j}{(1+\eta)^2\,2n_j\log\log n_j} = \frac{1}{2(1+\eta)^2\log(j\log2)},

which tends to 00 but is not summable in jj: Borel–Cantelli cannot conclude. The lower half of the law of the iterated logarithm applies the second Borel–Cantelli lemma to the independent increments Snj+1SnjS_{n_{j+1}} - S_{n_j}, using matching lower bounds for Gaussian-type tails. Both refinements are genuine Year 3 probability, one course further along; what this problem delivers unaided is the exact iterated-logarithm scale along geometric times.

24. Each p^i\hat p_i is an average of nn i.i.d. indicator variables with values in [0,1]\intcc01 and mean P(Ai)\P(A_i): Hoeffding gives P(p^iP(Ai)>ε)2e2nε2\P(\abs{\hat p_i - \P(A_i)} > \varepsilon) \leq 2\eu^{-2n\varepsilon^2}. The union bound multiplies by NN. Solving 2Ne2nε2δ2N\eu^{-2n\varepsilon^2} \leq \delta: nln(2N/δ)2ε2n \geq \frac{\ln(2N/\delta)}{2\varepsilon^2}. Numerically: ln21060.05=ln(4107)17.5\ln\frac{2\cdot10^6}{0.05} = \ln(4\cdot10^7) \approx 17.5, so n17.5210487600n \geq \frac{17.5}{2\cdot10^{-4}} \approx 87\,600: estimating one probability to ±1%\pm1\% takes about 1850018\,500 samples (ln(2/δ)/2ε2\ln(2/\delta)/2\varepsilon^2), and a million probabilities only 4.7\approx 4.7 times more — uniformity costs lnN\ln N, not NN: the observation that makes empirical risk minimization, and with it machine learning, statistically possible.

25. The variables Xn=εnnαX_n = \frac{\varepsilon_n} {n^\alpha} are independent, centered, bounded, with nV(Xn)=nn2α\sum_n\V(X_n) = \sum_nn^{-2\alpha}. If α>12\alpha > \frac12: the variance series converges, and the one-series theorem (Part VI) gives a.s. convergence of Xn\sum X_n. If α12\alpha \leq \frac12: the variance series diverges, and the converse half (Part VI’s Paley–Zygmund argument, applicable since the summands are bounded by 11) gives a.s. divergence. Absolute convergence asks nα<\sum n^{-\alpha} < \infty: α>1\alpha > 1. On (12,1]\intoc{\frac12}1, the series converges a.s. although Xn=\sum\abs{X_n} = \infty surely: the signs conspire to cancel, with probability one — convergence by cancellation, invisible to any absolute test, and (by the zero–one law) with a deterministic verdict all the same.