Erreichbare Korrelationen für logarithmische Zufallsvariablen


19

Betrachten Sie die logarithmischen Zufallsvariablen und mit und .X1X2Log(X1)N(0,1)Log(X2)N(0,σ2)

Ich versuche, und \ rho _ {\ min} für \ rho (X_1, X_2) zu berechnen . Ein Schritt in der gegebenen Lösung, die ich habe, ist:ρmaxρMindestρ(X1,X2)

ρmax=ρ(exp(Z),exp(σZ)) und ρmin=ρ(exp(Z),exp(σZ)) ,

aber sie haben einige Hinweise auf Komonotonie und Gegenkomonotonie gemacht. Ich hatte gehofft, jemand würde mir helfen zu verstehen, wie wichtig sie sind. (Ich weiß, wie man dies aus dem allgemeinen Ausdruck erhält, möchte aber genau wissen, was die Komonotonie-Teile sagten.)


8
Wer sind Sie"?
whuber

Antworten:


25

Ich beginne mit der Definition von Komonotonie und Gegenmonotonie . Dann werde ich erwähnen, warum dies relevant ist, um den minimal und maximal möglichen Korrelationskoeffizienten zwischen zwei Zufallsvariablen zu berechnen. Und schließlich berechne ich diese Grenzen für die logarithmischen Zufallsvariablen X1 und X2 .

Komonotonizität und Gegenmonotonizität
Die Zufallsvariablen gelten als komonoton, wenn ihre Copula die Fréchet-Obergrenze , die die stärkste ist Art der "positiven" Abhängigkeit. Es kann gezeigt werden, dass genau dann sind, wenn wobei eine Zufallsvariable ist, sind zunehmende Funktionen und M ( u 1 , , u d ) = min ( u 1 , , u d ) X 1 , , X d ( X 1 , , X d ) d = ( h 1 ( Z ) , , h d ( Z ) )X1,,Xd M(u1,,ud)=Mindest(u1,,ud)
X1,,XdZ h 1 , , h d d =

(X1,,Xd)=d(h1(Z),,hd(Z)),
Zh1,,hd=dbezeichnet Gleichheit in der Verteilung. Comonotone Zufallsvariablen sind also nur Funktionen einer einzelnen Zufallsvariablen.

Die Zufallsvariablen gelten als gegenmonoton, wenn ihre Copula die Fréchet - Untergrenze , die die stärkste Art der "negativen" Abhängigkeit in der ist bivariater Fall. Die Gegenmonotonizität verallgemeinert sich nicht auf höhere Dimensionen. Es kann gezeigt werden, dass genau dann sind, wenn wobei eine Zufallsvariable ist, und und sind jeweils eine zunehmende und eine abnehmende Funktion oder umgekehrt. W ( u 1 , u 2 ) = max ( 0 , u 1 + u 2 - 1 ) X 1 , X 2 ( X 1 , X 2 ) d = ( h 1 ( Z ) , h 2 ( Z ) ) , Z h 1 h 2X1,X2 W(u1,u2)=max(0,u1+u2-1)
X1,X2

(X1,X2)=d(h1(Z),h2(Z)),
Zh1h2

Erreichbare Korrelation
Sei und zwei Zufallsvariablen mit streng positiven und endlichen Varianzen und bezeichnen und den minimal und maximal möglichen Korrelationskoeffizienten zwischen und . Dann kann gezeigt werden, dassX 2 ρ min ρ max X 1 X 2X1X2ρMindestρmaxX1X2

  • ρ(X1,X2)=ρMindest genau dann, wenn und sind;X 2X1X2
  • ρ(X1,X2)=ρmax genau dann, wenn und comonoton sind.X 2X1X2

Erreichbare Korrelation für logarithmische Zufallsvariablen
Um wir die Tatsache, dass die maximale Korrelation nur dann erreicht wird, wenn und comonoton sind. Die Zufallsvariablen und wobei comonoton sind, da die Exponentialfunktion eine (streng) ansteigende Funktion ist und somit .ρmaxX1X2X1=eZX2=eσZZN(0,1)ρmax=corr(eZ,eσZ)

Unter Verwendung der Eigenschaften lognormaler Zufallsvariablen ergibt sich , , , und die Kovarianz ist Also E(eZ)=e1/2E(eσZ)=eσ2/2veinr(eZ)=e(e-1)veinr(eσZ)=eσ2(eσ2-1)

