Copula (Mathematik)

Eine Copula ist eine Funktion, die einen funktionalen Zusammenhang zwischen den Randverteilungsfunktionen verschiedener Zufallsvariablen und ihrer gemeinsamen Wahrscheinlichkeitsverteilung angeben kann.

Mit ihrer Hilfe kann man stochastische Abhängigkeit deutlich komplexer modellieren als beispielsweise mit dem Korrelationskoeffizienten.

Definition

Eine Copula ist eine Verteilungsfunktion

$C:[0,1]^{n}\rightarrow [0,1]$,

deren eindimensionale Randverteilungen gleichverteilt über dem Intervall $[0,1]$ sind. Formal ausgedrückt bedeutet dies:

$C(1,\ldots ,1,\underbrace {u_{j}} _{j{\text{-te Stelle}}},1,\ldots ,1)=u_{j}\quad \forall j\in \{1,\ldots ,n\}.$

Die Forderung an die Randverteilungen lässt sich wie folgt motivieren: Für $n\in \mathbb {N}$ beliebig verteilte Zufallsvariablen $X_{1},X_{2},\ldots ,X_{n}$ mit stetigen Verteilungen $F_{X_{i}},\;i=1,2,\ldots ,n$ ist die Zufallsvariable $F_{X_{i}}(X_{i})$ gleichverteilt über dem Intervall $[0,1]$. Zusammen mit dem folgenden Satz von Sklar wird die Trennung von Randverteilungen und Abhängigkeiten unter diesen möglich.

Satz von Sklar

Im Folgenden sei ${\overline {\mathbb {R} }}:=\mathbb {R} \cup \{-\infty ,+\infty \}$ eine Erweiterung der reellen Zahlen.
Sei $F:{\overline {\mathbb {R} }}^{n}\rightarrow [0,1]$ eine $n$-dimensionale Verteilungsfunktion mit eindimensionalen Randverteilungen $F_{1},\ldots ,F_{n}:{\overline {\mathbb {R} }}\rightarrow [0,1]$. Dann existiert eine $n$-dimensionale Copula $C$, so dass für alle $(x_{1},\ldots ,x_{n})\in {\overline {\mathbb {R} }}^{n}\$ gilt:

$F(x_{1},x_{2},\ldots ,x_{n})=C\left(F_{1}\left(x_{1}\right),\ldots ,F_{n}\left(x_{n}\right)\right)$

Sind alle $F_{i}$ stetig, so ist die Copula eindeutig.

Fréchet-Hoeffding-Schranken

Für jede $n$-variate Copula gilt:

• $W(u_{1},\ldots ,u_{n}):=\max \left\{1-n+\sum \limits _{i=1}^{n}{u_{i}},0\right\}\leq C(u_{1},\ldots ,u_{n})$ und
• $C(u_{1},\ldots ,u_{n})\leq \min _{j\in \{1,\ldots ,n\}}u_{j}=:M(u_{1},\ldots ,u_{n})$

$W$ wird als untere und $M$ als obere Fréchet-Hoeffding-Schranke bezeichnet.

Damit lässt sich unter anderem folgern:

• $C(u_{1},\ldots ,u_{j-1},0,u_{j+1},\ldots ,u_{n})=0$ für alle $j\in \{1,\ldots ,n\}$

Die untere Schranke, $W$, ist selbst auch eine Copula. Für $M$ gilt dies nicht.

Anwendung

Copulas werden eingesetzt, um Rückschlüsse auf die Art der stochastischen Abhängigkeit verschiedener Zufallsvariablen zu erzielen oder um Abhängigkeiten gezielt zu modellieren. Sie werden beispielsweise in der Kreditrisikoanalyse eingesetzt, um Aussagen über einen gehäuften Bankrott mehrerer Schuldner innerhalb eines Anleihenportfolios machen zu können. Analog sind Anwendungen im Versicherungsbereich üblich. Dort stellen gehäuft auftretende Schäden verschiedener Schadenarten ein finanzielles Problem dar. Beispiel hierfür ist ein zu beobachtender Zusammenhang zwischen Sturm- und Hochwasserschäden.

