Skip to content

4회 · 셈하는 분포와 섞인 분포 · 워크북

A. 계산 드릴 ​

pmf, 밀도, 적률, 사후분포를 손으로 계산한다. 답이 수나 식으로 떨어지므로 다 풀면 뒤의 최종답이나 노트북의 check()로 맞춰 본다. 계산기는 마지막 수치에만 쓴다.

A1 세 번째 6까지 ★ ​

정리 1.3 · 8분

공정한 주사위를 6이 세 번 나올 때까지 던진다. (a) 6이 아닌 눈이 정확히 5번 나올 확률을 구하라. (b) 던지는 횟수의 평균을 구하라.

힌트

6이 아닌 눈의 수는 NB(3,16)\mathrm{NB}(3,\tfrac16)이다. 던지는 횟수는 그것에 3을 더한 것이다.

풀이

(a) 6이 아닌 눈의 수 Y∼NB(3,16)Y\sim\mathrm{NB}(3,\tfrac16)이므로 정리 1.3으로

P(Y=5)=(5+3−15)(16)3(56)5=21⋅5568=21875559872≈0.0391.P(Y=5)=\binom{5+3-1}5\Bigl(\frac16\Bigr)^3\Bigl(\frac56\Bigr)^5=21\cdot\frac{5^5}{6^8}=\frac{21875}{559872}\approx0.0391 .

(b) 던지는 횟수는 Y+3Y+3이고 EY=rqp=3⋅5/61/6=15EY=\tfrac{rq}p=\tfrac{3\cdot5/6}{1/6}=15이므로 평균은 18이다. 6 하나를 보는 데 평균 6번이 걸리고, 그것이 세 번이다.

최종답 (a) (75)5568≈0.0391\tbinom75\tfrac{5^5}{6^8}\approx0.0391 (b) 1818

A2 브리지 패의 하트 ★ ​

정의 1.4, 정리 1.5 · 10분

잘 섞은 52장에서 13장을 받는다. 하트의 수를 XX라 하자. (a) EXEX, (b) Var⁡X\Var X, (c) P(X=0)P(X=0)을 구하라.

힌트

N=52N=52, K=13K=13, n=13n=13인 초기하분포다. 유한 모집단 수정 N−nN−1\tfrac{N-n}{N-1}을 잊지 않는다.

풀이

p=1352=14p=\tfrac{13}{52}=\tfrac14이다. (a) EX=13⋅14=134EX=13\cdot\tfrac14=\tfrac{13}4. (b) 정리 1.5로

Var⁡X=13⋅14⋅34⋅52−1352−1=3916⋅3951=507272≈1.864.\Var X=13\cdot\frac14\cdot\frac34\cdot\frac{52-13}{52-1}=\frac{39}{16}\cdot\frac{39}{51}=\frac{507}{272}\approx1.864 .

되돌려 넣는 뽑기라면 3916≈2.44\tfrac{39}{16}\approx2.44이다. (c) 하트가 아닌 39장에서 13장을 모두 받아야 하므로 P(X=0)=(3913)/(5213)≈0.0128P(X=0)=\tbinom{39}{13}/\tbinom{52}{13}\approx0.0128이다.

최종답 (a) 134\tfrac{13}4 (b) 507272≈1.864\tfrac{507}{272}\approx1.864 (c) (3913)/(5213)≈0.0128\tbinom{39}{13}\big/\tbinom{52}{13}\approx0.0128

A3 열 명의 혈액형 ★★ ​

정의 1.7, 정리 1.8 · 10분

혈액형 A, B, O, AB의 비율이 0.340.34, 0.270.27, 0.280.28, 0.110.11인 집단에서 열 명을 독립적으로 뽑는다. (a) A, B, O가 3명씩이고 AB가 1명일 확률을 구하라. (b) A형의 수와 B형의 수의 공분산을 구하라. (c) 두 수의 상관계수를 구하라.

