Mathematics · Book 5 · Bachelor Year 3

Matemáticas universitarias — Grado 3

Matemáticas universitarias — Grado 3 · Bachelor Year 3

22Probabilidad: fundamentos y la ley de los grandes Números

El año 2 construyó probabilidad en espacios contables; medida teoría ahora elimina todas las restricciones. Un espacio de probabilidad es un medir el espacio de masa total 11, variables aleatorias son mensurable mapas, expectativa es la integral de Lebesgue — y a la vez todo el arsenal analítico (Capítulos 9, 10 y 11) se aplica a oportunidad. Este capítulo instala el diccionario, construye secuencias infinitas de independiente variables aleatorias (en [0,1]\intcc01, de dígitos binarios: la aleatoriedad se esconde en su interior medida de lebesgue), demuestra los lemas de Borel-Cantelli y ley cero uno de Kolmogorov, clasifica los modos de convergencia, y demuestra el ley de grandes números — el teorema que hace que las frecuencias converjan a probabilidades y estadísticas posibles. El problema del fin de semana da la prueba de Etemadi del fuerte ley en su forma definitiva L1L^1.

22.1 el diccionario

Definición 22.1

Un espacio de probabilidad es un medir el espacio (Ω,A,P)(\Omega, \mathcal A, \P) con P(Ω)=1\P(\Omega) = 1; Los elementos deA\mathcal A son eventos y una propiedad. tiene casi seguramente (a.s.) si su evento tiene probabilidad 11. Un variable aleatoria es una aplicación mensurable X ⁣:ΩRX \colon \Omega \to \R (o Rd\R^d: un vector aleatorio); su ley es la probabilidad de avance medida PX=XP\P_X = X_*\PenR\R (Ejercicio 11.9), determinada por el función de distribuciónFX(t)=P(Xt)F_X(t) = \P(X \leq t) (Ejercicio 9.3). XX tiene densidad ff si PX=f ⁣dλ\P_X = f\,\dd\lambda; es discreto si PX\P_X es un combinación contable de masas de Dirac. el expectativa es

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)),

y el teorema de transferencia (Ejercicio 11.9) lo calcula en la ley: E[g(X)]=Rg ⁣dPX\E[g(X)] = \int_\R g\,\dd\P_X=g(xk)pk= \sum g(x_k)p_k en el caso discreto, =g(x)f(x) ⁣dx= \int g(x)f(x)\dd x en el caso densidad: fórmulas del año 2, ahora teoremas de una teoría. El diferencia es V(X)=E[(XEX)2]=E[X2](EX)2\V(X) = \E[(X - \E X)^2] = \E[X^2] - (\E X)^2paraXL2X \in L^2.

Ejemplo 22.2

El estándar leyes y sus transformaciones destacadas: Bernoulli B(p)\mathcal B(p), binomial B(n,p)\mathcal B(n, p), geométrico, Poisson P(λ)\mathcal P(\lambda) (discreto: tablas del año 2 siguen siendo válidos); uniforme en [0,1]\intcc01 (medida de lebesgue mismo); exponencial E(λ)\mathcal E(\lambda) (densidad λeλx1x>0\lambda\eu^{-\lambda x}\mathbf 1_{x>0}); el gaussianoN(m,σ2)\mathcal N(m, \sigma^2) con densidad 1σ2πexp((xm)22σ2)\frac1{\sigma\sqrt{2\pi}}\exp\bigl(-\frac{(x - m)^2}{2\sigma^2}\bigr) — una probabilidad densidad por Problema 10.1, con media mm y varianza σ2\sigma^2 (momentos gaussianos, Ejercicio 11.10).

Proposición 22.3 (Markov y Chebyshev)

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

Demostración. Ejercicio 10.5(a); Chebyshev es el aplicado a Markov (XEX)2(X - \E X)^2.

22.2 Independencia

Definición 22.4

Las sub-σ\sigma-álgebras A1,,AnA\mathcal A_1, \dots, \mathcal A_n \subseteq \mathcal A son independiente si P(A1An)=P(Ai)\P(A_1\cap\dots\cap A_n) = \prod\P(A_i) para todos los AiAiA_i \in \mathcal A_i; los eventos son independientes si las álgebrasσ\sigma {,Ai,Aic,Ω}\{\varnothing, A_i, A_i^c, \Omega\} son; variables aleatorias X1,,XnX_1, \dots, X_n si las σ\sigma-álgebras σ(Xi)=Xi1(B(R))\sigma(X_i) = X_i^{-1}(\mathcal B(\R)) son. Una familia infinita es independiente si cada subfamilia finita lo es.

Teorema 22.5

X1,,XnX_1, \dots, X_n son independiente si el ley del vector (X1,,Xn)(X_1, \dots, X_n) es el medida del producto PX1PXn\P_{X_1}\otimes\cdots\otimes\P_{X_n}. En ese caso, para gi0g_i \geq 0 (o tal que los productos sean integrable):

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

en particular E[XY]=EXEY\E[XY] = \E X\,\E Y y V(X1++Xn)=V(Xi)\V(X_1 + \dots + X_n) = \sum\V(X_i)para las variables independienteL2L^2.

Demostración. Si XiX_i es independiente, las dos probabilidades medidas P(X1,,Xn)\P_{(X_1,\dots,X_n)} y PXi\bigotimes\P_{X_i} están de acuerdo en todo productos B1××BnB_1\times\dots\times B_n de conjuntos Borel — a π\pi-sistema generando B(Rn)\mathcal B(\R^n) (Proposición 11.2(b)) — por lo tanto en todas partes (Teorema 9.7). Por el contrario, un producto ley factoriza todos los eventos iXi1(Bi)\bigcap_iX_i^{-1}(B_i): independencia. La fórmula expectativa es entonces Tonelli/Fubini (Teorema 11.5) a través del teorema de transferencia; E[XY]=EXEY\E[XY] = \E X\E Y es el caso gi=idg_i = \mathrm{id}, y al expandir el cuadrado se obtiene la aditividad de variaciones (términos cruzados E[(XiEXi)(XjEXj)]=0\E[(X_i - \E X_i)(X_j - \E X_j)] = 0).

Teorema 22.6 (Existencia de secuencias independientes.)

En ([0,1],L,λ)\bigl(\intcc01, \mathcal L, \lambda\bigr) existe un secuencia (Un)n1(U_n)_{n\geq1} de independiente variables aleatorias, cada uniforme en [0,1]\intcc01. En consecuencia, para cualquier prescripción leyes (μn)(\mu_n) en R\R existe independiente (Xn)(X_n) con PXn=μn\P_{X_n} = \mu_n.

Demostración. Dígitos. Para ω[0,1]\omega \in \intcc01, sea (bk(ω))(b_k(\omega)) ser sus dígitos binarios (ω=bk2k\omega = \sum b_k2^{-k}; elija el la expansión no termina en todos los 11 — la ambigüedad solo se refiere un conjunto contable, por lo tanto nulo). Cada bkb_k es un aleatorio variable ({bk=1}\{b_k = 1\} es una unión finita de diádicos intervalos) y el vector (b1,,bm)(b_1, \dots, b_m) toma cada valor en {0,1}m\{0,1\}^m en un intervalo diádico de longitud 2m2^{-m}: los bkb_k son independiente Bernoulli(12)(\frac12).

Reagrupación. Dividir N\N^* en infinitos desunidos conjuntos infinitos (In)(I_n) (por ejemplo, por potencias primas o diagonales); deje que (kjn)j(k^n_j)_j enumere InI_n y establezca

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

Cada UnU_n es uniforme: sus dígitos binarios son independiente justos bits, entonces P(Un[l2m,(l+1)2m))=2m\P(U_n \in [l2^{-m}, (l+1)2^{-m})) = 2^{-m} para cada intervalo diádico, y los intervalos diádicos determinan el ley (Teorema 9.7). Los UnU_n son independiente: son funciones de bloques disjuntos del independiente familia (bk)(b_k) — formalmente, eventos {UnDn}\{U_n \in D_n\} para diádico DnD_n depende de un número finito de dígitos disjuntos conjuntos y factorizar; el argumento del sistema π\pi se actualiza a todos los conjuntos de Borel.

Arbitrary leyes. Dejemos que Gn(u)=inf{t:Fμn(t)u}G_n(u) = \inf\{t : F_{\mu_n}(t) \geq u\} (el función cuantil de la distribución función FμnF_{\mu_n}); la equivalencia de clave Gn(u)t    uFμn(t)G_n(u) \leq t \iff u \leq F_{\mu_n}(t)(continuidad derecha deFF, monotonicidad) muestra Xn=Gn(Un)X_n = G_n(U_n) es mensurable con 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): leyμn\mu_n; independencia se hereda (funciones de variables independiente, Ejercicio 22.3).

Ejemplo 22.7 (El problema del cumpleaños, sinceramente.)

Entre nn personas con independiente, cumpleaños uniformes terminados N=365N = 365 días, la probabilidad de que todos los cumpleaños difieran es

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

por condicionamiento iterado (o directamente: el favorable N(N1)(Nn+1)N(N-1)\cdots(N - n + 1) sobre el total NnN^n, un conteo argumento que la fórmula del producto independencia hace riguroso). Tomando logaritmos y usando 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} .

