# Multivariate Verteilung

Eine multivariate Verteilung ist in der Wahrscheinlichkeitsrechnung und in der Statistik die Verteilung eines Zufallsvektors – also einer Zufallsvariablen, deren Werte Vektoren im $\mathbb {R} ^{n}$ sind. Im zweidimensionalen Fall $n=2$ spricht man auch von einer bivariaten Verteilung. Die multivariate Verteilung eines Zufallsvektors $X=(X_{1},\dotsc ,X_{n})$ ist somit ein Wahrscheinlichkeitsmaß auf $\mathbb {R} ^{n}$, das messbaren Teilmengen $A\subseteq \mathbb {R} ^{n}$ die Wahrscheinlichkeit zuordnet, dass $X$ einen Wert aus $A$ annimmt. Eine multivariate Verteilung kann durch eine multivariate Verteilungsfunktion charakterisiert werden. Die Verteilungen der einzelnen Komponenten $X_{i}$ eines Zufallsvektors werden die Randverteilungen von $X$ genannt. Beispiele für multivariate Verteilungen sind die Multinomialverteilung oder die multivariate Normalverteilung, weitere finden sich in der Liste multivariater und matrixvariater Wahrscheinlichkeitsverteilungen.

## Einführendes Beispiel

Wir betrachten zwei Zufallsexperimente:
* Zweimaliges Würfeln mit einem idealen Würfel. Dies ist äquivalent zu einem Urnenexperiment mit sechs unterscheidbaren Kugeln, wobei zweimal mit Zurücklegen gezogen wird. Es gibt 36 mögliche Ergebnispaare (da wir die Reihenfolge des Würfelns bzw. der Ziehung berücksichtigen), und alle 36 Möglichkeiten sind gleich wahrscheinlich, haben also eine Wahrscheinlichkeit von 1/36.
* Ein ähnliches Urnenexperiment, aber ohne Zurücklegen. In diesem Fall kommen die Ergebnisse (1,1), (2,2), …, (6,6) nicht vor, da die $i$-te Kugel beim zweiten Ziehen nicht vorkommen kann, wenn sie bereits bei der ersten Ziehung herausgenommen wurde. Die übrigen 30 Paare sind gleich wahrscheinlich und haben daher die Wahrscheinlichkeit 1/30.
Diese beiden Experimente ergeben nun zweidimensionale diskrete Zufallsvariablen $Z_{1}$ und $Z_{2}$, welche die gleichen Randverteilungen haben (jede Zahl von 1 bis 6 ist bei beiden Experimenten in beiden Ziehungen gleich wahrscheinlich und tritt mit Wahrscheinlichkeit 1/6 auf).
Jedoch sind die beiden Ziehungen im ersten Experiment unabhängig, da die gezogene Kugel zurückgelegt wird, während sie im zweiten Experiment nicht unabhängig sind. Das wird am deutlichsten, wenn man sich klarmacht, dass die Paare (1,1), (2,2), …, (6,6) bei einem unabhängigen Experiment jeweils mit Wahrscheinlichkeit 1/36 vorkommen müssen (Produkt der Randwahrscheinlichkeiten 1/6), sie aber beim zweiten Experiment überhaupt nicht auftreten können (Wahrscheinlichkeit 0 haben), da die Kugel nicht zurückgelegt wird.
Die Verteilungen von $Z_{1}$ und $Z_{2}$ sind daher verschieden; es handelt sich also um ein Beispiel zweier unterschiedlicher diskreter multivariater Verteilungen mit gleichen Randverteilungen.

## Bivariate Verteilung

Die Realisationen eines bivariaten Zufallsvektors (einer zweidimensionalen Zufallsvariablen) $Z=(X,Y)$ sind Vektoren in $\mathbb {R} ^{2}=\mathbb {R} \times \mathbb {R}$. Die bivariate Verteilung des Zufallsvektors $Z$ liegt durch die Angabe der Wahrscheinlichkeiten
    $P(Z\in B)\quad {\text{für alle }}B\in {\mathcal {B}}(\mathbb {R} ^{2})$
fest, wobei ${\mathcal {B}}(\mathbb {R} ^{2})$ die Borelsche σ-Algebra auf $\mathbb {R} ^{2}$ bezeichnet. Durch $P_{Z}(B):=P(Z\in B)$ für alle $B\in {\mathcal {B}}(\mathbb {R} ^{2})$ ist ein Wahrscheinlichkeitsmaß $P_{Z}\colon {\mathcal {B}}(\mathbb {R} ^{2})\to [0,1]$ definiert, das $(\mathbb {R} ^{2},{\mathcal {B}}(\mathbb {R} ^{2}),P_{Z})$ zu einem Wahrscheinlichkeitsraum macht. $P_{Z}$ ist die bivariate Wahrscheinlichkeitsverteilung oder kurz bivariate Verteilung von $Z$.