힌트

(a) 다항계수 10!3! 3! 3! 1!\tfrac{10!}{3!\,3!\,3!\,1!}. (c) 상관계수는 공분산을 두 표준편차로 나눈 것이고, 각 수의 주변분포는 이항이다.

풀이

(a) 10!3! 3! 3! 1!=16800\tfrac{10!}{3!\,3!\,3!\,1!}=16800이므로 확률은 16800⋅0.343⋅0.273⋅0.283⋅0.11≈0.031416800\cdot0.34^3\cdot0.27^3\cdot0.28^3\cdot0.11\approx0.0314이다. (b) 정리 1.8로 Cov⁡(XA,XB)=−10⋅0.34⋅0.27=−0.918\Cov(X_A,\allowbreak X_B)=-10\cdot0.34\cdot0.27=-0.918이다. (c) Var⁡XA=10⋅0.34⋅0.66\Var X_A=10\cdot0.34\cdot0.66, Var⁡XB=10⋅0.27⋅0.73\Var X_B=10\cdot0.27\cdot0.73이므로

ρ=−0.9182.244⋅1.971=−0.34⋅0.270.66⋅0.73≈−0.437.\rho=\frac{-0.918}{\sqrt{2.244\cdot1.971}}=-\sqrt{\frac{0.34\cdot0.27}{0.66\cdot0.73}}\approx-0.437 .

nn이 상관계수에서 사라진다. 상관은 몇 명을 뽑든 범주의 비율로만 정해진다.

최종답 (a) ≈0.0314\approx0.0314 (b) −0.918-0.918 (c) ≈−0.437\approx-0.437

A4 불량품 셋 이상 ★ ​

정리 2.3 · 10분

불량률이 0.0020.002인 공정에서 1000개를 만든다(각각 독립). 불량품이 3개 이상일 확률을 (a) Poisson 근사로, (b) 이항분포로 정확히 구하라. (c) 2.1절의 Le Cam 한계 np2np^2은 얼마인가? (a)와 (b)의 차이와 비교하라.

힌트

λ=np=2\lambda=np=2. (b)는 1−∑k=02(1000k)(0.002)k(0.998)1000−k1-\sum_{k=0}^2\tbinom{1000}k(0.002)^k(0.998)^{1000-k}이다.

풀이

(a) λ=2\lambda=2이므로 1−e−2(1+2+2)=1−5e−2≈0.3233241-e^{-2}(1+2+2)=1-5e^{-2}\approx0.323324이다. (b) 0.9981000≈0.1350650.998^{1000}\approx0.135065, 1000⋅0.002⋅0.998999≈0.2706711000\cdot0.002\cdot0.998^{999}\approx0.270671, (10002)(0.002)2(0.998)998≈0.270942\tbinom{1000}2(0.002)^2(0.998)^{998}\approx0.270942이므로 1−0.676677=0.3233231-0.676677=0.323323이다. (c) np2=0.004np^2=0.004이다. 실제 차이는 10−610^{-6} 정도로 한계보다 훨씬 작다. Le Cam의 한계는 모든 사건에 대한 최대 차이의 상한이다.

최종답 (a) 1−5e−2≈0.32331-5e^{-2}\approx0.3233 (b) ≈0.323323\approx0.323323 (c) 0.0040.004

A5 콜센터 ★ ​

정의 3.1, 정리 3.2 · 10분

전화가 시간당 평균 6통씩 Poisson 과정으로 걸려 온다. 시간은 분 단위로 쓴다. (a) 10분 동안 전화가 한 통도 없을 확률을 구하라. (b) 세 번째 전화가 올 때까지의 시간의 평균과 (c) 분산을 구하라. (d) 첫 전화가 10분과 20분 사이에 올 확률을 구하라.

힌트

