Skip to content

6회 · Student의 t · 워크북

A. 계산 드릴 ​

표본분산, χ2\chi^2, tt, FF와 구간을 손으로 계산한다. 분위수가 필요한 문항은 문제에 값을 적어 두었다. 답이 수나 식으로 떨어지므로 다 풀면 뒤의 최종답이나 노트북의 check()로 맞춰 본다.

A1 표본분산을 세 가지로 ★ ​

정리 1.1, 연습 1.1 · 6분

자료 3, 7, 5, 9, 63,\ 7,\ 5,\ 9,\ 6의 표본분산 s2s^2를 (i) 편차제곱합, (ii) ∑xi2−nxˉ2\sum x_i^2-n\bar x^2, (iii) 12n∑i,j(xi−xj)2\tfrac1{2n}\sum_{i,j}(x_i-x_j)^2의 세 가지로 계산해 같은 값이 나오는지 확인하라.

힌트

(iii)에서 i<ji<j인 쌍 10개의 차의 제곱을 더하면 12∑i,j\tfrac12\sum_{i,j}이다.

풀이

xˉ=6\bar x=6이고 편차는 −3,1,−1,3,0-3,1,-1,3,0이므로 (i) 편차제곱합은 9+1+1+9+0=209+1+1+9+0=20이다. (ii) ∑xi2=9+49+25+81+36=200\sum x_i^2=9+49+25+81+36=200이고 200−5⋅36=20200-5\cdot36=20이다. (iii) i<ji<j인 쌍의 차는 4,2,6,3, 2,2,1, 4,1, 34,2,6,3,\ 2,2,1,\ 4,1,\ 3이고 제곱의 합은 16+4+36+9+4+4+1+16+1+9=10016+4+36+9+4+4+1+16+1+9=100이다. ∑i,j=2⋅100=200\sum_{i,j}=2\cdot100=200이므로 110⋅200=20\tfrac1{10}\cdot200=20이다. 세 경우 모두 s2=20/4=5s^2=20/4=5이다.

최종답 s2=5s^2=5

A2 n=4n=4의 Helmert 회전 ★ ​

보조정리 2.1, 정리 2.2 · 8분

n=4n=4의 Helmert 행렬 HH의 네 번째 행을 쓰고, z=(1,2,3,6)′\z=(1,2,3,6)\T에 대해 Y=Hz\Y=H\z의 성분 Y1,…,Y4Y_1,\dots,Y_4를 구하라. (a) Y4Y_4는 얼마인가? (b) Y22+Y32+Y42Y_2^2+Y_3^2+Y_4^2을 구하고 ∑i(zi−zˉ)2\sum_i(z_i-\bar z)^2과 같음을 확인하라.

힌트

kk번째 행은 1k(k−1)(1,…,1,−(k−1),0,…,0)\tfrac1{\sqrt{k(k-1)}}(1,\allowbreak \dots,\allowbreak 1,\allowbreak -(k-1),\allowbreak 0,\allowbreak \dots,\allowbreak 0)이다. Y1=n zˉY_1=\sqrt n\,\bar z.

풀이

네 번째 행은 112(1,1,1,−3)\tfrac1{\sqrt{12}}(1,1,1,-3)이다. Y1=122=6Y_1=\tfrac{12}2=6, Y2=1−22=−12Y_2=\tfrac{1-2}{\sqrt2}=-\tfrac1{\sqrt2}, Y3=1+2−66=−36Y_3=\tfrac{1+2-6}{\sqrt6}=-\tfrac3{\sqrt6}, Y4=1+2+3−1812=−1212=−23≈−3.464Y_4=\tfrac{1+2+3-18}{\sqrt{12}}=-\tfrac{12}{\sqrt{12}}=-2\sqrt3\approx-3.464이다. 제곱의 합은 12+32+12=14\tfrac12+\tfrac32+12=14이고, zˉ=3\bar z=3이므로 ∑(zi−zˉ)2=4+1+0+9=14\sum(z_i-\bar z)^2=4+1+0+9=14로 같다. Y12=36Y_1^2=36을 더하면 ∥z∥2=50\|\z\|^2=50이다.

최종답 (a) −23≈−3.464-2\sqrt3\approx-3.464 (b) 1414

A3 표본분산이 두 배 넘게 나올 확률 ★ ​

정리 2.2(c) · 8분

X1,…,X5∼iidN(μ,4)X_1,\dots,X_5\iid\N(\mu,4)이다. P(S2>8)P(S^2>8)을 정확히 구하라. (힌트: χ2(4)\chi^2(4)는 모양모수 2인 감마분포다.)

힌트

