Die Maximum-a-posteriori-Methode (= MAP) ist in der mathematischen Statistik ein Schätzverfahren, genauer ein spezieller Bayes-Schätzer. Das Verfahren schätzt einen unbekannten Parameter durch den Modalwert der A-posteriori-Verteilung. Somit besteht eine gewisse Ähnlichkeit zur Maximum-Likelihood-Methode.

Beschreibung

Folgende Situation ist gegeben: ${\displaystyle \theta }$ ist ein unbekannter Parameter einer Grundgesamtheit, der auf der Basis von Beobachtungen ${\displaystyle x}$ geschätzt werden soll. Weiterhin sei ${\displaystyle f}$ die Stichprobenverteilung von ${\displaystyle x}$, also ${\displaystyle f(x|\theta )}$ die Wahrscheinlichkeit von ${\displaystyle x}$, wenn der (wahre) Parameter der Grundgesamtheit den Wert ${\displaystyle \theta }$ hat.

Die Funktion

${\displaystyle \theta \mapsto f(x|\theta )\!}$

ist als Likelihood-Funktion bekannt, und der Schätzwert

${\displaystyle {\hat {\theta }}_{\mathrm {ML} }(x)={\underset {\theta }{\operatorname {arg\,max} }}\ f(x|\theta )\!}$

als Maximum-Likelihood-Schätzer von ${\displaystyle \theta }$.

Nun stehe aber zusätzlich eine A-priori-Verteilung ${\displaystyle g}$ von ${\displaystyle \theta }$ zur Verfügung. Dadurch lässt sich ${\displaystyle \theta }$ als eine Zufallsvariable betrachten, wie es in der bayesschen Statistik üblich ist. Die A-posteriori-Verteilung von ${\displaystyle \theta }$ erhält man mit Hilfe des Satzes von Bayes durch:

${\displaystyle \theta \mapsto f(\theta |x)={\frac {f(x|\theta )\,g(\theta )}{\displaystyle \int _{\Theta }f(x|\theta ')\,g(\theta ')\,d\theta '}}}$

Die Maximum-a-Posteriori-Methode verwendet jetzt den Modalwert der A-posteriori-Verteilung als Schätzwert für ${\displaystyle \theta }$. Da das Integral im Nenner nicht von ${\displaystyle \theta }$ abhängt, braucht es bei der Maximierung nicht berücksichtigt zu werden:

${\displaystyle {\hat {\theta }}_{\mathrm {MAP} }(x)={\underset {\theta \in \Theta }{\operatorname {arg\,max} }}\left[{\frac {f(x|\theta )\,g(\theta )}{\displaystyle \int _{\Theta }f(x|\theta ')\,g(\theta ')\,d\theta '}}\right]={\underset {\theta \in \Theta }{\operatorname {arg\,max} }}\left[f(x|\theta )\,g(\theta )\right]}$.

Der MAP-Schätzer von ${\displaystyle \theta }$ ist identisch zum Maximum-Likelihood-Schätzer (= ML-Schätzer), wenn eine nichtinformative A-priori-Verteilung (z. B. Gleichverteilung) verwendet wird.

Vergleich mit anderen bayesschen (Punkt-)Schätzmethoden

In der Literatur wird der MAP-Schätzer als das Äquivalent der bayesschen Statistik zum ML-Schätzer beschrieben.

Allerdings spielt die MAP-Schätzung in der bayesschen Statistik nicht die gleiche Rolle wie der ML-Schätzer in der frequentistischen Statistik:

• Bayessche Statistiker drücken in der Regel die (A-posteriori-)Information über einen unbekannten Parameter in einer Wahrscheinlichkeitsverteilung aus, und nicht in einem Punktschätzer.
• Der Erwartungswert der A-posteriori-Verteilung ist dem MAP-Schätzer überlegen, wenn man, wie in der bayesschen Statistik üblich, die A-posteriori-Varianz eines Schätzers als Gütemaß verwendet.
• In vielen Fällen ist zudem auch der Median ein besserer Schätzer als der MAP-Schätzer.

In der bayesschen Entscheidungstheorie sind für die gebräuchlichsten Verlustfunktionen andere Schätzer als der MAP-Schätzer optimal:

• Bei einer quadratischen Verlustfunktion ist der Erwartungswert der A-posteriori-Verteilung die optimale Schätzfunktion.
• Wenn man den Absolutbetrag des Schätzfehlers als Verlustfunktion verwendet (also: ${\displaystyle L(\theta ,a)=|\theta -a|}$, mit a als beliebiger Schätzfunktion), ist der Median der A-posteriori-Verteilung die optimale Schätzfunktion.

Vergleich mit der Maximum-Likelihood-Methode

Im Unterschied zur ML-Methode wird bei der MAP-Methode Vorwissen in Form von A-priori-Wahrscheinlichkeiten berücksichtigt. Diese A-priori-Wahrscheinlichkeiten ergeben zusammen mit der Stichprobe nach dem Satz von Bayes die A-posteriori-Wahrscheinlichkeit. Die MAP-Methode verwendet den wahrscheinlichsten Parameterwert unter der A-posteriori-Verteilung, während die ML-Methode den Parameter mit der höchsten Likelihood (d. h. ohne Vorwissen) verwendet. Allerdings ist für einen streng frequentistischen Statistiker die Verwendung einer A-priori-Verteilung inakzeptabel. Daher wird in der klassischen Statistik die ML-Methode statt der MAP-Methode verwendet.

Der Maximum-Likelihood-Schätzer kann als Spezialfall eines Maximum-a-Posteriori-Schätzers angesehen werden, bei dem die A-priori-Verteilung ${\displaystyle g\left(\theta \right)}$ uninformativ ist (für beschränkte Wertebereiche von ${\displaystyle \theta }$ etwa eine Gleichverteilung). Umgekehrt ist jeder Maximum-a-Posteriori-Schätzer für Stichprobenverteilung ${\displaystyle f\left(x\mid \theta \right)}$ und A-priori-Verteilung ${\displaystyle g\left(\theta \right)}$ ein Maximum-Likelihood-Schätzer für die Stichprobenverteilung ${\displaystyle h_{\theta }\left(x\right)}$ mit

${\displaystyle h_{\theta }\left(x\right)={\frac {f\left(x\mid \theta \right)\cdot g\left(\theta \right)}{\displaystyle \int _{\Theta }f\left(x\mid \theta '\right)\cdot g\left(\theta '\right)\,\mathrm {d} \theta '}}}$.

Beide Verfahren können sich also gegenseitig simulieren und sind in diesem Sinne „gleichmächtig“.

Beispiel

In einer Urne befinden sich rote und schwarze Kugeln. Durch Ziehen mit Zurücklegen soll ermittelt werden, wie hoch der (wahre) Anteil der roten Kugeln in der Urne ist. Die Anzahl der roten Kugeln kann dann durch eine Binomialverteilung mit Stichprobengröße N=10 und unbekanntem Parameter ${\displaystyle \pi }$ (d. h. einer ${\displaystyle \operatorname {Bin} \left(10,\pi \right)}$-Verteilung) beschrieben werden. Im Folgenden nehmen wir an, dass bei solch einem Test 7 rote Kugeln gezogen wurden.

Siehe auch

• posterior predictive distribution

Literatur

• Bernhard Rüger: Induktive Statistik. Einführung für Wirtschafts- und Sozialwissenschaftler. R. Oldenbourg Verlag, München Wien 1988. ISBN 3-486-20535-8
• James O. Berger: Statistical decision theory and Bayesian analysis. Springer Series in Statistics, Springer-Verlag, New York Berlin Heidelberg 1985. ISBN 0-387-96098-8

Einzelnachweise

1. ↑ Bernhard Rüger: Induktive Statistik. Einführung für Wirtschafts- und Sozialwissenschaftler, S. 161f
2. ↑ James O. Berger: Statistical decision theory and Bayesian analysis, S. 133
3. ↑ James O. Berger: Statistical decision theory and Bayesian analysis, S. 136
4. ↑ James O. Berger: Statistical decision theory and Bayesian analysis, S. 134
5. ↑ James O. Berger: Statistical decision theory and Bayesian analysis, S. 161f.