비율은 분당 λ=0.1\lambda=0.1이다. 세 번째 전화의 시각은 Γ(3,1/λ)\Gamma(3,1/\lambda)이고, 첫 전화의 시각은 Exp(λ)\mathrm{Exp}(\lambda)이다.

풀이

분당 비율 λ=0.1\lambda=0.1이다. (a) N(10)∼Poi(1)N(10)\sim\mathrm{Poi}(1)이므로 e−1≈0.368e^{-1}\approx0.368이다. (b)(c) 정리 3.2로 T3∼Γ(3,10)T_3\sim\Gamma(3,10)이므로 평균 3⋅10=303\cdot10=30분, 분산 3⋅102=3003\cdot10^2=300분2^2이다. (d) T1∼Exp(0.1)T_1\sim\mathrm{Exp}(0.1)이므로 P(10<T1≤20)=e−1−e−2≈0.233P(10<T_1\le20)=e^{-1}-e^{-2}\approx0.233이다.

최종답 (a) e−1≈0.368e^{-1}\approx0.368 (b) 30분30\text{분} (c) 300분2300\text{분}^2 (d) e−1−e−2≈0.233e^{-1}-e^{-2}\approx0.233

A6 방문자와 구매자 ★★ ​

정리 2.5 · 10분

하루 방문자 수가 Poi(10)\mathrm{Poi}(10)이고, 방문자는 각자 독립적으로 확률 0.30.3으로 물건을 산다. (a) 구매자가 한 명도 없을 확률을 구하라. (b) 사지 않은 사람이 7명이었다는 것을 알 때 구매자 수의 조건부분산을 구하라. (c) 방문자가 10명이었다는 것을 알 때 구매자 수의 조건부분산을 구하라.

힌트

정리 2.5로 구매자 수와 비구매자 수는 독립이다. (c)에서는 전체가 고정되어 있으므로 이항이다.

풀이

(a) 정리 2.5로 구매자 수 N1∼Poi(3)N_1\sim\mathrm{Poi}(3)이므로 e−3≈0.0498e^{-3}\approx0.0498이다. (b) N1N_1과 비구매자 수 N2N_2는 독립이므로 N2=7N_2=7을 알아도 N1N_1의 분포는 그대로 Poi(3)\mathrm{Poi}(3)이고 분산은 3이다. (c) N=10N=10이 주어지면 N1∼b(10,0.3)N_1\sim\mathrm b(10,0.3)이므로 분산은 10⋅0.3⋅0.7=2.110\cdot0.3\cdot0.7=2.1이다. 전체를 알면 구매자 수에 대한 정보가 생기고, 비구매자 수만 알면 정보가 없다.

최종답 (a) e−3≈0.0498e^{-3}\approx0.0498 (b) 33 (c) 2.12.1

A7 감마함수로 적분하기 ★ ​

정의 3.4, 정리 3.5 · 8분

다음을 감마함수로 계산하라. (a) Γ(72)\Gamma(\tfrac72)  (b) ∫0∞x4e−2x dx\int_0^\infty x^4e^{-2x}\,dx  (c) ∫0∞x e−x dx\int_0^\infty\sqrt x\,e^{-x}\,dx

힌트

(a) Γ(α+1)=αΓ(α)\Gamma(\alpha+1)=\alpha\Gamma(\alpha)를 세 번. (b) y=2xy=2x로 치환한다.

풀이

(a) Γ(72)=52⋅32⋅12Γ(12)=15π8\Gamma(\tfrac72)=\tfrac52\cdot\tfrac32\cdot\tfrac12\Gamma(\tfrac12)=\tfrac{15\sqrt\pi}8. (b) y=2xy=2x이면 ∫0∞y416e−ydy2=Γ(5)32=2432=34\int_0^\infty\tfrac{y^4}{16}e^{-y}\tfrac{dy}2=\tfrac{\Gamma(5)}{32}=\tfrac{24}{32}=\tfrac34. (c) Γ(32)=12Γ(12)=π2\Gamma(\tfrac32)=\tfrac12\Gamma(\tfrac12)=\tfrac{\sqrt\pi}2.