El punto de inflexión pn=12p_n = \frac12 se sitúa en n2Nln21.18Nn \approx \sqrt{2N\ln2} \approx 1.18\sqrt N: paraN=365N = 365,n=23n = 23 (p23=0.4927p_{23} = 0.4927). Dos moralejas. Primero, colisiones entre nn Los elementos en los cuadros NN aparecen en la escala nNn \sim \sqrt N, no nNn \sim N — el escala de cumpleaños que gobierna el hash colisiones y el costo N\sqrt N de los ataques de cumpleaños en criptografía. En segundo lugar, el cálculo es una plantilla: el (n2)\binom n2 los eventos de colisión de pares aún no son independiente la respuesta se comporta como si lo fueran (e(n2)/N\eu^{-\binom n2/N} son exactamente los pares independientes heurística) — una primera instancia de la teoría de Poisson aproximación hecha rigurosa en el fin de semana de Capítulo 23 problema (desigualdad de Le Cam).

22.3 Borel–Cantelli y la ley cero-uno

Teorema 22.8 (Borel–Cantelli)

Sea (An)(A_n) eventos y lim supAn=NnNAn\limsup A_n = \bigcap_N \bigcup_{n\geq N}A_n(“AnA_n ocurre infinitamente a menudo”).

  1. Si P(An)<\sum\P(A_n) < \infty, entonces P(lim supAn)=0\P(\limsup A_n) = 0.
  2. Si P(An)=\sum\P(A_n) = \infty and the AnA_n are independiente, entonces P(lim supAn)=1\P(\limsup A_n) = 1.

Demostración. (1) es Ejercicio 9.4. (2): para NMN \leq M, independencia de complementos (Ejercicio 22.3) da

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}; la serie diverge). entonces P(nNAn)=1\P\bigl(\bigcup_{n\geq N}A_n\bigr) = 1 por cada NN, y la intersección decreciente sobre NN todavía tiene probabilidad 11 (continuidad desde arriba, Proposición 9.6).

Teorema 22.9 (Ley cero-uno de Kolmogorov)

Sea (Xn)(X_n) independiente y T=Nσ(XN,XN+1,)\mathcal T = \bigcap_N\sigma(X_N, X_{N+1}, \dots) tail σ\sigma-algebra (eventos insensible a cualquier número finito del XnX_n: convergencia de Xn\sum X_n, de Snn\frac{S_n}n, valores de lim sup\limsup’s, …). Entonces cada TTT \in \mathcal T tiene P(T){0,1}\P(T) \in \{0, 1\}.

Demostración. Reparar NN. Las álgebras σ\sigma σ(X1,,XN)\sigma(X_1, \dots, X_N) y σ(XN+1,)\sigma(X_{N+1}, \dots) son independiente: eventos dependiendo de los bloques disjuntos se factorizan en los sistemas generadores π\pi (cilindros iN{XiBi}\bigcap_{i\leq N}\{X_i \in B_i\}, resp. finito condiciones en variables posteriores), y Dynkin (Teorema 9.4, aplicado dos veces, un lado a la vez tiempo) extiende la factorización. Un evento de cola TT se encuentra en σ(XN+1,)\sigma(X_{N+1}, \dots) por cada NN: TT es independiente de cada σ(X1,,XN)\sigma(X_1, \dots, X_N), por lo tanto del álgebra σ\sigma generan, σ(X1,X2,)\sigma(X_1, X_2, \dots) (Dynkin una vez más: la unión del σ(X1,,XN)\sigma(X_1,\dots,X_N) es un sistema π\pi generarlo). Pero Tσ(X1,X2,)T \in \sigma(X_1, X_2, \dots) también: TT es independiente de si mismo, 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 Modos de convergencia

Definición 22.10

XnXX_n \to X casi seguramente si P(XnX)=1\P(X_n \to X) = 1; en probabilidad si P(XnXε)0\P(\abs{X_n - X} \geq \varepsilon) \to 0por cadaε>0\varepsilon > 0; inLpL^p si EXnXp0\E\abs{X_n - X}^p \to 0.

Proposición 22.11

(a) la convergencia a.s. implica convergencia en probabilidad; (b) La convergencia LpL^p implica convergencia en probabilidad; (c) la convergencia en probabilidad implica a.s. convergencia a lo largo de una subsecuencia; (d) ninguna otra implicación es válida en general.

Demostración. (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 bajo a.s. convergencia (continuidad desde arriba; el evento limsup excluye la convergencia). (b) Markov: P(XnXε)εpEXnXp\P(\abs{X_n - X} \geq \varepsilon) \leq \varepsilon^{-p}\,\E\abs{X_n - X}^p. (c) Elijankn_k con P(XnkX2k)2k\P(\abs{X_{n_k} - X} \geq 2^{-k}) \leq 2^{-k}; Borel–Cantelli (1) marca XnkX<2k\abs{X_{n_k} - X} < 2^{-k} finalmente, como. (d) La máquina de escribir (Ejercicio 12.3) en ([0,1],λ)(\intcc01, \lambda) converge en L1L^1 y en probabilidad pero en ninguna parte puntualmente; n1(0,1/n)0n\mathbf 1_{\intoo0{1/n}} \to 0 a.s. pero no en L1L^1; detalles y el resto contraejemplos en Ejercicio 22.6.

22.5 La ley de los grandes números.

En todo momento, (Xn)(X_n) son independiente con el mismo ley. (i.i.d.), Sn=X1++XnS_n = X_1 + \dots + X_n.

Teorema 22.12 (Ley débil de los grandes números)

Si X1L2X_1 \in L^2, con 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 en probabilidad (y en L2L^2).

Demostración. ESnn=m\E\frac{S_n}n = m y V(Snn)=nV(X1)n2\V\bigl(\frac{S_n}n\bigr) = \frac{n\V(X_1)}{n^2} (Teorema 22.5); Chebyshev.

Teorema 22.13 (Ley fuerte de los grandes números.)

Si X1L1X_1 \in L^1, entonces

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

Lo demostramos aquí bajo la hipótesis más fuerte X1L4X_1 \in L^4; el caso general (L1L^1: prueba de Etemadi) es el problema del fin de semana.

Prueba bajo EX14<\E X_1^4 < \infty. Centrado (XiXimX_i \mapsto X_i - m), suponga m=0m = 0. Ampliar:

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 ,

desde independencia y centrando eliminan todos los términos que contienen un factor aislado (E[XiXjXkXl]=E[Xi]E[]=0\E[X_iX_jX_kX_l] = \E[X_i]\E[\cdots] = 0 a menos que los índices se emparejen: el único Los supervivientes son los términos nn i=j=k=li=j=k=l y 3n(n1)3n(n-1). términos con dos pares distintos). 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},

sumable: Borel-Cantelli (1) da, para cada racional ε\varepsilon, que Sn/n<ε\abs{S_n/n} < \varepsilon eventualmente, como; intersección sobre εQ+\varepsilon \in \Q_+^* (contablemente muchos eventos de probabilidad-11): Sn/n0S_n/n \to 0 a.s.

Ejemplo 22.14 (Lo que compra la ley fuerte)

(a) Frecuencias: para i.i.d. lanzamientos de moneda, el observado frecuencia de cabezas converge a.s. a pp — el empírico justificación de la probabilidad misma. (b) Montecarlo: para gL1([0,1])g \in L^1(\intcc01)y(Un)(U_n) i.i.d. uniforme (Teorema 22.6), 1nkng(Uk)01g\frac1n\sum_{k\leq n}g(U_k) \to \int_0^1g a.s.: integrales por muestreo, en cualquier dimensión, en el independiente de la dimensión tarifa n1/2\sim n^{-1/2} precisada en Capítulo 23. (c) Números normales: casi todo número real tiene, en su expansión binaria, frecuencia asintótica 12\frac12 de unos (Aplique el fuerte ley a las variables de dígitos de Teorema 22.6) — Teorema de Borel, a declaración sobre los números de cadaday probados por medida: Problema 22.1 completa en todas las bases.

Método 22.15

El orden de funcionamiento de enunciados asintóticos sobre el azar. secuencias: (1) ¿Es el evento un evento de cola? Entonces es la probabilidad es 00 o 11 (Teorema 22.9) y solo hay que decidir cual. (2) To prove a.s. statements: Borel–Cantelli — probabilidades sumables para el Eventos "malos", a través de límites tipo Markov/Chebyshev en cualesquiera momentos que existan; independencia sólo es necesario para el dirección inversa. (3) Subsecuencia + sándwich: probar convergencia a lo largo de una subsecuencia manejable, controlar la oscilación intermedia por monotonicidad o máxima Desigualdades: el esqueleto de la prueba de Etemadi. (4) Para límites de distribución, espere Capítulo 23.

22.6 Ceremonias

Ejercicio 22.1

