Mathematics · Buku 4 · Bachelor Year 2

Matematika Universitas — Tahun 2

Matematika Universitas — Tahun 2 · Bachelor Year 2

22Peubah Acak Diskret

Peubah acak menata perhitungan peluang di sekitar fungsi alih-alih kejadian. Pada ruang terbilang teorinya digerakkan keluarga terjumlahkan pada Bab 7: sebab nilai harapan adalah jumlah sebuah keluarga yang diindeks ruang sampelnya, dan semua sifatnya — kelinearan, pemindahan, rumus hasil kali bagi peubah yang saling bebas — merupakan teorema tentang keluarga yang terjumlahkan. Bab ini membuktikan ketaksamaan kunci Markov, Chebyshev, Cauchy–Schwarz dan Jensen, lalu berakhir dengan hukum klasiknya beserta hukum bilangan besar yang lemah, yang buktinya dua baris begitu Chebyshev tersedia.

22.1 Peubah acak dan hukumnya

Definisi 22.1 (Peubah acak diskret; hukumnya)

Misalkan (Ω,P)(\Omega, \P) ruang peluang yang terbilang. Sebuah peubah acak adalah pemetaan X ⁣:ΩEX \colon \Omega \to E (dengan EE sembarang himpunan; dan disebut real bila E=RE = \R). Adapun hukum (atau distribusi) peubahnya adalah ukuran peluang PX\P_X pada himpunan terbilang X(Ω)X(\Omega) yang didefinisikan lewat

PX({x})=P(X=x)=P({ω:X(ω)=x}).\P_X(\{x\}) = \P(X = x) = \P\bigl(\{\omega : X(\omega) = x\}\bigr) .

Contoh 22.2 (Hukum klasiknya)

  • Bernoulli B(p)\mathcal{B}(p): dengan X{0,1}X \in \{0, 1\} dan P(X=1)=p\P(X = 1) = p. Yakni indikator sebuah kejadian.
  • Binomial B(n,p)\mathcal{B}(n, p): dengan P(X=k)=(nk)pk(1p)nk\P(X = k) = \binom nk p^k(1-p)^{n-k} untuk 0kn0 \leq k \leq n: yakni banyaknya keberhasilan pada nn percobaan Bernoulli yang saling bebas (jilid Kelas 10–12; yang dibuktikan ulang di bawah lewat jumlah peubah yang saling bebas).
  • Geometri G(p)\mathcal{G}(p): dengan P(X=k)=(1p)k1p\P(X = k) = (1-p)^{k-1}p untuk kNk \in \N^*: yakni pangkat keberhasilan pertamanya (Contoh 21.5).
  • Poisson P(λ)\mathcal{P}(\lambda): dengan P(X=k)=eλλkk!\P(X = k) = e^{-\lambda}\frac{\lambda^k}{k!} untuk kNk \in \N — yakni sebuah ukuran peluang berkat deret eksponensialnya. Inilah hukum kejadian langka (Bab 23).

Catatan 22.3 (Hukum mana memodelkan apa)

Keempat hukumnya menjawab empat pertanyaan purba: Bernoulli, “apakah ia terjadi?”; binomial, “berapa kali dalam nn percobaan?”; geometri, “berapa lama sampai kali pertamanya?”; dan Poisson, “berapa banyak kejadian pada laju tertentu, ketika percobaannya banyak dan masing-masing tak mungkin?”. Jadi mengenali pertanyaannya adalah sembilan persepuluh pemodelannya: sebab jumlah indikator menunjuk ke binomialnya, waktu tunggu ke geometrinya, dan cacah kejadian langka ke Poissonnya — dengan peralihan dari binomial ke Poisson dipersis oleh hukum kejadian langka pada Bab 23.

Proposisi 22.4 (Sifat tanpa ingatan pada hukum geometri)

Bila XG(p)X \sim \mathcal{G}(p), maka untuk setiap m,nNm, n \in \N:

P(X>m+nX>m)=P(X>n),\P(X > m + n \mid X > m) = \P(X > n) ,

dan hukum geometri adalah satu-satunya hukum pada N\N^* yang bersifat demikian.

Bukti. Menjumlahkan bobot geometrinya memberi P(X>n)=(1p)n\P(X > n) = (1-p)^n. Karena itu

P(X>m+nX>m)=P(X>m+n)P(X>m)=(1p)m+n(1p)m=(1p)n=P(X>n).\P(X > m + n \mid X > m) = \frac{\P(X > m + n)}{\P(X > m)} = \frac{(1-p)^{m+n}}{(1-p)^m} = (1-p)^n = \P(X > n).

Sebaliknya, bila G(n)=P(X>n)G(n) = \P(X > n) memenuhi G(m+n)=G(m)G(n)G(m + n) = G(m)G(n) dengan G(0)=1G(0) = 1, maka G(n)=G(1)nG(n) = G(1)^n secara induktif; lalu q=G(1)[0,1)q = G(1) \in \intco{0}{1}, dan q=0q = 0 atau hukumnya adalah G(1q)\mathcal{G}(1 - q): sebab P(X=k)=G(k1)G(k)=qk1(1q)\P(X = k) = G(k-1) - G(k) = q^{k-1}(1 - q).

Contoh 22.5 (Tak ada angka yang pernah “jatuh tempo”)

Lemparkan sebuah dadu sambil menunggu enam: maka waktu tunggunya XG(1/6)X \sim \mathcal G(1/6). Adapun sifat tanpa ingatannya mengatakan bahwa sesudah 1010 lemparan yang sia-sia, tunggu sisanya X10X - 10, jika diketahui X>10X > 10, kembali berhukum G(1/6)\mathcal G(1/6): jadi nilai harapan tunggu bersyaratnya masih 66 lemparan, persis seperti pada awalnya. Dadunya tak mengingat, dan tak ada enam yang pernah “jatuh tempo” — jadi kekeliruan penjudi adalah kepercayaan bahwa hukum bersyaratnya sepatutnya bergeser. Sebaliknya, paruh ketunggalan proposisinya mengatakan bahwa ketidakpedulian ini mencirikan waktu tunggu geometri: sebab sembarang waktu tunggu yang ramalannya tak pernah diperbarui bersifat geometri. Adapun antrean dan umur pakai yang nyata biasanya memperbarui, dan persis begitulah kita mendeteksi bahwa keduanya tak geometri.

22.2 Nilai harapan

Definisi 22.6 (Nilai harapan)

Sebuah peubah acak real XX pada (Ω,P)(\Omega, \P) disebut punya nilai harapan bila keluarga (X(ω)P({ω}))ωΩ\bigl(X(\omega)\,\P(\{\omega\})\bigr)_{\omega \in \Omega} bersifat terjumlahkan (Bab 7); dan nilai harapannya lalu

E(X)=ωΩX(ω)P({ω}).\E(X) = \sum_{\omega \in \Omega} X(\omega)\,\P(\{\omega\}) .

Teorema 22.7 (Teorema pemindahan)

Peubah XX punya nilai harapan bila dan hanya bila keluarga (xP(X=x))xX(Ω)\bigl(x\,\P(X = x)\bigr)_{x \in X(\Omega)} terjumlahkan, dan lalu

E(X)=xX(Ω)xP(X=x).\E(X) = \sum_{x \in X(\Omega)} x\,\P(X = x) .

Lebih umum, untuk f ⁣:X(Ω)Rf \colon X(\Omega) \to \R, peubah f(X)f(X) punya nilai harapan bila dan hanya bila xf(x)P(X=x)<\sum_x \abs{f(x)}\,\P(X = x) < \infty, dan lalu E(f(X))=xf(x)P(X=x)\E(f(X)) = \sum_x f(x)\,\P(X = x).

Bukti. Partisikanlah Ω\Omega menjadi himpunan arasnya Ωx={X=x}\Omega_x = \{X = x\} dengan xX(Ω)x \in X(\Omega). Menurut teorema penjumlahan lewat paket bagi keluarga yang terjumlahkan (Bab 7), keluarga (X(ω)P({ω}))ω(X(\omega)\P(\{\omega\}))_\omega terjumlahkan bila dan hanya bila tiap paketnya terjumlahkan (yang otomatis: sebab ωΩxxP({ω})=xP(X=x)\sum_{\omega \in \Omega_x}\abs{x}\P(\{\omega\}) = \abs x\,\P(X = x)) dan keluarga jumlah paketnya (xP(X=x))x\bigl(x\,\P(X = x)\bigr)_x terjumlahkan — dan lalu jumlah totalnya sepakat. Adapun untuk f(X)f(X): terapkanlah pernyataan yang terbukti itu pada peubah Y=fXY = f \circ X, yang himpunan arasnya {Y=y}=x:f(x)=y{X=x}\{Y = y\} = \bigsqcup_{x : f(x) = y}\{X = x\}; lalu penjumlahan lewat paket yang kedua mengubah yyP(Y=y)\sum_y y\,\P(Y = y) menjadi xf(x)P(X=x)\sum_x f(x)\,\P(X = x), dengan paketnya kini mengelompokkan nilai xx menurut petanya f(x)f(x), dan keterjumlahan mutlak satu keluarganya setara dengan keterjumlahan yang lain.

Teorema 22.8 (Sifat nilai harapan)

Pada himpunan peubah acak yang punya nilai harapan:

  1. (Kelinearan) E(aX+bY)=aE(X)+bE(Y)\E(aX + bY) = a\,\E(X) + b\,\E(Y).
  2. (Kepositifan dan kemonotonan) X0E(X)0X \geq 0 \Rightarrow \E(X) \geq 0; XYE(X)E(Y)X \leq Y \Rightarrow \E(X) \leq \E(Y); dan E(X)E(X)\abs{\E(X)} \leq \E(\abs X).
  3. (Dominasi) Bila XZ\abs X \leq Z dan ZZ punya nilai harapan, maka XX juga.

Bukti. Semuanya merupakan sifat jumlah keluarga yang terjumlahkan (Bab 7): yakni kelinearan jumlahnya, kepositifan suku demi suku, dan kriteria dominasi bagi keterjumlahannya. (Perhatikan bahwa kelinearannya seketika pada definisinya atas Ω\Omega, sedangkan ia akan janggal pada rumus pemindahannya — yakni satu keuntungan mendefinisikan E\E di hulu.)

Contoh 22.9

Untuk XB(n,p)X \sim \mathcal{B}(n, p): dengan menulis X=X1++XnX = X_1 + \dots + X_n sebagai jumlah indikator Bernoulli lalu memakai kelinearannya, E(X)=np\E(X) = np — tanpa perlu koefisien binomial. Untuk XG(p)X \sim \mathcal{G}(p): E(X)=k1k(1p)k1p=p1(1(1p))2=1p\E(X) = \sum_{k\geq1}k(1-p)^{k-1}p = p\cdot\frac{1}{(1 - (1-p))^2} = \frac1p, lewat menurunkan deret geometrinya di dalam cakramnya (Bab 11). Untuk XP(λ)X \sim \mathcal{P}(\lambda): E(X)=k1keλλkk!=λeλj0λjj!=λ\E(X) = \sum_{k\geq1}k e^{-\lambda}\frac{\lambda^k}{k!} = \lambda e^{-\lambda}\sum_{j\geq0}\frac{\lambda^j}{j!} = \lambda.

Contoh 22.10 (Pemindahan beraksi)

Untuk XP(λ)X \sim \mathcal P(\lambda), hitunglah E(11+X)\E\bigl(\frac1{1+X}\bigr) — sebab hukum 11+X\frac1{1+X} sendiri janggal, tetapi pemindahannya tak pernah memintanya:

E(11+X)=k01k+1eλλkk!=eλλk0λk+1(k+1)!=eλλ(eλ1)=1eλλ.\E\Bigl(\frac1{1+X}\Bigr) = \sum_{k\geq0}\frac{1}{k+1}\,\eu^{-\lambda} \frac{\lambda^k}{k!} = \frac{\eu^{-\lambda}}{\lambda}\sum_{k\geq0} \frac{\lambda^{k+1}}{(k+1)!} = \frac{\eu^{-\lambda}}{\lambda}\bigl(\eu^\lambda - 1\bigr) = \frac{1 - \eu^{-\lambda}}{\lambda} .

Ada dua pelajaran. Secara hitungan: mengenali deret eksponensial yang tergeser adalah seluruh kerjanya — sebab pemindahan menyusutkan nilai harapan f(X)f(X) menjadi manipulasi deret. Secara struktural: nilai colok naifnya adalah 11+EX=11+λ\frac1{1 + \E X} = \frac1{1 + \lambda}, sedangkan jawaban benarnya lebih besar,

1eλλ11+λ,\frac{1 - \eu^{-\lambda}}{\lambda} \geq \frac{1}{1 + \lambda},

persis seperti dituntut ketaksamaan Jensen bagi fungsi cembung t11+tt \mapsto \frac1{1+t}. Jadi nilai harapan peta cembung duduk di atas nilai colok naifnya, dan pemindahan ditambah pemeriksaan deret membuat ketaksamaan abstraknya menjadi konkret.

Teorema 22.11 (Kebebasan dan hasil kali)

