- 다양한 주제에 대해 자유롭게 글을 작성하는 게시판입니다.
Date 15/11/02 11:16:29
Name   Beer Inside
Subject   베이즈 정리, 몬티홀의 문제, 삶과 죽음의 확률


잊어 먹을만 하면 한번 씩 나와서 답을 아는 사람과 모르는 사람 모두 암걸리게 하는 문제입니다.

1970년대 방송인 몬티 홀이 진행하는 퀴즈쇼에서 일어난 상황인데 1990년 칼럼니스트 마릴린 사반트가
잡지 ‘퍼레이드’에서 이 문제를 질문한 독자의 편지에 대해 답을 하면서 유명해졌습니다.

한 곳에는 스포츠카가 나머지 두 곳에는 염소가 들어있습니다.
문을 열었을 때 스포츠카가 있을 확률은 3분의 1. 퀴즈 참가자가 1번 문을 찍었다.
이때 ‘스포츠카가 어디에 있는지 알고 있는’ 홀이 3번 문을 활짝 열었고 염소가 모습을 드러냈다.

그리고 홀이 참가자에게 물었습니다.

  “선택을 바꾸시겠습니까?”

몬티홀 문제의 핵심은 참가자가 새로 얻게 된 정보(3번 방에는 스포츠카가 없다)를 어떻게 추론에 반영하느냐 하는 것입니다. 

이 때 사용하는 것이 베이즈 정리입니다. 

베이즈는 영국의 목사로 살아 생전 단 두편의 논문만 출간했는데, 한편은 종교 다른 한편은 수학이였다고 합니다.

이 베이즈가 죽고난 후 친구인 리차드 프라이스가 정리해서 발표했습니다.

베이즈(싸이코의 노먼 베이츠 아닙니다.) 정리는 이전의 경험과 현재의 증거를 토대로 어떤 사건의 확률을 추론하는 알고리즘으로, 
극단적으로 단순화하면 먼저 일어나는 사건 A가 있고 시간상 다음에 일어나는 사건 B가 있을 때 
사건 B가 일어난 것을 전제로 한 사건 A의 조건부 확률을 구하는 것이라고 생각할 수 있습니다.



(이해 하기 쉬우라고 수식을 쓰지만 수식만 보면 토할 것 같습니다.)

직관적으로 생각하면 각 방에 스포츠카(C)가 있을 확률은 3분의 1로 똑 같고 (P(C1)=P(C2)=P(C3)=1/3)
참가자가 일단 1번 방을 선택한 뒤 진행자가 3번 방을 열었기 때문에(따라서 3번 방은 아니다)
1번 방과 2번 방에 스포츠카가 있을 확률이 1/2로 똑 같을 것 같습니다.
따라서 굳이 선택을 바꿀 필요는 없을 것 같습니다.
바꿔도 기대 확률은 마찬가지이기 때문입니다.

하지만, 마음 속 깊은 곳이나 뒷골 어디인가에서 바꾸고 싶다.
바꾼다고 특별히 달라지는 것은 없을 것 같지만 격렬하게 바꾸고 싶다는 생각이 들지 않습니까?

베이즈 정리에 따르면 선택을 2번으로 바꾸는 것이 답입니다.

이 경우 맞출 확률이 3분의 2로 2배나 높아지기 때문입니다.

이 칼럼이 나가고 설명을 이해할 수 없다는 독자 편지가 쇄도했고 전문가들 사이에서도 논란이 일었습니다.
심지어 폴 에르되시 같은 일급 수학자조차 “왜 선택을 바꿔야하는지 이해하지 못 하겠다”는 반응을 보였습니다.
(그러니 유머란에 이 문제가 다시 나와도 암걸릴 것 같은  사람은 또 생길겁니다.)

몬티홀 문제의 핵심은 참가자가 새로 얻게 된 정보(3번 방에는 스포츠카가 없다)를 어떻게 추론에 반영하느냐 하는 것입니다.

이때 베이즈 정리를 쓰면 선택을 바꿔야 하는 깔끔하게 이유가 설명됩니다.

새로운 정보, 즉 진행자가 3번 방을 열었을 때(O3) 1번 방에 스포츠카(C1)가 있을 확률은 다음의 베이즈 정리로 나타낼 수 있다.

  P(C1|O3)=P(O3|C1)P(C1)/P(O3)=(1/2)×(1/3)/(1/2)=1/3

