Wszystkie zagadnienia: Prawdopodobieństwo i statystyka
zaawansowany matura rozszerzona z matematykiwzór Bayesaprawdopodobieństwo warunkowedrzewo stochastyczne

Wzór Bayesa

01

Opis

Wzór Bayesa pozwala odwrócić warunek. Znasz prawdopodobieństwo, że test wyjdzie dodatni u osoby chorej, i chcesz wiedzieć, jak bardzo prawdopodobne jest, że osoba z dodatnim wynikiem jest chora. To dwa różne pytania: P(+∣C)P(+\mid C) i P(C∣+)P(C\mid +). Zwykle znamy to pierwsze („od przyczyny do skutku"), a pytają o drugie („od skutku do przyczyny"). Na maturze rozszerzonej poznasz to po zdaniu typu „oblicz prawdopodobieństwo, że bochenek pochodzi z piekarni P1P_1, jeśli wiadomo, że jest razowy".

Wzór składa się z elementów, które już znasz z haseł prawdopodobieństwo warunkowe i prawdopodobieństwo całkowite. Licznik to prawdopodobieństwo jednej ścieżki na drzewie (przyczyna BkB_k, potem skutek AA). Mianownik to suma wszystkich ścieżek prowadzących do skutku AA, czyli prawdopodobieństwo całkowite P(A)P(A). Wynik mówi, jaką część wszystkich ścieżek do AA stanowi ścieżka przez BkB_k.

Zasady, których pilnujesz:

  • zdarzenia B1,…,BnB_1, \ldots, B_n (przyczyny) muszą być parami rozłączne i razem dawać całą przestrzeń Ω\Omega, więc ich prawdopodobieństwa sumują się do 11,
  • skutek AA musi mieć P(A)>0P(A) > 0, bo przez to prawdopodobieństwo dzielimy,
  • prawdopodobieństwo wstępne P(Bk)P(B_k) (np. odsetek chorych w populacji) wchodzi do wzoru na równi z P(A∣Bk)P(A \mid B_k) – pominięcie go to najczęstszy błąd,
  • wzór jest wymaganiem wyłącznie rozszerzenia; na poziomie podstawowym go nie ma.
02

Wzory

Wzór Bayesa karta CKE
P(Bk∣A)=P(Bk)⋅P(A∣Bk)P(A∣B1)⋅P(B1)+…+P(A∣Bn)⋅P(Bn)P(B_k \mid A) = \frac{P(B_k) \cdot P(A \mid B_k)}{P(A \mid B_1) \cdot P(B_1) + \ldots + P(A \mid B_n) \cdot P(B_n)}
Dla 1≤k≤n1 \leq k \leq n, przy BiB_i rozłącznych, B1∪…∪Bn=ΩB_1 \cup \ldots \cup B_n = \Omega, P(Bi)>0P(B_i) > 0 i P(A)>0P(A) > 0.
Dwie przyczyny
P(B∣A)=P(B)⋅P(A∣B)P(B)⋅P(A∣B)+P(B′)⋅P(A∣B′)P(B \mid A) = \frac{P(B) \cdot P(A \mid B)}{P(B) \cdot P(A \mid B) + P(B') \cdot P(A \mid B')}
Najczęstszy przypadek, gdzie B′B' to zdarzenie przeciwne do BB (więcej w haśle zdarzenie przeciwne i suma zdarzeń): choroba i zdrowie, wyższe wykształcenie i brak, urna I i urna II. Wynika ze wzoru wyżej dla n=2n = 2.
Mianownik karta CKE
P(A)=P(A∣B1)⋅P(B1)+…+P(A∣Bn)⋅P(Bn)P(A) = P(A \mid B_1) \cdot P(B_1) + \ldots + P(A \mid B_n) \cdot P(B_n)
To twierdzenie o prawdopodobieństwie całkowitym; w Bayesie jest mianownikiem.
Definicja warunkowego karta CKE
P(Bk∣A)=P(Bk∩A)P(A)P(B_k \mid A) = \frac{P(B_k \cap A)}{P(A)}
Z niej wynika cały wzór, bo P(Bk∩A)=P(Bk)⋅P(A∣Bk)P(B_k \cap A) = P(B_k) \cdot P(A \mid B_k).
03

