Skip to content

11회 · 충분성 · 워크북

D. 퀄형 ​

여러 절의 도구를 잇는 종합 문제다. 머리의 제한 시간을 재고, 노트를 덮고 푼다. UMVUE를 구하는 문제는 “완비충분통계량 찾기 → 불편추정량 하나 → 조건부기댓값”의 순서를 답안에 드러낸다.

D1 지수분포의 생존확률 ★★★ ​

정리 3.3, 5.2, 4회 정리 4.2 · 제한 40분

X1,…,Xn∼iidX_1,\dots,X_n\iid Exp(평균 θ\theta), n≥2n\ge2이고 c>0c>0은 알려진 상수다. τ(θ)=Pθ(X1>c)=e−c/θ\tau(\theta)=P_\theta(X_1>c)=e^{-c/\theta}의 UMVUE를 구한다. S=∑XiS=\sum X_i이다.

  1. (a)

    SS가 완비충분통계량임을 보여라.

  2. (b)

    X1/S∼X_1/S\sim 베타(1,n−1)(1,n-1)이고 X1/SX_1/S가 SS와 독립임을 보여라.

  3. (c)

    P(X1>c∣S=s)=(1−c/s)n−1P(X_1>c\mid S=s)=(1-c/s)^{n-1} (s>cs>c), 00 (s≤cs\le c)임을 보이고 UMVUE를 적어라.

  4. (d)

    n=5n=5, ∑xi=10\sum x_i=10, c=1c=1이면 UMVUE의 값은 얼마인가? MLE e−c/xˉe^{-c/\bar x}와 비교하라.

힌트

(b) X1X_1은 감마(1,θ)(1,\theta), ∑i≥2Xi\sum_{i\ge2}X_i는 감마(n−1,θ)(n-1,\theta)이고 독립이다. 4회의 베타–감마 관계. 독립은 Basu로도 나온다(X1/SX_1/S는 척도에 불변). (c) 베타(1,n−1)(1,n-1)의 꼬리는 P(B>u)=(1−u)n−1P(B>u)=(1-u)^{n-1}.

풀이

(a) 밀도 1θe−x/θ\tfrac1\theta e^{-x/\theta}는 자연모수 η=−1/θ∈(−∞,0)\eta=-1/\theta\in(-\infty,0), T(x)=xT(x)=x인 지수족이고 자연모수 공간이 열린 구간이므로 정리 4.5로 SS는 완비충분이다.

(b) X1∼X_1\sim 감마(1,θ)(1,\theta)와 X2+⋯+Xn∼X_2+\cdots+X_n\sim 감마(n−1,θ)(n-1,\theta)는 독립이므로 4회 정리 4.2(베타와 감마)로 X1/S∼X_1/S\sim 베타(1,n−1)(1,n-1)이고 SS와 독립이다. Basu로도 같다. X1/SX_1/S는 Xi↦bXiX_i\mapsto bX_i에 대해 변하지 않아 보조이고 SS는 완비충분이다.

(c) W=1{X1>c}W=1\{X_1>c\}는 τ\tau의 불편추정량이다. (b)의 독립으로 P(X1>c∣S=s)=P(X1/S>c/s∣S=s)=P(B>c/s)P(X_1>c\mid S=s)=P(X_1/S>c/s\mid S=s)=P(B>c/s)이고, B∼B\sim 베타(1,n−1)(1,n-1)의 밀도 (n−1)(1−u)n−2(n-1)(1-u)^{n-2}를 적분하면 (1−c/s)n−1(1-c/s)^{n-1} (c<sc<s)이다. s≤cs\le c이면 X1≤S≤cX_1\le S\le c이므로 0이다. 정리 3.3으로 UMVUE는 φ(S)=(1−c/S)n−11{S>c}\varphi(S)=(1-c/S)^{n-1}1\{S>c\}이다.

(d) (1−0.1)4=0.6561(1-0.1)^4=0.6561. MLE는 e−1/2=0.6065e^{-1/2}=0.6065이다. 표본이 커지면 (1−c/S)n−1≈e−c(n−1)/S→e−c/θ(1-c/S)^{n-1}\approx e^{-c(n-1)/S}\to e^{-c/\theta}로 두 추정량이 가까워진다.

최종답 (c) (1−c/S)n−1 (S>c)(1-c/S)^{n-1}\ (S>c) (d) 0.94=0.65610.9^4=0.6561

D2 Poisson에서 “정확히 한 번”의 확률 ★★★ ​

정리 2.1, 3.3, 4.5 · 제한 30분

