σ11σ22Σρ=σ12/σ11σ22−−−−−√
En supposant une normalité bivariée, alors selon l'analyse à https://stats.stackexchange.com/a/71303, nous pouvons changer les variables en
X1=X, X2=ρX+(1−ρ2−−−−−√)Y
(X,Y)
E(X2(ρX+(1−ρ2−−−−−√)Y)2)=E(ρ2X4+(1−ρ2)X2Y2+cX3Y)
cYX2X1
E(X4)=3, E(X2)=E(Y2)=1, EY=0
XY
E(ρ2X4+(1−ρ2)X2Y2+cX3Y)=3ρ2+(1−ρ2)+0=1+2ρ2.
σ11σ22
E(X21X22)=σ11σ22+2σ212.
(X1,X2)(X,ρX+(1−ρ2−−−−−√)Y)XY
E(X2k)=E(Y2k)=(2k)!k!2k=π−1/22kΓ(k+12)
k≥0
E(X2p1X2q2)=(2q)!2−p−q∑i=0qρ2i(1−ρ2)q−i(2p+2i)!(2i)!(p+i)!(q−i)!
(avec toutes les autres attentes de monômes égales à zéro). Ceci est proportionnel à une fonction hypergéométrique (presque par définition: les manipulations impliquées ne sont ni profondes ni instructives),
1π2p+q(1−ρ2)qΓ(p+12)Γ(q+12)2F1(p+12,−q;12;ρ2ρ2−1).
(1−ρ2)qρ