Poziom zaawansowany

Prawdopodobieństwo warunkowe i twierdzenie Bayesa

Informacja zmienia szanse. Gdy wiesz, że zaszło zdarzenie B, przestrzeń możliwości kurczy się do samego B — i wszystkie prawdopodobieństwa trzeba przeliczyć w środku. Ta lekcja pokazuje, jak to policzyć, jak odwrócić warunek twierdzeniem Bayesa i dlaczego dodatni wynik testu przesiewowego prawie nigdy nie znaczy tego, co się wydaje.

Zanim zaczniesz

Ten temat opiera się na wcześniejszych zagadnieniach. Zanim zaczniesz, warto przerobić poniższe lekcje — dzięki nim wszystko pójdzie gładko:

Gdzie się tego używa

Konkretne sytuacje, w których liczysz dokładnie tak, jak uczy ta lekcja:

  • Wynik testu przesiewowego
    Test wykrywa chorobę u 99 na 100 chorych i myli się u 2 na 100 zdrowych, a choruje 1 osoba na 100. Wynik dodatni dostaje 0,01 · 0,99 = 0,0099 chorych i 0,99 · 0,02 = 0,0198 zdrowych, czyli razem 2,97 procenta badanych. Prawdopodobieństwo choroby po dodatnim wyniku wynosi 0,0099 / 0,0297, czyli dokładnie jedną trzecią — nie 99 procent. Dlatego dodatni wynik badania przesiewowego zawsze powtarza się testem potwierdzającym.
  • Filtr antyspamowy
    Filtr uczy się, jak często słowo pada w spamie, a jak często w zwykłej poczcie. Jeśli słowo „faktura” występuje w 40 procentach spamu i w 2 procentach normalnych listów, a spamem jest 30 procent poczty, to po zobaczeniu tego słowa prawdopodobieństwo spamu wynosi 0,3 · 0,4 / (0,3 · 0,4 + 0,7 · 0,02) = 0,12 / 0,134, czyli około 90 procent. Cały klasyfikator naiwny Bayesa to ten sam rachunek powtórzony dla każdego słowa wiadomości.
  • Kontrola jakości na dwóch liniach
    Zakład produkuje 70 procent sztuk na linii A, gdzie wadliwa jest 1 na 100, i 30 procent na linii B, gdzie wadliwe są 4 na 100. Wadliwość całej produkcji to 0,7 · 0,01 + 0,3 · 0,04 = 0,019, czyli 1,9 procenta. Gdy reklamacja dotyczy sztuki wadliwej, prawdopodobieństwo, że wyszła z linii B, wynosi 0,012 / 0,019, czyli około 63 procent — i to jest liczba, od której zaczyna się audyt.
  • Prognoza pogody na jutro
    Serwis zapowiada deszcz na 30 procent dni. Gdy zapowiada, pada w 80 procentach przypadków; gdy nie zapowiada, pada w 10 procentach. Deszcz w losowo wybrany dzień ma więc prawdopodobieństwo 0,3 · 0,8 + 0,7 · 0,1 = 0,31. Jeśli spadł deszcz, a nie pamiętasz prognozy, szansa, że była zapowiadana, wynosi 0,24 / 0,31, czyli około 77 procent.

Wszystkie wzory

  • Prawdopodobieństwo warunkowe

    P(AB)=P(AB)P(B)P(A|B) = \frac{P(A \cap B)}{P(B)}

    określone tylko dla P(B) > 0

  • Reguła mnożenia

    P(AB)=P(B)P(AB)P(A \cap B) = P(B) \cdot P(A|B)

    definicja przemnożona przez P(B)

  • Wersja licznościowa

    P(AB)=ABBP(A|B) = \frac{|A \cap B|}{|B|}

    gdy zdarzenia elementarne są równo prawdopodobne

  • Prawdopodobieństwo całkowite

    P(B)=i=1nP(Ai)P(BAi)P(B) = \sum_{i=1}^{n} P(A_i) \cdot P(B|A_i)

    A₁, …, Aₙ rozbijają Ω na rozłączne przypadki

  • Twierdzenie Bayesa

    P(AkB)=P(Ak)P(BAk)i=1nP(Ai)P(BAi)P(A_k|B) = \frac{P(A_k) \cdot P(B|A_k)}{\sum_{i=1}^{n} P(A_i) \cdot P(B|A_i)}

    odwraca warunek: z P(B|A) robi P(A|B)

  • Niezależność

    P(AB)=P(A)    P(AB)=P(A)P(B)P(A|B) = P(A) \iff P(A \cap B) = P(A) \cdot P(B)

    warunek nie zmienia szans zdarzenia A