4S24=S2∼χ2(4)\tfrac{4S^2}{4}=S^2\sim\chi^2(4)이다. χ2(4)\chi^2(4)의 밀도 14xe−x/2\tfrac14xe^{-x/2}를 부분적분한다.

풀이

정리 2.2(c)로 (n−1)S2σ2=4S24=S2∼χ2(4)\tfrac{(n-1)S^2}{\sigma^2}=\tfrac{4S^2}4=S^2\sim\chi^2(4)이다. χ2(4)\chi^2(4)의 밀도는 xe−x/2Γ(2)22=14xe−x/2\tfrac{xe^{-x/2}}{\Gamma(2)2^2}=\tfrac14xe^{-x/2}이고, 부분적분으로 P(χ2(4)>c)=e−c/2(1+c2)P(\chi^2(4)>c)=e^{-c/2}\bigl(1+\tfrac c2\bigr)이다. c=8c=8이면 5e−4≈0.09165e^{-4}\approx0.0916이다. 참분산의 두 배가 넘는 표본분산이 열에 한 번쯤 나온다.

최종답 5e−4≈0.09165e^{-4}\approx0.0916

A4 표본분산의 흔들림 ★ ​

따름정리 2.3 · 6분

X1,…,X9∼iidN(μ,2)X_1,\dots,X_9\iid\N(\mu,2)이다. (a) Var⁡S2\Var S^2를 구하라. (b) S2S^2의 상대 표준편차 Var⁡S2/ES2\sqrt{\Var S^2}/ES^2를 구하라.

힌트

따름정리 2.3의 Var⁡S2=2σ4n−1\Var S^2=\tfrac{2\sigma^4}{n-1}.

풀이

(a) 2⋅228=1\tfrac{2\cdot2^2}{8}=1. (b) ES2=2ES^2=2이므로 12\tfrac{1}{2}이다. 일반적으로 2/(n−1)\sqrt{2/(n-1)}이고 n=9n=9이면 12\tfrac12이다. 아홉 관측으로 잰 분산도 참값의 ±50%\pm50\%가 한 표준편차다.

최종답 (a) 11 (b) 12\tfrac12

A5 세 관측의 표본표준편차 ★ ​

연습 2.4 · 8분

n=3n=3인 정규 표본에서 ES=c4(3) σES=c_4(3)\,\sigma의 c4(3)c_4(3)을 감마함수로 구하고 소수로 나타내라.

힌트

c4(n)=2n−1 Γ(n/2)Γ((n−1)/2)c_4(n)=\sqrt{\tfrac2{n-1}}\,\tfrac{\Gamma(n/2)}{\Gamma((n-1)/2)}, Γ(32)=π2\Gamma(\tfrac32)=\tfrac{\sqrt\pi}2.

풀이

c4(3)=22 Γ(3/2)Γ(1)=π2≈0.886c_4(3)=\sqrt{\tfrac22}\,\tfrac{\Gamma(3/2)}{\Gamma(1)}=\tfrac{\sqrt\pi}2\approx0.886이다. 세 관측의 SS는 σ\sigma를 평균 11% 작게 어림한다.

최종답 π2≈0.886\tfrac{\sqrt\pi}2\approx0.886

A6 t(3)t(3)의 봉우리 ★ ​

정리 3.2 · 6분

t(3)t(3)의 밀도의 0에서의 값 f3(0)f_3(0)을 정확히 구하고, 표준정규의 ϕ(0)=12π≈0.399\phi(0)=\tfrac1{\sqrt{2\pi}}\approx0.399와 비교하라.

힌트

fr(0)=Γ((r+1)/2)rπ Γ(r/2)f_r(0)=\tfrac{\Gamma((r+1)/2)}{\sqrt{r\pi}\,\Gamma(r/2)}이고 Γ(2)=1\Gamma(2)=1.

풀이

f3(0)=Γ(2)3π Γ(3/2)=13π⋅π/2=23 π≈0.368f_3(0)=\tfrac{\Gamma(2)}{\sqrt{3\pi}\,\Gamma(3/2)}=\tfrac1{\sqrt{3\pi}\cdot\sqrt\pi/2}=\tfrac2{\sqrt3\,\pi}\approx0.368이다. 정규보다 낮다. 밀도의 넓이가 같으므로 가운데가 낮은 만큼 꼬리가 두껍다.

최종답 23 π≈0.368\tfrac2{\sqrt3\,\pi}\approx0.368

A7 자유도 1의 tt ★ ​

정리 3.3(c) · 8분

T∼t(1)T\sim t(1)이다. (a) P(T>3)P(T>\sqrt3), (b) P(∣T∣<1)P(|T|<1)을 정확히 구하라.