최종답 (a) 15π8\tfrac{15\sqrt\pi}8 (b) 34\tfrac34 (c) π2\tfrac{\sqrt\pi}2

A8 감마분포의 적률과 확률 ★★ ​

정리 3.2, 3.7 · 10분

X∼Γ(3,2)X\sim\Gamma(3,2)이다. (a) EX2EX^2, (b) E(1/X)E(1/X), (c) P(X≤2)P(X\le2)를 구하라.

힌트

(a)(b)는 EXk=βkΓ(α+k)/Γ(α)EX^k=\beta^k\Gamma(\alpha+k)/\Gamma(\alpha) (노트 연습 3.4). (c) XX를 비율 12\tfrac12인 Poisson 과정의 세 번째 사건의 시각으로 읽는다.

풀이

(a) EX2=22Γ(5)Γ(3)=4⋅12=48EX^2=2^2\tfrac{\Gamma(5)}{\Gamma(3)}=4\cdot12=48이다. 분산 αβ2=12\alpha\beta^2=12와 평균의 제곱 36을 더한 값과 같다. (b) E(1/X)=2−1Γ(2)Γ(3)=12⋅12=14E(1/X)=2^{-1}\tfrac{\Gamma(2)}{\Gamma(3)}=\tfrac12\cdot\tfrac12=\tfrac14이다. 1/EX=161/EX=\tfrac16보다 크다. (c) Γ(3,2)\Gamma(3,2)는 비율 λ=12\lambda=\tfrac12인 Poisson 과정의 T3T_3이므로 정리 3.2로 P(X≤2)=P(N(2)≥3)P(X\le2)=P(N(2)\ge3), N(2)∼Poi(1)N(2)\sim\mathrm{Poi}(1)이다. 1−e−1(1+1+12)=1−52e−1≈0.08031-e^{-1}(1+1+\tfrac12)=1-\tfrac52e^{-1}\approx0.0803이다.

최종답 (a) 4848 (b) 14\tfrac14 (c) 1−52e−1≈0.08031-\tfrac52e^{-1}\approx0.0803

A9 표준정규의 제곱합 ★ ​

정의 3.8, 정리 3.9 · 8분

Z1,Z2,Z3∼iidN(0,1)Z_1,Z_2,Z_3\iid\N(0,1)이다. (a) P(Z12+Z22>4)P(Z_1^2+Z_2^2>4), (b) Var⁡(Z12+Z22+Z32)\Var(Z_1^2+Z_2^2+Z_3^2), (c) Ee(Z12+Z22+Z32)/4Ee^{(Z_1^2+Z_2^2+Z_3^2)/4}를 구하라.

힌트

χ2(2)=Γ(1,2)\chi^2(2)=\Gamma(1,2)는 평균 2인 지수분포다. (c)는 χ2(3)\chi^2(3)의 mgf를 t=14t=\tfrac14에서.

풀이

(a) 정리 3.9로 Z12+Z22∼χ2(2)=Γ(1,2)Z_1^2+Z_2^2\sim\chi^2(2)=\Gamma(1,\allowbreak 2), 평균 2인 지수분포이므로 P(>4)=e−2≈0.135P(>4)=e^{-2}\approx0.135. (b) χ2(3)\chi^2(3)의 분산 2⋅3=62\cdot3=6. (c) mgf (1−2t)−3/2(1-2t)^{-3/2}에 t=14t=\tfrac14를 넣어 (12)−3/2=22(\tfrac12)^{-3/2}=2\sqrt2.

최종답 (a) e−2≈0.135e^{-2}\approx0.135 (b) 66 (c) 222\sqrt2

A10 베타분포 ★ ​

정의 4.1, 따름정리 4.3 · 8분