W prawdopodobieństwie klasycznym padło zdanie, że losowanie bez zwracania nie jest niezależne, bo po pierwszym losowaniu skład urny się zmienia. Zapisaliśmy to wtedy symbolem P(AB)P(A|B) i zostawiliśmy bez rachunku. Ta lekcja go dostarcza.

Punkt wyjścia jest prosty: informacja zmienia szanse. Gdy wiesz, że zaszło zdarzenie BB, wyniki spoza BB przestały być możliwe — przestrzeń zdarzeń elementarnych skurczyła się z całego Ω\Omega do samego BB i wszystko trzeba przeliczyć w środku.

Definicja

Prawdopodobieństwem warunkowym zdarzenia AA pod warunkiem BB nazywamy liczbę

P(AB)=P(AB)P(B),P(B)>0P(A|B) = \frac{P(A \cap B)}{P(B)}, \qquad P(B) > 0

W liczniku stoi to, co sprzyja AA i leży w nowej przestrzeni, w mianowniku — cała nowa przestrzeń. To dokładnie ta sama konstrukcja co wzór klasyczny, tylko rolę Ω\Omega przejęło BB.

Zastrzeżenie P(B)>0P(B) > 0 nie jest formalnością: warunkowanie przez zdarzenie niemożliwe nie ma sensu, bo nie ma czym dzielić.

Gdy wszystkie zdarzenia elementarne są jednakowo prawdopodobne, wolno skrócić Ω|\Omega| w liczniku i mianowniku i liczyć wprost na licznościach:

P(AB)=ABBP(A|B) = \frac{|A \cap B|}{|B|}

Tabela dwudzielcza

Najwygodniejszym narzędziem do liczenia prawdopodobieństw warunkowych jest tabela dwudzielcza: jedno kryterium w wierszach, drugie w kolumnach, liczności w komórkach.

Stu uczniów przed egzaminem — połowa się uczyła, połowa nie:

zdałnie zdałrazem
uczył się45550
nie uczył się153550
razem6040100

Każde prawdopodobieństwo warunkowe czyta się z tej tabeli jako iloraz komórki przez sumę wiersza albo kolumny — bez ani jednego wzoru.

Z tabeli powyżej odczytaj P(zdał | uczył się) oraz P(uczył się | zdał).

P(AB)P(A|B) to nie jest P(BA)P(B|A)

Przykład wyżej nie był ciekawostką, tylko najważniejszym ostrzeżeniem całej lekcji. Oba ilorazy dzielą to samo P(AB)P(A \cap B), ale przez co innego, więc równe są tylko w szczególnym przypadku P(A)=P(B)P(A) = P(B).

Codzienne wersje tej pomyłki brzmią zupełnie niewinnie:

  • „Większość chorych na grypę ma gorączkę" nie znaczy „większość osób z gorączką ma grypę".
  • „Prawie każdy wypadek ma w tle nadmierną prędkość" nie znaczy „prawie każda jazda za szybko kończy się wypadkiem".

W obu zdaniach zamieniono warunek z wnioskiem. Narzędzie, które robi tę zamianę poprawnie, poznasz dwie sekcje niżej.

Reguła mnożenia

Definicję warto raz przemnożyć przez mianownik, bo w tej postaci używa się jej najczęściej:

P(AB)=P(B)P(AB)P(A \cap B) = P(B) \cdot P(A|B)

Czyta się to jak przepis na dwa kroki doświadczenia: najpierw zaszło BB, potem, już w zawężonym świecie, zaszło AA.

Z talii 52 kart losujemy dwie karty bez zwracania. Jakie jest prawdopodobieństwo, że obie będą asami?

Prawdopodobieństwo całkowite

Często jest odwrotnie: nie znamy P(B)P(B), za to znamy je osobno w każdym przypadku. Wtedy rozbijamy Ω\Omega na rozłączne zdarzenia A1,,AnA_1, \dots, A_n pokrywające całość i sumujemy:

P(B)=i=1nP(Ai)P(BAi)P(B) = \sum_{i=1}^{n} P(A_i) \cdot P(B|A_i)

Każdy składnik to jedna reguła mnożenia. Na rysunku odpowiada mu jedna ścieżka: prawdopodobieństwa mnoży się wzdłuż niej, a gotowe ścieżki dodaje.

