Dirichlet-Verteilung

Beispiele einer Dirichlet-Verteilung mit K=3 für verschiedene Parametervektoren α. Im Uhrzeigersinn von oben links: α=(6, 2, 2), (3, 7, 5), (6, 2, 6), (2, 3, 4).

Die Dirichletverteilung (nach Peter Gustav Lejeune Dirichlet) ist eine Familie von stetigen, multivariaten Wahrscheinlichkeitsverteilungen.

Sie ist die multivariate Erweiterung der Beta-Verteilung und die konjugierte A-priori-Verteilung der Multinomialverteilung in der bayesschen Statistik. Ihre Dichtefunktion beschreibt die Verteilung der Wahrscheinlichkeiten für K verschiedene, exklusive Ereignisse. Sie wird durch einen Parametervektor $\alpha =({\alpha }_1,{\alpha }_2,...,{\alpha }_K)$ gesteuert, wobei jeder Parameter ${\alpha }_i$ das Vorwissen über die Häufigkeit des $i$-ten Ereignisses widerspiegelt. Konkret entspricht der Parameter ${\alpha }_i-1$ der Anzahl der angenommenen „Beobachtungen“ oder „Erfolge“ für das $i$-te Ereignis. Höhere Werte von ${\alpha }_i$ deuten auf eine größere Zuversicht in die Wahrscheinlichkeit des entsprechenden Ereignisses hin.

Veranschaulichung

Die Multinomialverteilung gibt die Wahrscheinlichkeiten $p_1$ bis $p_K$ für K unterschiedliche Ereignisse an, also z. B. wie wahrscheinlich es ist, in einem Wurf eine Eins, Zwei, Drei, Vier, Fünf oder Sechs zu würfeln. Im Gegensatz dazu gibt die Dirichlet-Verteilung an, wie wahrscheinlich eine solche Verteilung auftritt. Im Falle einer Würfelfabrik könnte die Dirichlet-Verteilung also angeben, wie wahrscheinlich die Verteilungen der Würfelergebnisse bei den fabrizierten Würfeln sind. Funktionieren die Maschinen der Würfelfabrik korrekt, wäre die Wahrscheinlichkeit für alles andere als die uniforme Verteilung (alle Augenzahlen sind gleich wahrscheinlich) sehr gering. Das entspräche einem Parametervektor $\alpha$ mit gleichen und sehr hohen Elementen wie etwa $(1000,1000,1000,1000,1000,1000)$. Hingegen würde $\alpha =(1000,500,500,500,500,500)$ bedeuten, dass die Maschinen Würfel fabrizieren, bei denen die Augenzahl Eins doppelt so häufig vorkommt wie jede andere Augenzahl. Und dies fast ausnahmslos, da die Werte wiederum sehr hoch sind und damit die Varianz niedrig. Wären die Werte in $\alpha$ aber z. B. alle $0,1$, dann würden Würfel hergestellt werden, die eine starke Tendenz zu einer Augenzahl haben. Welche die bevorzugte Augenzahl eines Würfels ist, wäre dabei zufällig, da alle Werte in $\alpha$ gleich sind. Je kleiner die Werte, desto ausgeprägter wäre die Unfairness der meisten Würfel, und desto seltener wären Würfel ohne eine bevorzugte Augenzahl.

Dichtefunktion

Die Dirichletverteilung der Ordnung K ≥ 2 mit den Parametern ${\alpha }_1,...,{\alpha }_K>0$ und dem Parametervektor $\alpha =({\alpha }_1,\dots ,{\alpha }_K)$ hat folgende Dichtefunktion:

$f(x_1,\dots ,x_K;{\alpha }_1,\dots ,{\alpha }_K)=\{\begin{array}{cc}\frac{1}{\mathrm{B}(\alpha )}\prod _{i=1}^K{x}_i^{{\alpha }_i-1}& \text{für alle }{x}_1>0,\dots ,x_K>0\text{ mit }\sum _{i=1}^K{x}_i=1\\ 0& \text{sonst}\end{array}.$

Die normierende Konstante $\mathrm{B}(\alpha )$ ist die multivariate Betafunktion an der Stelle $\alpha$, welche durch Werte der Gammafunktion an den Stellen ${\alpha }_1,\dots ,{\alpha }_K$ dargestellt werden kann:

$\mathrm{B}(\alpha )=\frac{\prod _{i=1}^K\mathrm{\Gamma }({\alpha }_i)}{\mathrm{\Gamma }\left(\sum _{i=1}^K{\alpha }_i\right)}.$

Die Dichtefunktion ist keine Dichte bezüglich des K-dimensionalen Lebesgue-Maßes, sondern eine Dichtefunktion bezüglich des (K-1)-dimensionalen Lebesgue-Maßes im durch die Restriktion ${\sum }_{i=1}^K{x}_i=1$ definierten (K-1)-dimensionalen Teilraum. Durch die Ersetzung $x_K=1-{\sum }_{i=1}^{K-1}{x}_i$ erhält man die (K-1)-dimensionale Dichtefunktion

$f(x_1,\dots ,x_{K-1};{\alpha }_1,\dots ,{\alpha }_K)=\{\begin{array}{cc}\frac{1}{\mathrm{B}(\alpha )}\prod _{i=1}^{K-1}{x}_i^{{\alpha }_i-1}(1-{\sum }_{i=1}^{K-1}{x}_i{)}^{{\alpha }_K-1}& \text{für alle }{x}_1>0,\dots ,x_{K-1}>0\text{ mit }\sum _{i=1}^{K-1}{x}_i<1\\ 0& \text{sonst}\end{array}.$

Bei einer Verwendung der Dirichlet-Verteilung als A-priori-Verteilung für eine Multinomialverteilung sind die Vektoren $(x_1,\dots ,x_K)$ mit positiver Dichte alternative Werte für den Parametervektor einer Multinomialverteilung.

Weblinks

• Eintrag in der Encyclopedia of Mathematics (Springer)