(a) Sea XX tener continuo distribución estrictamente creciente función FF. Demuestre que F(X)F(X) es uniforme en [0,1]\intcc01, y que G(U)FG(U) \sim F para UU uniforme, G=F1G = F^{-1}: simulación por inversión. (b) Calcule la función de distribución y densidad de X2X^2 para XX uniforme en [1,1]\intcc{-1}1, y de 1λlnU-\frac1\lambda\ln UparaUUuniforme en(0,1)\intoo01.

Solución

Solución de Ejercicio 22.1.

(a) Para 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 (continuidad y estricta La monotonicidad hace de FF una biyección sobre (0,1)\intoo01 con {F(X)u}={XF1(u)}\{F(X) \leq u\} = \{X \leq F^{-1}(u)\}): F(X)F(X) es uniforme. Por el contrario P(G(U)t)=P(UF(t))=F(t)\P(G(U) \leq t) = \P(U \leq F(t)) = F(t): a simular un ley, aplicar la función de distribución inversa a un muestra uniforme.

(b) Y=X2Y = X^2, XX uniforme en [1,1]\intcc{-1}1: para t[0,1]t \in \intcc01,FY(t)=P(tXt)=tF_Y(t) = \P(-\sqrt t \leq X \leq \sqrt t) = \sqrt t: densidad12t1(0,1)\frac1{2\sqrt t}\mathbf 1_{\intoo01}. y 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}: el exponencial E(λ)\mathcal E(\lambda) — inversión en acción.

Ejercicio 22.2

(a) Calcule la media y la varianza de Poisson P(λ)\mathcal P(\lambda) y geométrico leyes mediante el teorema de transferencia. (b) Demuestre que un variable aleatoria TT positivo con P(T>t)>0\P(T > t) > 0para todos losttsatisface la propiedad sin memoriaP(T>t+sT>t)=P(T>s)\P(T > t + s \mid T > t) = \P(T > s)para todos loss,t0s, t \geq 0siTT es exponencial. (The survival function satisfies Cauchy’s functional equation; monotonicity replaces continuidad.)

Solución

Solución de Ejercicio 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, entonces V=λ2+λλ2=λ\V = \lambda^2 + \lambda - \lambda^2 = \lambda. Geométrico (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} (diferenciar la serie geométrica dos veces).

(b) G(t)=P(T>t)G(t) = \P(T > t) no aumenta con G(0+)G(0^+)\dots G ⁣:[0,)(0,1]G \colon \intco0\infty \to \intoc01; la falta de memoria lee G(t+s)=G(t)G(s)G(t + s) = G(t)G(s). Entonces G(nt)=G(t)nG(n t) = G(t)^n y G(t/n)=G(t)1/nG(t/n) = G(t)^{1/n}:G(q)=G(1)qG(q) = G(1)^qparaq0q \geq 0 racional; escribiendo G(1)=eλG(1) = \eu^{-\lambda} ((0,1)\in \intoo01: G(1)=1G(1) = 1 sería fuerza G1G \equiv 1, imposible para un variable aleatoria finito; G(1)=0G(1) = 0 está excluido por hipótesis) y apretando un tt arbitrario entre racionales (monotonicidad): G(t)=eλtG(t) = \eu^{-\lambda t} — el exponencial ley. Lo contrario es un cálculo.

Ejercicio 22.3 ★★

(a) Demuestre que si X1,,XnX_1, \dots, X_n son independiente y fif_i son funciones de Borel, las fi(Xi)f_i(X_i) son independiente. (b) Demuestre que los eventos A1,,AnA_1, \dots, A_n son independiente si y así sus complementos son, si los indicadores 1Ai\mathbf 1_{A_i} son independiente variables aleatorias. (c) (Por parejas es más débil) Dos monedas justas: A=A = la primera es cara, B=B = el segundo es cara, C=C = los dos están de acuerdo. Mostrar A,B,CA, B, C son independiente en pares pero no independiente.

Solución

Solución de Ejercicio 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), y sub-σ\sigma-álgebras de independiente σ\sigma-álgebras son independiente (la identidad definitoria es válida 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}): los tres declaraciones afirman independencia del mismo σ\sigma-álgebras. (Esa factorización sobre el AiA_i se propaga a complementos es el argumento del sistema λ\lambda dentro de la equivalencia de Definición 22.4 — o inclusión-exclusión directa.)

(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 en pares: cada intersección es “ambas caras” o análogo, de probabilidad 14\frac14: por pares independiente. Pero P(ABC)=P(HH)=1418\P(A\cap B\cap C) = \P(\text{HH}) = \frac14 \neq \frac18: no independienteCCestá determinado porAA y BB.

Ejercicio 22.4 ★★

(a) (Mono infinito) Una secuencia i.i.d. de uniforme pulsaciones de teclas en un alfabeto finito como contiene cada finito texto infinitamente a menudo: pruébalo con Borel–Cantelli (2) en bloques disjuntos. (b) (Ejecuciones) Para i.i.d. bits justos, sea RnR_n la longitud de la serie de unos que comienzan en la posición nn. Muestra ese culo Rn(1+ε)log2nR_n \geq (1+\varepsilon)\log_2n con frecuencia finita y Rnlog2nR_n \geq \log_2 n con frecuencia infinita (both halves of Borel–Cantelli; for the second, pass to disjoint blocks to gain independencia): la ejecución más larga en los primeros dígitos nn crece como log2n\log_2n.

Solución

Solución de Ejercicio 22.4.

(a) Deje que el texto TT tenga una longitud LL y q=aLq = a^{-L} (aa el tamaño del alfabeto). Los eventos Ek={E_k = \{ posiciones kL+1,,(k+1)LkL+1, \dots, (k+1)LhechizoT}T\} son independiente (bloques disjuntos de letras i.i.d., cada una de probabilidad q>0q > 0: P(Ek)=\sum\P(E_k) = \infty, y Borel–Cantelli (2) da infinitas ocurrencias como

(b) Superior: 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)}, sumable: por Borel–Cantelli (1), a.s. sólo un número finito tal nn. Inferior: empaquetar bloques separados — el jj-ésimo de longitud j=log2sj\ell_j = \lceil\log_2s_j\rceil a partir de sj=i<jis_j = \sum_{i<j}\ell_i; los eventos “el bloquejj son todos unos” son independiente con probabilidad 2j1sj1jlog2j2^{-\ell_j} \asymp \frac1{s_j} \asymp \frac1{j\log_2 j}, cuya suma diverge: Borel–Cantelli (2) da infinitos bloques todos unos, es decir, Rsjlog2sjR_{s_j} \geq \log_2 s_j infinitamente frecuente. Juntos: la longitud máxima de ejecución en los primeros dígitos nn es (1+o(1))log2n(1 + o(1))\log_2n como

Ejercicio 22.5 ★★

Sea (Xn)(X_n) independiente. (a) Demuestre que el radio de convergencia de Xnzn\sum X_n z^n es una constante a.s. (posiblemente 00 o \infty). (b) Demuestre que P(Xn converges){0,1}\P(\sum X_n \text{ converges}) \in \{0, 1\} y P(Sn/nm){0,1}\P(S_n/n \to m) \in \{0,1\}. (c) Darle una cola a un evento sobre (Xn)(X_n) que sea no evento y verifique que el ley cero uno pueda fallar.

Solución

Solución de Ejercicio 22.5.

(a) R=(lim supXn1/n)1R = \bigl(\limsup\abs{X_n}^{1/n}\bigr)^{-1} es sin cambios si se modifican un número finito de XnX_n: por cada NN, RR es medible σ(XN,XN+1,)\sigma(X_N, X_{N+1}, \dots), es decir, cola medible. Entonces cada evento {Rc}\{R \leq c\} tiene probabilidad 00 o 11 (Teorema 22.9), entonces la función de distribución de RR toma solo los valores 0,10, 1: salta en un solo punto c0[0,+]c_0 \in \intcc0{+\infty}, yR=c0R = c_0 a.s.

(b) La convergencia de Xn\sum X_n y de Snn\frac{S_n}n son insensible a cambiar un número finito de términos (para el segundo: los términos modificados contribuyen O(1/n)0O(1/n) \to 0): eventos de cola; ley cero uno.

(c) {X1>0}\{X_1 > 0\} depende de X1X_1: para signos i.i.d. (P(X1=±1)=12\P(X_1 = \pm1) = \frac12), su probabilidad es 12{0,1}\frac12 \notin \{0,1\} — no hay contradicción, no es cola evento.

Ejercicio 22.6 ★★

En ([0,1],λ)(\intcc01, \lambda), presente — con pruebas — aleatorio variables tal que: (a) Xn0X_n \to 0 en probabilidad y en cada LpL^p, pero en ninguna parte como; (b) Xn0X_n \to 0 a.s. pero en el número LpL^p; (c) Xn0X_n \to 0 en L1L^1 pero no en L2L^2; (d) y muestre: si XnXX_n \to X en probabilidad y XnYL1\abs{X_n} \leq Y \in L^1, entoncesXnXX_n \to XenL1L^1 (subsecuencias + convergencia dominada + el truco de la subsecuencia).

Solución

Solución de Ejercicio 22.6.

