연습문제 6.23

(문제).  동전을 두 번 던졌을 때 가능한 iid 표본 x는 00, 01, 10, 11의 4 경우가 있다. 이 모든 경우 각각에 대하여 가능도, 주변확률 및 사후확률을 구하여라. 그리고 각 경우에 MLE 추정량, MAP 추정량 및 베이즈 기댓값 추정량을 구하여라.

이 문제는 [예제 6.20]의 연장이므로 문제가 요구하는 것이 무엇인지를 정확히 알기 위해서는 [예제 6.20]을 읽어 보아야만 한다.

모수 θ는 동전을 던졌을 때 앞면이 나올 확률이다. 앞면을 1, 뒷면을 0으로 나타낸다. 그러면 정의에 의해서 π(1|θ)=1, π(0|θ)=1θ 이 된다. 동전을 던지기 전에는 아무런 정보가 없으므로 사전확률을 π(θ)=1,(0θ1)로 둔다.

이 문제에서 사용하는 베이즈 규칙은 다음과 같다.

(1)π(θ|x)=f(x|θ)×π(θ)m(x), where m(x)=f(x|θ)×π(θ)dθ:주변확률, marginal probabilityπ(θ|x):사후확률π(θ):사전확률f(x|θ):가능도

(풀이).   [예제 6.20]에서는 x=01일 때의 사후확률 π(θ|01)=6θ(1θ)과 MAP 추정량 12를 구해 놓았다.

가능도 f(01|θ)는 첫 번째 관측과 두 번째 관측이 독립이므로 f(0|θ)f(1|θ)=(1θ)θ가 된다.

주변확률 m(01)을 구하기 위하여 (1)의 양변을 정적분하면 좌변은 항상 1이 되고 우변은 01(1θ)θdθ/m(01) 이므로 m(01)=01(1θ)θdθ=16을 얻는다.

MLE 추정량은 argmaxθθ(1θ)=12로 MAP 추정량과 일치한다.

베이즈 기댓값 추정량은 01θ6θ(1θ)dθ=12로 역시 MAP 추정량과 같다.

x가 00, 10, 11인 경우들에 대해서 알아보자.

주변확률은 항상, 아까와 마찬가지로 (1)의 양변을 정적분하여 얻는다. 좌변은 1이 되고 우변은 가능도의 적분값을 주변확률로 나눈 값이므로 다음과 같은 결과를 얻는다.

m(00)=01f(0|θ)2dθ=01(1θ)2dθ=13,L(θ)=(1θ)2m(10)=01f(1|θ)f(0|θ)dθ=01θ(1θ)dθ=16,L(θ)=θ(1θ)m(11)=01f(1|θ)2dθ=01θ2dθ=13,L(θ)=θ2

사후확률은 단순히 가능도를 주변확률로 나눈 것이다. 왜냐하면 사전확률이 1이므로.

π(θ|00)=3(1θ)2,θ^=0π(θ|10)=6θ(1θ),θ^=12π(θ|11)=3θ2,θ^=1

MLE 추정량은 항상 MAP 추정량과 일치한다. 왜냐하면 사전확률과 사후확률은 정비례하고 있으므로.

베이즈 기댓값 추정량은 x=10인 경우는 이전과 같이 12 임을 금방 알 수 있다. 나머지 경우에는 다음과 같이 계산된다.

x=00:013θ(1θ)2dθ=14x=11:013θ3dθ=34 [홈으로]