힌트

t(1)t(1)은 Cauchy이고 분포함수는 12+1πarctan⁡t\tfrac12+\tfrac1\pi\arctan t이다.

풀이

정리 3.3(c)로 밀도는 1π(1+t2)\tfrac1{\pi(1+t^2)}이고 분포함수는 F(t)=12+1πarctan⁡tF(t)=\tfrac12+\tfrac1\pi\arctan t이다. (a) 1−F(3)=12−1π⋅π3=161-F(\sqrt3)=\tfrac12-\tfrac1\pi\cdot\tfrac\pi3=\tfrac16. (b) F(1)−F(−1)=2π⋅π4=12F(1)-F(-1)=\tfrac2\pi\cdot\tfrac\pi4=\tfrac12.

최종답 (a) 16\tfrac16 (b) 12\tfrac12

A8 tt와 FF의 적률 ★ ​

정리 3.3(b), 4.2(b)(c) · 6분

(a) T∼t(6)T\sim t(6)의 분산을 구하라. (b) F∼F(3,8)F\sim F(3,8)의 평균을 구하라. (c) T2T^2의 분포를 쓰고 (a)와 ET2ET^2이 맞는지 확인하라.

힌트

Var⁡T=rr−2\Var T=\tfrac r{r-2}, EF=r2r2−2EF=\tfrac{r_2}{r_2-2}, T2∼F(1,r)T^2\sim F(1,r).

풀이

(a) 66−2=32\tfrac6{6-2}=\tfrac32. (b) 88−2=43\tfrac8{8-2}=\tfrac43. (c) T2∼F(1,6)T^2\sim F(1,6)이고 EF(1,6)=64=32=ET2EF(1,\allowbreak 6)=\tfrac64=\tfrac32=ET^2로 (a)와 같다(ET=0ET=0이므로 Var⁡T=ET2\Var T=ET^2).

최종답 (a) 32\tfrac32 (b) 43\tfrac43

A9 다섯 관측의 tt-구간 ★ ​

정리 5.2 · 10분

정규 모집단에서 10.2, 9.6, 10.8, 10.1, 9.810.2,\allowbreak \ 9.6,\allowbreak \ 10.8,\allowbreak \ 10.1,\allowbreak \ 9.8을 얻었다. t4, 0.025=2.776t_{4,\,0.025}=2.776을 써서 μ\mu의 95% 신뢰구간의 (a) 아래 끝과 (b) 위 끝을 구하라.

힌트

xˉ=10.1\bar x=10.1이다. 편차제곱합을 먼저 구한다.

풀이

xˉ=10.1\bar x=10.1, 편차는 0.1,−0.5,0.7,0,−0.30.1,-0.5,0.7,0,-0.3이고 편차제곱합은 0.840.84이므로 s2=0.21s^2=0.21, s=0.4583s=0.4583이다. 반폭은 2.776⋅0.4583/5=0.5692.776\cdot0.4583/\sqrt5=0.569이므로 구간은 [9.531, 10.669][9.531,\ 10.669]이다.

최종답 (a) 9.5319.531 (b) 10.66910.669

A10 대응 tt 통계량 ★ ​

정리 3.4, 예제 5.3 · 10분

환자 6명의 처리 전후 차이가 2, 4, −1, 3, 5, 12,\ 4,\ -1,\ 3,\ 5,\ 1이다. (a) 대응 tt 통계량 t=dˉ/(sd/n)t=\bar d/(s_d/\sqrt n)을 정확한 수로 구하라. (b) 귀무가설 μD=0\mu_D=0 아래의 자유도는?

힌트

dˉ=73\bar d=\tfrac73, ∑di2=56\sum d_i^2=56. 분수로 계산하면 근호 하나가 남는다.

풀이

dˉ=146=73\bar d=\tfrac{14}6=\tfrac73이고 편차제곱합은 ∑di2−ndˉ2=56−6⋅499=703\sum d_i^2-n\bar d^2=56-6\cdot\tfrac{49}9=\tfrac{70}3이므로 sd2=143s_d^2=\tfrac{14}3이다. sd2/n=79s_d^2/n=\tfrac79이므로 t=7/37/3=7≈2.646t=\tfrac{7/3}{\sqrt7/3}=\sqrt7\approx2.646이다. (b) 자유도는 n−1=5n-1=5이다. t5,0.025=2.571t_{5,0.025}=2.571보다 크므로 양쪽 5% 수준에서 차이가 있다.

최종답 (a) 7≈2.646\sqrt7\approx2.646 (b) 55

A11 합동분산의 두 표본 tt ★ ​

정리 5.4 · 10분