cOv(eZ,eσZ)=E(e(σ+1)Z)-E(eσZ)E(eZ)=e(σ+1)2/2-e(σ2+1)/2=e(σ2+1)/2(eσ-1).
ρmax=e(σ2+1)/2(eσ-1)e(e-1)eσ2(eσ2-1)=(eσ-1)(e-1)(eσ2-1).

Ähnliche Berechnungen mit ergeben X2=e-σZ

ρMindest=(e-σ-1)(e-1)(eσ2-1).

Kommentar
Dieses Beispiel zeigt, dass es möglich ist, ein Paar von Zufallsvariablen zu haben, die stark abhängig sind - Komonotonie und Gegenmonotonie sind die stärkste Art der Abhängigkeit -, die jedoch eine sehr geringe Korrelation aufweisen. Die folgende Tabelle zeigt diese Grenzen als Funktion von .σ

Bildbeschreibung hier eingeben

Dies ist der R-Code, mit dem ich die obige Tabelle erstellt habe.

curve((exp(x)-1)/sqrt((exp(1) - 1)*(exp(x^2) - 1)), from = 0, to = 5,
      ylim = c(-1, 1), col = 2, lwd = 2, main = "Lognormal attainable correlation",
      xlab = expression(sigma), ylab = "Correlation", cex.lab = 1.2)
curve((exp(-x)-1)/sqrt((exp(1) - 1)*(exp(x^2) - 1)), col = 4, lwd = 2, add = TRUE)
legend(x = "bottomright", col = c(2, 4), lwd = c(2, 2), inset = 0.02,
       legend = c("Correlation upper bound", "Correlation lower bound"))
abline(h = 0, lty = 2)

7
(+6) Schöne ausführliche Darstellung und gut illustriert. Es ist interessant, dass Versuche, Ihr Diagramm durch Simulation zu bestätigen, zum Scheitern verurteilt sind, wenn viel größer als da der Beispielkorrelationskoeffizient extrem variabel ist (aufgrund der Möglichkeit, einen extrem hohen Wert von , der eine hohe Hebelwirkung hat). . Das legt einen höheren Wert als üblich auf eine solide theoretische Analyse. σ3X2
Whuber

5
Diese Darlegung ist eine Anpassung von Beispiel 2.1 (pg. 23) von M. J. und Denuit Dhaene (2003), einfachen Charakterisierungen comonotonicity und countermonotonicity durch extremal Korrelationen , Belgian Actuarial Bulletin , Vol. 3, 22 & ndash; 27.
Kardinal

3
@ Cardinal Mir war dieser Artikel nicht bekannt, danke. Weitere potenzielle Referenzen sind ebooks.cambridge.org/… oder McNeil, AJ, Frey, R. und Embrechts, P. (2005). Quantitatives Risikomanagement: Konzepte, Techniken und Tools. Princeton: Princeton University Press.
QuantIbex

2
Das Beispiel geht mindestens auf RD De Veaux (1976), Enge obere und untere Schranken für die Korrelation der in Luftverschmutzungsmodellen auftretenden bivariaten Verteilungen , Tech. Bericht 5, Abteilung Statistik, Stanford University. Siehe Kapitel 3 ab Seite 6. Die zugrunde liegenden Werkzeuge waren Höffding bekannt.
Kardinal

@QuantIbex in Ihrem Beweis gibt es etwas Unklares für mich. Sie behaupten zunächst, dass und genau dann comonoton sind, wenn ihre gemeinsame Verteilung gleich , für usw., aber wenn Sie dieses Ergebnis auf den logarithmischen Zufall anwenden Variablen, Sie sagen, dass dies impliziert, dass die Zufallsvariablen selbst so sind, dass und , dh, Sie scheinen die Behauptung auf die Zufallsvariablen selbst anzuwenden, nicht nur auf deren Verteilungen. Wie ist das? X1X2(h1(Z),h2(Z))h1,h2X1=eZX1=eσZ
RandomGuy
Durch die Nutzung unserer Website bestätigen Sie, dass Sie unsere Cookie-Richtlinie und Datenschutzrichtlinie gelesen und verstanden haben.
Licensed under cc by-sa 3.0 with attribution required.