X1,…,Xn∼iidPoi(λ)X_1,\allowbreak \dots,\allowbreak X_n\iid\mathrm{Poi}(\lambda)이고 τ(λ)=P(X=1)=λe−λ\tau(\lambda)=P(X=1)=\lambda e^{-\lambda}이다. T=∑XiT=\sum X_i이다.

  1. (a)

    τ\tau의 불편추정량을 하나 찾아라.

  2. (b)

    (a)에 TT로 조건을 걸어 UMVUE φ(T)\varphi(T)를 구하라(답의 변수는 tt와 nn).

  3. (c)

    φ\varphi가 UMVUE인 근거를 정리 번호로 적어라.

  4. (d)

    n=10n=10, T=7T=7이면 UMVUE의 값은 얼마인가?

  5. (e)

    같은 자료에서 MLE τ(xˉ)\tau(\bar x)의 값을 구하고 비교하라.

힌트

W=1{X1=1}W=1\{X_1=1\}. T=tT=t가 주어지면 X1∼b(t,1n)X_1\sim\mathrm b(t,\tfrac1n)이다.

풀이

(a) W=1{X1=1}W=1\{X_1=1\}은 EλW=λe−λE_\lambda W=\lambda e^{-\lambda}이므로 불편이다. (b) T=tT=t가 주어지면 X1∼b(t,1n)X_1\sim\mathrm b(t,\tfrac1n)이므로(노트 연습 1.1) φ(t)=P(X1=1∣T=t)=t⋅1n(1−1n)t−1\varphi(t)=P(X_1=1\mid T=t)=t\cdot\tfrac1n\bigl(1-\tfrac1n\bigr)^{t-1}이다. (c) TT는 1-모수 지수족의 완비충분통계량이고(정리 4.5), φ(T)\varphi(T)는 불편추정량에 조건을 건 것이므로 불편이다(정리 2.1(b)). 정리 3.3으로 유일한 UMVUE다. (d) 0.7⋅0.96=0.37200.7\cdot0.9^6=0.3720이다. (e) xˉ=0.7\bar x=0.7이므로 0.7e−0.7=0.34760.7e^{-0.7}=0.3476이다. 두 값의 차는 0.024이고, λ=0.7\lambda=0.7에서 UMVUE의 표준편차 0.051보다 작다.

최종답 (b) Tn(1−1n)T−1\tfrac Tn(1-\tfrac1n)^{T-1} (d) 0.7⋅0.96≈0.37200.7\cdot0.9^6\approx0.3720 (e) 0.7e−0.7≈0.34760.7e^{-0.7}\approx0.3476

D3 정규분포의 σ\sigma와 μ2\mu^2 ★★★ ​

정리 3.3, 4.5, 6회 정리 2.2 · 제한 35분

X1,…,Xn∼iidN(μ,σ2)X_1,\allowbreak \dots,\allowbreak X_n\iid\N(\mu,\allowbreak \sigma^2), n≥2n\ge2이다.

  1. (a)

    ES=σ2n−1 Γ(n/2)Γ((n−1)/2)ES=\sigma\sqrt{\tfrac2{n-1}}\,\dfrac{\Gamma(n/2)}{\Gamma((n-1)/2)}임을 보여라.

  2. (b)

    σ\sigma의 UMVUE가 cnSc_nS, cn=n−12 Γ((n−1)/2)Γ(n/2)c_n=\sqrt{\tfrac{n-1}2}\,\dfrac{\Gamma((n-1)/2)}{\Gamma(n/2)}임을 보여라.

  3. (c)

    c5c_5를 닫힌 꼴과 소수로 구하라.

  4. (d)

    μ2\mu^2의 UMVUE를 구하고, n=5n=5, xˉ=3\bar x=3, s2=2s^2=2일 때 값을 구하라.

  5. (e)

    (d)의 UMVUE가 음수가 될 수 있는가? 7절의 터무니없는 UMVUE와 비교하라.

힌트

(a) (n−1)S2/σ2∼χ2(n−1)(n-1)S^2/\sigma^2\sim\chi^2(n-1)이므로 EVE\sqrt{V}, V∼χ2(k)V\sim\chi^2(k)를 감마 적분으로 구한다. (d) EXˉ2=μ2+σ2/nE\bar X^2=\mu^2+\sigma^2/n.

풀이

(a) V=(n−1)S2/σ2∼χ2(n−1)V=(n-1)S^2/\sigma^2\sim\chi^2(n-1)이다(6회 정리 2.2). k=n−1k=n-1이면

EV=∫0∞v1/2vk/2−1e−v/2Γ(k/2)2k/2dv=Γ((k+1)/2)2(k+1)/2Γ(k/2)2k/2=2 Γ(n/2)Γ((n−1)/2)E\sqrt V=\int_0^\infty v^{1/2}\frac{v^{k/2-1}e^{-v/2}}{\Gamma(k/2)2^{k/2}}dv=\frac{\Gamma((k+1)/2)2^{(k+1)/2}}{\Gamma(k/2)2^{k/2}}=\sqrt2\,\frac{\Gamma(n/2)}{\Gamma((n-1)/2)}

