Die Gumbel-Verteilung

Stetige Verteilungen bzw. Verteilungen von stetigen Zufallsvariablen lassen sich durch Dichtefunktionen beschreiben, die oft von einem oder mehreren reellen Parametern abhängen. Verschiedene Parameter können sich unterschiedlich auf die Dichtefunktion auswirken:

Die Gumbel-Verteilung \(\text{Gum}(\mu,\beta)\) (nach Emil Julius Gumbel, 1891 - 1966) wird beschrieben durch einen Lageparameter \(\mu\in\mathbb{R}\) und einen Skalenparameter \(\beta>0\). Die Dichtefunktion ist \(f_{\mu,\beta}:\mathbb{R}\longrightarrow\mathbb{R}\) mit \[f_{\mu,\beta}(x)=\frac{1}{\beta}\exp\left(-\frac{x-\mu}{\beta}\right)\exp\left(-\exp\left(-\frac{x-\mu}{\beta}\right)\right).\]

Da Lagemaße wie der Erwartungswert und Streuungsmaße wie die Standardabweichung mithilfe der Dichtefunktion berechnet werden können, hängen auch diese direkt von den Parametern ab. Für \(X\sim\text{Gum}(\mu,\beta)\) ist der Erwartungswert \[\mu_X=\text{E}(X)=\mu+\beta\gamma\] (hierbei ist \(\gamma\approx0.57721\) die Euler-Mascheroni-Konstante) und die Standardabweichung ist \[\sigma_X=\sqrt{\text{Var}(X)}=\pi\beta/\sqrt{6}.\]

Die Funktion \(f_{\mu,\beta}\) ist eine Dichtefunktion, denn \(f_{\mu,\beta}(x)\geq0\) für alle \(x\in\mathbb{R}\) und mit der Substitution \(y=\exp(-(x-\mu)/\beta)\), d.h. \(x=\mu-\beta\log(y)\) und \(\text{d}x=-\beta y^{-1}\text{ d}y\), gilt \begin{align*} \int\limits_{-\infty}^\infty f_{\mu,\beta}(x)\text{ d}x &=\dfrac{1}{\beta}\int\limits_{-\infty}^\infty\exp\left(-\frac{x-\mu}{\beta}\right)\exp\left(-\exp\left(-\frac{x-\mu}{\beta}\right)\right)\text{ d}x\\ &=\dfrac{1}{\beta}\int\limits_{\infty}^0ye^{-y}\dfrac{-\beta}{y}\text{ d}y\\ &=\int\limits_0^\infty e^{-y}\text{ d}y\\ &=\Bigg[-e^{-y}\Bigg]_{y=0}^\infty\\ &=1. \end{align*} Bei der Berechnung des Erwartungswerts und der Varianz der Gumbel-Verteilung werden die folgenden Identitäten für die Euler-Mascheroni-Konstante \(\gamma\) benötigt: \begin{equation*} \int\limits_0^\infty\log(y)e^{-y}\text{ d}y=-\gamma \quad\text{und}\quad \int\limits_0^\infty\log(y)^2e^{-y}\text{ d}y=\gamma^2+\frac{\pi^2}{6}. \end{equation*} Sei \(X\sim\text{Gum}(\mu,\beta)\). Wieder mit der Substitution \(y=\exp(-(x-\mu)/\beta)\) ist der Erwartungswert von \(X\) \begin{align*} \text{E}(X) &=\int\limits_{-\infty}^\infty xf_{\mu,\beta}(x)\text{ d}x\\ &=\frac{1}{\beta}\int\limits_{-\infty}^\infty x\exp\left(-\frac{x-\mu}{\beta}\right)\exp\left(-\exp\left(-\frac{x-\mu}{\beta}\right)\right)\text{ d}x\\ &=\frac{1}{\beta}\int\limits_\infty^0\big(\mu-\beta\log(y)\big)ye^{-y}\frac{-\beta}{y}\text{ d}y\\ &=\mu\int\limits_0^\infty e^{-y}\text{ d}y-\beta\int\limits_0^\infty\log(y)e^{-y}\text{ d}y\\ &=\mu+\beta\gamma. \end{align*} Mit derselben Substitution ist das zweite Moment von \(X\) \begin{align*} \text{E}(X^2) &=\int\limits_{-\infty}^\infty x^2f_{\mu,\beta}(x)\text{ d}x\\ &=\frac{1}{\beta}\int\limits_{-\infty}^\infty x^2\exp\left(-\frac{x-\mu}{\beta}\right)\exp\left(-\exp\left(-\frac{x-\mu}{\beta}\right)\right)\text{ d}x\\ &=\frac{1}{\beta}\int\limits_\infty^0\big(\mu-\beta\log(y)\big)^2ye^{-y}\frac{-\beta}{y}\text{ d}y\\ &=\mu^2\int\limits_0^\infty e^{-y}\text{ d}y-2\mu\beta\int\limits_0^\infty\log(y)e^{-y}\text{ d}y+\beta^2\int\limits_0^\infty\log(y)^2e^{-y}\text{ d}y\\ &=\mu^2+2\mu\beta\gamma+\beta^2\left(\gamma^2+\frac{\pi^2}{6}\right)\\ &=(\mu+\beta\gamma)^2+\frac{(\pi\beta)^2}{6}. \end{align*} Somit ist die Varianz \begin{align*} \text{Var}(X) &=\text{E}(X^2)-\text{E}(X)^2\\ &=(\mu+\beta\gamma)^2+\frac{(\pi\beta)^2}{6}-(\mu+\beta\gamma)^2\\ &=\frac{(\pi\beta)^2}{6} \end{align*} und wir erhalten schließlich die Standardabweichung \(\sigma_X=\sqrt{\text{Var}(X)}=\pi\beta/\sqrt{6}\).