Tabela

Test przesiewowy z zadania CKE, przeliczony na 100 000 osób (choroba 0,2%; test dodatni u 99% chorych i u 2% zdrowych)
GrupaLiczba osóbWynik dodatniWynik ujemny
chorzy2001982
zdrowi99 8001 99697 804
razem100 0002 19497 806
chorzy wśród dodatnich–1982194≈0,09\frac{198}{2194} \approx 0{,}09–
04

Przykład

Przykład 1 – test przesiewowy. Pewna choroba dotyczy 0,2%0{,}2\% populacji. Test daje wynik dodatni u 99%99\% chorych i wynik ujemny u 98%98\% zdrowych. Pan X ma wynik dodatni. Oblicz prawdopodobieństwo, że jest chory.

Najpierw nazywamy zdarzenia: CC – osoba chora, ZZ – osoba zdrowa, ++ – wynik dodatni. Z treści: P(C)=0,002P(C) = 0{,}002, więc P(Z)=0,998P(Z) = 0{,}998, oraz P(+∣C)=0,99P(+ \mid C) = 0{,}99. Uwaga na drugi warunek: zdrowi mają wynik ujemny w 98%98\% przypadków, więc dodatni w 2%2\%, czyli P(+∣Z)=0,02P(+ \mid Z) = 0{,}02. Szukamy P(C∣+)P(C \mid +):

P(C∣+)=0,99⋅0,0020,99⋅0,002+0,02⋅0,998=0,001980,02194=991097≈0,09P(C \mid +) = \frac{0{,}99 \cdot 0{,}002}{0{,}99 \cdot 0{,}002 + 0{,}02 \cdot 0{,}998} = \frac{0{,}00198}{0{,}02194} = \frac{99}{1097} \approx 0{,}09

Mimo że test „trafia" w 99%99\% chorych, dodatni wynik oznacza chorobę tylko w około 9%9\% przypadków. Choroba jest rzadka, więc zdrowi z fałszywym alarmem (1 996 osób na 100 000) mocno przeważają nad chorymi z prawdziwym (198 osób) – widać to w tabeli wyżej.

Przykład 2 – trzy przyczyny. Chleb w supermarkecie pochodzi z piekarni P1,P2,P3P_1, P_2, P_3 w stosunku 2:3:42 : 3 : 4. Chleb razowy stanowi 10%10\% dostaw z P1P_1, 12%12\% z P2P_2 i 15%15\% z P3P_3. Oblicz prawdopodobieństwo, że losowo wybrany bochenek pochodzi z P1P_1, jeśli jest razowy.

Stosunek 2:3:42 : 3 : 4 daje razem 99 części, więc P(A1)=29P(A_1) = \frac{2}{9}, P(A2)=39P(A_2) = \frac{3}{9}, P(A3)=49P(A_3) = \frac{4}{9} (suma równa 11 – kontrola). Niech BB oznacza „chleb jest razowy". Mianownik, czyli P(B)P(B):

P(B)=29⋅0,10+39⋅0,12+49⋅0,15=20+36+60900=116900P(B) = \frac{2}{9} \cdot 0{,}10 + \frac{3}{9} \cdot 0{,}12 + \frac{4}{9} \cdot 0{,}15 = \frac{20 + 36 + 60}{900} = \frac{116}{900}

Licznik to ścieżka przez P1P_1: 29⋅0,10=20900\frac{2}{9} \cdot 0{,}10 = \frac{20}{900}. Stąd

P(A1∣B)=20/900116/900=20116=529P(A_1 \mid B) = \frac{20/900}{116/900} = \frac{20}{116} = \frac{5}{29}

Przy okazji: pozostałe dwie piekarnie dają odpowiednio 929\frac{9}{29} i 1529\frac{15}{29}. Razem 5+9+1529=1\frac{5 + 9 + 15}{29} = 1, co jest dobrą kontrolą – skoro bochenek jest razowy, musi skądś pochodzić.

