Distribuzione del massimo di due variabili normali correlate


Risposte:


22

Secondo Nadarajah e Kotz, 2008 , distribuzione esatta del massimo / minimo di due variabili casuali gaussiane , il PDF di X=max(X1,X2) sembra essere

f(X)=2φ(X)Φ(1-r1-r2X),

dove φ è il PDF e Φ è il CDF della distribuzione normale standard.

inserisci qui la descrizione dell'immagine


Come sarà se (nessuna correlazione)? Ho problemi a visualizzarlo. r=0
Mitch,

3
Ho aggiunto una figura che visualizza la distribuzione. Sembra un gaussiano schiacciato leggermente inclinato a destra.
Lucas,

22

Sia il PDF normale bivariato per con marginali standard e correlazione . Il CDF del massimo è, per definizione, ( X , Y ) ρfρ(X,Y)ρ

Pr(max(X,Y)z)=Pr(Xz, Yz)=zzfρ(x,y)dydx.

Il PDF normale bivariato è simmetrico (tramite riflessione) attorno alla diagonale. Pertanto, aumentando da a aggiungono due strisce di probabilità equivalente al quadrato semi-infinito originale: quello superiore infinitesimamente spesso è mentre la sua controparte riflessa, il striscia di destra, è .z + d z ( - , z ] × ( z , z + d z ] ( z , z + d z ] × ( - , z ]zz+dz(-,z]×(z,z+dz](z,z+dz]×(-,z]

figura

La densità di probabilità della striscia di destra è la densità di a volte la probabilità condizionata totale che sia nella striscia, . La distribuzione condizionale di è sempre normale, quindi per trovare questa probabilità condizionale totale abbiamo solo bisogno della media e della varianza. La media condizionale di in è la previsione di regressione e la varianza condizionale è la varianza "inspiegabile" .z Y Pr ( Y zXzYY Y X ρ X var ( Y ) - var ( ρ X ) = 1 - ρ 2Pr(Yz|X=z)YYXρXvar(Y)-var(ρX)=1-ρ2

Ora che conosciamo la media condizionale e la varianza, il CDF condizionale di dato può essere ottenuto standardizzando Y e applicando il normale CDF Φ :XYXYΦ

Pr(Yy|X)=Φ(y-ρX1-ρ2).

Valutandolo in e X = z e moltiplicandolo per la densità di X in z (un PDF normale normale ϕ ) si ottiene la densità di probabilità della seconda striscia (a destra)y=zX=zXzφ

φ(z)Φ(z-ρz1-ρ2)=φ(z)Φ(1-ρ1-ρ2z).

Raddoppiando questo si spiega la striscia superiore equi-probabile, dando il PDF del massimo come

ddzPr(max(X,Y)z)=2φ(z)Φ(1-ρ1-ρ2z).

Ricapitolazione

2φ(z)Φ()1-ρ1-ρ2zY=zX=z


Questo può essere esteso a più di due variabili normali standard con una data matrice di correlazione?
A. Donda,

1
@ A.Donda Sì, ma l'espressione diventa più complicata. Ad ogni nuova dimensione nasce la necessità di integrarsi ancora una volta.
whuber
Utilizzando il nostro sito, riconosci di aver letto e compreso le nostre Informativa sui cookie e Informativa sulla privacy.
Licensed under cc by-sa 3.0 with attribution required.