Trabajar en ([0,1],λ)(\intcc01, \lambda). (a) La máquina de escribir 1In\mathbf 1_{I_n} (Ejercicio 12.3): Xnpp=λ(In)0\norm{X_n}_p^p = \lambda(I_n) \to 0 (todos p<p < \infty), por lo tanto también en probabilidad; en cada ω\omega los valores 00 y 11 se repiten: no puntualmente convergencia en cualquier lugar. (b) Xn=n1(0,1/n)0X_n = n\mathbf 1_{\intoo0{1/n}} \to 0 de 00, pero 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) De cualquier extracto de subsecuencia (convergencia en probabilidad) otra subsecuencia convergente a.s. (Proposición 22.11(c)); convergencia dominada da convergencia L1L^1 a lo largo de él, con el límite mismo XX. Así, cada subsecuencia de la secuencia numérica EXnX\E\abs{X_n - X} tiene una subsecuencia que tiende a 00: la toda la secuencia tiende a 00.

Ejercicio 22.7 ★★

Una encuesta de opinión estima una proporción desconocida pp según el Frecuencia empírica p^n\hat p_n de sorteos nn independiente. (a) Chebyshev: muestre P(p^npε)14nε2\P(\abs{\hat p_n - p} \geq \varepsilon) \leq \frac1{4n\varepsilon^2}(usep(1p)14p(1-p) \leq \frac14). (b) ¿Cuántos sorteos garantizan un error 3%\leq 3\% con probabilidad 95%\geq 95\% por este límite? (La verdadera respuesta, vía Capítulo 23, se trata de 10701070: Chebyshev es honesto pero crudo.)

Solución

Solución de Ejercicio 22.7.

(a) p^n=Snn\hat p_n = \frac{S_n}n con el binomio SnS_n: V(p^n)=p(1p)n14n\V(\hat p_n) = \frac{p(1-p)}n \leq \frac1{4n} y Chebyshev (Proposición 22.3) da el límite. (b) Resuelva 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. el central El teorema del límite justificará n1070n \approx 1070 para el mismo garantía: Chebyshev paga su generalidad con un factor 5\approx 5.

Ejercicio 22.8 ★★★

(Bernstein) Para fC([0,1])f \in \mathcal C(\intcc01) defina el Polinomio de Bernstein 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) Reconozca Bnf(x)=E[f(Snn)]B_nf(x) = \E\bigl[f\bigl(\frac {S_n}n\bigr)\bigr]para el binomioSnS_n B(n,x)\mathcal B(n, x). (b) Demostrar BnffB_nf \to f uniformemente en [0,1]\intcc01: división en {Snnxδ}\{\abs{\frac{S_n}n - x} \leq \delta\} y su complemento, utilizando el uniforme continuidad y Chebyshev con el cota uniformemente V(Snn)14n\V(\frac{S_n}n) \leq \frac1{4n}. (c) Concluir: una segunda prueba probabilística de la Teorema de aproximación de Weierstrass (Corolario 7.16), con la tarifa explícita Bnff32ωf(n1/2)\norm{B_nf - f}_\infty \leq \frac32\,\omega_f(n^{-1/2}) para el módulo de continuidad ωf\omega_f — demostrar al menos el formulario O(ωf(n1/2))O(\omega_f(n^{-1/2})).

Solución

Solución de Ejercicio 22.8.

(a) Si SnB(n,x)S_n \sim \mathcal B(n, x), el teorema de transferencia da 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) Sea ω=ωf\omega = \omega_f el módulo de continuidad (f(u)f(v)ω(uv)\abs{f(u) - f(v)} \leq \omega(\abs{u - v}) y ω(cδ)(1+c)ω(δ)\omega(c \delta) \leq (1 + c)\,\omega(\delta) encadenando pasos). Luego, para cualquier δ>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)

(si uxδ\abs{u - x} \leq \delta, claro; en caso contrario ω(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)). Tome esperanzas de heredar enu=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) ;

con δ=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 (uniforme continuidad en compacto): un Weierstrass probabilístico teorema, con una tasa explícita y uniforme.

Ejercicio 22.9 ★★★

(Coleccionista de cupones) Las tarjetas del tipo nn se extraen de manera uniforme con reemplazo; sea TnT_n el número de sorteos hasta que todos Se ven tipos. (a) Escriba Tn=k=1nτkT_n = \sum_{k=1}^{n}\tau_k con τk\tau_k geométrico del parámetro nk+1n\frac{n - k + 1}n, el τk\tau_k independiente, y deducir ETn=nHnnlnn\E T_n = n\,H_n \sim n\ln n (HnH_n el número de armónico) y V(Tn)π26n2\V(T_n) \leq \frac{\pi^2}6n^2. (b) Chebyshev: Tnnlnn1\frac{T_n}{n\ln n} \to 1 en probabilidad. (c) Refinar con Borel–Cantelli: mostrar directamente P(Tn>βnlnn)n1β\P(T_n > \beta n\ln n) \leq n^{1 - \beta}paraβ>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}), y deducir eso a lo largo de n=2mn = 2^m, a.s. TnβnlnnT_n \leq \beta n\ln neventualmente, por cadaβ>2\beta > 2.

Solución

Solución de Ejercicio 22.9.

(a) Después de recopilar los tipos k1k - 1, cada sorteo es nuevo con probabilidad pk=nk+1np_k = \frac{n-k+1}n: τk\tau_k es geométrica (pk)(p_k) y τk\tau_k son independiente (los sorteos son). Sumas: 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 0ynHnnlnn1\frac{nH_n}{n\ln n} \to 1:Tnnlnn1\frac{T_n}{n\ln n} \to 1 en probabilidad.

(c) Con destino a la unión: Tn>tT_n > t significa que algún tipo no se ve después t\lceil t\rceil empata, por lo que P(Tn>t)n(11n)tnet/n\P(T_n > t) \leq n(1 - \frac1n)^{t} \leq n\,\eu^{-t/n}; ent=βnlnnt = \beta n\ln n: n1β\leq n^{1 - \beta}. Para β>1\beta > 1, m2m(1β)<\sum_m 2^{m(1-\beta)} < \infty: Borel–Cantelli da, junto conn=2mn = 2^m, a.s.TnβnlnnT_n \leq \beta n\ln n eventualmente — en particular para cada β>2\beta > 2 como se indica (cualquier β>1\beta > 1 funciona a lo largo de la subsecuencia).

Ejercicio 22.10 ★★

Usando la construcción de dígitos (Teorema 22.6): (a) verificar mediante cálculo directo que U=b2k2kU = \sum b_{2k}2^{-k} (dígitos pares indexados de un ω\omega uniforme) es uniforme y independiente de V=b2k12kV = \sum b_{2k-1}2^{-k}; (b) deducir una biyección mensurable hasta conjuntos nulos entre [0,1]\intcc01 y [0,1]2\intcc01^2 conservando medida y comenta: un número aleatorio uniforme contiene dos (y muchos contables) independiente unos — comparar con la curva de Peano (Problema 6.1), que logró sobreyectividad pero no preservación de la medida o inyectividad.

Solución

Solución de Ejercicio 22.10.

(a) Los dígitos pares (b2k)k(b_{2k})_k son i.i.d. justos bits (una subfamilia de la familia de dígitos independiente), por lo que U=kb2k2kU = \sum_kb_{2k}2^{-k} le da a cada intervalo diádico su valor correcto probabilidad (como en Teorema 22.6): uniforme; igualmente VV; y (U,V)(U, V) dependen de bloques de dígitos separados: independiente (factorización en rectángulos diádicos, luego Dynkin).

(b) Φ(ω)=(U(ω),V(ω))\Phi(\omega) = (U(\omega), V(\omega)) es mensurable con Φλ=λλ=λ2\Phi_*\lambda = \lambda\otimes\lambda = \lambda_2 (acuerdo sobre rectángulos diádicos + unicidad). Intercalado dígitos define un inverso definido a partir del conjunto (nulo) de racionales diádicos en cualquier factor: una medida que preserva biyección entre medida completa subconjuntos de [0,1]\intcc01 y [0,1]2\intcc01^2. Contraste con Peano (Problema 6.1): continuidad sobreyectividad forzada sin inyectividad; dejando caer continuidad por mera mensurabilidad compra un isomorfismo de medida — la dimensión es invisible para la teoría medida, visible para topología.

Ejercicio 22.11 ★★

(Registros) Sea (Xn)n1(X_n)_{n\geq1} i.i.d. con continuo función de distribución, y digamos que ocurre un registro en tiempo nn si Xn>max(X1,,Xn1)X_n > \max(X_1, \dots, X_{n-1}) (el tiempo 11 es un registro). Sea RnR_n el indicador de registro. (a) Mostrar 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) Demuestre que los RnR_n son independiente (count orderings compatible with prescribed record positions, or argue that the relative order of X1,,Xn1X_1, \dots, X_{n-1} is independiente of the rank of XnX_n among them). (c) Deducir de Borel–Cantelli (Teorema 22.8, ambas mitades) que infinitos registros ocurren como, pero los registros en veces consecutivas n,n+1n, n+1 ocurren infinitamente con probabilidad — ¡decide cuál! — y calcular nP(Rn=1,Rn+1=1)\sum_n\P(R_n = 1, R_{n+1} = 1).

Solución

Solución de Ejercicio 22.11.