Peubah acak X,YX, Y disebut saling bebas bila P(X=x,Y=y)=P(X=x)P(Y=y)\P(X = x, Y = y) = \P(X = x)\P(Y = y) untuk setiap x,yx, y — setara dengan itu, bila kejadian {XA}\{X \in A\} dan {YB}\{Y \in B\} saling bebas untuk setiap A,BA, B. Bila XX dan YY merupakan peubah real yang saling bebas dan punya nilai harapan, maka XYXY punya nilai harapan dan

E(XY)=E(X)E(Y).\E(XY) = \E(X)\,\E(Y) .

Bukti. Kesetaraan kedua perumusannya menyusul dengan menjumlahkan kesamaan titik demi titiknya atas (x,y)A×B(x, y) \in A \times B (lewat keaditifan-σ\sigma dua kali). Adapun untuk hasil kalinya: keluarga gandanya (xyP(X=x)P(Y=y))(x,y)\bigl(xy\,\P(X = x)\P(Y = y)\bigr)_{(x,y)} terjumlahkan, sebab menurut Fubini bagi keluarga (Bab 7)

x,yxyP(X=x)P(Y=y)=(xxP(X=x))(yyP(Y=y))<;\sum_{x, y}\abs x \abs y\,\P(X{=}x)\P(Y{=}y) = \Bigl(\sum_x \abs x \P(X{=}x)\Bigr) \Bigl(\sum_y \abs y \P(Y{=}y)\Bigr) < \infty ;

lalu berkat kebebasannya keluarga ini persis (xyP(X=x,Y=y))\bigl(xy\,\P(X = x, Y = y)\bigr), yang jumlahnya adalah E(XY)\E(XY) menurut pemindahan yang diterapkan pada peubah (X,Y)xy(X, Y) \mapsto xy; dan Fubini lagi menilai jumlah tak bertandanya sebagai hasil kali E(X)E(Y)\E(X)\E(Y).

Contoh 22.12 (Hasil kali, dengan dan tanpa kebebasan)

Lemparkan dua dadu setimbang. Bila YY adalah dadu kedua (yang bebas dari yang pertama), maka E(XY)=E(X)E(Y)=3.52=12.25\E(XY) = \E(X)\E(Y) = 3.5^2 = 12.25. Sedangkan bila Y=XY = X (yakni “hasil kali” sebuah dadu dengan dirinya sendiri),

E(X2)=1+4+9+16+25+366=91615.1712.25:\E(X^2) = \frac{1 + 4 + 9 + 16 + 25 + 36}{6} = \frac{91}{6} \approx 15.17 \neq 12.25 :

jadi hukum marginal yang sama pada kedua skenarionya, hukum bersama yang berbeda, dan nilai harapan hasil kali yang berbeda. Adapun pelajarannya, yang layak dipahat: E(XY)\E(XY) adalah fungsional pasangannya, bukan kedua marginalnya — dan senjangnya E(X2)E(X)22.92\E(X^2) - \E(X)^2 \approx 2.92 persis, menurut König–Huygens, merupakan ragam 3512\frac{35}{12} dadunya.

22.3 Ragam, kovariansi, dan ketaksamaan klasiknya

Definisi 22.13 (Momen, ragam)

Peubah XX disebut punya momen orde 2 bila X2X^2 punya nilai harapan (dan lalu XX juga, berkat dominasi: sebab X1+X22\abs X \leq \frac{1 + X^2}{2}). Adapun ragam dan simpangan baku peubahnya lalu

V(X)=E((XE(X))2)=E(X2)E(X)2,σ(X)=V(X),V(X) = \E\bigl((X - \E(X))^2\bigr) = \E(X^2) - \E(X)^2 , \qquad \sigma(X) = \sqrt{V(X)} ,

(yakni bentuk keduanya — rumus König–Huygens — lewat menguraikan kuadratnya lalu memakai kelinearannya:

E((XEX)2)=E(X22XEX+E(X)2)=E(X2)2E(X)2+E(X)2,\E\bigl((X - \E X)^2\bigr) = \E\bigl(X^2 - 2X\,\E X + \E(X)^2\bigr) = \E(X^2) - 2\,\E(X)^2 + \E(X)^2 ,

dengan suku tengahnya memakai bahwa EX\E X sebuah konstanta). Adapun untuk X,YX, Y yang bermomen kedua, kovariansi keduanya adalah

Cov(X,Y)=E((XEX)(YEY))=E(XY)E(X)E(Y).\operatorname{Cov}(X, Y) = \E\bigl((X - \E X)(Y - \E Y)\bigr) = \E(XY) - \E(X)\E(Y) .

Teorema 22.14 (Perkakas ragam)

Untuk peubah yang bermomen kedua:

  1. V(aX+b)=a2V(X)V(aX + b) = a^2\,V(X);
  2. V(X+Y)=V(X)+V(Y)+2Cov(X,Y)V(X + Y) = V(X) + V(Y) + 2\operatorname{Cov}(X, Y), dan lebih umum

    V(i=1nXi)=i=1nV(Xi)+2i<jCov(Xi,Xj);V\Bigl(\sum_{i=1}^n X_i\Bigr) = \sum_{i=1}^n V(X_i) + 2\sum_{i < j}\operatorname{Cov}(X_i, X_j) ;
  3. bila X,YX, Y saling bebas, maka Cov(X,Y)=0\operatorname{Cov}(X, Y) = 0 (sedangkan konversnya salah), jadi ragam peubah yang saling bebas bertambah.

Bukti. Butir 1 dan 2 merupakan penguraian kuadrat ditambah kelinearannya; sedangkan hasil kali XiXjX_iX_j punya nilai harapan berkat Cauchy–Schwarz di bawah (atau berkat XiXjXi2+Xj22\abs{X_iX_j} \leq \frac{X_i^2 + X_j^2}{2}). Adapun butir 3 adalah Teorema 22.11 yang diterapkan pada peubah yang dipusatkan. Contoh tandingan baku bagi konversnya: XX yang berhukum seragam pada {1,0,1}\{-1, 0, 1\} dan Y=X2Y = X^2 tak berkorelasi (sebab E(XY)=E(X3)=0=EXEY\E(XY) = \E(X^3) = 0 = \E X \cdot \E Y) tetapi jelas bergantungan.

Teorema 22.15 (Ketaksamaan Markov dan Chebyshev)

  1. (Markov) Bila X0X \geq 0 punya nilai harapan, maka untuk setiap a>0a > 0:

    P(Xa)E(X)a.\P(X \geq a) \leq \frac{\E(X)}{a} .
  2. (Chebyshev) Bila XX bermomen kedua, maka untuk setiap ε>0\varepsilon > 0:

    P(XE(X)ε)V(X)ε2.\P\bigl(\abs{X - \E(X)} \geq \varepsilon\bigr) \leq \frac{V(X)}{\varepsilon^2} .

Bukti. 1. Titik demi titik berlaku a1XaXa\,\mathbf{1}_{X \geq a} \leq X (sebab pada kejadiannya ruas kirinya aXa \leq X; sedangkan di luarnya 0X0 \leq X). Lalu ambillah nilai harapannya: aP(Xa)E(X)a\,\P(X \geq a) \leq \E(X) berkat kemonotonannya dan E(1A)=P(A)\E(\mathbf{1}_A) = \P(A). 2. Terapkan Markov pada peubah taknegatif (XEX)2(X - \E X)^2 di aras a=ε2a = \varepsilon^2: sebab kejadian {(XEX)2ε2}\{(X - \E X)^2 \geq \varepsilon^2\} persis merupakan {XEXε}\{\abs{X - \E X} \geq \varepsilon\}.

Contoh 22.16 (Tak berkorelasi tetapi terlem bersama)

Lemparkan dua dadu setimbang, dengan XX dan YY saling bebas, lalu tetapkan S=X+YS = X + Y dan D=XYD = X - Y. Menurut kebilinearan kovariansinya,

Cov(S,D)=V(X)V(Y)+Cov(Y,X)Cov(X,Y)=V(X)V(Y)=0:\operatorname{Cov}(S, D) = V(X) - V(Y) + \operatorname{Cov}(Y, X) - \operatorname{Cov}(X, Y) = V(X) - V(Y) = 0 :

jadi jumlah dan selisihnya tak berkorelasi. Saling bebas? Tentu tidak: sebab S=12S = 12 memaksa D=0D = 0, sedangkan P(D=0)=16\P(D = 0) = \frac16 tanpa syarat. Jadi korelasi hanya menguji bagian linear sebuah kebergantungan; sebab di sini kebergantungannya diusung kendala bahwa SS dan DD berparitas sama, yang tak terlihat oleh kovariansinya. (Adapun bagi pasangan ini, kovariansi nolnya menuntut V(X)=V(Y)V(X) = V(Y): jadi distribusi yang identik, bukan kebebasan, yang mengerjakannya.)

Contoh 22.17 (Ketika Markov eksak)

Ketaksamaan Markov menjadi kesamaan persis ketika tak ada yang terbuang pada batas a1XaXa\,\mathbf 1_{X\geq a} \leq X: yakni peubahnya hanya boleh mengambil nilai 00 dan aa. Secara konkret, bila P(X=a)=π\P(X = a) = \pi dan P(X=0)=1π\P(X = 0) = 1 - \pi, maka E(X)=aπ\E(X) = a\pi dan

P(Xa)=π=E(X)a.\P(X \geq a) = \pi = \frac{\E(X)}{a} .

Adapun pembacaan realistisnya: pada sebuah populasi yang kekayaan rata-ratanya 100100 dan kekayaannya entah 00 entah 10610^6, proporsi jutawannya persis 10410^{-4} — yakni batas Markov, yang tercapai persis lewat ketaksamaan maksimalnya. Jadi setiap kali XX menyebar atas nilai antaranya batasnya menjadi tegas, kerap secara liar; tetapi seperti ditunjukkan kasus ekstremnya, tak ada ketaksamaan yang lebih baik yang dapat disarikan dari rerataannya sendirian.

Contoh 22.18 (Chebyshev bersifat tajam — tanpa hipotesis lebih lanjut)

Tetapkan ε>0\varepsilon > 0 dan q(0,1]q \in \intoc01, lalu biarkanlah XX mengambil nilai ±ε\pm\varepsilon masing-masing berpeluang q2\frac q2 dan 00 berpeluang 1q1 - q. Maka E(X)=0\E(X) = 0, V(X)=qε2V(X) = q\varepsilon^2, dan

P(XEXε)=q=V(X)ε2:\P\bigl(\abs{X - \E X} \geq \varepsilon\bigr) = q = \frac{V(X)}{\varepsilon^2} :

yakni kesamaan pada Chebyshev. Jadi ketaksamaannya tak dapat diperbaiki dengan memakai ragamnya saja — sebab peluruhan 1/ε21/\varepsilon^2-nya adalah harga persis informasi momen keduanya. Adapun peluruhan yang lebih cepat menuntut hipotesis yang lebih kuat: sebab keterbatasan peubahnya membeli pemusatan eksponensial, sebagaimana dipratinjau Latihan 22.7 dan dikembangkan secara sistematis oleh soal akhir pekan bab ini.

Teorema 22.19 (Cauchy–Schwarz dan Jensen)

  1. (Cauchy–Schwarz) Bila X,YX, Y bermomen kedua, maka XYXY punya nilai harapan dan E(XY)2E(X2)E(Y2)\E(XY)^2 \leq \E(X^2)\,\E(Y^2); sehingga Cov(X,Y)2V(X)V(Y)\operatorname{Cov}(X,Y)^2 \leq V(X)V(Y).
  2. (Jensen) Bila φ ⁣:IR\varphi \colon I \to \R cembung pada sebuah interval yang memuat X(Ω)X(\Omega), dan XX serta φ(X)\varphi(X) punya nilai harapan, maka

    φ(E(X))E(φ(X)).\varphi\bigl(\E(X)\bigr) \leq \E\bigl(\varphi(X)\bigr) .

Bukti. 1. Keterjumlahan XYXY: sebab XYX2+Y22\abs{XY} \leq \frac{X^2 + Y^2}2. Adapun pemetaan (X,Y)E(XY)(X, Y) \mapsto \E(XY) merupakan bentuk bilinear simetrik yang positif pada ruang peubah yang bermomen kedua, jadi ketaksamaan Cauchy–Schwarz abstrak pada Bab 12 berlaku (sebab semidefinit positif sudah cukup bagi ketaksamaannya). Lalu menerapkannya pada peubah yang dipusatkan memberi batas kovariansinya.

2. Pertama, m=E(X)m = \E(X) terletak di II: sebab II merupakan interval yang memuat semua nilai XX, dan nilai harapan bersifat monoton, jadi mm berada di antara infX(Ω)\inf X(\Omega) dan supX(Ω)\sup X(\Omega). Lalu menurut teorema garis penyangga bagi fungsi cembung (Bab 8), ada α,β\alpha, \beta dengan φ(t)αt+β\varphi(t) \geq \alpha t + \beta untuk setiap tIt \in I dan φ(m)=αm+β\varphi(m) = \alpha m + \beta. Maka, titik demi titik pada Ω\Omega, φ(X)αX+β\varphi(X) \geq \alpha X + \beta; jadi setelah nilai harapannya diambil,

