- 회원들이 추천해주신 좋은 글들을 따로 모아놓는 공간입니다.
- 추천글은 매주 자문단의 투표로 선정됩니다.
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)을 추측하는 것이 타당한가에 대한 입장이 엇갈리고 있습니다. 

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



* 수박이두통에게보린님에 의해서 자유 게시판으로부터 게시물 복사되었습니다 (2015-11-16 21:21)
* 관리사유 :



8


    목록
    번호 제목 이름 날짜 조회 추천
    186 음악홍차넷 지상파 입성 기념 뮤직비디오 241 Toby 16/04/20 13793 9
    795 의료/건강오늘 받은 정관수술 후기 21 미스터주 19/04/17 13822 37
    44 요리/음식이탈리안 식당 주방에서의 일년(4) - 토마토소스만들기 29 뤼야 15/07/09 13851 0
    45 꿀팁/강좌캐나다 영주권에 대한 이런저런 이야기 16 이젠늙었어 15/07/10 13914 0
    1123 기타남자바지3종(청바지,검은바지,베이지면바지) 입는방법에 대한 연구 22 흑마법사 21/08/31 14283 17
    781 여행타베로그 이용 팁 8 温泉卵 19/03/18 14285 12
    127 의료/건강의전은 어떻게 실패했는가 ? 41 Zel 15/12/09 14289 2
    274 IT/컴퓨터컴퓨터는 어떻게 빠르게 검색을 할까 - 보이어-무어-호스풀 알고리즘 18 April_fool 16/10/04 14561 1
    473 기타필름포장지 이야기 24 헬리제의우울 17/07/14 14602 3
    263 게임[삼국지 영걸전] 1599 클리어 기념 팁 + 후기와 기타 등등 이야기 37 조홍 16/09/09 14607 8
    412 과학음수X음수는 왜 양수인가요? 62 캡틴아메리카 17/04/15 14910 14
    205 요리/음식덴뿌라와 튀김의 기원 29 마르코폴로 16/05/14 15307 11
    800 일상/생각불안 애착 유형과 회피 애착 유형의 연애/이별기 4 자일리톨 19/05/01 15400 17
    41 기타하고 싶은 일이 무언지 모르겠다는 그대에게 32 ArcanumToss 15/07/08 18978 3
    37 게임독수리의 눈으로 입문자를 노리는 HOMM3 소개 (1) 기초 14 조홍 15/06/22 19230 0
    107 기타베이즈 정리, 몬티홀의 문제, 삶과 죽음의 확률 27 Beer Inside 15/11/02 20418 8
    236 IT/컴퓨터어느 게임 회사 이야기 (1) 26 NULLPointer 16/07/19 22092 29
    415 요리/음식집에서 스테이크 구우면 맛이가 없는 이유 37 Morpheus 17/04/19 22618 18
    15 꿀팁/강좌책장에서 책을 치우자! (북스캐너 + 스캔이북 써먹기) 21 damianhwang 15/06/08 31011 1
    69 IT/컴퓨터거장의 그림을 30초만에 만들다: DeepStyle 33 Azurespace 15/09/08 31901 10
    1237 일상/생각만년필 덕후가 인정하는 찰스 3세의 착한 빡침 95 SCV 22/09/13 32515 49
    11 체육/스포츠남성의 정력을 증강시키는 운동 69 스타-로드 15/06/05 43253 0
    695 정치/사회강제추행으로 법정구속되었다는 판결문 감상 - 랴 리건.... 30 烏鳳 18/09/07 50951 85
    737 꿀팁/강좌지금 쓰는 안경에만 돈을 75만원씩 퍼부은 사람이 알려주는 안경 렌즈 선택의 거의 모든 것 22 April_fool 18/11/28 61833 52
    869 일상/생각따뜻함에 대해서 22 19/09/29 246008 27
    목록

    + : 최근 6시간내에 달린 댓글
    + : 최근 12시간내에 달린 댓글

    댓글