(a) Continuidad de la distribución hace que los empates sean eventos nulos (como en los argumentos de las estadísticas de orden del capítulo), y el Los pedidos relativos de n!n! de (X1,,Xn)(X_1, \dots, X_n) son intercambiables y, por tanto, igualmente probables. Rn=1R_n = 1 significa el el máximo se sitúa en la última posición: probabilidad (n1)!n!=1n\frac{(n-1)!}{n!} = \frac1n.

(b) Corrija nn y condición sobre el orden relativo de X1,,Xn1X_1, \dots, X_{n-1}: insertarXnX_nen el rango posiblenn las ranuras son uniformes y independiente de ese orden (intercambiabilidad de la tupla nn). Por lo tanto RnR_n (el evento “XnX_n ocupa el primer lugar”) es independiente del total registrar el historial (R1,,Rn1)(R_1, \dots, R_{n-1}), que es una función del orden relativo de las primeras variables n1n - 1. La inducción proporciona independencia completo con P(Rn=1)=1n\P(R_n = 1) = \frac1n.

(c) P(Rn=1)=1n=\sum\P(R_n = 1) = \sum\frac1n = \infty con independencia: la segunda parte Borel–Cantelli da récords infinitamente a menudo a.s. (los registros nunca se detienen — pero adelgazar logarítmicamente: E[#recordsn]=Hnlnn\E[\#\text{records} \leq n] = H_n \approx \ln n). Registros consecutivos:P(Rn=Rn+1=1)=1n(n+1)\P(R_n = R_{n+1} = 1) = \frac1{n(n+1)} (independencia), y

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

la primera mitad Borel-Cantelli se aplica — sólo de forma limitada Se producen muchos pares de registros consecutivos, como.

Ejercicio 22.12 ★★

(Carrera más larga) Lanza una moneda justa infinitas veces y sea LnL_n la longitud de la ejecución más larga de secuencias consecutivas cabezas dentro de los primeros giros nn. (a) Demuestre que para cada ε>0\varepsilon > 0, a.s. Ln(1+ε)log2nL_n \leq (1 + \varepsilon)\log_2n eventualmente (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) Demuestre que a.s. Ln(1ε)log2nL_n \geq (1 - \varepsilon)\log_2n eventualmente (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 independiente, 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) Concluir Lnlog2n1\frac{L_n}{\log_2n} \to 1 a.s.: en un millón En los lanzamientos justos se debe esperar una serie de aproximadamente 2020 cabezas — y probablemente se haya fabricado un conjunto de datos sin uno.

Solución

Solución de Ejercicio 22.12.

(a) Un tramo de longitud \ell que comienza en la posición ini \leq n tiene probabilidad 22^{-\ell}; vinculado a la unión: P(Ln)n2\P(L_n \geq \ell) \leq n2^{-\ell}. Conn=(1+ε)log2n\ell_n = (1 + \varepsilon)\log_2n:P(Lnn)nε\P(L_n \geq \ell_n) \leq n^{-\varepsilon}. A lo largo den=2kn = 2^k: k2kε<\sum_k2^{-k\varepsilon} < \infty, entonces a.s. L2k<(1+ε)kL_{2^k} < (1+\varepsilon)k eventualmente (Borel–Cantelli); para generales nn elige 2k1<n2k2^{k-1} < n \leq 2^k y usa la monotonicidad de LnL_n más 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, y el factor extra se absorbe ampliando ε\varepsilon ligeramente.

(b) Con =(1ε)log2n\ell = \lceil(1 - \varepsilon)\log_2n\rceil y m=n/m = \lfloor n/\ell\rfloor bloques disjuntos: los bloques son independiente, cada cara con probabilidad 2n(1ε)/22^{-\ell} \geq n^{-(1-\varepsilon)}/2, por lo que

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)

para una constante c>0c > 0 y nn grande. Estas probabilidades son sumables a lo largo de n=2kn = 2^k (de hecho, a lo largo de todo nn): Borel–Cantelli finalmente da a.s. Ln(1ε)log2nL_n \geq (1 - \varepsilon)\log_2n (la monotonicidad completa entre el 2k2^k como en (a), sin causar daño).

(c) Ambos límites a lo largo de una secuencia ε=1j\varepsilon = \frac1j, intersectando contablemente muchos eventos medida completa: Lnlog2n1\frac{L_n}{\log_2n} \to 1 como Para n=106n = 10^6: log2n19.9\log_2n \approx 19.9— una serie de cabezales20\approx 20 no es un anomalía sospechosa sino una certeza matemática, y su La ausencia es evidencia de un ser humano fingiendo "aleatoriedad" (los humanos rara vez se atreve a escribir más de 55 o 66 cabezas seguidas).

22.7 Problema: la prueba de la ley fuerte de Etemadi

Problema 22.1

Problema del fin de semana — la ley fuerte de los grandes números para i.i.d. variables integrables

El fuerte ley de Kolmogorov — SnnEX1\frac{S_n}n \to \E X_1 a.s. para i.i.d. XnL1X_n \in L^1 — durante mucho tiempo sólo tuvo pruebas intrincadas; en 1981 N. Etemadi encontró uno de sorprendente economía, utilizando nada más allá de este capítulo (e incluso debilitando independencia a pares independencia). Lo seguimos. Sea (Xn)(X_n) independiente por pares, distribuido idénticamente, integrable; m=EX1m = \E X_1, Sn=X1++XnS_n = X_1 + \dots + X_n.

Parte I — Reductions.

  1. Demostrar que es suficiente tratar Xn0X_n \geq 0 (split Xn=Xn+XnX_n = X_n^+ - X_n^-: check the two halves are again pairwise independiente i.i.d. integrable). Supongamos en adelante Xn0X_n \geq 0.
  2. (truncamiento) Deje Yn=Xn1XnnY_n = X_n\,\mathbf 1_{X_n \leq n} y Sn=Y1++YnS_n^* = Y_1 + \dots + Y_n. Mostrar

    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

    (Ejercicio 11.3), y deducir vía Borel–Cantelli que SnSnn0\frac{S_n - S_n^*}{n} \to 0 a.s.: basta con probar Snnm\frac{S^*_n}n \to m como

  3. Mostrar EYn=E[X11X1n]m\E Y_n = \E\bigl[X_1\mathbf 1_{X_1\leq n}\bigr] \to m (convergencia monótona), por lo tanto 1nknEYkm\frac1n\sum_{k\leq n}\E Y_k \to m (Cesàro): es basta con probar SnESnn0\frac{S_n^* - \E S_n^*}{n} \to 0 como

Parte II — The variance estimate.

  1. Mostrar

    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]

    y, usando el pastel de capas (Proposición 11.8), la llave encuadernada

    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 expectativa — 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).

Parte III — Convergence along geometric subsequences. Repare α>1\alpha > 1 y deje que kj=αjk_j = \lfloor\alpha^j\rfloor.

  1. Usando independencia por pares (las variaciones se suman, Teorema 22.5 — comprueba que la aditividad de las varianzas solo necesita por pares independencia) y Chebyshev, espectáculo para todos ε>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. Mostrar 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) y concluir con pregunta 4 y 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.}

Parte IV — Sandwich and conclusion.

  1. Para kjnkj+1k_j \leq n \leq k_{j+1}, use la monotonicidad de SnS^*_n (¡resúmenes no negativos!) para mostrar

    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}},

    y deducir, como:

    mαlim infSnnlim supSnnαm.\frac m\alpha \leq \liminf\frac{S^*_n}n \leq \limsup\frac{S^*_n}n \leq \alpha\,m .
  2. Deje que α1\alpha \downarrow 1 siga una secuencia y concluir Snnm\frac{S_n^*}n \to m a.s., por lo tanto (Parte I) el strong ley of large numbers:

     Snnna.s.E[X1]. \boxed{\ \frac{S_n}{n} \xrightarrow[n\to\infty]{\text{a.s.}} \E[X_1].\ }
  3. ¿Dónde estaba exactamente independencia por pares (en lugar de completo independencia) ¿suficiente? Enumera los tres lugares donde se invocaron hipótesis de tipo independentista.

Part V — Dividends.

  1. (Números normales de Borel) Demuestre que λ\lambda: casi todos los x[0,1]x \in \intcc01 son normal en cada base b2b \geq 2: cada dígito 0,,b10, \dots, b-1 aparece con frecuencia asintótica 1b\frac1b(fix bb and a digit, apply the strong ley 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 Teorema 22.6 — then intersect the countably many probability-one events). Exhibe un número explícito no normal, y reflexione: el teorema afirma la normalidad de casi todos los números, pero demostrando la normalidad de 2\sqrt2 o π\pi permanece abierto.
  2. (Montecarlo, garantizado) Justificar completamente el método de Ejemplo 22.14(b) para gL1([0,1]d)g \in L^1(\intcc01^d): construir el i.i.d. muestra uniforme en [0,1]d\intcc01^d de Teorema 22.6 y Ejercicio 22.10, e indique cuál es el fuerte ley cumple.