E(φ(X))αE(X)+β=φ(E(X)).\E\bigl(\varphi(X)\bigr) \geq \alpha\,\E(X) + \beta = \varphi\bigl(\E(X)\bigr). \qedhere

Contoh 22.20

Jensen dengan φ(t)=t2\varphi(t) = t^2 memberi E(X)2E(X2)\E(X)^2 \leq \E(X^2) — yakni kepositifan ragamnya; sedangkan dengan φ(t)=1/t\varphi(t) = 1/t pada (0,)\intoo{0}{\infty}: 1EXE(1X)\frac{1}{\E X} \leq \E\bigl(\frac1X\bigr) — yakni rerata harmoniknya di bawah rerata aritmetikanya, kini dalam bentuk acak.

Catatan 22.21 (Jebakan yang sering muncul)

(i) Rumus E(XY)=E(X)E(Y)\E(XY) = \E(X)\E(Y) menuntut kebebasan (atau sekurangnya kovariansi nol): sebab mengambil Y=XY = X memberi E(X2)E(X)2\E(X^2) \neq \E(X)^2 setiap kali V(X)>0V(X) > 0. (ii) Demikian pula V(X+X)=4V(X)V(X + X) = 4V(X), bukan 2V(X)2V(X): sebab ragam hanya bertambah antar suku yang saling bebas (atau yang tak berkorelasi). (iii) Nilai E(f(X))\E(f(X)) bukanlah f(E(X))f(\E(X)); dan untuk ff yang cembung Jensen bahkan memberitahukanmu arah galatnya, seperti pada Contoh 22.10. (iv) Keberadaannya sebuah hipotesis yang sejati: sebab untuk peubah St Petersburg X=2KX = 2^K dengan P(K=k)=2k\P(K = k) = 2^{-k} (untuk k1k \geq 1),

k12k2k=k11=:\sum_{k\geq1}2^k\cdot2^{-k} = \sum_{k\geq1}1 = \infty :

jadi XX berhingga hampir pasti namun tak punya nilai harapan, sehingga tak ada harga masuk yang adil bagi permainannya. Jadi keterjumlahan pada definisi E\E bukanlah kerewelan tata buku — melainkan tempat ekor yang berat terdeteksi. (v) Akhirnya, teorema pemindahannya menuntut keterjumlahan mutlak sebelum sembarang penataan ulang jumlahnya atas nilainya menjadi sah (Bab 7).

Contoh 22.22 (Chebyshev pada seratus lemparan)

Untuk XB(100,12)X \sim \mathcal B(100, \frac12): berlaku EX=50\E X = 50 dan V(X)=25V(X) = 25. Lalu Chebyshev dengan ε=6\varepsilon = 6:

P(45X55)=P(X50<6)125360.31,\P(45 \leq X \leq 55) = \P(\abs{X - 50} < 6) \geq 1 - \frac{25}{36} \approx 0.31 ,

sedangkan jumlah binomial eksaknya memberi 0.73\approx 0.73. Jadi jaminan 31%31\%-nya jauh dari kebenarannya, tetapi ia hanya menuntut rerataannya dan ragamnya saja — sebab sertifikat yang sama berlaku kata demi kata bagi sembarang peubah dengan E=50\E = 50 dan V=25V = 25, seeksotis apa pun, dan Contoh 22.18 menunjukkan bahwa suatu peubah semacam itu menjenuhkannya. Jadi keuniversalannya berharga; sedangkan ketika distribusinya sungguh binomial, perkakas eksponensial pada soal akhir pekan menutup sebagian besar senjangnya.

Contoh 22.23 (Korelasi sebuah bagian dengan keseluruhannya)

Untuk X,YX, Y yang saling bebas dan berdistribusi identik dengan ragam σ2>0\sigma^2 > 0, seberapa berkorelasikah satu sukunya dengan jumlahnya S=X+YS = X + Y? Hitunglah

Cov(X,S)=Cov(X,X)+Cov(X,Y)=σ2+0=σ2,V(S)=2σ2,\operatorname{Cov}(X, S) = \operatorname{Cov}(X, X) + \operatorname{Cov}(X, Y) = \sigma^2 + 0 = \sigma^2, \qquad V(S) = 2\sigma^2,

jadi koefisien korelasinya adalah

ρ(X,S)=Cov(X,S)σ(X)σ(S)=σ2σσ2=120.707,\rho(X, S) = \frac{\operatorname{Cov}(X, S)}{\sigma(X)\,\sigma(S)} = \frac{\sigma^2}{\sigma\cdot\sigma\sqrt2} = \frac{1}{\sqrt2} \approx 0.707 ,

apa pun hukum bersamanya — dadu, koin, atau cacah Poisson. Adapun dengan nn suku perhitungan yang sama memberi ρ(X1,Sn)=1/n\rho(X_1, S_n) = 1/\sqrt n: jadi pengaruh tiap suku individualnya pada totalnya mengencer seperti akar kuadrat, dan itulah bayangan korelasional bagi skala fluktuasi n\sqrt n. Adapun Cauchy–Schwarz menjamin ρ1\abs\rho \leq 1 selalu; dan di sini batasnya tercapai persis pada kasus merosot n=1n = 1 lalu meluruh secara terduga sesudahnya.

Contoh 22.24 (AM–GM terbobot dari Jensen)

Misalkan YY mengambil nilai positif a1,,aka_1, \dots, a_k dengan peluang λ1,,λk\lambda_1, \dots, \lambda_k. Adapun fungsi ln-\ln cembung pada (0,)\intoo0\infty, jadi Jensen memberi lnE(Y)E(lnY)-\ln\E(Y) \leq \E(-\ln Y), yakni

a1λ1a2λ2akλk    λ1a1+λ2a2++λkak:a_1^{\lambda_1}a_2^{\lambda_2}\cdots a_k^{\lambda_k} \;\leq\; \lambda_1a_1 + \lambda_2a_2 + \dots + \lambda_ka_k :

yaitu ketaksamaan aritmetika–geometri yang terbobot, dengan kesamaan bila dan hanya bila YY tetap. Adapun bobot yang sama λi=1k\lambda_i = \frac1k memulihkan AM–GM klasiknya. Jadi peluang diam-diam telah membuktikan sebuah teorema yang murni aljabar: sebab memilih sebuah hukum peluang tak lain alat tata buku bagi kombinasi cembung — yakni sudut pandang barisentrik Bab 17 sekali lagi, kini dengan Jensen sebagai mesinnya.

22.4 Hukum bilangan besar yang lemah

Teorema 22.25 (Hukum bilangan besar yang lemah)

Misalkan (Xk)k1(X_k)_{k \geq 1} peubah acak yang saling bebas berpasangan dan berhukum sama, serta bermomen kedua; lalu tulislah m=E(X1)m = \E(X_1) dan Sn=X1++XnS_n = X_1 + \dots + X_n. Maka untuk setiap ε>0\varepsilon > 0:

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 .

Bukti. Berkat kelinearannya E(Sn/n)=m\E(S_n/n) = m; lalu menurut Teorema 22.14 (sebab kebebasan berpasangannya membunuh kovariansinya) V(Sn)=nV(X1)V(S_n) = n\,V(X_1), jadi V(Sn/n)=V(X1)/nV(S_n/n) = V(X_1)/n. Sehingga ketaksamaan Chebyshev yang diterapkan pada Sn/nS_n/n memberi batasnya.

Catatan 22.26

Inilah teorema yang menghubungkan peluang dengan frekuensi: sebab untuk XkX_k berupa indikator sebuah kejadian AA pada pengulangan yang saling bebas, Sn/nS_n/n merupakan frekuensi AA yang teramati, dan hukum bilangan besar mengatakan bahwa ia memusat di sekitar P(A)\P(A) dengan laju p(1p)nε2\frac{p(1-p)}{n\varepsilon^2}. Adapun hukum yang kuat (yakni Sn/nmS_n/n \to m hampir pasti) merupakan teorema Tahun ke-3 — namun buktinya untuk momen keempat sudah terjangkau: lihat Latihan 22.9, yang menjalankan Borel–Cantelli pada batas bertipe Chebyshev. Adapun taksiran Chebyshev yang sama menggerakkan bukti polinomial Bernstein bagi teorema hampiran Weierstrass pada Bab 10 — sebab lema pencacahan di sana adalah hukum bilangan besar yang lemah yang menyamar.

Contoh 22.27 (Mengumpulkan lima puluh kupon)

Pengumpul kupon pada Latihan 22.3 dengan n=50n = 50 mainan berbeda: maka total yang diharapkan adalah

E(T50)=50H50=50k=1501k50×4.499225\E(T_{50}) = 50\,H_{50} = 50\sum_{k=1}^{50}\frac1k \approx 50 \times 4.499 \approx 225

kotak — yakni empat setengah kali tebakan naifnya 5050. Adapun pertumbuhan harmoniknya seluruh ceritanya: sebab 2525 mainan pertamanya tiba dalam sekitar 50ln23550\ln2 \approx 35 kotak, sedangkan mainan terakhirnya sendirian berharga 5050 kotak secara rata-rata (yakni tunggu geometri berparameter 150\frac1{50}). Jadi masalah penyelesaian didominasi babak akhirnya, dan itulah sebabnya Latihan 22.12 menemukan fluktuasi berorde nn — yakni besarnya tunggu geometri terakhir itu — di sekitar rerataannya nlnnn\ln n.

Contoh 22.28 (Seberapa besar nn harus?)

Untuk memaku frekuensi teramatinya dalam ε=0.01\varepsilon = 0.01 dari P(A)\P(A) dengan kepercayaan 95%95\%, batas Chebyshev menuntut

p(1p)nε214nε20.05,yaknin140.05(0.01)2=50000.\frac{p(1-p)}{n\varepsilon^2} \leq \frac{1}{4n\varepsilon^2} \leq 0.05, \qquad\text{yakni}\qquad n \geq \frac{1}{4\cdot0.05\cdot(0.01)^2} = 50\,000 .

Kebergantungannya brutal dalam ε\varepsilon (yakni kuadratik) dan lunak dalam kepercayaannya (yakni linear dalam 1/α1/\alpha). Adapun kedua ciri itu merupakan sifat batasnya, bukan sifat kebenarannya: sebab ketaksamaan eksponensial pada soal akhir pekan menurunkan harga kepercayaannya dari 1/α1/\alpha menjadi ln(1/α)\ln(1/\alpha) — jadi spesifikasi yang sama kelak berharga sekitar 1850018\,500 cuplikan di sana — sedangkan skala 1/ε21/\varepsilon^2-nya sejati dan tak terperbaiki. Jadi mengetahui bagian mana sebuah batas yang longgar sama bergunanya dengan batasnya sendiri.

Hukum bilangan besar sebagai gambar: sebab hukum S_n/n (yang digambar secara skematis) mempertahankan pusatnya m tetapi menyempit ketika n tumbuh, jadi peluang di luar pita [m- , m+ ] — yakni kedua ekornya — menciut ke nol. Adapun Chebyshev membatasi ekornya oleh V(X_1)/(n 2); sedangkan soal akhir pekan menunjukkan bahwa keduanya sebenarnya kecil secara eksponensial.
Hukum bilangan besar sebagai gambar: sebab hukum Sn/nS_n/n (yang digambar secara skematis) mempertahankan pusatnya mm tetapi menyempit ketika nn tumbuh, jadi peluang di luar pita [mε,m+ε]\intcc{m-\varepsilon}{m+\varepsilon} — yakni kedua ekornya — menciut ke nol. Adapun Chebyshev membatasi ekornya oleh V(X1)/(nε2)V(X_1)/(n\varepsilon^2); sedangkan soal akhir pekan menunjukkan bahwa keduanya sebenarnya kecil secara eksponensial.

Catatan 22.29 (Pandangan ke depan di dalam jilid ini)

Ke depan, segalanya di sini menyuapi Bab 23: sebab nilai harapan E(tX)\E(t^X) bagi satu fungsi XX yang cerdik memampatkan seluruh hukumnya menjadi sebuah deret pangkat, momennya menjadi turunan di 11, dan kesamaan bertipe Wald bagi jumlah acak mengusung teori proses bercabangnya; sedangkan teorema hasil kali bagi peubah yang saling bebas menjadi kemultiplikatifan fungsi pembangkitnya. Ke belakang, nilai harapan adalah barisentrum dengan bobot peluang (Bab 17), ketaksamaan Jensen adalah geometri garis penyangga bagi fungsi cembung (Bab 8), dan metode momen eksponensial pada soal akhir pekan bab ini adalah Markov yang diterapkan pada etX\eu^{tX} — jadi satu ketaksamaan, yang ditingkatkan oleh satu penggantian peubah yang baik, membentang tiga bab.

22.5 Latihan

Latihan 22.1