Beispiele für Copulas

• Die einfachste Form der Copula ist die Unabhängigkeits-Copula

$C(u_{1},\ldots ,u_{n})=\prod \limits _{i=1}^{n}u_{i}=u_{1}\cdot \ldots \cdot u_{n}$.

Sie steht für stochastisch unabhängige Zufallsvariablen $U_{1},\ldots ,U_{n}$, die gemäß der Copula C verteilt sind. In Zeichen: $(U_{1},\ldots ,U_{n})\sim C$

• Die obere Fréchet-Hoeffding-Schranke, ebenfalls eine Copula, ist gegeben durch

$C(u_{1},\ldots ,u_{n})=\min _{i=1,\ldots ,n}u_{i}$

Sie beschreibt perfekte positive stochastische Abhängigkeit (totale positive Korrelation).

• Die Normal- oder auch Gauß-Copula wird mit Hilfe der Verteilungsfunktion der Normalverteilung $F(\cdot )$ definiert. So ist

$C(u_{1},u_{2})=F_{2}(F^{-1}(u_{1}),F^{-1}(u_{2}),\rho )\,$

eine Copula, wobei $F_{2}(\cdot ,\cdot ,\rho )$ die bivariate Verteilungsfunktion zweier standard-normalverteilter Zufallsvariablen mit dem Korrelationskoeffizienten $\rho$ ist.

Erzeugt man Punkte, die gemäß der Normal-Copula mit Parameter $\rho =0.5$ verteilt sind, ergibt sich bereits eine leichte Konzentration dieser entlang der Winkelhalbierenden.

Simulation der bivariaten Normal-Copula, rho = 0.5, 1500 Punkte

• Die Gumbel-Copula wird mit Hilfe der Exponentialfunktion und dem natürlichen Logarithmus definiert

$C_{\lambda }(u_{1},u_{2})=\exp \left(-\left(\left(-\ln u_{1}\right)^{\lambda }+\left(-\ln u_{2}\right)^{\lambda }\right)^{1/\lambda }\right)$,

wobei $\lambda \geq 1$ als Parameter fest zu wählen ist.

Erzeugt man Punkte, die gemäß der Gumbel-Copula mit Parameter $\lambda >1$ verteilt sind, ergibt sich insbesondere eine Punkthäufung in der Nähe des Punktes $(1,1)$.

Simulation der bivariaten Gumbel-Copula, lambda = 2, 1500 Punkte

Archimedische Copulas

Archimedische Copulas stellen eine Klasse von Copulas dar. Diese lassen sich wie folgt beschreiben:

Sei $\ \varphi :[0,1]\rightarrow [0,\infty ]\$ eine stetige, streng monoton fallende Funktion mit $\ \varphi (1)=0$. Bezeichne $\ \varphi ^{[-1]}:[0,\infty ]\rightarrow [0,1]\$ die Pseudo-Inverse von $\ \varphi$, d.h. $\varphi ^{[-1]}(t)\ :=\ {\begin{cases}\varphi ^{-1}(t),&{\textrm {falls}}\ 0\leq t\leq \varphi (0)\\0,&{\textrm {sonst}}\end{cases}}$

Mit Hilfe von $\ \varphi \$ und $\ \varphi ^{[-1]}\$ lässt sich nun eine bivariate Funktion definieren: $C:[0,1]^{2}\rightarrow [0,1],\quad C(u,v):=\varphi ^{[-1]}\left(\varphi \left(u\right)+\varphi \left(v\right)\right)$

Die Funktion $\ C\$ ist genau dann eine Copula, wenn $\ \varphi \$ konvex ist. In diesem Fall heißt $\ \varphi \$ Erzeuger der Copula. Offensichtlich ist $\ C\$ symmetrisch, d.h. $\ C(u,v)=C(v,u)\$ für alle $\ u,v\in [0,1]$.

