6장. 큰 수의 법칙과 중심극한정리

확률과 통계 · 2부. 극한정리

심화 · 약 10분 · 연습문제 5 · 그림 1

요약

마르코프·체비쇼프 부등식, 표본평균의 분포, 큰 수의 법칙, 중심극한정리, 정규근사와 오차 막대

동전을 많이 던지면 앞면의 비율은 12\frac12에 가까워지고(큰 수의 법칙), 그 비율의 흔들림은 어떤 분포에서 시작하든 종 모양(정규분포)이 됩니다(중심극한정리). 이 두 정리가 “데이터가 많으면 평균을 믿을 수 있다”, “오차 막대는 σ/n\sigma/\sqrt n”이라는 실험 과학과 머신러닝의 상식을 정당화합니다. 정규분포가 어디에나 나오는 이유도 여기에 있습니다.

6.1 확률 부등식

분포를 몰라도 평균과 분산만으로 꼬리 확률을 막을 수 있습니다.

정리 6.1 (마르코프 부등식). X≥0X \ge 0이고 a>0a > 0이면 P(X≥a)≤E[X]aP(X \ge a) \le \dfrac{E[X]}{a}.

증명. X≥a⋅[X≥a]X \ge a\cdot[X \ge a](지시변수)이므로 양변의 기댓값을 취하면 E[X]≥aP(X≥a)E[X] \ge aP(X \ge a). □

정리 6.2 (체비쇼프 부등식). P(∣X−μ∣≥kσ)≤1k2P(|X - \mu| \ge k\sigma) \le \dfrac{1}{k^2}.

증명. (X−μ)2(X - \mu)^2에 마르코프 부등식을 씁니다. □

어떤 분포든 평균에서 2σ2\sigma 이상 벗어날 확률은 25% 이하입니다(정규분포라면 약 5%). 분포를 모를 때 쓰는 보수적인 상한입니다.

6.2 표본평균

독립이고 같은 분포(independent and identically distributed, i.i.d.)를 따르는 X1,…,XnX_1, \dots, X_n(평균 μ\mu, 분산 σ2\sigma^2)의 표본평균 Xˉn=1n∑Xi\bar X_n = \frac1n\sum X_i에 대해

E[Xˉn]=μ,Var(Xˉn)=σ2n.E[\bar X_n] = \mu, \qquad \mathrm{Var}(\bar X_n) = \frac{\sigma^2}{n}.

표준편차는 σ/n\sigma/\sqrt n입니다. 데이터를 4배 모으면 오차가 절반이 됩니다. 이 σ/n\sigma/\sqrt n을 표준오차(standard error)라 합니다.

6.3 큰 수의 법칙

정리 6.3 (큰 수의 약법칙). 모든 ε>0\varepsilon > 0에 대해 lim⁡n→∞P(∣Xˉn−μ∣≥ε)=0\displaystyle\lim_{n\to\infty}P(|\bar X_n - \mu| \ge \varepsilon) = 0.

증명. 체비쇼프 부등식으로 P(∣Xˉn−μ∣≥ε)≤σ2nε2→0P(|\bar X_n - \mu| \ge \varepsilon) \le \frac{\sigma^2}{n\varepsilon^2} \to 0. □

몬테카를로 방법(12장)이 정당한 이유입니다: E[g(X)]E[g(X)]를 표본 평균 1n∑g(Xi)\frac1n\sum g(X_i)로 어림하면 nn이 클수록 맞습니다.

6.4 중심극한정리

정리 6.4 (중심극한정리). i.i.d. XiX_i의 평균이 μ\mu, 분산이 σ2<∞\sigma^2 < \infty이면

Xˉn−μσ/n→ d N(0,1),\frac{\bar X_n - \mu}{\sigma/\sqrt n} \xrightarrow{\ d\ } \mathcal{N}(0, 1),

곧 모든 zz에 대해 P ⁣(Xˉn−μσ/n≤z)→Φ(z)P\!\left(\frac{\bar X_n - \mu}{\sigma/\sqrt n} \le z\right) \to \Phi(z)입니다.

증명 개요. 표준화한 합의 적률생성함수(3장)는 (M(t/n))n\left(M(t/\sqrt n)\right)^n 꼴이고, M(s)=1+s22+O(s3)M(s) = 1 + \frac{s^2}2 + O(s^3)(평균 0, 분산 1)를 넣으면 (1+t22n+⋯ )n→et2/2\left(1 + \frac{t^2}{2n} + \cdots\right)^n \to e^{t^2/2}, 곧 표준정규분포의 적률생성함수로 갑니다. □

원래 분포가 균등이든, 치우친 지수분포든, 0과 1만 나오는 동전이든 상관없습니다. 많은 작은 독립 효과의 합은 정규분포가 됩니다. 측정 오차, 키, 시험 점수가 정규분포에 가까운 이유입니다.

