Satz von Bayes

Aus testwiki
Zur Navigation springen Zur Suche springen
Illustration des Satzes von Bayes durch Überlagerung der beiden ihm zugrundeliegenden Entscheidungsbäume bzw. Baumdiagramme

Der Satz von Bayes (IPA: [[[:Vorlage:IPA]]], Vorlage:Audio) ist ein mathematischer Satz aus der Wahrscheinlichkeitstheorie, der die Berechnung bedingter Wahrscheinlichkeiten beschreibt. Er ist nach dem englischen Mathematiker Thomas Bayes benannt, der ihn erstmals in einem Spezialfall in der 1763 posthum veröffentlichten Abhandlung An Essay Towards Solving a Problem in the Doctrine of Chances beschrieb. Er wird auch Formel von Bayes oder als Lehnübersetzung Bayes-Theorem genannt.

Formel

Für zwei Ereignisse A und B mit P(B)>0 lässt sich die Wahrscheinlichkeit von A unter der Bedingung, dass B eingetreten ist, durch die Wahrscheinlichkeit von B unter der Bedingung, dass A eingetreten ist, errechnen:

P(AB)=P(BA)P(A)P(B).

Hierbei ist

P(AB) die (bedingte) Wahrscheinlichkeit des Ereignisses A unter der Bedingung, dass B eingetreten ist,
P(BA) die (bedingte) Wahrscheinlichkeit des Ereignisses B unter der Bedingung, dass A eingetreten ist,
P(A) die A-priori-Wahrscheinlichkeit des Ereignisses A und
P(B) die A-priori-Wahrscheinlichkeit des Ereignisses B.

Bei endlich vielen Ereignissen lautet der Satz von Bayes:

Wenn Ai,i=1,,N eine Zerlegung der Ergebnismenge in disjunkte Ereignisse ist, gilt für die A-posteriori-Wahrscheinlichkeit P(AiB)

P(AiB)=P(BAi)P(Ai)P(B)=P(BAi)P(Ai)j=1NP(BAj)P(Aj).

Den letzten Umformungsschritt bezeichnet man auch als Marginalisierung.

Da ein Ereignis A und sein Komplement Ac stets eine Zerlegung der Ergebnismenge darstellen, gilt insbesondere

P(AB)=P(BA)P(A)P(BA)P(A)+P(BAc)P(Ac).

Des Weiteren gilt der Satz auch für eine Zerlegung des Grundraumes Ω in abzählbar viele paarweise disjunkte Ereignisse.

Beweis

Der Wahrscheinlichkeitsbaum illustriert P(AB)P(B)=P(AB).

Der Satz folgt unmittelbar aus der Definition der bedingten Wahrscheinlichkeit:

P(AB)=P(AB)P(B)=P(AB)P(A)P(A)P(B)=P(BA)P(A)P(B).

Die Beziehung

P(B)=j=1NP(AjB)=j=1NP(BAj)P(Aj)

ist eine Anwendung des Gesetzes der totalen Wahrscheinlichkeit.

Interpretation

Vorlage:Siehe auch Der Satz von Bayes erlaubt in gewissem Sinn das Umkehren von Schlussfolgerungen: Man geht von einem bekannten Wert P(BA) aus, ist aber eigentlich an dem Wert P(AB) interessiert. Beispielsweise ist es von Interesse, wie groß die Wahrscheinlichkeit ist, dass jemand eine bestimmte Krankheit hat, wenn ein dafür entwickelter Schnelltest ein positives Ergebnis zeigt. Aus empirischen Studien kennt man in der Regel die Wahrscheinlichkeit dafür, mit der der Test bei einer von dieser Krankheit befallenen Person zu einem positiven Ergebnis führt. Die gewünschte Umrechnung ist nur dann möglich, wenn man die Prävalenz der Krankheit kennt, das heißt die (absolute) Wahrscheinlichkeit, mit der die betreffende Krankheit in der Gesamtpopulation auftritt (siehe Rechenbeispiel 2).

Für das Verständnis kann ein Entscheidungsbaum oder eine Vierfeldertafel helfen. Das Verfahren ist auch als Rückwärtsinduktion bekannt.

Mitunter begegnet man dem Fehlschluss, direkt von P(BA) auf P(AB) schließen zu wollen, ohne die A-priori-Wahrscheinlichkeit P(A) zu berücksichtigen, beispielsweise indem angenommen wird, die beiden bedingten Wahrscheinlichkeiten müssten ungefähr gleich groß sein (siehe Prävalenzfehler). Wie der Satz von Bayes zeigt, ist das aber nur dann der Fall, wenn auch P(A) und P(B) ungefähr gleich groß sind.