차가 1번에 있다면 진행자는 2번이나 3번 문을 열 수 있습니다.

따라서 P(O3|C1)는 2분의 1입니다.

한편 참가자의 관점에서 자기가 1번 방을 선택했기 때문에 진행자는 2번이나 3번 문을 열 수 밖에 없으므로 P(O3) 역시 2분의 1입니다.

이제 진행자가 3번 방을 열었을 때(O3) 2번 방에 스포츠카가 있을 확률을 베이즈 정리로 구해봅시다.

차가 2번에 있다면 진행자는 3번 문을 열 수 밖에 없습니다. 따라서 P(O3|C2)는 1입니다.

  P(C2|O3)=P(O3|C2)P(C2)/P(O3)=1×(1/3)/(1/2)=2/3

결국 홀이 3번 문을 열고 난 뒤, 즉 새로운 정보가 알려진 뒤 2번 방에 스포츠카가 있을 확률은 3분의 2로 2배 높아집니다.

따라서 참가자가 자신의 감을 믿지 않고 순전히 확률이 높은 쪽을 택하기로 했다면 무조건 2번으로 선택을 바꿔야 합니다.

그깟 몬티홀 문제 풀지 않으면 그만 아니냐구요.

그럼 이 문제는 어떤가요?

45세 골드미스 A씨는 자신의 가슴을 만져보다가 이상한 느낌이 들어서 병원에 갔습니다. 
(몬티홀 문제를 너무 오래 쳐다 보았나 봅니다.) 

유방암 검사를 했고 검사 결과는 양성이였습니다.

그리고 이 검사의 검사 정확도는 90%입니다.

A씨는 결혼도 못하고 아이도 없는데 유방암이라는 진단이 나와서 기가 막혀서 진료실에서 눈물을 흘렸습니다.

하지만 이 놈의 의사 자식은 웃으면서 이 상태에서는 유방암일 확률이 10%도 되지 않으니 추가 검사를 해 보자고 합니다.

자꾸 돈들게 말입니다.

이 싸이코패스 같은  의사는 유방암에 걸린 여성은 성인 여성의 1% 수준이고 (실제는 10만명당 321명이 발생합니다. 대략 0.3%)

검사 정확도가 90%이므로 정상인데도 검사에서 유방암에 걸린 것으로 나올 확률은 10%다.

따라서 설사 검사에서 양성으로 나왔더라도 진짜 유방암에 걸렸을 확률은 8%에 불과하다는 겁니다.
(암이라는데 확률강의를 하는 싸이코패스 의사야!)




A씨의 사례는 ‘양성반응일 때 유방암일 확률’을 구하는 것 입니다.

즉 양성반응이라는 조건에서 유방암일 확률을 ‘P(암|양성)’으로 나타내면,  

  P(암|양성)=P(양성|암)P(암)/P(양성) 
               
P(양성|암) = 0.9
P(암) = 0.01
P(양성) = 암일 때 양성일 확률 + 암이 아닐 때 양성일 확률  = 0.9*0.01+0.1*0.99= 0.108
                   
따라서 P(암|양성) = 0.9* 0.01/0.108 = 0.083 

즉 검사에서 양성일 경우 유방암일 확률은 8.3%가 됩니다. 

만약 이 검사에서 음성으로 나오면 음성예측치는 90%이므로 추가 검사할 필요가 없다고 보아야겠지요. 

이 환자는 8.3%의 확률을 확인하기 위해서 조직 검사를 했고 조직검사상 양성이 나왔습니다. 
(조직검서의 민감도는 100%라고 합시다. Borderline pathologisit가 아니라면...) 

이후 A씨는 잘 지내다가 5년 후인 50세에 다시 가슴에 무엇인가 느껴집니다. 

이번에는 별것 아니라는 생각에 병원에 방문해서 저번에 그 의사가 미친놈 같았지만 결과가 좋았으므로

다시 그 의사에게 유방암 검사를 합니다. 

이번에도 양성이 나왔는데, 의사의 표정이 좋지 않습니다. 

A씨는 이 넘이 오늘 기분이 좋지 않나 생각을 하는데, 이 싸이코패스가 또 통계 강의를 합니다.  

5년 동안 유방암 진단기술이 발전해서 민감도가 99%로 증가했다는 것 입니다.  






P(양성|암) = 0.99
P(암) = 0.01
P(양성) = 암일 때 양성일 확률 + 암이 아닐 때 양성일 확률  = 0.99*0.01+0.01*0.99= 0.0198