Hitunglah E(X)\E(X) dan V(X)V(X) untuk XB(n,p)X \sim \mathcal{B}(n, p) (lewat indikator), XP(λ)X \sim \mathcal{P}(\lambda) (tunjukkan V(X)=λV(X) = \lambda), dan XG(p)X \sim \mathcal{G}(p) (tunjukkan V(X)=1pp2V(X) = \frac{1-p}{p^2}; pakailah E(X(X1))\E(X(X-1)) beserta turunan kedua deret geometrinya).

Solusi

Solusi Latihan 22.1.

Binomial: X=i=1nXiX = \sum_{i=1}^n X_i dengan Bernoulli XiX_i yang saling bebas; V(Xi)=E(Xi2)E(Xi)2=pp2V(X_i) = \E(X_i^2) - \E(X_i)^2 = p - p^2, dan ragam peubah yang saling bebas itu menjumlah (Teorema 22.14):

E(X)=np,V(X)=np(1p).\E(X) = np, \qquad V(X) = np(1-p) .

Poisson: E(X(X1))=k2k(k1)eλλkk!=λ2eλj0λjj!=λ2\E\bigl(X(X-1)\bigr) = \sum_{k\geq2}k(k-1)e^{-\lambda}\frac{\lambda^k}{k!} = \lambda^2 e^{-\lambda}\sum_{j\geq0}\frac{\lambda^j}{j!} = \lambda^2, jadi

V(X)=E(X2)E(X)2=λ2+λλ2=λ.V(X) = \E(X^2) - \E(X)^2 = \lambda^2 + \lambda - \lambda^2 = \lambda .

Geometrik (q=1pq = 1 - p): dengan menurunkan k0qk=11q\sum_{k\geq0}q^k = \frac{1}{1-q} dua kali di dalam cakramnya (Bab 11), k2k(k1)qk2=2(1q)3\sum_{k\geq2}k(k-1)q^{k-2} = \frac{2}{(1-q)^3}, jadi

E(X(X1))=pqk2k(k1)qk2=2qp2,V(X)=2qp2+1p1p2=qp2=1pp2.\E\bigl(X(X-1)\bigr) = pq\sum_{k\geq2}k(k-1)q^{k-2} = \frac{2q}{p^2}, \qquad V(X) = \frac{2q}{p^2} + \frac1p - \frac{1}{p^2} = \frac{q}{p^2} = \frac{1-p}{p^2} .

Latihan 22.2

Misalkan XP(λ)X \sim \mathcal{P}(\lambda) dan YP(μ)Y \sim \mathcal{P}(\mu) saling bebas. Tunjukkan bahwa X+YP(λ+μ)X + Y \sim \mathcal{P}(\lambda + \mu) (lewat konvolusi bobotnya dan teorema binomial), dan bahwa hukum bersyarat XX jika diketahui X+Y=nX + Y = n bersifat binomial B(n,λλ+μ)\mathcal{B}\bigl(n, \frac{\lambda}{\lambda + \mu}\bigr).

Solusi

Solusi Latihan 22.2.

Jumlah: untuk nNn \in \N, menurut kesalinglepasan dan kesalingbebasannya,

P(X+Y=n)=k=0nP(X=k)P(Y=nk)=e(λ+μ)1n!k=0n(nk)λkμnk=e(λ+μ)(λ+μ)nn!\P(X + Y = n) = \sum_{k=0}^n \P(X = k)\P(Y = n - k) = e^{-(\lambda + \mu)}\frac{1}{n!} \sum_{k=0}^n \binom nk \lambda^k\mu^{n-k} = e^{-(\lambda+\mu)}\frac{(\lambda + \mu)^n}{n!}

menurut teorema binomial: X+YP(λ+μ)X + Y \sim \mathcal{P}(\lambda + \mu). Hukum bersyarat: untuk 0kn0 \leq k \leq n,

P(X=kX+Y=n)=P(X=k)P(Y=nk)P(X+Y=n)=(nk)(λλ+μ)k(μλ+μ)nk,\P(X = k \mid X + Y = n) = \frac{\P(X = k)\P(Y = n - k)}{\P(X + Y = n)} = \binom nk \Bigl(\frac{\lambda}{\lambda+\mu}\Bigr)^{k} \Bigl(\frac{\mu}{\lambda+\mu}\Bigr)^{n-k} ,

yakni hukum binomial B(n,λλ+μ)\mathcal{B}\bigl(n, \frac{\lambda}{\lambda+\mu}\bigr): bila cacah totalnya diketahui, setiap kejadian secara saling bebas “memilih” sumber pertamanya dengan peluang yang sebanding dengan lajunya.

Latihan 22.3

(Pengumpul kupon, nilai harapannya) Sebuah merek sereal menyembunyikan satu dari nn mainan berbeda, secara seragam, di tiap kotaknya. Misalkan TnT_n banyaknya kotak yang diperlukan untuk mengumpulkan seluruh nn mainannya. Dengan menulis TnT_n sebagai jumlah peubah geometri yang saling bebas (yakni waktu melihat mainan baru ketika kk masih hilang), tunjukkanlah

E(Tn)=nk=1n1knlnn\E(T_n) = n\sum_{k=1}^{n}\frac{1}{k} \sim n\ln n

(yang setara lewat pembandingan deret–integral pada Bab 6).

Solusi

Solusi Latihan 22.3.

Ketika kk mainan masih kurang, setiap kotak baru membawa mainan baru dengan peluang kn\frac kn, bebas dari masa lalunya: waktu tunggu WkW_k sampai mainan baru berikutnya berhukum geometrik G(kn)\mathcal{G}\bigl(\frac kn\bigr), dengan E(Wk)=nk\E(W_k) = \frac nk, dan Tn=Wn+Wn1++W1T_n = W_n + W_{n-1} + \dots + W_1 (kotak pertama selalu memberi mainan baru: Wn=1W_n = 1, sesuai dengan E=n/n\E = n/n). Menurut kelinearannya,

E(Tn)=k=1nnk=nk=1n1knlnn,\E(T_n) = \sum_{k=1}^n \frac nk = n\sum_{k=1}^n\frac1k \sim n\ln n ,

dengan memakai kn1k=lnn+γ+o(1)\sum_{k\leq n}\frac1k = \ln n + \gamma + o(1) (Bab 6). Mengumpulkan segelintir mainan terakhirlah yang mahal: separuh kotaknya habis untuk sisa yang sedikit itu.

Latihan 22.4 ★★

Misalkan X0X \geq 0 bernilai bulat. Buktikan rumus ekornya

E(X)=n=1P(Xn)\E(X) = \sum_{n=1}^{\infty} \P(X \geq n)

(ketika salah satu ruasnya berhingga), dengan menulis X=n11XnX = \sum_{n\geq1}\mathbf{1}_{X \geq n} lalu menukarkan penjumlahannya (lewat Fubini bagi keluarga taknegatif). Lalu perolehlah E(X)=1p\E(X) = \frac1p bagi hukum geometrinya.

Solusi

Solusi Latihan 22.4.

Secara titik demi titik, X(ω)=#{n1:X(ω)n}=n11Xn(ω)X(\omega) = \#\{n \geq 1 : X(\omega) \geq n\} = \sum_{n\geq1}\mathbf{1}_{X \geq n}(\omega). Keluarga rangkap (1Xn(ω)P({ω}))n,ω\bigl(\mathbf{1}_{X \geq n}(\omega)\,\P(\{\omega\})\bigr)_{n, \omega} taknegatif, jadi Fubini untuk keluarga (Bab 7) berlaku tanpa syarat: menjumlah dahulu atas nn memberi E(X)\E(X), menjumlah dahulu atas ω\omega memberi nP(Xn)\sum_n \P(X \geq n); keduanya berhingga serentak dan sama. Untuk XG(p)X \sim \mathcal{G}(p): P(Xn)=qn1\P(X \geq n) = q^{n-1} (q=1pq = 1-p), jadi E(X)=n1qn1=11q=1p\E(X) = \sum_{n\geq1}q^{n-1} = \frac{1}{1 - q} = \frac1p.

Latihan 22.5 ★★

(Pencuplikan tanpa pengembalian lebih memusat) Sebuah guci berisi NN bola, dan MM di antaranya putih. Tariklah nNn \leq N tanpa pengembalian lalu misalkan XX mencacah yang putihnya (yakni hukum hipergeometrik). Dengan memakai indikator X=i=1nYiX = \sum_{i=1}^n Y_i dengan YiY_i menyatakan tarikan ke-ii: tunjukkanlah bahwa tiap YiY_i bersifat Bernoulli berparameter p=M/Np = M/N (berkat kesimetriannya!), lalu simpulkan E(X)=np\E(X) = np persis seperti dengan pengembalian, kemudian tunjukkan Cov(Yi,Yj)=p(1p)N1<0\operatorname{Cov}(Y_i, Y_j) = -\frac{p(1-p)}{N-1} < 0 untuk iji \neq j, sehingga V(X)=np(1p)NnN1np(1p)V(X) = np(1-p)\frac{N - n}{N - 1} \leq np(1-p).

Solusi

Solusi Latihan 22.5.

Kesetangkupan: bola ke-ii yang terambil merupakan bola seragam acak dari gucinya (setiap dari NN bolanya sama mungkin menempati kedudukan ii pada urutan pengambilannya), jadi P(Yi=1)=MN=p\P(Y_i = 1) = \frac MN = p dan E(X)=np\E(X) = np menurut kelinearannya — tanpa perlu kesalingbebasan.

Kovariansi: untuk iji \neq j, E(YiYj)=P(pengambilan i,j sama-sama putih)=M(M1)N(N1)\E(Y_iY_j) = \P(\text{pengambilan } i, j \text{ sama-sama putih}) = \frac{M(M-1)}{N(N-1)} (pasangan terurut kedudukan yang berbeda memperoleh pasangan terurut bola yang berbeda, secara seragam). Jadi

Cov(Yi,Yj)=M(M1)N(N1)M2N2=M(NM)N21N1=p(1p)N1<0:\operatorname{Cov}(Y_i, Y_j) = \frac{M(M-1)}{N(N-1)} - \frac{M^2}{N^2} = \frac{M(N - M)}{N^2}\cdot\frac{-1}{N-1} = -\frac{p(1-p)}{N-1} < 0 :

mengambil satu bola putih membuat yang putih makin langka bagi pengambilan lainnya. Menurut Teorema 22.14,

V(X)=np(1p)+n(n1)(p(1p)N1)=np(1p)NnN1np(1p):V(X) = np(1-p) + n(n-1)\Bigl(-\frac{p(1-p)}{N-1}\Bigr) = np(1-p)\,\frac{N - n}{N - 1} \leq np(1-p) :

pencuplikan tanpa pengembalian punya rerata yang sama tetapi ragam yang lebih kecil daripada dengan pengembalian (kesamaannya hanya untuk n=1n = 1), dengan korelasi negatifnya berperan sebagai penstabil. Untuk n=Nn = N ragamnya lenyap: cacahnya lalu bersifat deterministik.

Latihan 22.6 ★★

Misalkan XX bermomen kedua. Tunjukkan bahwa cE((Xc)2)c \mapsto \E\bigl((X - c)^2\bigr) minimum persis di c=E(X)c = \E(X), dengan minimum V(X)V(X). Lalu tunjukkan bahwa P(X=E(X))=1\P(X = \E(X)) = 1 bila dan hanya bila V(X)=0V(X) = 0. (Untuk butir keduanya: bila V(X)=0V(X) = 0, pakailah Chebyshev dengan ε=1/n\varepsilon = 1/n beserta kekontinuan monoton, Teorema 21.6.)

Solusi

Solusi Latihan 22.6.

Dengan menguraikan di sekitar m=E(X)m = \E(X):

E((Xc)2)=E((Xm)2)+2(mc)E(Xm)+(mc)2=V(X)+(mc)2,\E\bigl((X - c)^2\bigr) = \E\bigl((X - m)^2\bigr) + 2(m - c)\,\E(X - m) + (m - c)^2 = V(X) + (m - c)^2 ,

yang minimum persis di c=mc = m dengan nilai V(X)V(X)nilai harapan merupakan penduga konstan terbaik dalam rerata kuadrat.

Bila P(X=m)=1\P(X = m) = 1 maka (Xm)2(X - m)^2 lenyap dengan peluang 11, jadi V(X)=0V(X) = 0 (keluarga pendefinisinya bersuku nol kecuali pada himpunan nol). Sebaliknya, bila V(X)=0V(X) = 0, Chebyshev (Teorema 22.15) memberi P(Xm1n)n2V(X)=0\P\bigl(\abs{X - m} \geq \frac1n\bigr) \leq n^2\,V(X) = 0 untuk setiap nn; kejadian {Xm1n}\bigl\{\abs{X - m} \geq \frac1n\bigr\} naik menuju {Xm}\{X \neq m\}, jadi kekontinuan monoton (Teorema 21.6) memberikan P(Xm)=0\P(X \neq m) = 0.

Latihan 22.7 ★★★

(Pemusatan mengalahkan Markov) Misalkan SnB(n,12)S_n \sim \mathcal{B}(n, \frac12) (yakni banyaknya gambar pada nn lemparan setimbang). Bandingkanlah batas yang diberikan Markov (yakni P(Sn3n4)\P(S_n \geq \frac{3n}{4})), Chebyshev, dan metode eksponensial (Chernoff):