Part VI — What full independencia buys: maximal inequalities and random series. Etemadi gasta sólo por pares independencia; las partes restantes explotan el versión completa (mutua). Sea (Zn)(Z_n) independiente centrado variables de L2L^2 y Sk=Z1++ZkS_k = Z_1 + \dots + Z_k (una nueva notación, no relacionada con el XnX_n anterior).

  1. (desigualdad máxima de Kolmogorov) Para ε>0\varepsilon > 0 demuestre

    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) :

    El precio de Chebyshev compra el máximo (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 independencia of the coalitions (Z1,,Zk)(Z_1, \dots, Z_k) and (Zk+1,,Zn)(Z_{k+1}, \dots, Z_n), Teorema 22.5). señalar el paso donde independencia en pares no lo haría ya es suficiente.

  2. (Teorema de una serie de Khhinchin-Kolmogorov) Deducir: si nV(Zn)<\sum_n\V(Z_n) < \infty, entonces nZn\sum_nZ_n converge casi con seguridad (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. (serie Rademacher) Sea (εn)(\varepsilon_n) signos i.i.d. P(εn=±1)=12\P(\varepsilon_n = \pm1) = \frac12 (Teorema 22.6), y sean (xn)(x_n) números reales. mostrar eso nxnεn\sum_nx_n\varepsilon_n converge a.s. tan pronto como nxn2<\sum_nx_n^2 < \infty; mostrar también que, sea lo que sea (xn)(x_n), la probabilidad de que nxnεn\sum_nx_n\varepsilon_n converge es 00 o 11 (Teorema 22.9).
  4. Lo contrario, elementalmente. Configure Tn=knxkεkT_n = \sum_{k\leq n}x_k\varepsilon_kysn2=knxk2s_n^2 = \sum_{k\leq n}x_k^2y supongamossns_n \to \infty. (a) Demostrar el Paley–desigualdad de Zygmund: para Z0Z \geq 0 con EZ2<\E Z^2 < \infty y 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) Mostrar ETn43sn4\E T_n^4 \leq 3s_n^4. (c) Deduzca P(Tn>sn2)316\P\bigl(\abs{T_n} > \frac{s_n}2\bigr) \geq \frac3{16} y concluya que nxnεn\sum_nx_n\varepsilon_n diverge a.s.; de ahí el dicotomía

    nxnεn converges a.s.    nxn2<.\sum_nx_n\varepsilon_n\ \text{converges a.s.} \iff \sum_nx_n^2 < \infty .
  5. (Serie armónica aleatoria) Concluimos que nεnns\sum_n\frac{\varepsilon_n}{n^s} converge como si y sólo si s>12s > \frac12. Para 12<s1\frac12 < s \leq 1 la serie converge a.s. mientras nns=\sum_nn^{-s} = \infty: se producen signos aleatorios cancelación de fuerza de raíz cuadrada — comparar con la serie alterna n(1)nns\sum_n\frac{(-1)^n}{n^s}, que converge para cada s>0s > 0.

Part VII — Concentration: Hoeffding’s inequality. El fuerte ley dice Snnm\frac{S_n}n \to m; Las desigualdades de concentración dicen qué tan improbable es una desviación. en cada nn fijo.

  1. (lema de Hoeffding) (a) Mostrar coshλeλ2/2\cosh\lambda \leq \eu^{\lambda^2/2} para todos λR\lambda \in \R, comparando las dos series término a término. (b) Sea ZZ centrado con aZba \leq Z \leq b, a<ba < b. Mostrar

    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 esperanzas de heredar, and study φ(t)=pt+log(1p+pet)\varphi(t) = -pt + \log(1 - p + p\eu^t)withp=abap = \frac{-a}{b-a}andt=λ(ba)t = \lambda(b - a): show φ(0)=φ(0)=0\varphi(0) = \varphi'(0) = 0 and φ14\varphi'' \leq \frac14).

  2. (desigualdad de Hoeffding) Sea X1,,XnX_1, \dots, X_n independiente con aiXibia_i \leq X_i \leq b_i y Sn=X1++XnS_n = X_1 + \dots + X_n. Demuestre, parat>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),

    y lo mismo para la cola inferior (exponential Chebyshev: bound Eeλ(SnESn)\E\,\eu^{\lambda(S_n - \E S_n)} using independencia and question 17, then optimize over λ>0\lambda > 0).

  3. (El ley fuerte, caso acotado, con una tasa) Dejemos que XiX_i sea i.i.d. con valores en [a,b]\intcc ab y m=EX1m = \E X_1. Mostrar

    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)

    y recuperar Snnm\frac{S_n}n \to m a.s. por Borel–Cantelli: una segunda prueba de la fuerza ley para variables acotadas — sin truncamiento, un tasa exponencial en cada nn finito, pero acotado demandas y independencia completo. comparar el hipótesis con la de Etemadi.

  4. (Montecarlo, garantizado en nn fijo) Sea g ⁣:[0,1]d[0,1]g \colon \intcc01^d \to \intcc01 mensurable y (Uk)(U_k) la muestra uniforme i.i.d. de la pregunta 11. Dado ε,δ>0\varepsilon, \delta > 0, mostrar

    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,

    y evaluar el umbral para ε=δ=102\varepsilon = \delta = 10^{-2}. El límite no involucra add: comparar con la pregunta 11 y con determinista rejillas.

Part VIII — How big is a random walk? Toward the iterated logarithm. Sea Sn=ε1++εnS_n = \varepsilon_1 + \dots + \varepsilon_n el paseo aleatorio simple construido a partir de i.i.d. signos justos.

  1. (colas subgaussianas) Mostrar EeλSn=(coshλ)nenλ2/2\E\,\eu^{\lambda S_n} = (\cosh\lambda)^n \leq \eu^{n\lambda^2/2} y deducir, para 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. Deducir, vía 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). En particular, el paseo sigue vivo. la escala CLT n\sqrt n hasta un factor logarítmico — muy por debajo del límite del crudo Snn\abs{S_n} \leq n.

  3. A lo largo de la subsecuencia de duplicación nj=2jn_j = 2^j, muestre

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

    y reflejar: el ley of the iterated logarithm (Khinchin; Hartman–Wintner para general centrados L2L^2) establece que

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

    Explique con precisión qué separa la subsecuencia estimación recién demostrada de la mitad superior de este declaración (se debe controlar maxnjnnj+1Sn\max_{n_j \leq n \leq n_{j+1}}S_n dentro de cada bloque, lo que requiere una desigualdad máxima en la escala exponencial) y comprobar cuantitativamente que la pregunta 12 la desigualdad es demasiado débil para ese propósito. el inferior la mitad se basa en el segundo lema de Borel-Cantelli aplicado a bloques independiente; ambas mitades son Material honesto del tercer año para una probabilidad dedicada. curso.

  4. (Desviación uniforme sobre una clase finita) Sean A1,,ANA_1, \dots, A_N eventos en un experimento repetible, y estimar cada probabilidad por su valor empírico frecuencia p^i\hat p_i sobre nn i.i.d. repeticiones. Combinando la desigualdad de Hoeffding con un límite sindical, mostrar

    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},

    y deducir la regla del tamaño de la muestra: nln(2N/δ)2ε2n \geq \frac{\ln(2N/\delta)}{2\varepsilon^2} garantiza todas las estimaciones NN simultáneamente ε\varepsilon: exacto con probabilidad 1δ\geq 1 - \delta. CalculennparaN=106N = 10^6,ε=0.01\varepsilon = 0.01,δ=0.05\delta = 0.05: el precio logarítmico de uniformidad.

  5. (La ventana armónica aleatoria) Combinando los dos mitades de la teoría de series aleatorias, muestran que para i.i.d. firma (εn)(\varepsilon_n) la serie nεnnα\sum_n\frac{\varepsilon_n}{n^\alpha} converge a.s. si α>12\alpha > \frac12 y diverge a.s. si α12\alpha \leq \frac12; contraste con absoluto convergencia (que requiere α>1\alpha > 1): en el ventana α(12,1]\alpha \in \intoc{\frac12}1, convergencia es un fenómeno genuinamente probabilístico — cancelación, no tamaño.
Solución

Solución de Problema 22.1.

1. Xn±X_n^{\pm} son funciones Borel de XnX_n: permanecer en pares independiente (Ejercicio 22.3(a)) e idénticamente distribuido, integrable, con EX1=EX1+EX1\E X_1 = \E X_1^+ - \E X_1^-. Si el teorema se cumple para variables no negativas, aplícalo a ambas mitades y resta: 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) (idéntico leyes) y 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 (Ejercicio 11.3(a)). Borel–Cantelli (1): a.s. Xn=YnX_n = Y_n para todos los nn grandes, por lo que SnSnS_n - S_n^* es eventualmente constante en nn: SnSnn0\frac{S_n - S_n^*}n \to 0 a.s., y las dos sumas normalizadas comparten su asintótica comportamiento.

3. X11X1nX1X_1\mathbf 1_{X_1 \leq n} \nearrow X_1: MCT da EYnm\E Y_n \to m; Cesàro significa un convergente La secuencia converge al mismo límite: ESnn=1nknEYkm\frac{\E S_n^*}n = \frac1n\sum_{k\leq n}\E Y_k \to m. Por lo tanto basta con probar 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}]. Por Tonelli para serie,

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,