Beispiele für archimedische Copulas sind:

• Gumbel-Copula: Ihr Erzeuger ist die Funktion $\varphi (t)=(-\ln t)^{\lambda }$ mit Parameter $\lambda \geq 1$.

Damit ergibt sich $\varphi ^{[-1]}(t)=\exp \left(-t^{\frac {1}{\lambda }}\right)$ und damit die Gumbel-Copula $C_{\lambda }(u,v)$ wie oben.

• Clayton-Copula: Ihr Erzeuger ist die Funktion $\varphi (t)={\frac {1}{\Theta }}\left(t^{-\Theta }-1\right)$ mit Parameter $\Theta >0$.

Damit ist $\varphi ^{[-1]}(t)=\left(\Theta \cdot t+1\right)^{-{\frac {1}{\Theta }}}$ und die bivariate Clayton-Copula ergibt sich zu: $C(u,v)=\max \left\{0,\left(u^{-\Theta }+v^{-\Theta }-1\right)^{-{\frac {1}{\Theta }}}\right\}$

• Frank-Copula: Ihr Erzeuger ist die Funktion $\varphi (t)=-\ln \left({\frac {e^{-\Theta \cdot t}-1}{e^{-\Theta }-1}}\right)$ mit Parameter $\Theta >0$.

Extremwertcopula

Definition

Eine Copula $C$ heißt Extremwertcopula, wenn es die Copula einer multivariaten Extremwertverteilung ist, d. h. es exisitiert eine multivariate Extremwertverteilung $G$ mit univariaten Rändern $G_{1},\dots ,G_{n}$, dass gilt $C(u_{1},\dots ,u_{n})=G(G_{1}^{-1}(u_{1}),\dots ,G_{n}^{-1}(u_{n}))$.

Lemma

Eine Copula $C$ ist genau dann eine Extremwertcopula, wenn für $\mathbf {0} \leq \mathbf {u} =(u_{1},\dots ,u_{n})^{T}\leq \mathbf {1}$ und $t>0$ gilt $C(u_{1}^{t},\dots ,u_{n}^{t})=C^{t}(u_{1},\dots ,u_{n})$.

Ist $C$ eine Extremwertcopula und sind $G_{1},\dots ,G_{n}$ univariate Extremwertverteilungen, dann ist $G((x_{1},\dots ,x-n)^{T}):=C(G_{1}(x_{1}),\dots ,G_{n}(x_{n}))$ eine multivariate Extremwertverteilung.

Weblinks

• http://library.soa.org/library/naaj/1997-09/naaj9801_1.pdf - Understanding relationships using copulas, Frees, E. W., Valdez, E. A. (1998), North American Actuarial Journal, 2, pp. 1-25.

• http://www.math.ethz.ch/~baltes/ftp/copchapter.pdf - Modelling Dependence with Copulas and Applications to Risk Management, Embrechts, P., Lindskog, F., McNeil, A. (2003), Handbook of Heavy Tailed Distributions in Finance, ed. S. Rachev, Elsevier, Chapter 8, pp. 329-384.

• http://www.math.ethz.ch/%7Estrauman/preprints/pitfalls.pdf - Correlation and dependence in risk management: properties and pitfalls, Embrechts, P., McNeil, A., Straumann, D. (2002), Risk Management: Value at Risk and Beyond, ed. M.A.H. Dempster, Cambridge University Press, Cambridge, pp. 176-223

Literatur

• Nelsen, Roger B.: An Introduction to Copulas (Lecture Notes in Statistics), Springer Verlag, 2006, ISBN 0387286594

• Sklar, A.: Random variables, distribution functions, and copulas – a personal look backward and forward in Rüschendorf, L., Schweizer, B. und Taylor, M. (eds) Distributions With Fixed Marginals & Related Topics (Lecture Notes - Monograph Series Number 28), 1997, ISBN 0940600404