P(Sn3n4)E(etSn)e3nt/4=(1+et2)ne3nt/4(t>0),\P\Bigl(S_n \geq \frac{3n}4\Bigr) \leq \E\bigl(e^{tS_n}\bigr)e^{-3nt/4} = \Bigl(\frac{1 + e^t}{2}\Bigr)^n e^{-3nt/4} \quad (t > 0),

lalu optimumkan tt untuk memperoleh batas yang kecil secara eksponensial. (Di t=ln3t = \ln 3: batasnya (233/4)n(0.877)n\bigl(2\cdot 3^{-3/4}\bigr)^n \approx (0.877)^n.)

Solusi

Solusi Latihan 22.7.

E(Sn)=n2\E(S_n) = \frac n2 dan V(Sn)=n4V(S_n) = \frac n4. Markov: P(Sn3n4)n/23n/4=23\P\bigl(S_n \geq \frac{3n}4\bigr) \leq \frac{n/2}{3n/4} = \frac23 — sebuah batas konstan, tak berguna untuk nn yang besar. Chebyshev: kejadiannya mengakibatkan Snn2n4\abs{S_n - \frac n2} \geq \frac n4, jadi peluangnya n/4(n/4)2=4n\leq \frac{n/4}{(n/4)^2} = \frac4n — meluruh, tetapi hanya secara suku banyak. Chernoff: menurut kesalingbebasannya, E(etSn)=i=1nE(etXi)=(1+et2)n\E(e^{tS_n}) = \prod_{i=1}^n\E(e^{tX_i}) = \bigl(\frac{1 + e^t}{2}\bigr)^n, dan Markov yang diterapkan pada etSne3nt/4e^{tS_n} \geq e^{3nt/4} memberi, untuk setiap t>0t > 0,

P(Sn3n4)(1+et2)ne3nt/4=exp(n(ln1+et23t4)).\P\Bigl(S_n \geq \frac{3n}4\Bigr) \leq \Bigl(\frac{1 + e^t}{2}\Bigr)^n e^{-3nt/4} = \exp\Bigl(n\bigl(\ln\tfrac{1 + e^t}{2} - \tfrac{3t}4\bigr)\Bigr).

Minimumkan pangkatnya:  ⁣d ⁣dtln1+et2=et1+et=34\frac{\dd}{\dd t}\ln\frac{1+e^t}{2} = \frac{e^t}{1 + e^t} = \frac34 di et=3e^t = 3, yakni t=ln3t = \ln 3, sehingga

P(Sn3n4)(42)n33n/4=(233/4)n(0.877)n,\P\Bigl(S_n \geq \frac{3n}4\Bigr) \leq \Bigl(\frac{4}{2}\Bigr)^n 3^{-3n/4} = \bigl(2 \cdot 3^{-3/4}\bigr)^n \approx (0.877)^n ,

yang kecil secara eksponensial. Hierarki Markov \to Chebyshev \to Chernoff merupakan tangga bakunya: setiap anak tangganya menerapkan Markov pada fungsi peubahnya yang tumbuh lebih cepat.

Latihan 22.8 ★★★

(Weierstrass lagi, secara peluang) Misalkan f ⁣:[0,1]Rf \colon [0,1] \to \R kontinu dan SnB(n,x)S_n \sim \mathcal{B}(n, x). Tunjukkan bahwa polinomial Bernstein Bnf(x)=k=0nf(kn)(nk)xk(1x)nkB_nf(x) = \sum_{k=0}^n f\bigl(\frac kn\bigr)\binom nk x^k(1-x)^{n-k} sama dengan E[f(Snn)]\E\bigl[f\bigl(\frac{S_n}{n}\bigr)\bigr], lalu turunkan ulang taksiran Bnf(x)f(x)ωf(δ)+2f4nδ2\abs{B_nf(x) - f(x)} \leq \omega_f(\delta) + \frac{2\norm f_\infty}{4n\delta^2} pada Bab 10 dalam bahasa peluang ini (yakni pecahlah pada Snnxδ\bigl|\frac{S_n}{n} - x\bigr| \geq \delta lalu pakailah Chebyshev).

Solusi

Solusi Latihan 22.8.

Menurut teorema pemindahan (Teorema 22.7) yang diterapkan pada f(Snn)f\bigl(\frac{S_n}{n}\bigr) dengan SnB(n,x)S_n \sim \mathcal{B}(n, x):

E[f(Snn)]=k=0nf(kn)(nk)xk(1x)nk=Bnf(x).\E\Bigl[f\Bigl(\frac{S_n}{n}\Bigr)\Bigr] = \sum_{k=0}^n f\Bigl(\frac kn\Bigr)\binom nk x^k(1-x)^{n-k} = B_nf(x) .

Tetapkan δ>0\delta > 0 lalu pisahkan f(Sn/n)f(x)\abs{f(S_n/n) - f(x)} pada kejadian D={Snnxδ}D = \bigl\{\abs{\frac{S_n}{n} - x} \geq \delta\bigr\}: di luar DD, selisihnya paling besar sama dengan modulus kekontinuannya ωf(δ)=supstδf(s)f(t)\omega_f(\delta) = \sup_{\abs{s - t}\leq\delta}\abs{f(s) - f(t)}; pada DD, paling besar 2f2\norm f_\infty. Dengan mengambil nilai harapannya lalu memakai Chebyshev dengan V(Snn)=x(1x)n14nV\bigl(\frac{S_n}{n}\bigr) = \frac{x(1-x)}{n} \leq \frac{1}{4n}:

Bnf(x)f(x)Ef(Sn/n)f(x)ωf(δ)+2fP(D)ωf(δ)+2f4nδ2.\abs{B_nf(x) - f(x)} \leq \E\,\abs{f(S_n/n) - f(x)} \leq \omega_f(\delta) + 2\norm f_\infty\,\P(D) \leq \omega_f(\delta) + \frac{2\norm f_\infty}{4n\delta^2} .

Kekontinuan seragam ff pada [0,1][0, 1] membuat ωf(δ)0\omega_f(\delta) \to 0: pilih δ\delta lalu nn, maka BnffB_nf \to f secara seragam — itulah teorema hampiran Weierstrass pada Bab 10, yang “lema pencacahannya” kini terkenali sebagai ketaksamaan Chebyshev bagi hukum binomialnya.

Latihan 22.9 ★★★

(Hukum kuat di bawah momen keempat) Misalkan (Xk)(X_k) saling bebas, berdistribusi identik, terpusat (dengan EX1=0\E X_1 = 0), dan E(X14)<\E(X_1^4) < \infty. Dengan menguraikan E(Sn4)\E(S_n^4) lalu mencacah suku yang bertahan (yakni hanya suku E(Xi4)\E(X_i^4) dan E(Xi2Xj2)\E(X_i^2X_j^2) dengan iji \neq j), tunjukkanlah E(Sn4)Cn2\E(S_n^4) \leq C n^2 untuk sebuah konstanta CC. Lalu turunkan nP(Sn/nε)<\sum_n \P\bigl(\abs{S_n/n} \geq \varepsilon\bigr) < \infty untuk tiap ε>0\varepsilon > 0 (lewat Markov pada orde 4) kemudian rampungkan dengan Borel–Cantelli (Teorema 21.25) bahwa Sn/n0S_n/n \to 0 hampir pasti dalam perumusan yang sesuai: yakni bahwa kejadian jNnN{Sn/n<1j}\bigcap_{j}\bigcup_N\bigcap_{n \geq N}\{\abs{S_n/n} < \frac1j\} berpeluang 11.

Solusi

Solusi Latihan 22.9.

Uraikan Sn4=i,j,k,lXiXjXkXlS_n^4 = \sum_{i,j,k,l}X_iX_jX_kX_l lalu ambil nilai harapannya. Menurut kesalingbebasan dan pemusatannya, setiap suku yang memuat satu indeks yang muncul tepat sekali akan lenyap (E(Xi)=0\E(X_i) = 0 terfaktorkan keluar). Suku yang bertahan: nn suku diagonal E(Xi4)\E(X_i^4), dan suku yang memasangkan dua pasang indeks yang sama, E(Xi2Xj2)=E(X12)2\E(X_i^2X_j^2) = \E(X_1^2)^2 untuk iji \neq j, yang muncul 3n(n1)3n(n-1) kali: pilih pasangan takterurut nilainya ((n2)\binom n2 cara), lalu 4!2!2!=6\frac{4!}{2!\,2!} = 6 cara menempatkannya pada keempat lubangnya — 6(n2)=3n(n1)6\binom n2 = 3n(n-1). Jadi, dengan E(X12)2E(X14)\E(X_1^2)^2 \leq \E(X_1^4) (Jensen atau Cauchy–Schwarz),

E(Sn4)=nE(X14)+3n(n1)E(X12)2Cn2,C=4E(X14).\E(S_n^4) = n\,\E(X_1^4) + 3n(n-1)\,\E(X_1^2)^2 \leq C n^2, \qquad C = 4\,\E(X_1^4) .

Markov pada orde 4:

P(Snnε)=P(Sn4n4ε4)Cn2n4ε4=Cn2ε4,\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}{n^2\varepsilon^4} ,

yakni deret yang terjumlahkan. Menurut Borel–Cantelli 1 (Teorema 21.25), untuk setiap jj kejadian Bj=lim supn{Sn/n1j}B_j = \limsup_n\bigl\{\abs{S_n/n} \geq \frac1j\bigr\} berpeluang 00, jadi P(jBj)=0\P\bigl(\bigcup_j B_j\bigr) = 0 menurut kesubaditifan terbilang. Pada komplemennya — yang berpeluang 11 — untuk setiap jj ada NN dengan Sn/n<1j\abs{S_n/n} < \frac1j bagi semua nNn \geq N: itu persis Sn/n0S_n/n \to 0. Hukum kuat bilangan besar berlaku di bawah momen keempat; menyingkirkan hipotesis itu (teorema Kolmogorov) merupakan pekerjaan Tahun ke-3.

Latihan 22.10

Dua dadu setimbang dilempar; misalkan MM yang lebih besar di antara kedua hasilnya. Dengan memakai rumus ekor pada Latihan 22.4 (versi berhingganya), tunjukkanlah

E(M)=k=16P(Mk)=6j=05(j6)2=161364.47.\E(M) = \sum_{k=1}^{6}\P(M \geq k) = 6 - \sum_{j=0}^5\Bigl(\frac j6\Bigr)^2 = \frac{161}{36} \approx 4.47 .
Solusi

Solusi Latihan 22.10.

P(Mk)=(k6)2\P(M \leq k) = \bigl(\frac k6\bigr)^2 (kedua dadunya paling besar kk, secara saling bebas), jadi P(Mk)=1(k16)2\P(M \geq k) = 1 - \bigl(\frac{k-1}6\bigr)^2 dan

E(M)=k=16P(Mk)=60+1+4+9+16+2536=65536=161364.47,\E(M) = \sum_{k=1}^6\P(M \geq k) = 6 - \frac{0 + 1 + 4 + 9 + 16 + 25}{36} = 6 - \frac{55}{36} = \frac{161}{36} \approx 4.47 ,

yang cukup jauh di atas rerata 3.53.5 satu dadu tunggal, sebagaimana layaknya sebuah maksimum.

Latihan 22.11 ★★

Misalkan FnF_n banyaknya titik tetap sebuah permutasi acak seragam atas {1,,n}\{1, \dots, n\} (dengan n2n \geq 2). Dengan menulis Fn=i1σ(i)=iF_n = \sum_i\mathbf 1_{\sigma(i) = i}, hitunglah E(Fn)=1\E(F_n) = 1 dan Cov(1σ(i)=i,1σ(j)=j)=1n2(n1)\operatorname{Cov}(\mathbf 1_{\sigma(i)=i}, \mathbf 1_{\sigma(j)=j}) = \frac1{n^2(n-1)} untuk iji \neq j, lalu simpulkan V(Fn)=1V(F_n) = 1: jadi rata-rata satu huruf tetap, dengan ragam persis 11, berapa pun nn.

Solusi

Solusi Latihan 22.11.

Dengan Ii=1σ(i)=iI_i = \mathbf 1_{\sigma(i) = i}: P(σ(i)=i)=(n1)!n!=1n\P(\sigma(i) = i) = \frac{(n-1)!}{n!} = \frac1n, jadi E(Fn)=n1n=1\E(F_n) = n\cdot\frac1n = 1. Untuk iji \neq j: P(σ(i)=i,σ(j)=j)=(n2)!n!=1n(n1)\P(\sigma(i) = i, \sigma(j) = j) = \frac{(n-2)!}{n!} = \frac1{n(n-1)}, sehingga

Cov(Ii,Ij)=1n(n1)1n2=1n2(n1).\operatorname{Cov}(I_i, I_j) = \frac1{n(n-1)} - \frac1{n^2} = \frac{1}{n^2(n-1)} .

Menurut perkakas ragamnya (Teorema 22.14),

