Inizierò fornendo la definizione di comonotonicity e countermonotonicity . Quindi, menzionerò perché questo è rilevante per calcolare il coefficiente di correlazione minimo e massimo possibile tra due variabili casuali. E infine, calcolerò questi limiti per le variabili casuali lognormali e .X 2X1X2
Comonotonicità e Si dice che
le variabili casuali sono comonotoniche se la loro copula è il limite superiore di Fréchet , che è il più forte tipo di dipendenza "positiva".
Si può dimostrare che sono comonotonici se e solo se
dove è una variabile casuale, sono funzioni crescenti e
X1,…,Xd X 1 , … , X d ( X 1 , … , X d ) d = ( h 1 ( Z ) , … , h d ( Z ) ) , Z h 1 , … ,M(u1,…,ud) = min(u1, ...,ud)
X1, ... , Xd
( X1, ... , Xd) =d( h1( Z) , ... , hd(Z) ) ,
Zd =h1,…,hd=ddenota l'uguaglianza nella distribuzione. Quindi, le variabili casuali comonotoniche sono solo funzioni di una singola variabile casuale.
Le variabili aleatorie sono detti countermonotonic se la loro copula rappresenta l' Fréchet limite inferiore , che è il tipo più forte di dipendenza "negativo" nel caso bivariato. La contromonotonocità non si generalizza a dimensioni superiori.
Si può dimostrare che sono contrmonotonici se e solo se
dove è una variabile casuale, e e sono rispettivamente una funzione crescente e una diminuzione, o viceversa.X1,X2 X 1 , X 2 ( X 1 , X 2 ) d = ( h 1 ( Z ) , h 2 ( Z ) ) , Z h 1 h 2W(u1,u2)=max(0,u1+u2−1)
X1,X2
(X1,X2)=d(h1(Z),h2(Z)),
Zh1h2
Correlazione ottenibile
Sia e due variabili casuali con varianze rigorosamente positive e finite, e lascia che e denotino il coefficiente di correlazione minimo e massimo possibile tra e . Quindi, si può dimostrare cheX 2 ρ min ρ max X 1 X 2X1X2ρminρmaxX1X2
- ρ(X1,X2)=ρmin se e solo se e sono contronotonici;X 2X1X2
- ρ(X1,X2)=ρmax se e solo se e sono comonotonici.X 2X1X2
Correlazione
ottenibile per variabili casuali lognormali Per ottenere utilizziamo il fatto che si ottiene la massima correlazione se e solo se e sono comonotonici. Le variabili casuali e dove sono comonotoniche poiché la funzione esponenziale è una funzione (strettamente) crescente, e quindi .ρmaxX1X2X1=eZX2=eσZZ∼N(0,1)ρmax=corr(eZ,eσZ)
Usando le proprietà delle variabili casuali lognormali , abbiamo
,
,
,
e la covarianza è
Pertanto,
E(eZ)=e1/2E(eσZ)=eσ2/2Z)=e(e−1)var(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)−−−−−−−−−−−−√.
simili con yield
X2=e−σZ
ρmin=(e−σ−1)(e−1)(eσ2−1)−−−−−−−−−−−−√.
Commento
Questo esempio mostra che è possibile avere una coppia di variabili casuali che sono fortemente dipendenti - la comonotonicità e la contromanotonicità sono il tipo più forte di dipendenza - ma che hanno una correlazione molto bassa. La seguente tabella mostra questi limiti in funzione di .σ
Questo è il codice R che ho usato per produrre il grafico sopra.
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)