### Bivariate Verteilungsfunktion

Die Verteilung des Zufallsvektors $Z$ liegt bereits dann fest, wenn die Wahrscheinlichkeiten
    $P(Z\in (-\infty ,x]\times (-\infty ,y])\quad {\text{für alle }}x,y\in \mathbb {R}$
gegeben sind. Dies motiviert das Konzept der bivariaten Verteilungsfunktion des Zufallsvektors $Z$, die durch
    $F_{Z}(x,y)=P(X\leq x,Y\leq y)\quad {\text{für alle }}x,y\in \mathbb {R}$
definiert ist. Durch Angabe der Funktion $F_{Z}\colon \mathbb {R} ^{2}\to [0,1]$ liegt die Wahrscheinlichkeitsverteilung von $Z$ fest, da sich aus den durch die Verteilungsfunktion angegebenen Wahrscheinlichkeiten mit den Regeln der Wahrscheinlichkeitsrechnung die Wahrscheinlichkeiten aller anderen Ereignisse ergeben.

### Bivariate Dichtefunktion

Falls der Zufallsvektor $Z$ eine bivariate (oder zweidimensionale) Dichtefunktion $f_{X,Y}$ besitzt, dann besteht zwischen der bivariaten Verteilungsfunktion und der bivariaten Dichtefunktion der Zusammenhang
    $F_{Z}(x,y)=\int _{-\infty }^{y}\int _{-\infty }^{x}f_{X,Y}(u,v)\mathrm {d} u\,\mathrm {d} v$.
Somit liegt durch die Angabe einer bivariaten Dichtefunktion die bivariate Verteilungsfunktion und damit die bivariate Verteilung fest.

### Bivariate Verteilung und bedingte Wahrscheinlichkeiten

Wenn die Zufallsvariable diskret ist, dann kann man die gemeinsame Verteilung mit Hilfe bedingter Wahrscheinlichkeiten so schreiben:
    ${\begin{aligned}\mathrm {P} (X=x\ \mathrm {und} \ Y=y)&=\mathrm {P} (Y=y\mid X=x)\cdot \mathrm {P} (X=x)\\&=\mathrm {P} (X=x\mid Y=y)\cdot \mathrm {P} (Y=y)\end{aligned}}$
und im stetigen Fall entsprechend
    $f_{X,Y}(x,y)=f_{Y|X}(y|x)f_{X}(x)=f_{X|Y}(x|y)f_{Y}(y).$
Hier sind $f_{Y|X}(y|x)$ und $f_{X|Y}(x|y)$ die bedingten Dichten ($Y$ unter der Bedingung $X=x$, bzw. von $X$ unter der Bedingung $Y=y$) und $f_{X}(x),f_{Y}(y)$ die Dichten der Randverteilungen von $X$ und $Y$.

### Stochastische Unabhängigkeit

Die Komponenten des Zufallsvektors $Z=(X,Y)$ sind genau dann stochastisch unabhängig, wenn
    $P(X\leq x,Y\leq y)=P(X\leq x)\cdot P(Y\leq y)\quad {\text{für alle }}x,y\in \mathbb {R}$
gilt. Anderenfalls liegt stochastische Abhängigkeit vor. Wenn für diskrete Zufallsvariablen
    $P(X=x,Y=y)=P(X=x)\cdot P(Y=y)\quad {\text{für alle }}x,y\in \mathbb {R}$
gilt oder wenn für stetige Zufallsvariablen
    $f_{X,Y}(x,y)=f_{X}(x)\cdot f_{Y}(y)\quad {\text{für alle }}x,y\in \mathbb {R}$
gilt, dann sind die Zufallsvariablen $X$ und $Y$ unabhängig.

### Abhängigkeitsstruktur und Copula

Wenn die Komponenten des Zufallsvektors $(X,Y)$ nicht stochastisch unabhängig sind, kann die Abhängigkeitsstruktur durch die sogenannte Copula charakterisiert werden. In der Abbildung ist ein Beispiel für die Modellierung der Abhängigkeitsstruktur mit Hilfe einer speziellen Copula gezeigt. Insbesondere ist das ein Beispiel dafür, dass eine bivariate Zufallsvariable mit normalen Randverteilungen nicht bivariat normalverteilt sein muss.

## Der allgemeine mehrdimensionale Fall