V(Fn)=n1n(11n)+n(n1)1n2(n1)=11n+1n=1.V(F_n) = n\cdot\frac1n\Bigl(1 - \frac1n\Bigr) + n(n-1)\cdot\frac1{n^2(n-1)} = 1 - \frac1n + \frac1n = 1 .

Rerata 11, ragam 11, tak bergantung pada nn — sesuai dengan limit Poisson pada masalah pencocokannya (Latihan 21.5).

Latihan 22.12 ★★★

(Pengumpul kupon, pemusatannya) Pada latar Latihan 22.3, tunjukkanlah

V(Tn)=k=1n1k/n(k/n)2n2k=1n1k2π26n2,V(T_n) = \sum_{k=1}^n\frac{1 - k/n}{(k/n)^2} \leq n^2\sum_{k=1}^n\frac{1}{k^2} \leq \frac{\pi^2}{6}n^2,

dengan memakai kebebasan tahap geometrinya beserta V(G(p))=1pp2V(\mathcal G(p)) = \frac{1-p}{p^2} (Latihan 22.1; dan nilai π2/6\pi^2/6-nya adalah Contoh 14.12). Lalu turunkan dengan Chebyshev bahwa Tnnlnn1\dfrac{T_n}{n\ln n} \to 1 dalam peluang: jadi total waktu pengumpulnya adalah nlnnn\ln n hingga fluktuasi berorde nn.

Solusi

Solusi Latihan 22.12.

Tn=k=1nGkT_n = \sum_{k=1}^nG_k dengan GkG(k/n)G_k \sim \mathcal G(k/n) menyatakan waktu sampai terlihat mainan baru ketika kk mainan masih kurang, dan tahapannya saling bebas. Jadi

V(Tn)=k=1n1k/n(k/n)2k=1nn2k2π26n2,V(T_n) = \sum_{k=1}^n\frac{1 - k/n}{(k/n)^2} \leq \sum_{k=1}^n\frac{n^2}{k^2} \leq \frac{\pi^2}6\,n^2 ,

menurut Contoh 14.12. Dengan E(Tn)=nHn\E(T_n) = nH_n, Hn=1n1kH_n = \sum_1^n\frac1k (Latihan 22.3), Chebyshev memberi, untuk ε>0\varepsilon > 0,

P(TnnHnεnlnn)π2n2/6ε2n2ln2n=π26ε2ln2nn0.\P\bigl(\abs{T_n - nH_n} \geq \varepsilon\,n\ln n\bigr) \leq \frac{\pi^2n^2/6}{\varepsilon^2n^2\ln^2 n} = \frac{\pi^2}{6\,\varepsilon^2\ln^2n} \xrightarrow[n\to\infty]{} 0 .

Karena HnlnnH_n \sim \ln n, membaginya dengan nlnnn\ln n menunjukkan Tn/(nlnn)1T_n/(n\ln n) \to 1 dalam peluang: fluktuasi TnT_n berorde nn, yang dapat diabaikan terhadap reratanya nlnnn\ln n.

22.6 Soal: kotak perkakas pemusatan, dari Markov ke Hoeffding

Soal 22.1

Soal akhir pekan — pemusatan eksponensial dengan tangan, dan berapa orang yang harus ditanya sebuah jajak pendapat

Ketaksamaan Markov berharga satu momen dan membeli peluruhan 1/a1/a; sedangkan Chebyshev berharga dua momen dan membeli 1/ε21/\varepsilon^2 — dan Contoh 22.18 menunjukkan bahwa itulah semua yang dapat dibeli momen tersebut. Soal ini mendaki sisa tangganya: yakni metode eksponensial (Chernoff) dengan lajunya yang eksak bagi lemparan koin, ketaksamaan Hoeffding bagi semua peubah yang terbatas, dan panennya — yakni ukuran cuplikan yang gamblang dan jujur bagi jajak pendapat, penetapan pemenang pemilu, dan pengujian koin. Di sepanjang soal ini, SnB(n,p)S_n \sim \mathcal B(n, p) merupakan jumlah nn peubah Bernoulli yang saling bebas dan p^n=Sn/n\widehat p_n = S_n/n frekuensi empirisnya.

Bagian I — Penakaran pada koin yang setimbang. Di sini p=12p = \frac12 dan a(12,1)a \in \intoo{\frac12}{1}.

  1. Markov di aras anan: tunjukkan P(Snan)12a\P(S_n \geq an) \leq \frac1{2a}, yakni batas yang bahkan tak menuju 00. Lalu di mana Markov kehilangan sebanyak itu?
  2. Chebyshev: dengan memakai kesetangkupan binomial setimbangnya terhadap n/2n/2, tunjukkanlah

    P(Snan)=12P(Snn2n(a12))18n(a1/2)2,\P(S_n \geq an) = \tfrac12\, \P\bigl(\abs{S_n - \tfrac n2} \geq n(a - \tfrac12)\bigr) \leq \frac{1}{8n(a - 1/2)^2},

    yakni 2n\frac2n di a=34a = \frac34: jadi peluruhan polinomial pada akhirnya.

  3. (Chernoff, aras umum) Hitunglah E(etSn)=(1+et2)n\E(\eu^{tS_n}) = \bigl(\frac{1 + \eu^t}2\bigr)^n lalu optimumkan P(Snan)E(etSn)etan\P(S_n \geq an) \leq \E(\eu^{tS_n})\eu^{-tan} atas t>0t > 0: tunjukkan bahwa tt optimumnya adalah lna1a\ln\frac{a}{1-a} dan

    P(Snan)enI(a),I(a)=ln2+alna+(1a)ln(1a)>0.\P(S_n \geq an) \leq \eu^{-n\,I(a)}, \qquad I(a) = \ln 2 + a\ln a + (1-a)\ln(1-a) > 0 .

    Lalu periksalah bahwa a=34a = \frac34 memulihkan batas (233/4)n\bigl(2\cdot3^{-3/4}\bigr)^n pada Latihan 22.7.

  4. (Eksponennya eksak) Misalkan k=ank = an bilangan bulat. Dari fakta bahwa (nk)ak(1a)nk\binom nk a^k(1-a)^{n-k} merupakan yang terbesar di antara n+1n + 1 suku sebuah distribusi peluang, buktikanlah (nk)enH(a)n+1\binom nk \geq \frac{\eu^{nH(a)}}{n+1} dengan H(a)=alna(1a)ln(1a)H(a) = -a\ln a - (1-a)\ln(1-a), lalu turunkan batas bawah yang bersesuaian

    P(Snan)(nan)2nenI(a)n+1.\P(S_n \geq an) \geq \binom{n}{an}2^{-n} \geq \frac{\eu^{-n\,I(a)}}{n + 1} .
  5. Tabelkan ketiga batasnya di n=100n = 100 dan a=34a = \frac34: yakni Markov 23\frac23, Chebyshev 0.020.02, dan Chernoff 2.1106\approx 2.1\cdot10^{-6} (sedangkan nilai benarnya 2.8107\approx 2.8\cdot10^{-7}). Lalu pelajarannya, dalam satu kalimat?

Bagian II — Ketaksamaan Hoeffding.

  1. (Kasus Rademacher) Untuk ε=±1\varepsilon = \pm1 yang masing-masing berpeluang 12\frac12, buktikanlah

    E(etε)=coshtet2/2(tR)\E(\eu^{t\varepsilon}) = \cosh t \leq \eu^{t^2/2} \qquad (t \in \R)

    dengan membandingkan kedua deretnya suku demi suku (sebab (2k)!2kk!(2k)! \geq 2^kk!).

  2. Turunkan, untuk peubah Rademacher yang saling bebas ε1,,εn\varepsilon_1, \dots, \varepsilon_n dan setiap s>0s > 0:

    P(i=1nεis)es2/(2n).\P\Bigl(\sum_{i=1}^n\varepsilon_i \geq s\Bigr) \leq \eu^{-s^2/(2n)} .
  3. Terjemahkan ke koin setimbang (dengan Xi=1+εi2X_i = \frac{1+\varepsilon_i}2): P(p^n12δ)e2nδ2\P\bigl(\widehat p_n - \tfrac12 \geq \delta\bigr) \leq \eu^{-2n\delta^2}, beserta versi dua sisinya dengan faktor 22.
  4. (Lema Hoeffding) Misalkan X[0,1]X \in \intcc01 dengan EX=p\E X = p, dan ψ(t)=lnE(etX)\psi(t) = \ln\E(\eu^{tX}). Benarkanlah bahwa ψ\psi terdiferensialkan dua kali dengan

    ψ(t)=Et(X2)Et(X)2,Et(Y):=E(YetX)E(etX),\psi''(t) = \E_t(X^2) - \E_t(X)^2, \qquad \E_t(Y) := \frac{\E(Y\eu^{tX})}{\E(\eu^{tX})},

    yakni sebuah ragam bagi peubah terbobot ulang yang masih bernilai di [0,1]\intcc01; lalu batasilah ia dengan 14\frac14 (lewat hujah keminimalan Latihan 22.6) kemudian rampungkan lewat Taylor:

    E(et(Xp))et2/8.\E\bigl(\eu^{t(X - p)}\bigr) \leq \eu^{t^2/8} .
  5. (Ketaksamaan Hoeffding) Untuk Xi[0,1]X_i \in \intcc01 yang saling bebas dengan rerataan bersama pp, turunkanlah

    P(p^npδ)2e2nδ2(δ>0).\P\bigl(\abs{\widehat p_n - p} \geq \delta\bigr) \leq 2\,\eu^{-2n\delta^2} \qquad (\delta > 0).
  6. Bandingkan laju Chebyshev p(1p)nδ2\frac{p(1-p)}{n\delta^2} dengan laju Hoeffding 2e2nδ22\eu^{-2n\delta^2}: hipotesis mana yang dituntut masing-masingnya, dan mulai dari nn yang mana (secara kasar) batas eksponensialnya menang di δ=0.03\delta = 0.03 dan p=12p = \frac12?

Bagian III — Berapa orang yang harus ditanya sebuah jajak pendapat? Sebuah jajak pendapat menanyai nn pemilih yang saling bebas dan dipilih seragam; masing-masing menjawab jujur; dengan pp menyatakan skor benarnya dan p^n\widehat p_n angka jajak pendapatnya.

  1. Tunjukkan bahwa jajak pendapatnya akurat sampai ±δ\pm\delta dengan kepercayaan 1α1 - \alpha (yakni P(p^npδ)α\P(\abs{\widehat p_n - p} \geq \delta) \leq \alpha) segera setelah

    n    ln(2/α)2δ2.n \;\geq\; \frac{\ln(2/\alpha)}{2\,\delta^2} .
  2. Hitunglah nn yang dituntut bagi spesifikasi baku “tiga poin, sembilan puluh lima persen” (dengan δ=0.03\delta = 0.03 dan α=0.05\alpha = 0.05): yakni n2050n \geq 2050; dan bagi satu poin: n18445n \geq 18\,445. Amatilah — lalu jelaskanlah — fakta yang mencolok bahwa jawabannya tak melibatkan besarnya populasi.
  3. Kerjakan ulang pertanyaan 13 lewat Chebyshev (dengan V(X1)=p(1p)14V(X_1) = p(1-p) \leq \frac14): yakni n14αδ2=5556n \geq \frac1{4\alpha\delta^2} = 5556 pada tiga poin. Perhatikan bahwa pencuplikan tanpa pengembalian hanya membantu (Latihan 22.5: sebab ragamnya menciut sebesar NnN1\frac{N-n}{N-1}).
  4. (Menetapkan pemenang pemilu) Skor benar seorang kandidat adalah p=0.52p = 0.52. Berapa banyak pemilih yang harus dijajaki agar P(p^n12)0.01\P(\widehat p_n \leq \tfrac12) \leq 0.01? Tunjukkanlah nln1002(0.02)25757n \geq \frac{\ln 100}{2\cdot(0.02)^2} \approx 5757 — jadi menetapkan pemenang perlombaan yang ketat jauh lebih mahal daripada menaksir sebuah skor.
  5. Apa yang tak diliput matematikanya: daftarkanlah asumsi pemodelan yang dipakai (yakni pencuplikan seragam yang saling bebas, jawaban yang jujur, dan pp yang tetap), lalu jelaskanlah dalam satu paragraf pendek mengapa galat jajak pendapat yang nyata didominasi bias (yakni pencuplikan yang tak seragam dan ketakresponsan), yang tak dikurangi oleh kenaikan nn mana pun.