이고 S=σV/(n−1)S=\sigma\sqrt{V/(n-1)}이므로 주어진 식이 나온다. (b) (a)로 cnSc_nS는 불편이고 (Xˉ,S2)(\bar X,S^2)는 완비충분이므로(정리 4.5) 정리 3.3으로 UMVUE다. (c) c5=2 Γ(2)Γ(5/2)=2⋅43π=423π≈1.0638c_5=\sqrt2\,\tfrac{\Gamma(2)}{\Gamma(5/2)}=\sqrt2\cdot\tfrac{4}{3\sqrt\pi}=\tfrac{4\sqrt2}{3\sqrt\pi}\approx1.0638이다. SS는 σ\sigma를 약 6% 작게 추정한다. (d) EXˉ2=μ2+σ2/nE\bar X^2=\mu^2+\sigma^2/n, ES2=σ2ES^2=\sigma^2이므로 Xˉ2−S2/n\bar X^2-S^2/n은 불편이고 UMVUE다. 값은 9−0.4=8.69-0.4=8.6이다. (e) Xˉ2<S2/n\bar X^2<S^2/n이면 음수다. μ2≥0\mu^2\ge0인데 음수로 추정할 수 있다. 노트 예제 7.2의 (−1)X(-1)^X처럼 불편이라는 조건이 값의 범위를 지키지 못한다. max⁡(Xˉ2−S2/n,0)\max(\bar X^2-S^2/n,0)은 평균제곱오차가 더 작지만 불편이 아니다.

최종답 (c) c5=423π≈1.0638c_5=\tfrac{4\sqrt2}{3\sqrt\pi}\approx1.0638 (d) 8.68.6

D4 균등분포의 분포함수 값 ★★★ ​

정리 2.1, 3.3 · 제한 35분

X1,…,Xn∼iidU(0,θ)X_1,\allowbreak \dots,\allowbreak X_n\iid\mathrm U(0,\allowbreak \theta)이고 c>0c>0은 알려진 상수다. τ(θ)=Pθ(X1≤c)=min⁡(c/θ,1)\tau(\theta)=P_\theta(X_1\le c)=\min(c/\theta,\allowbreak 1)의 UMVUE를 구한다. M=X(n)M=X_{(n)}이다.

  1. (a)

    M=mM=m이 주어지면 확률 1n\tfrac1n로 X1=mX_1=m이고, 나머지 확률로 X1∼U(0,m)X_1\sim\mathrm U(0,m)임을 설명하라.

  2. (b)

    φ(m)=E[1{X1≤c}∣M=m]\varphi(m)=E\bigl[1\{X_1\le c\}\bigm|M=m\bigr]을 구하라.

  3. (c)

    φ(M)\varphi(M)이 UMVUE인 근거를 적어라.

  4. (d)

    n=4n=4, M=5M=5, c=2c=2이면 UMVUE의 값은 얼마인가?

힌트

(a) 대칭으로 최댓값이 X1X_1일 확률은 1n\tfrac1n이고, 최댓값이 아닌 관측은 (0,m)(0,m)에서 균등하다(8회의 순서통계량). (b) m≤cm\le c와 m>cm>c를 나눈다.

풀이

(a) 표본은 교환가능하므로 nn개 가운데 X1X_1이 최댓값일 확률은 1n\tfrac1n이다. X1X_1이 최댓값이 아니면, M=mM=m이 주어질 때 최댓값이 아닌 n−1n-1개는 (0,m)(0,m)에서 독립인 균등분포다(독일 전차의 연속판, 노트 연습 2.2와 같은 논리). (b) m≤cm\le c이면 X1≤M≤cX_1\le M\le c이므로 φ(m)=1\varphi(m)=1이다. m>cm>c이면 φ(m)=1n⋅0+n−1n⋅cm=(n−1)cnm\varphi(m)=\tfrac1n\cdot0+\tfrac{n-1}n\cdot\tfrac cm=\tfrac{(n-1)c}{nm}이다. (c) 1{X1≤c}1\{X_1\le c\}는 불편이고 MM은 완비충분이므로(예제 3.2(c)) 정리 2.1과 3.3으로 φ(M)\varphi(M)이 UMVUE다. 직접 확인: θ>c\theta>c이면 ∫0cnmn−1θndm+∫cθ(n−1)cnmnmn−1θndm=cnθn+c(θn−1−cn−1)θn=cθ\int_0^c\tfrac{nm^{n-1}}{\theta^n}dm+\int_c^\theta\tfrac{(n-1)c}{nm}\tfrac{nm^{n-1}}{\theta^n}dm=\tfrac{c^n}{\theta^n}+\tfrac{c(\theta^{n-1}-c^{n-1})}{\theta^n}=\tfrac c\theta이다. (d) 3⋅24⋅5=310\tfrac{3\cdot2}{4\cdot5}=\tfrac3{10}이다. MLE c/M=0.4c/M=0.4보다 작다.

최종답 (d) 310\tfrac3{10}