P(암|양성) = 0.99*0.01/0.0198 = 0.5 

따라서 암일 확률은 50%라고 합니다. 

그러니 의사를 잘 만나야 합니다. 

이런 설명을 자세하게 해 주는 의사는 높은 확률로 싸이코패스일겁니다. 

하지만 스텐포드의 에프론 교수는 무시무시한 Impact Factor를 가진 Science에서 

오늘날 응용통계를 이용한 논문의 25%는 베이즈 정리를 이용하고 있다고 합니다. 


베이즈 정리가 나온지 250년이 넘었고,  최근들어 여러 문제를 해결하는데 자주 쓰이고 있지만 이에 대한 비판도 있습니다. 

베이즈 정리는 사전 정보(prior)가 확실한 것일 때만 성립하는 것인데 실제 상황에서는 이 정보가 100% 확실한 경우가 별로 없기 때문입니다. 

결국 불확실한 사전 정보를 토대로 사후 확률(posterior)을 추측하는 것이 타당한가에 대한 입장이 엇갈리고 있습니다. 

이글은 아래 글을 참고해서 작성했습니다.  




7


    목록
    번호 제목 이름 날짜 조회 추천
    60 기타살아오며 가장 부끄러운 일에 대하여 29 DEICIDE 15/05/30 11476 14
    1323 IT/컴퓨터인공지능, 고용 없는 성장 19 kpark 15/10/23 11474 1
    3708 의료/건강니코틴과 히로뽕 이야기 5 모모스 16/09/15 11465 6
    667 요리/음식연인 혹은 아내에게 선물하기 좋은 의미를 가진 와인(~3만원이하) 24 마르코폴로 15/07/27 11452 0
    277 기타한인 천재소녀는 진짜인가 37 Azurespace 15/06/09 11450 0
    192 기타어제 있었던 "메르스 확산, 어떻게 막을 것인가?" 에 대한 토론 전문입니다.(스압) 2 아나키 15/06/03 11449 0
    1027 일상/생각분탕질이란 단어에 대한 생각 20 王天君 15/09/18 11433 0
    262 기타제가 생각하는 최고의 게임 Best 5 12 김치찌개 15/06/08 11431 0
    709 IT/컴퓨터컴퓨터 파일 관리 어떻게 하시나요? 12 한아 15/08/03 11416 0
    2808 의료/건강정맥주사라인 잡기 이야기 21 원하 16/05/14 11413 0
    1736 IT/컴퓨터막귀 입장에서 써보는 5만원 이하 스피커 사용기 21 전크리넥스만써요 15/12/09 11408 2
    1488 음악Janis Ian -Stars 6 새의선물 15/11/07 11380 0
    1426 꿀팁/강좌통계, 그 바이어스의 가면을 벗겨보자 28 눈부심 15/11/02 11374 1
    562 역사피맛골, 사라진 골목에 대한 아쉬움 18 마르코폴로 15/07/10 11332 0
    5950 기타필름포장지 이야기 18 헬리제의우울 17/07/14 11331 15
    335 기타쥬라기 월드의 흥행에 적어보는 공룡이야기....왜 새가 공룡의 후예인가 13 개평3냥 15/06/15 11321 0
    1376 요리/음식라면 41 헬리제의우울 15/10/29 11309 25
    1375 영화우리 2017년에 만나요...가능하면 말이죠... 17 Neandertal 15/10/29 11303 0
    7035 IT/컴퓨터딥러닝으로 만든 유인나 오디오북 10 Toby 18/02/02 11299 1
    1669 역사러시아에서 푸틴의 인기에 대하여... 9 조홍 15/11/30 11292 6
    1428 기타베이즈 정리, 몬티홀의 문제, 삶과 죽음의 확률 25 Beer Inside 15/11/02 11292 7
    1414 음악Halloween을 맞이하여 Helloween을 소개합니다. 24 맥주만땅 15/11/01 11286 0
    1712 문화/예술니 꺼는 내 꺼, 땡큐 18 눈부심 15/12/05 11285 1
    7213 기타여성의 날을 맞이하여 써 보는 한 가상의(?) 사례 20 烏鳳 18/03/08 11277 13
    1133 경제한국형 시장경제체제 6 MANAGYST 15/09/30 11277 1
    목록

    + : 최근 2시간내에 달린 댓글
    + : 최근 4시간내에 달린 댓글

    댓글