usando nxn24x\sum_{n\geq x}n^{-2} \leq \frac4x para x1x \geq 1 (para x2x \geq 2: 1x12x\leq \frac1{x-1} \leq \frac2x; para 1x<21 \leq x < 2:π264x\leq \frac{\pi^2}6 \leq \frac4x desde 4x>2\frac4x > 2) y X12/max(X1,1)X1X_1^2/\max(X_1, 1) \leq X_1 en ambos casos X11X_1 \gtrless 1.

5. Por pares independencia da E[(YiEYi)(YjEYj)]=0\E[(Y_i - \E Y_i)(Y_j - \E Y_j)] = 0poriji \neq j (el producto fórmula para dos variables), por lo que las variaciones suman: V(Sk)=nkV(Yn)\V(S^*_k) = \sum_{n\leq k}\V(Y_n). Chebyshev en cadakjk_j y sumando:

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 para la serie doble no negativa).

6. kj=αjαj2k_j = \lfloor\alpha^j\rfloor \geq \frac{\alpha^j}2(válido una vezαj1\alpha^j \geq 1, es decir todos j0j \geq 0: xx2\lfloor x\rfloor \geq \frac x2 para x1x \geq 1). Por lo tanto

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},

(serie geométrica del primer jj con αjn\alpha^j \geq n). Combinando con las preguntas 4 a 5, la suma doble es finito; Borel–Cantelli (1), aplicado para cada racional ε\varepsilon y se cruza, da SkjESkjkj0\frac{S^*_{k_j} - \E S^*_{k_j}}{k_j} \to 0 a.s., y con pregunta 3: Skjkjm\frac{S^*_{k_j}}{k_j} \to m a.s.

7. Yn0Y_n \geq 0 hace nSnn \mapsto S^*_n no decreciente: para 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},

cuál es el sándwich que se muestra después de insertar kjkj+1\frac{k_j}{k_{j+1}} y kj+1kj\frac{k_{j+1}}{k_j}. desde kj+1kjα\frac{k_{j+1}}{k_j} \to \alpha, la pregunta 6 da 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. Aplique la pregunta 7 para α=1+1p\alpha = 1 + \frac1p, pNp \in \N^*: contablemente muchos eventos a.s. en su intersección, dejando pp \to \infty: limSnn=m\lim\frac{S^*_n}n = ma.s. Con preguntas 1–3,SnnEX1\frac{S_n}n \to \E X_1 a.s.: el fuerte ley de grandes números, bajo pares independencia.

9. Aparecieron tres hipótesis de tipo independencia veces: (i) aditividad de varianzas (pregunta 5) — por pares es suficiente; (ii) distribución idéntica, en el sumas de truncamiento (pregunta 2) y el cálculo de la media (pregunta 3) — no independencia en absoluto; (iii) Borel–Cantelli (1) (preguntas 2 y 6) — válido sin cualquier independencia. independencia mutuo completo nunca fue invocado: observación de Etemadi.

10. Fija una base bb y un dígito rr. La base-bb Los dígitos (dk)(d_k) de un uniforme ω\omega son i.i.d. uniforme en {0,,b1}\{0, \dots, b-1\} (cada valor de vector de dígito ocupa un intervalo de longitud bmb^{-m}: el argumento de Teorema 22.6 textualmente). el fuerte ley aplicado a las variables acotadas i.i.d. 1dk=r\mathbf 1_{d_k = r}da: a.s., la frecuencia del dígitorr tiende a 1b\frac1b. Intersección sobre los contablemente muchos pares (b,r)(b, r): casi todos los números son simplemente normal en cada base. Una anormalidad explícita número: x=0.1001001002x = 0.100100100\ldots_2 (frecuencia de unos 1312\frac13 \neq \frac12). El contraste es humillante: casi todos los números son normales, pero para 2\sqrt2, e\eu o π\pi la normalidad sigue sin demostrarse — medida la teoría cuenta sin exponer.

11. Por Ejercicio 22.10 iterado, un una sola variable uniforme produce una secuencia de i.i.d. uniforme vectores UkU_k en [0,1]d\intcc01^d (dividir el conjunto de dígitos de cada UnU_n de Teorema 22.6 en subfamilias dd). Para gL1([0,1]d)g \in L^1(\intcc01^d), las variables g(Uk)g(U_k) son i.i.d. integrable con media g ⁣dλd\int g\,\dd\lambda_d (transferencia): el fuerte ley da

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 :

La integración Montecarlo converge casi con seguridad, en cada dimensión — el tamaño del error es asunto de la central teorema del límite (Capítulo 23).

12. Vamos Ak={Skε}j<k{Sj<ε}A_k = \{\abs{S_k} \geq \varepsilon\} \cap \bigcap_{j<k}\{\abs{S_j} < \varepsilon\}: losAkA_k son disjunto con la unión A={maxknSkε}A = \{\max_{k\leq n}\abs{S_k} \geq \varepsilon\}. entonces

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],

porque el término cruzado desaparece: Sk1AkS_k\mathbf 1_{A_k} es un Función Borel de la coalición (Z1,,Zk)(Z_1, \dots, Z_k), que es independiente de SnSkS_n - S_k, una función de (Zk+1,,Zn)(Z_{k+1}, \dots, Z_n) (Teorema 22.5), por lo que 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. EnAkA_k,Sk2ε2S_k^2 \geq \varepsilon^2, de dondeESn2ε2kP(Ak)=ε2P(A)\E S_n^2 \geq \varepsilon^2\sum_k\P(A_k) = \varepsilon^2\P(A); yESn2=knV(Zk)\E S_n^2 = \sum_{k\leq n}\V(Z_k) (se agregan variaciones). el El paso decisivo es la factorización: Sk1AkS_k\mathbf 1_{A_k} es una función no lineal de todo el primer bloque, y su independencia del segundo bloque es coalición independenciaindependencia por pares del ZiZ_i únicamente descorrelaciona pares y no lo justificaría.

13. Corrija NN y aplique la pregunta 12 a 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) .

Los eventos aumentan con mm; continuidad desde abajo da P(supk>NSkSN>ε)rN/ε2\P(\sup_{k>N}\abs{S_k - S_N} > \varepsilon) \leq r_N/\varepsilon^2yrN0r_N \to 0 por hipótesis. Por lo tanto para cada 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: casi seguramente, por cada pp hay NN con supk>NSkSN1p\sup_{k>N}\abs{S_k - S_N} \leq \frac1p (cruza el contablemente muchos eventos a.s. sobre pp), de modo que SkSl2p\abs{S_k - S_l} \leq \frac2ppara todos losk,l>Nk, l > N: las sumas parciales son a.s. Cauchy, por lo tanto a.s. convergentes.

14. Las variables Zn=xnεnZ_n = x_n\varepsilon_n son independiente (Funciones Borel de variables independiente, Ejercicio 22.3(a)), centrado, con V(Zn)=xn2\V(Z_n) = x_n^2: la pregunta 13 aplica cuandonxn2<\sum_nx_n^2 < \infty y da a.s. convergencia. En general, para cada NN el La convergencia de nxnεn\sum_nx_n\varepsilon_n no se ve afectada por los valores de ε1,,εN\varepsilon_1, \dots, \varepsilon_N: el evento de convergencia se encuentra en el álgebra σ\sigma de cola del independiente secuencia (εn)(\varepsilon_n), por lo que Kolmogorov Fuerzas ley cero uno (Teorema 22.9) su probabilidad es 00 o 11.

15. (a) División en el nivel θEZ\theta\E Z y usando Cauchy–Schwarz en la pieza superior,

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)} ,

entonces (1θ)EZEZ2P(Z>θEZ)(1 - \theta)\E Z \leq \sqrt{\E Z^2\,\P(Z > \theta\E Z)}; cuadrado. (b) AmpliarTn4=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]: el expectativa es 11 cuando los índices se emparejan (los cuatro iguales, o dos pares distintos, este último en 33 arreglos) y 00 en caso contrario (un signo no apareado tiene cero media y factorizado por independencia). Por lo tanto

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 con 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} .

Si la serie convergiera con probabilidad positiva, convergería a.s. (pregunta 14), por lo supnTn<\sup_n\abs{T_n} < \inftya.s., y algo deMM satisfaría P(supnTn>M)<316\P(\sup_n\abs{T_n} > M) < \frac3{16}; pero en cuanto 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}: contradicción. Entonces la divergencia es casi claro, y con la pregunta 14 la dicotomía es completo.

16. Aquí xn=nsx_n = n^{-s} y nn2s<\sum_nn^{-2s} < \inftyexactamente cuandos>12s > \frac12: mediante las preguntas 14–15, nεnns\sum_n\frac{\varepsilon_n}{n^s} converge a.s. si y sólo si s>12s > \frac12 (para s12s \leq \frac12, a.s. divergencia). Para 12<s1\frac12 < s \leq 1 la convergencia es nunca absoluto. La comparación es instructiva: perfectamente Los signos alternos se cancelan con la fuerza nsn^{-s} por cada s>0s > 0, mientras que los signos aleatorios típicos se cancelan solo con la fuerzansn^{-s}. fuerza de la raíz cuadrada — el paseo aleatorio de la pregunta 21 crece como n\sqrt n, y la suma de Abel convierte exactamente ese crecimiento hacia la convergencia de εnns\sum\varepsilon_nn^{-s} para s>12s > \frac12.