Przykład 3 – dwie urny i moneta. Rzucamy symetryczną monetą. Jeśli wypadnie orzeł, losujemy kulę z urny I (3 białe, 2 czarne), jeśli reszka – z urny II (1 biała, 4 czarne). Wylosowano kulę białą. Z jakim prawdopodobieństwem pochodzi z urny I?

Ścieżka przez urnę I: 12⋅35=310\frac{1}{2} \cdot \frac{3}{5} = \frac{3}{10}. Ścieżka przez urnę II: 12⋅15=110\frac{1}{2} \cdot \frac{1}{5} = \frac{1}{10}. Białą kulę dają obie ścieżki, razem 410\frac{4}{10}, więc

P(urna I∣biała)=3/104/10=34P(\text{urna I} \mid \text{biała}) = \frac{3/10}{4/10} = \frac{3}{4}

Wynik jest większy niż 12\frac{1}{2}, bo urna I ma więcej białych kul – zgodnie z intuicją, że biała kula „wskazuje" na urnę I.

05

Typowe błędy

Odwrócony warunek. W przykładzie z testem uczeń podaje P(+∣C)=0,99P(+ \mid C) = 0{,}99 jako odpowiedź na pytanie „czy pan X jest chory". To prawdopodobieństwo wyniku przy założeniu choroby, a pytanie dotyczy choroby przy założeniu wyniku. Poprawna odpowiedź to około 0,090{,}09. Przed obliczeniami zapisz, co jest szukane, i sprawdź, po której stronie kreski stoi warunek.

Mianownik bez sumy. Zapis P(C∣+)=0,99⋅0,0020,99⋅0,002P(C \mid +) = \frac{0{,}99 \cdot 0{,}002}{0{,}99 \cdot 0{,}002} daje 11 – wynik niemożliwy w zadaniu, w którym ktoś ma fałszywy alarm. W mianowniku jest suma po wszystkich przyczynach (tu chorzy i zdrowi), a nie tylko składnik z licznika. Kontrola: wynik musi mieścić się między 00 a 11 i nie może być równy 11, gdy istnieje druga przyczyna.

Pominięcie prawdopodobieństwa wstępnego. W zadaniu o chlebie ktoś liczy tylko udziały razowego w piekarniach: 1010+12+15=1037≈0,27\frac{10}{10 + 12 + 15} = \frac{10}{37} \approx 0{,}27. To ignoruje, że z P1P_1 pochodzi tylko 29\frac{2}{9} chleba. Poprawnie wychodzi 529≈0,17\frac{5}{29} \approx 0{,}17. Prawdopodobieństwo wstępne P(Bk)P(B_k) zawsze mnoży odpowiedni składnik.

Zła interpretacja „specyficzności". Zdanie „wynik ujemny u 98% zdrowych" to P(−∣Z)=0,98P(- \mid Z) = 0{,}98, a w obliczeniach potrzeba P(+∣Z)=0,02P(+ \mid Z) = 0{,}02. Wstawienie 0,980{,}98 do mianownika zawyża go prawie pięćdziesięciokrotnie. W informatorze CKE osobny punkt przyznano właśnie za obliczenie prawdopodobieństw P(−∣C)P(- \mid C) i P(+∣Z)P(+ \mid Z) – to znak, że ten krok się liczy.

Przyczyny, które nie sumują się do jedynki. Stosunek 2:3:42 : 3 : 4 to udziały 29,39,49\frac{2}{9}, \frac{3}{9}, \frac{4}{9}, a nie 24,34,44\frac{2}{4}, \frac{3}{4}, \frac{4}{4}. Zawsze sprawdź, czy P(B1)+…+P(Bn)=1P(B_1) + \ldots + P(B_n) = 1, zanim użyjesz wzoru.

06

Sprawdź się

W fabryce maszyna M1M_1 wytwarza 60%60\% produkcji, z czego 3%3\% to braki, a maszyna M2M_2 wytwarza 40%40\% produkcji, z czego 5%5\% to braki. Uzupełnij prawdopodobieństwa (jako ułamki dziesiętne).