Ebenso ist zu beachten, dass bedingte Wahrscheinlichkeiten für sich allein nicht dazu geeignet sind, eine bestimmte Kausalbeziehung nachzuweisen.

Anwendungsgebiete

Rechenbeispiel 1

Urnenversuch

In den beiden Urnen A und B befinden sich jeweils zehn Kugeln. In A sind sieben rote und drei weiße Kugeln, in B eine rote und neun weiße. Es wird nun eine beliebige Kugel aus einer zufällig gewählten Urne gezogen. Anders ausgedrückt: Ob aus Urne A oder B gezogen wird, ist a priori gleich wahrscheinlich. Das Ergebnis der Ziehung ist: Die Kugel ist rot. Gesucht ist die Wahrscheinlichkeit, dass diese rote Kugel aus Urne A stammt.

Es sei:

  • A das Ereignis „Die Kugel stammt aus Urne A“,
  • B das Ereignis „Die Kugel stammt aus Urne B“ und
  • R das Ereignis „Die Kugel ist rot“.

Dann gilt: P(A)=P(B)=12  (beide Urnen sind a priori gleich wahrscheinlich)

P(RA)=710  (in Urne A sind 10 Kugeln, davon 7 rote)

P(RB)=110  (in Urne B sind 10 Kugeln, davon 1 rote)

P(R)=P(RA)P(A)+P(RB)P(B)=71012+11012=25  (totale Wahrscheinlichkeit, eine rote Kugel zu ziehen)

Damit ist P(AR)=P(RA)P(A)P(R)=7101225=78 .

Die bedingte Wahrscheinlichkeit, dass die gezogene rote Kugel aus der Urne A gezogen wurde, beträgt also 78=87,5%.

Das Ergebnis der Bayes-Formel in diesem einfachen Beispiel kann leicht anschaulich eingesehen werden: Da beide Urnen a priori mit der gleichen Wahrscheinlichkeit ausgewählt werden und sich in beiden Urnen gleich viele Kugeln befinden, haben alle Kugeln – und damit auch alle acht roten Kugeln – die gleiche Wahrscheinlichkeit, gezogen zu werden. Wenn man wiederholt eine Kugel aus einer zufälligen Urne zieht und wieder in dieselbe Urne zurücklegt, wird man im Durchschnitt in acht von 20 Fällen eine rote und in zwölf von 20 Fällen eine weiße Kugel ziehen (deshalb ist auch die totale Wahrscheinlichkeit, eine rote Kugel zu ziehen, gleich 820=25). Von diesen acht roten Kugeln kommen im Mittel sieben aus Urne A und eine aus Urne B. Die Wahrscheinlichkeit, dass eine gezogene rote Kugel aus Urne A stammt, ist daher gleich 78.

Rechenbeispiel 2

Eine bestimmte Krankheit tritt mit einer Prävalenz von 20 pro 100 000 Personen auf. Der Sachverhalt K, dass ein Mensch diese Krankheit in sich trägt, hat also die Wahrscheinlichkeit P(K)=0,0002.

Ist ein Screening der Gesamtbevölkerung ohne Rücksicht auf Risikofaktoren oder Symptome geeignet, Träger dieser Krankheit zu ermitteln? Es würden dabei weit überwiegend Personen aus dem Komplement Kc von K getestet, also Personen, die diese Krankheit nicht in sich tragen: Die Wahrscheinlichkeit, dass eine zu testende Person nicht Träger der Krankheit ist, beträgt P(Kc)=1P(K)=0,9998.

T bezeichne die Tatsache, dass der Test bei einer Person „positiv“ ausgefallen ist, also die Krankheit anzeigt. Es sei bekannt, dass der Test K mit 95 % Wahrscheinlichkeit anzeigt (Sensitivität =P(TK)=0,95), aber manchmal auch bei Gesunden anspricht, d. h. ein falsch positives Testergebnis liefert, und zwar mit einer Wahrscheinlichkeit von P(TKc)=1% (Spezifität =1P(TKc)=10,01=0,99).

Nicht nur für die Eingangsfrage, sondern in jedem Einzelfall T, insbesondere vor dem Ergebnis weiterer Untersuchungen, interessiert die positiver prädiktiver Wert genannte bedingte Wahrscheinlichkeit P(KT), dass positiv Getestete tatsächlich Träger der Krankheit sind.

Berechnung mit dem Satz von Bayes

P(KT)=P(TK)P(K)P(TK)P(K)+P(TKc)P(Kc)=0,950,00020,950,0002+0,010,99980,0186.

Berechnung mittels Baumdiagramm

Ereignisbaum zum Beispiel