Bagian IV — Lebih tajam dan lebih murah.

  1. (Median rerata: peluruhan eksponensial dari dua momen) Pecahlah sebuah anggaran berisi kmkm cuplikan menjadi kk kelompok saling bebas berisi mm; lalu misalkan p^(1),,p^(k)\widehat p^{(1)}, \dots, \widehat p^{(k)} rerata kelompoknya dan MM mediannya. Pilihlah mm sedemikian sehingga tiap kelompoknya memenuhi P(p^(i)pδ)18\P(\abs{\widehat p^{(i)} - p} \geq \delta) \leq \frac18 (lewat Chebyshev: sebab m2δ2m \geq \frac2{\delta^2} sudah cukup). Tunjukkanlah bahwa bila Mpδ\abs{M - p} \geq \delta maka sekurangnya k/2k/2 kelompoknya keliru, lalu turunkan

    P(Mpδ)(kk/2)(18)k/22k8k/2=2k/2:\P(\abs{M - p} \geq \delta) \leq \binom{k}{\lceil k/2\rceil}\Bigl(\frac18 \Bigr)^{k/2} \leq 2^k\cdot 8^{-k/2} = 2^{-k/2} :

    jadi pemusatan eksponensial dengan memakai apa pun yang tak melampaui ragamnya.

  2. (Paley–Zygmund) Untuk X0X \geq 0 yang bermomen kedua, buktikanlah P(X>0)E(X)2E(X2)\P(X > 0) \geq \dfrac{\E(X)^2}{\E(X^2)} (lewat Cauchy–Schwarz pada X1X>0X\mathbf 1_{X>0}): yakni perkakas arah baliknya — sebab momen juga dapat memaksa kejadian untuk terjadi.
  3. (Pinsker ringan) Tunjukkan I(a)2(a12)2I(a) \geq 2\bigl(a - \tfrac12\bigr)^2 pada (12,1)\intoo{\frac12}1 (sebab selisihnya lenyap sampai orde dua di 12\frac12 dan turunan keduanya 1a(1a)40\frac1{a(1-a)} - 4 \geq 0): jadi eksponen eksak Chernoff selalu mengalahkan eksponen kuadratik Hoeffding.
  4. Uraikan I(12+δ)=2δ2+O(δ4)I\bigl(\tfrac12 + \delta\bigr) = 2\delta^2 + O(\delta^4) lalu gabungkan dengan pertanyaan 4: jadi untuk simpangan yang kecil eksponen Hoeffding 2nδ22n\delta^2 bersifat eksak secara asimtotik — sehingga tak ada metode yang dapat mengalahkannya lebih daripada faktor polinomial.
  5. Susunlah tabel kotak perkakasnya: yakni bagi Markov, Chebyshev, batas momen keempat pada Latihan 22.9, Hoeffding, dan Chernoff dengan eksponen II, nyatakanlah dalam satu baris masing-masing: hipotesis yang dituntut, peluruhan yang diperoleh, dan pertanyaan pada soal ini tempat ia paling tajam.

Bagian V — Panennya.

  1. (Menguji sebuah koin) Sebuah koin entah setimbang entah berat sebelah dengan p=0.55p = 0.55. Kamu melemparnya nn kali lalu mengumumkan “berat sebelah” ketika p^n>0.525\widehat p_n > 0.525. Tunjukkanlah bahwa kedua peluang galatnya paling banyak e2n(0.025)2\eu^{-2n(0.025)^2}, dan bahwa n3685n \geq 3685 lemparan menjamin keduanya di bawah 1%1\%.
  2. (Kejadian langka menuntut batas yang sadar ragam) Misalkan p=0.01p = 0.01 lalu ambillah spesifikasi relatifnya δ=p/2=0.005\delta = p/2 = 0.005 dan α=0.05\alpha = 0.05. Bandingkanlah ukuran cuplikan yang dituntut Hoeffding (yakni n74000n \approx 74\,000) dan yang dituntut Chebyshev dengan ragam benarnya p(1p)p(1-p) (yakni n7920n \approx 7920): jadi batas eksponensial yang buta ragam kalah oleh momen kedua yang sederhana. Nyatakanlah pelajarannya, dan dari mana perkakas yang hilangnya (yakni batas eksponensial yang sadar ragam; dan hampiran Poisson pada Bab 23) kelak datang.
  3. (Hukum kuat bagi koin) Dari n2e2nδ2<\sum_n 2\eu^{-2n\delta^2} < \infty dan Borel–Cantelli (Teorema 21.25), buktikanlah bahwa p^np\widehat p_n \to p hampir pasti bagi lemparan koin yang saling bebas: yakni rumuskanlah kejadian hampir pastinya sebagai jNnN{p^np<1j}\bigcap_j\bigcup_N\bigcap_{n\geq N} \{\abs{\widehat p_n - p} < \tfrac1j\} seperti pada Latihan 22.9, lalu simpulkan. (Adapun keterbatasannya menggantikan momen keempat yang dipakai di sana.)
  4. Rangkuman. Dalam lima kalimat: apa yang diharga dan dibeli tiap anak tangganya (yakni momen satu, dua, empat; eksponensial yang terbatas; dan eksponen yang eksak); mengapa menjajaki 20502050 orang sudah cukup bagi sebuah negeri sebesar apa pun; dan batas mana di antaranya yang kelak dipertajam jilid Tahun ke-3 menjadi konstanta eksak teorema limit pusat.
Solusi

Solusi Soal 22.1.

1. E(Sn)=n2\E(S_n) = \frac n2 dan Markov (Teorema 22.15) memberi P(Snan)n/2an=12a\P(S_n \geq an) \leq \frac{n/2}{an} = \frac1{2a}. Markov hanya mengenal reratanya: ia tak dapat membedakan peubah yang memusat di n/2n/2 dari peubah yang tersebar antara 00 dan nn, jadi ia menghargai ekornya seolah-olah seluruh massanya boleh berada di sana.

2. Binomial setimbangnya setangkup terhadap n/2n/2 (SnS_n dan nSnn - S_n berhukum sama), jadi dengan x=n(a12)>0x = n(a - \frac12) > 0 kedua kejadian {Snn2x}\{S_n - \frac n2 \geq x\} dan {Snn2x}\{S_n - \frac n2 \leq -x\} saling lepas dan sama peluangnya: P(Snan)=12P(Snn2x)\P(S_n \geq an) = \frac12\P(\abs{S_n - \frac n2} \geq x). Chebyshev dengan V(Sn)=n4V(S_n) = \frac n4:

P(Snan)12n/4n2(a1/2)2=18n(a1/2)2,\P(S_n \geq an) \leq \frac12\cdot\frac{n/4}{n^2(a - 1/2)^2} = \frac1{8n(a - 1/2)^2},

yang bernilai 2n\frac2n di a=34a = \frac34.

3. Menurut kesalingbebasan dan teorema hasil kalinya, E(etSn)=(EetX1)n=(1+et2)n\E(\eu^{tS_n}) = \bigl(\E \eu^{tX_1}\bigr)^n = \bigl(\frac{1 + \eu^t}2\bigr)^n. Markov yang diterapkan pada etSn\eu^{tS_n}:

P(Snan)etan(1+et2) ⁣n=exp(n(ln1+et2ta)).\P(S_n \geq an) \leq \eu^{-tan}\Bigl(\frac{1 + \eu^t}2\Bigr)^{\!n} = \exp\Bigl(n\bigl(\ln\tfrac{1 + \eu^t}2 - ta\bigr)\Bigr).

Turunan pangkatnya terhadap tt adalah et1+eta\frac{\eu^t}{1 + \eu^t} - a, yang lenyap di et=a1a\eu^t = \frac a{1-a}, yakni t=lna1a>0t^* = \ln\frac a{1-a} > 0; di sana 1+et2=12(1a)\frac{1 + \eu^{t^*}}2 = \frac1{2(1-a)} dan pangkatnya sama dengan

n(ln2ln(1a)alna1a)=n(ln2+alna+(1a)ln(1a))=nI(a),n\Bigl(-\ln 2 - \ln(1-a) - a\ln\frac a{1-a}\Bigr) = -n\bigl(\ln2 + a\ln a + (1-a)\ln(1-a)\bigr) = -n\,I(a),

dengan I(12)=0I(\frac12) = 0 dan I(a)=lna1a>0I'(a) = \ln\frac a{1-a} > 0 pada (12,1)\intoo{\frac12}1: I(a)>0I(a) > 0. Di a=34a = \frac34: eI(3/4)=12(34)3/4(14)1/4=233/4\eu^{-I(3/4)} = \frac12(\tfrac34)^{-3/4}(\tfrac14)^{-1/4} = 2\cdot3^{-3/4}, yakni batas pada Latihan 22.7.

4. Sebanyak n+1n + 1 bilangan (nj)aj(1a)nj\binom nja^j(1-a)^{n-j} berjumlah 11, dan yang terbesar adalah yang di j=k=anj = k = an (modus B(n,a)\mathcal B(n, a) di sini adalah (n+1)a=k\floor{(n+1)a} = k). Maksimum dari n+1n + 1 bilangan yang berjumlah 11 paling kecil adalah 1n+1\frac1{n+1}:

(nk)ak(1a)nk1n+1(nk)aan(1a)n(1a)n+1=enH(a)n+1.\binom nk a^k(1-a)^{n-k} \geq \frac1{n+1} \quad\Longrightarrow\quad \binom nk \geq \frac{a^{-an}(1-a)^{-n(1-a)}}{n+1} = \frac{\eu^{nH(a)}}{n+1}.

Jadi P(Snan)(nan)2nen(H(a)ln2)/(n+1)=enI(a)/(n+1)\P(S_n \geq an) \geq \binom{n}{an}2^{-n} \geq \eu^{n(H(a) - \ln2)}/(n+1) = \eu^{-nI(a)}/(n+1): sampai faktor suku banyaknya n+1n + 1, pangkat Chernoff itulah kebenarannya.

5. n=100n = 100, a=34a = \frac34: Markov 23\frac23; Chebyshev 2100=0.02\frac2{100} = 0.02; Chernoff (233/4)100=e100I(3/4)2.1106(2\cdot3^{-3/4})^{100} = \eu^{-100\,I(3/4)} \approx 2.1\cdot10^{-6}, terhadap nilai persisnya 2.81072.8\cdot10^{-7}. Pelajarannya: setiap momen informasi membagi batasnya secara suku banyak; adapun momen eksponensialnya mengubah sifatnya.

6. cosht=k0t2k(2k)!\cosh t = \sum_{k\geq0}\frac{t^{2k}}{(2k)!} dan et2/2=k0t2k2kk!\eu^{t^2/2} = \sum_{k\geq0}\frac{t^{2k}}{2^kk!}; klaimnya menyusul suku demi suku dari (2k)!2kk!(2k)! \geq 2^kk!, yang berlaku menurut induksi: (2k)!=2k(2k1)(2k2)!2k2k1(k1)!=2kk!(2k1)2kk!(2k)! = 2k(2k-1)\cdot(2k-2)! \geq 2k\cdot 2^{k-1}(k-1)! = 2^kk!\cdot(2k-1) \geq 2^kk!.

7. Menurut kesalingbebasannya, E(etεi)=(cosht)nent2/2\E\bigl(\eu^{t\sum\varepsilon_i} \bigr) = (\cosh t)^n \leq \eu^{nt^2/2}, jadi Markov memberi P(εis)ent2/2ts\P(\sum\varepsilon_i \geq s) \leq \eu^{nt^2/2 - ts}; dengan meminimumkannya di t=s/nt = s/n diperoleh es2/(2n)\eu^{-s^2/(2n)}.

8. Dengan Xi=1+εi2X_i = \frac{1 + \varepsilon_i}2, p^n12=12nεi\widehat p_n - \frac12 = \frac1{2n}\sum\varepsilon_i, jadi {p^n12δ}={εi2nδ}\{\widehat p_n - \frac12 \geq \delta\} = \{\sum\varepsilon_i \geq 2n\delta\} dan pertanyaan 7 memberi batas e(2nδ)2/(2n)=e2nδ2\eu^{-(2n\delta)^2/(2n)} = \eu^{-2n\delta^2}. Kejadian setangkupnya berbatas sama, dari situlah faktor 22 bagi p^n12δ\abs{\widehat p_n - \frac12} \geq \delta.

9. E(etX)=xetxP(X=x)\E(\eu^{tX}) = \sum_x\eu^{tx}\P(X = x) merupakan deret fungsi mulus dari tt yang turunan suku demi sukunya terdominasi, pada setiap selang-tt kompak, oleh etP(X=x)\eu^{\abs t}\P(X = x) (sebab 0x10 \leq x \leq 1): menurut teorema penurunan bagi deret yang konvergen normal (Teorema 10.7) ia dua kali terdiferensialkan, dan aturan hasil baginya memberi ψ=Et(X)\psi' = \E_t(X) dan ψ=Et(X2)Et(X)2\psi'' = \E_t(X^2) - \E_t(X)^2, dengan Et\E_t menyatakan nilai harapan bagi bobot yang ditimbang ulang etxP(X=x)/E(etX)\eu^{tx}\P(X{=}x)/\E(\eu^{tX}) — taknegatif, berjumlah 11, ditopang oleh nilai yang sama x[0,1]x \in \intcc01. Ragam sebuah peubah bernilai-[0,1]\intcc01 paling besar 14\frac14: menurut Latihan 22.6, ia sama dengan mincEt((Xc)2)Et((X12)2)14\min_c\E_t((X - c)^2) \leq \E_t\bigl((X - \tfrac12)^2\bigr) \leq \tfrac14. Taylor dengan sisa integral, memakai ψ(0)=0\psi(0) = 0, ψ(0)=p\psi'(0) = p:

ψ(t)=tp+0t(ts)ψ(s) ⁣dstp+t2214,\psi(t) = tp + \int_0^t(t - s)\,\psi''(s)\,\dd s \leq tp + \frac{t^2}2\cdot\frac14,

yakni E(et(Xp))et2/8\E(\eu^{t(X - p)}) \leq \eu^{t^2/8} untuk semua tt real.

10. Menurut kesalingbebasannya, E(et(Snnp))ent2/8\E\bigl(\eu^{t(S_n - np)}\bigr) \leq \eu^{nt^2/8}; Markov dan pengoptimuman t=4δt = 4\delta memberi

P(p^npδ)ent2/8tnδt=4δ=e2nδ2;\P(\widehat p_n - p \geq \delta) \leq \eu^{nt^2/8 - tn\delta}\Big|_{t = 4\delta} = \eu^{-2n\delta^2};

adapun penerapannya pada peubah 1Xi1 - X_i (yang juga di [0,1]\intcc01) membatasi ekor lainnya, dari situlah batas dua sisi 2e2nδ22\eu^{-2n\delta^2}.

11. Chebyshev hanya perlu momen kedua lalu memberi p(1p)nδ2\frac{p(1-p)}{n\delta^2}; Hoeffding perlu keterbatasan lalu memberi 2e2nδ22\eu^{-2n\delta^2}. Di p=12p = \frac12, δ=0.03\delta = 0.03: batasnya adalah 278n\frac{278}{n} (kira-kira) berbanding 2e0.0018n2\eu^{-0.0018n}; keduanya berpotongan dekat n1200n \approx 1200, dan setelah itu batas eksponensialnya menang, bahkan menang telak (n=5000n = 5000: 0.0560.056 berbanding 2.51042.5\cdot10^{-4}).

12. Menurut Hoeffding (pertanyaan 10), P(p^npδ)2e2nδ2α\P(\abs{\widehat p_n - p} \geq \delta) \leq 2\eu^{-2n\delta^2} \leq \alpha segera setelah 2nδ2ln2α2n\delta^2 \geq \ln\frac2\alpha, yakni nln(2/α)2δ2n \geq \frac{\ln(2/\alpha)}{2\delta^2}.

13. δ=0.03\delta = 0.03, α=0.05\alpha = 0.05: nln4020.00092049.4n \geq \frac{\ln 40}{2\cdot0.0009} \approx 2049.4: 20502050 orang. Untuk δ=0.01\delta = 0.01: nln400.000218445n \geq \frac{\ln40}{0.0002} \approx 18\,445. Ukuran populasinya tak pernah muncul sebab setiap pemilih yang tercuplik dimodelkan sebagai pengambilan Bernoulli(p)(p) yang baru: kesulitan jajak pendapatnya terletak pada ragam sekeping koin, bukan pada besarnya negerinya. Memotong separuh galat batasnya berharga empat kali lipat cuplikannya — itulah hukum 1/δ21/\delta^2.

14. Chebyshev: P(p^npδ)p(1p)nδ214nδ2α\P(\abs{\widehat p_n - p} \geq \delta) \leq \frac{p(1-p)}{n\delta^2} \leq \frac1{4n\delta^2} \leq \alpha untuk n14αδ2n \geq \frac1{4\alpha\delta^2}, yakni 55565556 pada tiga poin — sekitar 2.72.7 kali kebutuhan Hoeffding. Tanpa pengembalian, ragamnya terkalikan NnN1<1\frac{N - n}{N-1} < 1 (Latihan 22.5), jadi nn yang sama hanya bisa lebih baik: perhitungan dengan pengembaliannyalah yang konservatif.

15. {p^n12}{p^n0.520.02}\{\widehat p_n \leq \frac12\} \subseteq \{\widehat p_n - 0.52 \leq -0.02\}, jadi menurut batas Hoeffding satu sisi P(p^n12)e2n(0.02)20.01\P(\widehat p_n \leq \tfrac12) \leq \eu^{-2n(0.02)^2} \leq 0.01 segera setelah nln10020.00045756.5n \geq \frac{\ln 100}{2\cdot0.0004} \approx 5756.5: 57575757 pemilih. Biayanya berskala seperti kuadrat kebalikan keunggulannya, bukan kebalikan ketelitian yang diinginkan: pemilihan yang ketat itu mahal.

16. Yang dipakai: cuplikannya diambil secara seragam dan saling bebas dari kumpulan pemilihnya; setiap orang yang tercuplik menjawab, dengan jujur, dan pp tak bergerak selama penjajakannya. Jajak pendapat nyata melanggar ketiganya: responden yang terjangkau dan bersedia bukanlah cuplikan seragam (bias pemilihan dan bias tanpa jawaban), dan jawabannya bisa tak jujur atau tak stabil. Semuanya galat bias: ia menggeser E(p^n)\E(\widehat p_n) menjauhi pp sebesar suatu jumlah yang tak bergantung pada nn, jadi tak ada ukuran cuplikan yang menguranginya — matematika pada Bagian ini hanya mengendalikan suku fluktuasinya.

17. Chebyshev untuk satu kelompok berukuran mm: P(p^(i)pδ)14mδ218\P(\abs{ \widehat p^{(i)} - p} \geq \delta) \leq \frac{1}{4m\delta^2} \leq \frac18 untuk m2δ2m \geq \frac2{\delta^2}. Bila kelompok yang keliru kurang dari k/2k/2, maka lebih dari k/2k/2 di antara nilai p^(i)\widehat p^{(i)} terletak pada selang buka (pδ,p+δ)\intoo{p - \delta}{p + \delta}, begitu pula mediannya; jadi {Mpδ}\{\abs{M - p} \geq \delta\} memaksa sedikitnya k/2\lceil k/2\rceil kekeliruan di antara kk kelompok yang saling bebas. Batas gabungan atas (kk/2)\binom k{\lceil k/2\rceil} himpunan kelompok keliru yang mungkin memberi

P(Mpδ)(kk/2)(18)k/22k8k/2=2k/2:\P(\abs{M - p} \geq \delta) \leq \binom{k}{\lceil k/2\rceil} \Bigl(\frac18\Bigr)^{k/2} \leq 2^k\,8^{-k/2} = 2^{-k/2} :

yakni peluruhan eksponensial terhadap banyaknya kelompok, yang ditebus dengan ragam belaka — berguna justru ketika sukunya takterbatas dan Hoeffding tak tersedia.

18. Cauchy–Schwarz (Teorema 22.19):

E(X)=E(X1X>0)E(X2)E(1X>02)=E(X2)P(X>0);\E(X) = \E(X\,\mathbf 1_{X>0}) \leq \sqrt{\E(X^2)}\sqrt{\E(\mathbf 1_{X>0}^2)} = \sqrt{\E(X^2)\,\P(X > 0)} ;

kuadratkan lalu bagi.

19. Ambil h(a)=I(a)2(a12)2h(a) = I(a) - 2(a - \tfrac12)^2. Maka h(12)=0h(\tfrac12) = 0, h(a)=lna1a4(a12)h'(a) = \ln\frac a{1-a} - 4(a - \tfrac12) lenyap di 12\tfrac12, dan

h(a)=1a+11a4=1a(1a)40h''(a) = \frac1a + \frac1{1-a} - 4 = \frac{1}{a(1-a)} - 4 \geq 0

sebab a(1a)14a(1-a) \leq \frac14. Jadi hh' naik dari 00 pada [12,1)\intco{\frac12}1, sehingga h0h' \geq 0 dan h0h \geq 0: I(a)2(a12)2I(a) \geq 2(a - \tfrac12)^2.

20. I(12)=I(12)=0I(\tfrac12) = I'(\tfrac12) = 0, I(a)=1a(1a)I''(a) = \frac1{a(1-a)} memberi I(12)=4I''(\tfrac12) = 4, dan I(12)=0I'''(\tfrac12) = 0 (fungsinya setangkup terhadap 12\tfrac12), jadi I(12+δ)=2δ2+O(δ4)I(\tfrac12 + \delta) = 2\delta^2 + O(\delta^4). Pertanyaan 4 lalu membatasi ekor sejatinya dari bawah oleh en(2δ2+O(δ4))/(n+1)\eu^{-n(2\delta^2 + O(\delta^4))}/(n+1): untuk δ\delta yang kecil pangkat Hoeffding 2nδ22n\delta^2 persis secara asimtotik — hanya perbaikan suku-banyak-dalam-nn yang mungkin.

21. Markov: satu momen, peluruhan 1/a1/a, berguna hanya sebagai mesin di balik yang lain (pertanyaan 1 menunjukkannya datar). Chebyshev: dua momen, peluruhan Vnδ2\frac{V}{n\delta^2}, tajam tanpa hipotesis tambahan (Contoh 22.18), sekaligus perkakas terbaik pada pertanyaan 23. Momen keempat (Latihan 22.9): peluruhan C/n2C/n^2, keterjumlahannya pas cukup untuk sebuah hukum kuat. Hoeffding: peubah terbatas, peluruhan 2e2nδ22\eu^{-2n\delta^2}, kuda beban Bagian III. Chernoff dengan laju persisnya I(a)I(a): momen eksponensial penuh, pangkat yang tak terkalahkan (pertanyaan 4, 20), titik acuan bagi segala yang lain.

22. Bila koinnya setimbang: P(p^n>0.525)P(p^n120.025)e2n(0.025)2\P(\widehat p_n > 0.525) \leq \P(\widehat p_n - \tfrac12 \geq 0.025) \leq \eu^{-2n(0.025)^2}. Bila p=0.55p = 0.55: P(p^n0.525)P(p^n0.550.025)e2n(0.025)2\P(\widehat p_n \leq 0.525) \leq \P(\widehat p_n - 0.55 \leq -0.025) \leq \eu^{-2n(0.025)^2}. Kedua galatnya di bawah 0.010.01 ketika 2n(0.025)2ln1002n(0.025)^2 \geq \ln 100, yakni n3684.2n \geq 3684.2: 36853685 lemparan. (Membedakan hipotesis yang terpaut 2.52.5 poin berharga sama dengan menaksir sampai ketelitian ±2.5\pm2.5 poin.)

23. Hoeffding: nln402(0.005)273778n \geq \frac{\ln 40}{2(0.005)^2} \approx 73\,778. Chebyshev dengan ragam sejatinya p(1p)=0.0099p(1-p) = 0.0099: n0.00990.05(0.005)2=7920n \geq \frac{0.0099}{0.05\cdot(0.005)^2} = 7920 — sembilan kali lebih murah. Pangkat Hoeffding 2nδ22n\delta^2 menghargai ragamnya pada kasus terburuknya 14\frac14, yang pesimistis secara ganjil ketika p=0.01p = 0.01; adapun momen kedua yang sederhana itu lebih tahu. Perkakas yang hilang adalah batas eksponensial yang sadar-ragam (ketaksamaan Bernstein, Tahun ke-3) — atau, untuk kejadian langka, hampiran Poisson yang dibuktikan pada Bab 23, yang bekerja pada skala nisbi yang alami.

24. Tetapkan δ>0\delta > 0: n2e2nδ2<\sum_n 2\eu^{-2n\delta^2} < \infty (deret bertipe geometrik), jadi Borel–Cantelli 1 (Teorema 21.25) memberi P(p^npδ tak hingga kali)=0\P(\abs{\widehat p_n - p} \geq \delta \text{ tak hingga kali}) = 0, yakni kejadian Ej=NnN{p^np<1j}E_j = \bigcup_N\bigcap_{n\geq N}\{\abs{\widehat p_n - p} < \tfrac1j\} berpeluang 11 untuk setiap jj. Irisan terbilang jEj\bigcap_jE_j tetap berpeluang 11 (kesubaditifan pada komplemennya), dan padanya p^np\widehat p_n \to p: itulah hukum kuat bilangan besar untuk lemparan koin, dengan keterbatasan memainkan peran yang dimainkan momen keempat pada Latihan 22.9.

25. Satu momen membeli batas yang datar; dua membeli 1/(nδ2)1/(n\delta^2), tak lebih (contoh ketajamannya); empat membeli 1/n21/n^2, cukup untuk meneleskopkannya menjadi hukum hampir pasti; keterbatasan membeli e2nδ2\eu^{-2n\delta^2}; dan momen eksponensial penuh membeli laju persisnya II, yang tak terkalahkan oleh metode mana pun. Menjajaki 20502050 orang memadai bagi negeri mana pun sebab fluktuasi cuplikannya dikendalikan oleh ragam koinnya, bukan oleh besarnya populasi — banderol 1/δ21/\delta^2 dan ln(1/α)\ln(1/\alpha) berlaku semesta. Teorema limit pusat pada jilid Tahun ke-3 menggantikan ketaksamaan ini, pada skala n\sqrt n, dengan hukum limit persis berkonstanta eksplisit — yang mengubah setiap batas pada masalah ini menjadi kesamaan asimtotik.

Istilah yang didefinisikan dalam bab ini

Lihat semua 395 istilah di glosarium