Erkunden Sie in der folgenden interaktiven Grafik, wie sich die Parameter der Gumbel-Verteilung auf ihre Dichte- und Verteilungsfunktion sowie auf den Erwartungswert und die Standardabweichung auswirken. Außerdem können Sie Intervallwahrscheinlichkeiten und Quantile grafisch berechnen lassen.

Wechseln Sie mit den Radiobuttons zwischen den Ansichten der Funktionsgraphen, Intervallwahrscheinlichkeiten und Quantilen. Blenden Sie mit den Checkboxen die Dichtefunktion (durchgezogene Kurve) bzw. die Verteilungsfunktion (gestrichelte Kurve) ein oder aus. Verändern Sie die Werte der Parameter über die Schieberegler oder die Eingabefelder.



Zum Nachdenken

(a) Stellen Sie in der interaktiven Grafik verschiedene Werte von \(\mu\) und \(\beta\) ein und beschreiben Sie anhand Ihrer Beobachtungen qualitativ die Abhängigkeit des Erwartungswerts und der Standardabweichung der Gumbel-Verteilung von den Parametern \(\mu\) und \(\beta\).

(b) Beschreiben Sie, wie sich die Parameter \(\mu\) und \(\beta\) qualitativ auf die Form der Dichtefunktion und der Verteilungsfunktion auswirken.

(a) Der Erwartungswert der Gumbel-Verteilung hängt sowohl von \(\mu\) als auch von \(\beta\) ab und ist in beiden Parametern streng monoton steigend. Die Standardabweichung der Gumbel-Verteilung hängt nur von \(\beta\) ab (nicht von \(\mu\)) und ist streng monoton steigend in \(\beta\).

(b) Der Lageparameter \(\mu\) bestimmt die Lage des Hochpunkts der Dichtefunktion sowie die Lage des Wendepunkts der Verteilungsfunktion. Bei steigendem \(\mu\) verschieben sie sich nach rechts, bei fallendem \(\mu\) nach links. Der Skalenparameter \(\beta\) bestimmt die Breite der Funktionsgraphen. Je größer bzw. kleiner der Wert von \(\beta\), desto flacher bzw. steiler sind die Funktionsgraphen.