Probleme mit wenigen Klassen und einfachen Verteilungen lassen sich übersichtlich im Baumdiagramm für die Aufteilung der Häufigkeiten darstellen. Geht man von den Häufigkeiten über auf relative Häufigkeiten bzw. auf (bedingte) Wahrscheinlichkeiten, wird aus dem Baumdiagramm ein Ereignisbaum, ein Sonderfall des Entscheidungsbaums.

Den obigen Angaben folgend ergeben sich als absolute Häufigkeit bei 100 000 Personen 20 tatsächlich erkrankte Personen, 99 980 Personen sind gesund. Der Test diagnostiziert bei den 20 kranken Personen in 19 Fällen (95 Prozent Sensitivität) korrekt die Erkrankung; aber in einem Fall versagt der Test und zeigt die vorliegende Krankheit nicht an (falsch negativ). Bei 99 Prozent der 99 980 gesunden Personen (99 Prozent Spezifität) diagnostiziert der Test korrekt; aber bei 1 Prozent, also etwa 1000 der 99 980 gesunden Personen zeigt der Test fälschlicherweise eine Erkrankung an. Von den insgesamt etwa 1019 positiv getesteten Personen sind also nur 19 tatsächlich krank (denn P(KT)1910190,0186).

Bedeutung des Ergebnisses

Der Preis, 19 Träger der Krankheit zu finden, möglicherweise rechtzeitig genug für eine Behandlung oder Isolation, besteht nicht nur in den Kosten für 100 000 Tests, sondern auch in den unnötigen Ängsten und womöglich Behandlungen von 1000 falsch positiv Getesteten. Die Ausgangsfrage, ob bei diesen Zahlenwerten ein Massenscreening sinnvoll ist, ist daher wohl zu verneinen.

Die intuitive Annahme, dass eine – auf den ersten Blick beeindruckende – Sensitivität von 95 % bedeutet, dass eine positiv getestete Person auch tatsächlich mit hoher Wahrscheinlichkeit krank ist, ist also falsch. Dieses Problem tritt immer dann auf, wenn die tatsächliche Rate, mit der ein Merkmal in der untersuchten Gesamtmenge vorkommt (im Beispiel 0,0002), klein ist gegenüber der Rate der falsch positiven Ergebnisse (im Beispiel 0,01).

Ohne Training in der Interpretation statistischer Aussagen werden Risiken oft falsch eingeschätzt oder vermittelt. Der Psychologe Gerd Gigerenzer spricht von Zahlenanalphabetismus im Umgang mit Unsicherheit und plädiert für eine breit angelegte didaktische Offensive.[1]

Bayessche Statistik

Vorlage:Hauptartikel

Die Bayessche Statistik verwendet den Satz von Bayes im Rahmen der induktiven Statistik zur Schätzung von Parametern und zum Testen von Hypothesen.

Problemstellung

Folgende Situation sei gegeben: ϑ ist ein unbekannter Umweltzustand (z. B. ein Parameter einer Wahrscheinlichkeitsverteilung), der auf der Basis einer Beobachtung x einer Zufallsvariable X geschätzt werden soll. Weiterhin ist Vorwissen in Form einer A-priori-Wahrscheinlichkeitsverteilung des unbekannten Parameters ϑ gegeben. Diese A-priori-Verteilung enthält die gesamte Information über den Umweltzustand ϑ, die vor der Beobachtung der Stichprobe gegeben ist.

Je nach Kontext und philosophischer Schule wird die A-priori-Verteilung verstanden

  • als mathematische Modellierung des subjektiven degrees of belief (subjektiver Wahrscheinlichkeitsbegriff),
  • als adäquate Darstellung des allgemeinen Vorwissens (wobei Wahrscheinlichkeiten als natürliche Erweiterung der aristotelischen Logik in Bezug auf Unsicherheit verstanden werden – Cox’ Postulate),
  • als aus Voruntersuchungen bekannte Wahrscheinlichkeitsverteilung eines tatsächlich zufälligen Parameters oder
  • als eine spezifisch gewählte Verteilung, die auf ideale Weise mit Unwissen über den Parameter korrespondiert (objektive A-priori-Verteilungen, zum Beispiel mithilfe der Maximum-Entropie-Methode).

Die bedingte Verteilung von X unter der Bedingung, dass ϑ den Wert ϑ0 annimmt, wird im Folgenden mit f(xϑ0) bezeichnet. Diese Wahrscheinlichkeitsverteilung kann nach Beobachtung der Stichprobe bestimmt werden und wird auch als Likelihood des Parameterwerts ϑ0 bezeichnet.

Die A-posteriori-Wahrscheinlichkeit P(ϑ=ϑ0x) kann mit Hilfe des Satzes von Bayes berechnet werden. Im Spezialfall einer diskreten A-priori-Verteilung erhält man:[2]

P(ϑ=ϑ0x)=f(xϑ0)P(ϑ=ϑ0)ϑΘf(xϑ)P(ϑ=ϑ)