urnaU₁ 1/3biała 2/52/15czarna 3/53/15U₂ 1/3biała 3/53/15czarna 2/52/15U₃ 1/3biała 1/51/15czarna 4/54/15
Trzy urny wybierane po równo, w każdej inny skład. Liczba przy końcu ścieżki to iloczyn prawdopodobieństw na jej gałęziach. Białą kulę dają trzy ścieżki — ich suma jest prawdopodobieństwem całkowitym.
Trzy urny wybieramy z równym prawdopodobieństwem. W pierwszej są 2 kule białe na 5, w drugiej 3 na 5, w trzeciej 1 na 5. Jakie jest prawdopodobieństwo wylosowania kuli białej?

Twierdzenie Bayesa

Teraz odwrócenie warunku. Skoro P(AkB)P(A_k \cap B) da się rozpisać na dwa sposoby, to

P(AkB)=P(AkB)P(B)=P(Ak)P(BAk)i=1nP(Ai)P(BAi)P(A_k|B) = \frac{P(A_k \cap B)}{P(B)} = \frac{P(A_k) \cdot P(B|A_k)}{\sum_{i=1}^{n} P(A_i) \cdot P(B|A_i)}

To jest twierdzenie Bayesa, a w mianowniku stoi dokładnie prawdopodobieństwo całkowite z poprzedniej sekcji. Na drzewie ma ono krótki opis: jedna ścieżka podzielona przez sumę wszystkich ścieżek kończących się tym samym wynikiem.

badanychory 0,01wynik + 0,990,0099wynik − 0,010,0001zdrowy 0,99wynik + 0,020,0198wynik − 0,980,9702
Dwie ścieżki kończą się wynikiem dodatnim: chory wykryty (0,0099) i zdrowy pomylony (0,0198). Druga jest dwa razy grubsza od pierwszej — mimo że test myli się rzadko, zdrowych jest po prostu bardzo dużo.
Test wykrywa chorobę u 99% chorych i daje fałszywy alarm u 2% zdrowych. Choruje 1% populacji. Jakie jest prawdopodobieństwo, że osoba z wynikiem dodatnim jest chora?

Dlaczego tak mało? Bo zdrowych jest 9999 razy więcej niż chorych, więc nawet dwuprocentowy odsetek pomyłek daje dwa razy więcej fałszywych alarmów niż wszystkich trafień. Widać to na drzewie: cienka gałąź 0,00990{,}0099 przeciwko grubej 0,01980{,}0198.

Częstość bazowa i błąd prokuratora

Liczba P(C)=0,01P(C) = 0{,}01 z przykładu ma własną nazwę: częstość bazowa — prawdopodobieństwo hipotezy, zanim zobaczyliśmy jakikolwiek dowód. Wchodzi i do licznika, i do mianownika wzoru Bayesa, więc rządzi wynikiem.

Zobaczmy, co się dzieje, gdy zmienia się tylko ona, a test zostaje ten sam:

częstość bazowaprawdopodobieństwo choroby po wyniku dodatnim
0,1%0{,}1\%4,7%4{,}7\%
1%1\%33%33\%
10%10\%85%85\%
50%50\%98%98\%

Ten sam test, ta sama czułość, cztery zupełnie różne odpowiedzi. Pominięcie częstości bazowej — czyli czytanie czułości testu tak, jakby była odpowiedzią na pytanie o pacjenta — jest właśnie tym błędem, który tabela pokazuje w liczbach.

Jego wersja sądowa nazywa się błędem prokuratora. Brzmi tak: ślad pasuje do jednej osoby na tysiąc, więc oskarżony jest winny z prawdopodobieństwem 99,9%99{,}9\%. To nieprawda, bo pomylono warunek z wnioskiem.

Ślad pasuje do 1 osoby na 1000. Zostawić go mógł każdy z miliona mieszkańców, a innych dowodów brak. Co można powiedzieć o winie osoby dopasowanej?

Niezależność jako przypadek szczególny

Skoro warunkowanie ma zmieniać szanse, warto nazwać sytuację, w której nic nie zmienia:

P(AB)=P(A)P(A|B) = P(A)

Podstawienie tego do definicji daje natychmiast wzór, który znasz z lekcji o prawdopodobieństwie klasycznym:

P(AB)=P(A)P(B)P(A \cap B) = P(A) \cdot P(B)

Niezależność nie jest więc osobnym pojęciem, tylko szczególnym przypadkiem prawdopodobieństwa warunkowego — tym, w którym informacja o BB okazuje się bezużyteczna. W tabeli z uczniami niezależności nie ma: P(zdał)=60100P(\text{zdał}) = \tfrac{60}{100}, a P(zdałuczył się)=4550P(\text{zdał}|\text{uczył się}) = \tfrac{45}{50}, więc uczenie się jednak pomaga.

Zadania