17. (a) coshλ=kλ2k(2k)!\cosh\lambda = \sum_k\frac{\lambda^{2k}}{(2k)!}yeλ2/2=kλ2k2kk!\eu^{\lambda^2/2} = \sum_k\frac{\lambda^{2k}}{2^kk!}; y(2k)!2kk!(2k)! \geq 2^kk! se mantiene a lo largo de los términos, porque (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! (cada factor satisface k+i2ik + i \geq 2i para iki \leq k), de modo que de hecho (2k)!2k(k!)22kk!(2k)! \geq 2^k(k!)^2 \geq 2^kk!. (b) Nota a0ba \leq 0 \leq b (ZZ es centrado), y por la convexidad de zeλzz \mapsto \eu^{\lambda z}, para 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)}

con 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). entonces φ(0)=0\varphi(0) = 0, φ(t)=p+pet1p+pet\varphi'(t) = -p + \frac{p\eu^t}{1 - p + p\eu^t}desaparece en00yφ(t)=ρ(1ρ)14\varphi''(t) = \rho(1 - \rho) \leq \frac14paraρ=pet1p+pet[0,1]\rho = \frac{p\eu^t}{1 - p + p\eu^t} \in \intcc01: Taylor en la orden22 da φ(t)t28=λ2(ba)28\varphi(t) \leq \frac{t^2}8 = \frac{\lambda^2(b-a)^2}8.

18. Para λ>0\lambda > 0, Markov solicitó al variable positiva eλ(SnESn)\eu^{\lambda(S_n - \E S_n)} (Proposición 22.3) y el producto la fórmula para las variables independiente da

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),

por la pregunta 17(b) aplicada a cada XiEXi[aiEXi,biEXi]X_i - \E X_i \in \intcc{a_i - \E X_i}{b_i - \E X_i} centrado (mismo ancho). Minimizar el exponente en λ=4tD\lambda = \frac{4t}{D}, D=i(biai)2D = \sum_i(b_i - a_i)^2produce2t2D-\frac{2t^2}D. el inferior tail sigue aplicando el resultado a (Xi)(-X_i).

19. Tome t=nεt = n\varepsilon y 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),

que es sumable en nn (una serie de tipo geométrico): Borel–Cantelli (Teorema 22.8) da ese culo Snnm<ε\abs{\frac{S_n}n - m} < \varepsilon eventualmente; intersección sobre ε=1p\varepsilon = \frac1p produce Snnm\frac{S_n}n \to m a.s. Comparación: Etemadi pregunta solo X1L1X_1 \in L^1 y independencia por pares, y entrega sin tarifa; Hoeffding pide limitación y completa independencia, y ofrece una garantía exponencial explícita en cada finito nn — los dos teoremas responden preguntas diferentes aproximadamente el mismo límite.

20. Los g(Uk)g(U_k) son i.i.d. con valores en [0,1]\intcc01 y significa g ⁣dλd\int g\,\dd\lambda_d (transferencia), entonces pregunta 18 con biai=1b_i - a_i = 1, t=nεt = n\varepsilon da el cota de dos caras 2e2nε2δ2\eu^{-2n\varepsilon^2} \leq \delta tan pronto como e2nε22δ\eu^{2n\varepsilon^2} \geq \frac2\delta, es decir nlog(2/δ)2ε2n \geq \frac{\log(2/\delta)}{2\varepsilon^2}. Para ε=δ=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 :

Aproximadamente las muestras 2650026\,500 garantizan una precisión 1%1\% con 99%99\% confianza — en todas las dimensiones dd, para cada mensurable integrando con valores en [0,1]\intcc01. Pregunta El fuerte ley de 11 prometió convergencia sin finito-nn garantía; una cuadrícula determinista con puntos kk por eje cuesta evaluaciones kdk^d, exponencial en dd. Concentración es lo que hace que Montecarlo sea un método en lugar de un esperanza.

21. Independencia y la fórmula del producto: 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} por pregunta 17(a). Markov en 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,

y el límite simétrico para Sn-S_n (mismo ley) duplica el constante para Sn\abs{S_n}.

22. Repare η>0\eta > 0 y configure xn=(1+η)2nlognx_n = (1 + \eta)\sqrt{2n\log n}paran2n \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}},

sumable desde (1+η)2>1(1 + \eta)^2 > 1. Borel–Cantelli: a.s. Sn<(1+η)2nlogn\abs{S_n} < (1 + \eta)\sqrt{2n\log n} para todos los grandes nn, entonces lim supnSn2nlogn1+η\limsup_n\frac{\abs{S_n}}{\sqrt{2n\log n}} \leq 1 + \etaa.s.; al cruzar los eventos a.s. paraη=1p\eta = \frac1p,pNp \in \N^*, se obtiene el reclamo. El paseo del tamaño nn tiene la amplitud típica n\sqrt n (su varianza), y incluso sus peores excursiones superan esa escala como máximo 2logn\sqrt{2\log n}.

23. Con nj=2jn_j = 2^j y x=(1+η)2njloglognjx = (1 + \eta)\sqrt{2n_j\log\log n_j}(definido paraj2j \geq 2), la pregunta 21 da

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},

sumable en jj desde (1+η)2>1(1 + \eta)^2 > 1: Borel–Cantelli y η=1p\eta = \frac1p le dan a lim supjSnj/2njloglognj1\limsup_jS_{n_j}/\sqrt{2n_j \log\log n_j} \leq 1 un.s. Lo que falta para el completo La mitad superior es el puente entre los puntos de control: hay que mostrar maxnjnnj+1Sn\max_{n_j\leq n\leq n_{j+1}}S_n excede (1+η)2njloglognj(1+\eta)\sqrt{2n_j\log\log n_j} sólo un número limitado de veces, lo que exige una desigualdad máxima con gaussiano colas (reflejo de la desigualdad de Lévy o Ottaviani) desigualdad, no demostrada aquí). La pregunta 12 es cuantitativa demasiado débil: limita la probabilidad por

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)},

que tiende a 00 pero es no sumable en jj: Borel-Cantelli no puede concluir. La mitad inferior del ley del logaritmo iterado se aplica el segundo Borel–Lema de Cantelli a los incrementos independiente Snj+1SnjS_{n_{j+1}} - S_{n_j}, utilizando límites inferiores coincidentes para Colas de tipo gaussiano. Ambos refinamientos son genuinos del Año 3. probabilidad, un curso más adelante; que problema este entrega sin ayuda es la escala exacta de logaritmo iterado a lo largo de tiempos geométricos.

24. Cada p^i\hat p_i es un promedio de nn i.i.d. variables indicadoras con valores en [0,1]\intcc01 y media P(Ai)\P(A_i): Hoeffding da P(p^iP(Ai)>ε)2e2nε2\P(\abs{\hat p_i - \P(A_i)} > \varepsilon) \leq 2\eu^{-2n\varepsilon^2}. El sindicato obligado se multiplica por NN. Resolviendo 2Ne2nε2δ2N\eu^{-2n\varepsilon^2} \leq \delta:nln(2N/δ)2ε2n \geq \frac{\ln(2N/\delta)}{2\varepsilon^2}. Numéricamente: ln21060.05=ln(4107)17.5\ln\frac{2\cdot10^6}{0.05} = \ln(4\cdot10^7) \approx 17.5, entoncesn17.5210487600n \geq \frac{17.5}{2\cdot10^{-4}} \approx 87\,600: estimación uno probabilidad de ±1%\pm1\% toma aproximadamente 1850018\,500 muestras (ln(2/δ)/2ε2\ln(2/\delta)/2\varepsilon^2), y un millón probabilidades solamente 4.7\approx 4.7 veces más — la uniformidad cuesta lnN\ln N, no NN: la observación de que hace la minimización empírica del riesgo, y con ello la máquina aprendizaje, estadísticamente posible.

25. Las variables Xn=εnnαX_n = \frac{\varepsilon_n} {n^\alpha} son independiente, centradas, acotadas, con nV(Xn)=nn2α\sum_n\V(X_n) = \sum_nn^{-2\alpha}. Si α>12\alpha > \frac12: la serie de varianzas converge y la El teorema de una serie (Parte VI) da la convergencia a.s. de Xn\sum X_n. Si α12\alpha \leq \frac12: la serie de varianza diverge, y la mitad inversa (Parte VI Argumento Paley-Zygmund, aplicable ya que las demandas son delimitado por 11) da una divergencia. absoluto convergencia pregunta nα<\sum n^{-\alpha} < \infty: α>1\alpha > 1. En (12,1]\intoc{\frac12}1, la serie converge a.s. aunque Xn=\sum\abs{X_n} = \infty seguramente: los signos conspiran para cancelar, con probabilidad uno — convergencia por cancelación, invisible a cualquier prueba absoluta, y (por el ley cero uno) con un veredicto determinista de todos modos.