Falls die Menge aller möglichen Umweltzustände endlich ist, lässt sich die A-posteriori-Verteilung im Wert ϑ0 als die Wahrscheinlichkeit interpretieren, mit der man nach Beobachtung der Stichprobe und unter Einbeziehung des Vorwissens den Umweltzustand ϑ0 erwartet.

Als Schätzwert verwendet ein Anhänger der subjektivistischen Schule der Statistik in der Regel den Erwartungswert der A-posteriori-Verteilung, in manchen Fällen auch den Modalwert.

Beispiel

Ähnlich wie oben werde wieder eine Urne betrachtet, die mit zehn Kugeln gefüllt ist, aber nun sei unbekannt, wie viele davon rot sind. Die Anzahl ϑ der roten Kugeln ist hier der unbekannte Umweltzustand und als dessen A-priori-Verteilung soll angenommen werden, dass alle möglichen Werte von null bis zehn gleich wahrscheinlich sein sollen, d. h., es gilt P(ϑ=ϑ0)=111 für alle ϑ0{0,1,,10}.

Nun werde fünfmal mit Zurücklegen eine Kugel aus der Urne gezogen und X bezeichne die Zufallsvariable, die angibt, wie viele davon rot sind. Unter der Annahme ϑ=ϑ0 ist dann X binomialverteilt mit den Parametern n=5 und p=ϑ010, es gilt also

f(xϑ0)=(5x)(ϑ010)x(1ϑ010)5x

für x{0,1,,5}.

Beispielsweise für x=2, d. h., zwei der fünf gezogenen Kugeln waren rot, ergeben sich die folgenden Werte (auf drei Nachkommastellen gerundet)

ϑ0 0 1 2 3 4 5 6 7 8 9 10
P(ϑ=ϑ0) 0,091 0,091 0,091 0,091 0,091 0,091 0,091 0,091 0,091 0,091 0,091
P(ϑ=ϑ0X=2) 0,000 0,044 0,123 0,185 0,207 0,188 0,138 0,079 0,031 0,005 0,000

Man sieht, dass im Gegensatz zur A-priori-Verteilung in der zweiten Zeile, in der alle Werte von ϑ als gleich wahrscheinlich angenommen wurden, unter der A-posteriori-Verteilung in der dritten Zeile ϑ=4 die größte Wahrscheinlichkeit besitzt, das heißt, der A-posteriori-Modus ist 4.

Als Erwartungswert der A-posteriori-Verteilung ergibt sich hier:

ϑ0=010ϑ0P(ϑ=ϑ0X=2)4,287.

Erweiterung

Um bedingte Wahrscheinlichkeitsdichten zu definieren, kann auf die Reguläre bedingte Verteilung zurückgegriffen werden. Dann ist die reguläre bedingte Verteilung von Y gegeben X gegeben durch die Dichte

fY|X(y|x):=f(x,y)fX(x).

Siehe auch

Literatur

  • Alan F. Chalmers: Wege der Wissenschaft: Einführung in die Wissenschaftstheorie. 6. Auflage. Springer, Berlin [u. a.], 2007, ISBN 3-540-49490-1, S. 141–154, Vorlage:DOI (Einführung in wissenschaftsgeschichtlicher Perspektive).
  • Sharon Bertsch McGrayne: Die Theorie, die nicht sterben wollte. Wie der englische Pastor Thomas Bayes eine Regel entdeckte, die nach 150 Jahren voller Kontroversen heute aus Wissenschaft, Technik und Gesellschaft nicht mehr wegzudenken ist. Springer, Berlin/Heidelberg 2014, ISBN 978-3-642-37769-3, doi:10.1007/978-3-642-37770-9
  • F. Thomas Bruss: 250 years of ’An Essay towards solving a Problem in the Doctrine of Chance. By the late Rev. Mr. Bayes, communicated by Mr. Price, in a letter to John Canton, A. M. F. R. S.‘ . In: Jahresbericht der Deutschen Mathematiker-Vereinigung. Vol. 115, Issue 3–4, 2013, S. 129–133, doi:10.1365/s13291-013-0069-z.
  • Wolfgang Tschirk: Statistik: Klassisch oder Bayes. Zwei Wege im Vergleich. Springer Spektrum, 2014, ISBN 978-3-642-54384-5, doi:10.1007/978-3-642-54385-2.

Vorlage:Wikibooks Vorlage:Wikibooks Vorlage:Wikibooks

Einzelnachweise

  1. Gerd Gigerenzer: Das Einmaleins der Skepsis. Piper, Berlin 2014, ISBN 978-3-8270-7792-9 (Rezension des englischen Originals. In: NEJM).
  2. Vorlage:Literatur

Vorlage:Normdaten