Zestaw pyta o trzy rzeczy, w kolejności lekcji. Prompt P(AB)P(A|B) podaje prawdopodobieństwa P(AB)P(A \cap B) oraz P(B)P(B) — nie liczności — i prosi o iloraz. Prompt P(B)P(B) podaje rozbicie na przypadki: pary P(Ai)P(A_i) i P(BAi)P(B|A_i), z których trzeba złożyć prawdopodobieństwo całkowite. Prompt P(A1B)P(A_1|B) podaje te same dane i pyta o odwrócenie warunku, czyli o twierdzenie Bayesa.

Odpowiedzią w każdym z trzech przypadków jest ułamek nieskracalny — konwencja z lekcji o prawdopodobieństwie klasycznym obowiązuje dalej.

Ćwiczenia

Rozwiąż zestaw zadań — trudność rośnie z każdym kolejnym. Na końcu zobaczysz wynik i listę pomyłek do powtórki.

Zadanie 1 z 8Punkty: 0
P(A|B): P(A∩B) = 1/10, P(B) = 4/5

Generator nie pyta o interpretację wyniku — o to, czy przy danej częstości bazowej dodatni wynik testu jest już powodem do niepokoju. To rozumowanie, a nie liczba, więc masz je wyłożone w przykładach powyżej.

Częste błędy

  • Mylenie P(AB)P(A|B) z P(BA)P(B|A) — ten sam licznik, inne mianowniki. Zawsze sprawdź, co jest warunkiem: to on stoi pod kreską ułamka.
  • Ignorowanie częstości bazowej — czułość testu nie jest odpowiedzią na pytanie o pacjenta. Przy chorobie rzadkiej większość alarmów jest fałszywa, choć test jest dobry.
  • Dzielenie przez P(A)P(A) zamiast przez P(B)P(B) — w mianowniku stoi prawdopodobieństwo warunku, czyli tego zdarzenia, o którym już wiadomo, że zaszło.
  • Rozbicie, które nie pokrywa całego Ω\Omega — wzór na prawdopodobieństwo całkowite wymaga przypadków rozłącznych i wyczerpujących; jeśli suma P(Ai)P(A_i) nie wynosi 11, rachunek jest bez sensu.
  • Zakładanie niezależności bez sprawdzenia — mnożenie P(A)P(B)P(A) \cdot P(B) jest legalne dopiero wtedy, gdy warunek naprawdę nic nie zmienia. Losowanie bez zwracania nigdy takie nie jest.
  • Warunkowanie przez zdarzenie o prawdopodobieństwie zero — nie ma czym dzielić, więc wyrażenie P(AB)P(A|B) po prostu nie jest określone.

Karta wzorów

Temat: Prawdopodobieństwo warunkowe

  • Prawdopodobieństwo warunkowe

    P(AB)=P(AB)P(B)P(A|B) = \frac{P(A \cap B)}{P(B)}

    określone tylko dla P(B) > 0

  • Reguła mnożenia

    P(AB)=P(B)P(AB)P(A \cap B) = P(B) \cdot P(A|B)

    definicja przemnożona przez P(B)

  • Wersja licznościowa

    P(AB)=ABBP(A|B) = \frac{|A \cap B|}{|B|}

    gdy zdarzenia elementarne są równo prawdopodobne

  • Prawdopodobieństwo całkowite

    P(B)=i=1nP(Ai)P(BAi)P(B) = \sum_{i=1}^{n} P(A_i) \cdot P(B|A_i)

    A₁, …, Aₙ rozbijają Ω na rozłączne przypadki

  • Twierdzenie Bayesa

    P(AkB)=P(Ak)P(BAk)i=1nP(Ai)P(BAi)P(A_k|B) = \frac{P(A_k) \cdot P(B|A_k)}{\sum_{i=1}^{n} P(A_i) \cdot P(B|A_i)}

    odwraca warunek: z P(B|A) robi P(A|B)

  • Niezależność

    P(AB)=P(A)    P(AB)=P(A)P(B)P(A|B) = P(A) \iff P(A \cap B) = P(A) \cdot P(B)

    warunek nie zmienia szans zdarzenia A

badanychory 0,01wynik + 0,990,0099wynik − 0,010,0001zdrowy 0,99wynik + 0,020,0198wynik − 0,980,9702
Test przesiewowy w populacji, w której choruje 1 osoba na 100. Cztery ścieżki wyczerpują wszystkie możliwości, a ich prawdopodobieństwa sumują się do jedynki. Wynik dodatni dają dwie z nich — i to ta druga, grubsza, przeważa.

Najczęstsze pytania

Powiązane artykuły