ZdarzeniePrawdopodobieństwo
wyrób z M1M_1 i brak
wyrób z M2M_2 i brak
wyrób jest brakiem (razem)

Losowo wybrany wyrób okazał się brakiem. Prawdopodobieństwo, że pochodzi z M1M_1, wynosi (do części setnych) , a że pochodzi z M2M_2 – .

W teście z informatora na 100 000100\,000 osób wynik dodatni dostaje zdrowych osób. Prawdopodobieństwo, że osoba z wynikiem dodatnim jest chora, wynosi w zaokrągleniu do części setnych .

Gdyby chorowało 20%20\% populacji, a test miał te same parametry, to prawdopodobieństwo choroby przy wyniku dodatnim wynosiłoby (do części setnych) .

Przyczyny B1,…,BnB_1, \ldots, B_n we wzorze Bayesa muszą wykluczać się nawzajem, czyli być parami , a razem mają prawdopodobieństwo równe .

07

Najczęściej zadawane pytania

Czym wzór Bayesa różni się od prawdopodobieństwa warunkowego?

Prawdopodobieństwo warunkowe P(A∣B)=P(A∩B)P(B)P(A \mid B) = \frac{P(A \cap B)}{P(B)} odpowiada na pytanie „jak często zachodzi AA w sytuacji, gdy zaszło BB". Wzór Bayesa jest jego zastosowaniem do odwrócenia warunku: znając P(A∣Bk)P(A \mid B_k) i P(Bk)P(B_k), liczy P(Bk∣A)P(B_k \mid A). W mianowniku używa twierdzenia o prawdopodobieństwie całkowitym.

Skąd w mianowniku suma?

Mianownik to prawdopodobieństwo zdarzenia AA („wynik dodatni", „bochenek razowy") niezależnie od przyczyny. Rozbijamy je na przyczyny rozłączne: P(A)=∑P(A∣Bi)⋅P(Bi)P(A) = \sum P(A \mid B_i) \cdot P(B_i). Gdyby brać tylko jedną przyczynę, wynik byłby zawsze równy 11.

Dlaczego test, który trafia w 99% chorych, daje tylko 9% pewności?

Bo choroba jest rzadka. Na 100 000100\,000 osób jest 200200 chorych (z nich 198198 ma wynik dodatni) i 99 80099\,800 zdrowych, z których 2%2\% czyli 1 9961\,996 dostaje fałszywy alarm. Wśród 2 1942\,194 osób z wynikiem dodatnim chorych jest 198198, czyli około 9%9\%. Im rzadsze zjawisko, tym więcej fałszywych alarmów w stosunku do prawdziwych.

Czy muszę rysować drzewo?

Nie, ale drzewo albo tabela z liczebnościami porządkuje dane i pozwala uniknąć pomyłki w warunkach. Na drzewie licznik to jedna gałąź (iloczyn prawdopodobieństw wzdłuż niej), a mianownik to suma wszystkich gałęzi kończących się w AA. W rozwiązaniu pisanym musisz jednak zapisać prawdopodobieństwa z treści i wzór, bo z nich bierze się większość punktów.

Kiedy użyć wzoru Bayesa, a kiedy twierdzenia o prawdopodobieństwie całkowitym?

Gdy w zadaniu liczysz zwykłe ułamki z liczby wyników sprzyjających (hasło prawdopodobieństwo klasyczne), Bayes nie jest potrzebny. Gdy pytanie brzmi „jakie jest prawdopodobieństwo skutku" (wynik dodatni, wylosowana biała kula), wystarczy prawdopodobieństwo całkowite. Gdy pytanie brzmi „jakie jest prawdopodobieństwo przyczyny, jeśli znamy skutek" (osoba jest chora, jeśli test dodatni), potrzebny jest Bayes.

Czy wzór Bayesa jest wymagany na poziomie podstawowym?

Nie. Pojawia się w wymaganiach rozszerzonych i jest tam razem z twierdzeniem o prawdopodobieństwie całkowitym. Podstawowe zadania z „jeśli wiadomo, że" rozwiązuje się zwykle liczeniem elementów zbioru, bez tego wzoru.

↗

Powiązane zagadnienia