Multivariate Verteilungen ordnen einem geeigneten Teilmengensystem von $\mathbb {R} ^{n}$, dem sogenannten Ereignissystem, Wahrscheinlichkeiten zu. Typischerweise wird als Ereignissystem die Borelsche σ-Algebra auf $\mathbb {R} ^{n}$ gewählt, die mit ${\mathcal {B}}(\mathbb {R} ^{n})$ bezeichnet wird. Aus wahrscheinlichkeitstheoretischer Sicht gilt: „Multivariate Verteilungen sind Wahrscheinlichkeitsmaße auf ${\mathcal {B}}(\mathbb {R} ^{n})$.“
Bei Anwendungen, insbesondere im Bereich der Statistik, steht ein Zufallsvektor (eine $n$-dimensionale Zufallsvariable), $X:\Omega \to \mathbb {R} ^{n}$ im Vordergrund, dessen Komponenten Zufallsvariablen sind. Der Zufallsvektor $X=(X_{1},\dots ,X_{n})$ besitzt dann eine multivariate Verteilung, die das Bildmaß $P_{X}$ von $P$ unter $X$ ist. Dabei gilt
    $P_{X}(B)=P(X\in B)\quad {\text{für alle }}B\in {\mathcal {B}}(\mathbb {R} ^{n})\;.$

### Verteilungsfunktion

Die Verteilungsfunktion $F_{P}$ einer multivariaten Verteilung $P\colon {\mathcal {B}}(\mathbb {R} ^{n})\to [0,1]$ ist durch
    $F_{P}(x_{1},\dots ,x_{n})=P((-\infty ,x_{1}]\times \dots \times (-\infty ,x_{n}])\quad {\text{für alle }}(x_{1},\dots ,x_{n})\in \mathbb {R} ^{n}$
definiert.
Die Verteilungsfunktion $F_{X}$ eines Zufallsvektors $X\colon \Omega \to \mathbb {R} ^{n}$ ist durch
    $F_{X}(x_{1},\dots ,x_{n})=P(X_{1}\leq x_{1},\dots ,X_{n}\leq x_{n})=P(X\in (-\infty ,x_{1}]\times \dots \times (-\infty ,x_{n}])\quad {\text{für alle }}(x_{1},\dots ,x_{n})\in \mathbb {R} ^{n}$
definiert.
Besitzt der Zufallsvektor $X=(X_{1},\dots ,X_{n})$ eine Dichtefunktion $f_{X}\colon \mathbb {R} ^{n}\to [0,\infty )$, dann hängen die Verteilungsfunktion $F_{X}$ und die Dichte $f_{X}$ über
    $F_{X}(x_{1},\dots ,x_{n})=\int _{-\infty }^{x_{n}}\dots \int _{-\infty }^{x_{1}}f_{X}(u_{1},\dots ,u_{n})\mathrm {d} u_{1}\dots \mathrm {d} u_{n}$
zusammen.

### Randverteilungen

Es gibt für Randverteilungen mehr Möglichkeiten als im zweidimensionalen Fall, da nun Randverteilungen für jede niedrigere Dimension $1\leq k<n$ existieren und man ${n \choose k}$ Möglichkeiten hat, einen $k$-dimensionalen Teilvektor aus einem $n$-dimensionalen Vektor auszuwählen. Beispielsweise gibt es im dreidimensionalen Fall zur Verteilung des Zufallsvektors $(X_{1},X_{2},X_{3})$ die drei eindimensionale Randverteilungen der Komponenten $X_{1}$, $X_{2}$ und $X_{3}$ sowie die drei zweidimensionale Randverteilungen der Teilvektoren $(X_{1},X_{2})$, $(X_{1},X_{3})$ und $(X_{2},X_{3})$.

## Verteilungen von Zufallsmatrizen

Eine Zufallsmatrix $\mathbf {X}$ ist eine Matrix, deren Elemente Zufallsvariablen sind. Sie kann als eine Funktion $\mathbf {X} \colon \Omega \to \mathbb {R} ^{m\times n}$ aufgefasst werden. Die Verteilung einer Zufallsmatrix ordnet den (Borel-messbaren) Teilmengen von $\mathbb {R} ^{m\times n}$ Wahrscheinlichkeiten zu.
Im Zusammenhang mit der Schätzung der Kovarianzmatrix einer $d$-dimensionalen multivariaten Normalverteilung spielt eine $d\times d$-dimensionale Wishart-Verteilung eine entscheidende Rolle. Die Wishart-Verteilung ist eine Matrix-Verallgemeinerung der Chi-Quadrat-Verteilung.

## Spezielle multivariate Verteilungen
### Diskrete multivariate Verteilungen

Spezielle multivariate Verteilungen sind z. B. die Multinomialverteilung, die negative Multinomialverteilung, die multivariate Poissonverteilung, die multivariate hypergeometrische Verteilung, die multivariate Pólya-Eggenberger-Verteilung und die multivariate Ewens-Verteilung.

### Stetige multivariate Verteilungen

Spezielle stetige multivariate Verteilungen sind z. B. die multivariate Normalverteilung, verschiedene Konzepte multivariater Exponentialverteilungen, Gammaverteilungen, logistischer Verteilungen, die multivariate Dirichlet-Verteilung, die multivariate Liouville-Verteilung und die multivariaten Extremwertverteilungen.