U∼Beta(3,2)U\sim\mathrm{Beta}(3,2)이다. (a) EUEU, (b) Var⁡U\Var U, (c) P(U>12)P(U>\tfrac12)를 구하라.

힌트

밀도는 12u2(1−u)12u^2(1-u)이다. B(3,2)=Γ(3)Γ(2)Γ(5)=112B(3,\allowbreak 2)=\tfrac{\Gamma(3)\Gamma(2)}{\Gamma(5)}=\tfrac1{12}.

풀이

(a) 33+2=35\tfrac3{3+2}=\tfrac35. (b) 3⋅252⋅6=125\tfrac{3\cdot2}{5^2\cdot6}=\tfrac1{25}. (c) ∫1/2112u2(1−u) du=12[u33−u44]1/21=12(112−5192)=1116\int_{1/2}^112u^2(1-u)\,du=12\bigl[\tfrac{u^3}3-\tfrac{u^4}4\bigr]_{1/2}^1=12\bigl(\tfrac1{12}-\tfrac5{192}\bigr)=\tfrac{11}{16}.

최종답 (a) 35\tfrac35 (b) 125\tfrac1{25} (c) 1116\tfrac{11}{16}

A11 사후분포와 다음 시행 ★★ ​

정리 4.5, 따름정리 4.3 · 10분

성공 확률 pp의 사전분포가 Beta(2,2)\mathrm{Beta}(2,2)이고 10번 가운데 7번 성공했다. (a) 다음 시행이 성공할 확률, (b) 사후분산, (c) 다음 두 시행이 모두 성공할 확률을 구하라.

힌트

사후분포는 Beta(9,5)\mathrm{Beta}(9,5)다. (c)는 E(p2∣자료)E(p^2\mid\text{자료})이다. 두 시행은 pp가 주어지면 독립이지만, pp를 모르면 독립이 아니다.

풀이

정리 4.5로 사후분포는 Beta(2+7,2+3)=Beta(9,5)\mathrm{Beta}(2+7,\allowbreak 2+3)=\mathrm{Beta}(9,\allowbreak 5)이다. (a) 반복기댓값으로 E(p∣자료)=914E(p\mid\text{자료})=\tfrac9{14}. (b) 9⋅5142⋅15=3196\tfrac{9\cdot5}{14^2\cdot15}=\tfrac3{196}. (c) pp가 주어지면 두 시행이 모두 성공할 확률은 p2p^2이므로 E(p2∣자료)=9⋅1014⋅15=37≈0.429E(p^2\mid\text{자료})=\tfrac{9\cdot10}{14\cdot15}=\tfrac37\approx0.429이다. (914)2≈0.413(\tfrac9{14})^2\approx0.413보다 크다. 차이가 사후분산 3196\tfrac3{196}이다.

최종답 (a) 914\tfrac9{14} (b) 3196\tfrac3{196} (c) 37\tfrac37

A12 두 Poisson의 혼합 ★★ ​

정의 5.1, 정리 5.2 · 10분

확률 0.30.3으로 Poi(1)\mathrm{Poi}(1), 확률 0.70.7로 Poi(4)\mathrm{Poi}(4)에서 뽑은 XX가 있다. (a) EXEX, (b) Var⁡X\Var X를 구하라. (c) X=0X=0이었을 때 Poi(1)\mathrm{Poi}(1)에서 뽑았을 확률을 구하라.

힌트

(b) 집단 안의 분산의 평균과 집단 평균의 분산을 더한다. (c) 베이즈 정리.

풀이

(a) EX=0.3⋅1+0.7⋅4=3.1EX=0.3\cdot1+0.7\cdot4=3.1. (b) 집단 안: 0.3⋅1+0.7⋅4=3.10.3\cdot1+0.7\cdot4=3.1. 집단 사이: 0.3⋅12+0.7⋅42−3.12=11.5−9.61=1.890.3\cdot1^2+0.7\cdot4^2-3.1^2=11.5-9.61=1.89. 합은 4.994.99로 평균보다 크다. (c) 베이즈 정리로