분포와 표본 크기 n을 바꾸면, 표본평균의 분포(막대)가 n이 커질수록 표준정규분포(곡선)에 가까워집니다. 치우친 지수분포도 n이 30쯤이면 거의 종 모양입니다

이 위젯은 JavaScript가 필요합니다

6.5 정규근사와 오차 막대

예제 6.5. 공정한 동전을 100번 던져 앞면이 60번 이상 나올 확률: S≈N(50,25)S \approx \mathcal{N}(50, 25)이므로 P(S≥60)≈1−Φ ⁣(59.5−505)=1−Φ(1.9)≈0.0287P(S \ge 60) \approx 1 - \Phi\!\left(\frac{59.5 - 50}{5}\right) = 1 - \Phi(1.9) \approx 0.0287. (59.5는 이산을 연속으로 바꿀 때의 연속성 보정입니다. 정확한 이항 값은 0.02840.0284.)

예제 6.6 (오차 막대). 실험을 nn번 반복해 평균 xˉ\bar x와 표본 표준편차 ss를 얻으면, 결과를 xˉ±s/n\bar x \pm s/\sqrt n으로 적습니다. 중심극한정리로 참값이 xˉ±1.96 s/n\bar x \pm 1.96\,s/\sqrt n 안에 있을 확률이 약 95%입니다(8장의 신뢰구간).

6.6 이 장의 요약

  • 마르코프 P(X≥a)≤E[X]/aP(X \ge a) \le E[X]/a, 체비쇼프 P(∣X−μ∣≥kσ)≤1/k2P(|X - \mu| \ge k\sigma) \le 1/k^2.
  • 표본평균의 표준오차 σ/n\sigma/\sqrt n. 큰 수의 법칙: Xˉn→μ\bar X_n \to \mu.
  • 중심극한정리: 표준화한 평균은 분포와 상관없이 표준정규분포로.

쓰이는 곳. 미니배치 SGD가 작동하는 이유(미니배치 그래디언트는 평균이라 분산이 1/∣B∣1/|B|로 줄어듭니다, 다변수 미적분 5장), 실험 결과의 오차 막대와 A/B 테스트, 몬테카를로 적분의 오차 O(1/n)O(1/\sqrt n), 학습 이론의 일반화 상한(호에프딩 부등식은 체비쇼프의 강화판)이 이 장의 정리를 씁니다.

문제

6.1 [계산] · 확인 문제

평균 50, 표준편차 10인 분포에서 ∣X−50∣≥30|X - 50| \ge 30일 확률의 체비쇼프 상한은? (분수로 답할 수 있습니다)

정답
1/9
풀이

30=3σ30 = 3\sigma이므로 1/32=1/91/3^2 = 1/9입니다.

6.2 [계산] · 확인 문제

표준편차가 12인 분포에서 표본 36개의 평균의 표준오차는?

정답
2
풀이

12/36=212/\sqrt{36} = 2입니다.

6.3 [계산] · 확인 문제

주사위를 100번 던진 눈의 합 SS가 380 이상일 확률을 정규근사(연속성 보정 없이)로 구합니다. (소수 넷째 자리까지)

힌트 1

주사위 하나의 평균 3.5, 분산 35/12입니다.

정답
0.0395
풀이

E[S]=350E[S] = 350, Var(S)=100⋅3512≈291.67\mathrm{Var}(S) = 100\cdot\frac{35}{12} \approx 291.67이므로 σS≈17.08\sigma_S \approx 17.08입니다. P(S≥380)≈1−Φ ⁣(3017.08)=1−Φ(1.757)≈0.0395P(S \ge 380) \approx 1 - \Phi\!\left(\frac{30}{17.08}\right) = 1 - \Phi(1.757) \approx 0.0395입니다.

6.4 [계산]

오차를 지금의 1/5로 줄이려면 데이터를 몇 배로 늘려야 할까요?

정답
25
풀이

표준오차가 1/n1/\sqrt n에 비례하므로 25=5\sqrt{25} = 5, 25배입니다.

6.5 [증명]

i.i.d. XiX_i에 대해 Var(Xˉn)=σ2/n\mathrm{Var}(\bar X_n) = \sigma^2/n임을 보입니다.

풀이

독립이면 분산이 더해지므로 Var(∑Xi)=nσ2\mathrm{Var}(\sum X_i) = n\sigma^2이고, Var(aY)=a2Var(Y)\mathrm{Var}(aY) = a^2\mathrm{Var}(Y)이므로 Var(1n∑Xi)=nσ2n2=σ2n\mathrm{Var}(\frac1n\sum X_i) = \frac{n\sigma^2}{n^2} = \frac{\sigma^2}n입니다.