XX: 5,7,65,7,6과 YY: 3,4,2,33,4,2,3이 등분산 정규 모집단의 독립 표본이다. (a) μ1=μ2\mu_1=\mu_2 아래의 통계량 TT를 정확한 수로 구하라. (b) 자유도는?

힌트

두 편차제곱합이 모두 2다. Sp2=2+25S_p^2=\tfrac{2+2}{5}.

풀이

xˉ=6\bar x=6, yˉ=3\bar y=3이고 편차제곱합은 각각 22와 22이므로 Sp2=45S_p^2=\tfrac45이다. Sp2(13+14)=45⋅712=715S_p^2(\tfrac13+\tfrac14)=\tfrac45\cdot\tfrac7{12}=\tfrac7{15}이므로 T=37/15=315/7≈4.392T=\tfrac3{\sqrt{7/15}}=3\sqrt{15/7}\approx4.392이다. (b) m+n−2=5m+n-2=5.

최종답 (a) 315/7≈4.3923\sqrt{15/7}\approx4.392 (b) 55

A12 분산비의 구간 ★ ​

정리 4.3 · 10분

m=7m=7, n=9n=9인 독립 정규 표본에서 sX2=6s_X^2=6, sY2=2s_Y^2=2를 얻었다. F6,8, 0.025=4.652F_{6,8,\,0.025}=4.652, F8,6, 0.025=5.600F_{8,6,\,0.025}=5.600을 써서 σ12/σ22\sigma_1^2/\sigma_2^2의 95% 신뢰구간의 (a) 아래 끝과 (b) 위 끝을 구하라.

힌트

F6,8,0.975=1/F8,6,0.025F_{6,8,0.975}=1/F_{8,6,0.025} (노트 연습 4.2).

풀이

노트 연습 4.3의 구간 [sX2/sY2Fm−1,n−1,0.025, sX2/sY2Fm−1,n−1,0.975]\bigl[\tfrac{s_X^2/s_Y^2}{F_{m-1,n-1,0.025}},\allowbreak \ \tfrac{s_X^2/s_Y^2}{F_{m-1,n-1,0.975}}\bigr]에 비 33을 넣는다. F6,8,0.975=1/5.600F_{6,8,0.975}=1/5.600이므로 구간은 [3/4.652, 3⋅5.600]=[0.645, 16.80][3/4.652,\allowbreak \ 3\cdot5.600]=[0.645,\allowbreak \ 16.80]이다. 1을 포함한다.

최종답 (a) 0.6450.645 (b) 16.8016.80

A13 다음 한 관측의 구간 ★ ​

정리 8.1 · 8분

A9의 자료로 다음 관측 하나에 대한 95% 예측구간의 (a) 아래 끝과 (b) 위 끝을 구하라. A9의 신뢰구간보다 몇 배 넓은가?

힌트

반폭이 t4,0.025 s1+1/5t_{4,0.025}\,s\sqrt{1+1/5}로 바뀐다.

풀이

반폭은 2.776⋅0.4583⋅1.2=1.3942.776\cdot0.4583\cdot\sqrt{1.2}=1.394이므로 구간은 [8.706, 11.494][8.706,\ 11.494]이다. 신뢰구간의 반폭과의 비는 1.2⋅5=6≈2.45\sqrt{1.2}\cdot\sqrt5=\sqrt6\approx2.45배다.

최종답 (a) 8.7068.706 (b) 11.49411.494

A14 세 관측에 1.96을 쓰면 ★★ ​

정리 3.4, 노트 연습 3.3 · 8분

n=3n=3인 정규 표본에서 Xˉ±1.96 S/3\bar X\pm1.96\,S/\sqrt3의 실제 포함확률을 정확한 식으로 구하라. t(2)t(2)의 분포함수 F2(t)=12(1+t2+t2)F_2(t)=\tfrac12\bigl(1+\tfrac t{\sqrt{2+t^2}}\bigr)를 써도 된다.

힌트

포함되는 사건은 ∣T2∣≤1.96|T_2|\le1.96이다.

풀이

포함되는 사건은 ∣T∣≤1.96|T|\le1.96, T∼t(2)T\sim t(2)이다. P(∣T∣≤c)=F2(c)−F2(−c)=c2+c2P(|T|\le c)=F_2(c)-F_2(-c)=\tfrac c{\sqrt{2+c^2}}이므로 1.965.8416≈0.811\tfrac{1.96}{\sqrt{5.8416}}\approx0.811이다. “95%” 구간이 다섯 번에 한 번 가까이 틀린다.

최종답 1.962+1.962≈0.811\tfrac{1.96}{\sqrt{2+1.96^2}}\approx0.811