0.3e−10.3e−1+0.7e−4≈0.11040.1232≈0.896.\frac{0.3e^{-1}}{0.3e^{-1}+0.7e^{-4}}\approx\frac{0.1104}{0.1232}\approx0.896 .

0이 나온 것은 위험이 낮은 집단이었을 증거가 된다. 사전확률 0.30.3이 사후확률 0.8960.896이 된다.

최종답 (a) 3.13.1 (b) 4.994.99 (c) ≈0.896\approx0.896

A13 다섯 점의 중앙값 ★★ ​

예제 4.4, 따름정리 4.3 · 10분

U1,…,U5∼iidU(0,1)U_1,\allowbreak \dots,\allowbreak U_5\iid\mathrm U(0,\allowbreak 1)의 중앙값을 MM이라 하자. (a) P(M≤14)P(M\le\tfrac14)를 구하라. (b) Var⁡M\Var M을 구하라.

힌트

(a) M≤14M\le\tfrac14은 다섯 점 가운데 3개 이상이 14\tfrac14 이하라는 것과 같다. (b) M∼Beta(3,3)M\sim\mathrm{Beta}(3,3).

풀이

(a) 14\tfrac14 이하인 점의 수는 b(5,14)\mathrm b(5,\tfrac14)이므로

P(M≤14)=∑j=35(5j)(14)j(34)5−j=10⋅9+5⋅3+11024=1061024=53512≈0.104.P(M\le\tfrac14)=\sum_{j=3}^5\binom5j\Bigl(\frac14\Bigr)^j\Bigl(\frac34\Bigr)^{5-j}=\frac{10\cdot9+5\cdot3+1}{1024}=\frac{106}{1024}=\frac{53}{512}\approx0.104 .

(b) 예제 4.4로 M=U(3)∼Beta(3,3)M=U_{(3)}\sim\mathrm{Beta}(3,\allowbreak 3)이고 Var⁡M=936⋅7=128\Var M=\tfrac{9}{36\cdot7}=\tfrac1{28}이다. 표준편차 0.190.19는 한 점의 표준편차 0.290.29보다 작다.

최종답 (a) 53512≈0.104\tfrac{53}{512}\approx0.104 (b) 128\tfrac1{28}

A14 닳는 전구와 닳지 않는 전구 ★ ​

정리 3.3 · 8분

이미 500시간 켜진 전구가 앞으로 1000시간 이상 버틸 확률을 구하라. 수명이 (a) 평균 1000시간인 지수분포일 때, (b) Γ(2,500)\Gamma(2,500)일 때. 두 분포의 평균은 같다.

힌트

(b)에서 P(X>x)=e−x/500(1+x500)P(X>x)=e^{-x/500}(1+\tfrac x{500})이다. 이것은 비율 1500\tfrac1{500}인 Poisson 과정에서 xx까지 사건이 1개 이하일 확률이다.

풀이

(a) 무기억성(정리 3.3)으로 P(X>1000)=e−1≈0.368P(X>1000)=e^{-1}\approx0.368이다. (b) 정리 3.2로 P(X>x)=P(N(x)≤1)=e−x/500(1+x500)P(X>x)=P(N(x)\le1)=e^{-x/500}(1+\tfrac x{500})이므로

P(X>1500)P(X>500)=4e−32e−1=2e−2≈0.271.\frac{P(X>1500)}{P(X>500)}=\frac{4e^{-3}}{2e^{-1}}=2e^{-2}\approx0.271 .

같은 평균이라도 α=2>1\alpha=2>1인 감마는 이미 쓴 시간만큼 닳아서 남은 수명이 짧아진다.

최종답 (a) e−1≈0.368e^{-1}\approx0.368 (b) 2e−2≈0.2712e^{-2}\approx0.271