Erwartung an Produkte höherer Ordnung von Normaldistributionen

9

Ich habe zwei normalverteilte Variablen und mit mittlerer Null und Kovarianzmatrix . Ich bin daran interessiert, den Wert von anhand der Einträge von zu berechnen .X 2 Σ E [ X 2 1 X 2 2 ] ΣX1X2ΣE[X12X22]Σ

Ich habe das Gesetz der Gesamtwahrscheinlichkeit verwendet, um aber ich bin mir nicht sicher, auf was sich die innere Erwartung reduziert. Gibt es hier eine andere Methode?E[X12X22]=E[X12E[X22|X1]]

Vielen Dank.

Bearbeiten: Die Variablen sind auch multivariat normalverteilt.

AGK
quelle
5
Do und genießen Sie eine bivariate auch Normalverteilung? (Nur zu sagen, dass und mit der Kovarianzmatrix normal sind, reicht nicht aus, um zu dem Schluss zu kommen, dass die gemeinsame Verteilung bivariat normal ist). X 2X1X2X 2 ΣX1X2Σ
Dilip Sarwate
1
Für die spezifische Anwendung, an die ich denke, haben und eine bivariate Normalverteilung nach dem multivariaten zentralen Grenzwertsatz. Ich habe vergessen, dies in meinem ursprünglichen Beitrag zu erwähnen. X 2X1X2
AGK
1
@AGK Wenn Sie Ihren Beitrag klarstellen möchten, gibt es eine Schaltfläche "Bearbeiten", mit der Sie Änderungen vornehmen können. Dies ist besser für zukünftige Leser, die dann keine wichtigen Informationen in den Kommentaren unter der Frage nachschlagen müssen.
Silverfish

Antworten:

8

Die Erwartung ist eindeutig proportional zum Produkt der quadratischen Skalierungsfaktoren . Die Proportionalitätskonstante wird durch Standardisierung der Variablen erhalten, wodurch auf die Korrelationsmatrix mit der Korrelation reduziert wird . Σ ρ = σ 12 / σ11σ22Σρ=σ12/σ11σ22

Unter der Annahme einer bivariaten Normalität können wir gemäß der Analyse unter https://stats.stackexchange.com/a/71303 Variablen in ändern

X1=X, X2=ρX+(1ρ2)Y

Dabei hat eine standardmäßige (unkorrelierte) bivariate Normalverteilung, und wir müssen nur berechnen(X,Y)

E(X2(ρX+(1ρ2)Y)2)=E(ρ2X4+(1ρ2)X2Y2+cX3Y)

wobei der genaue Wert der Konstanten keine Rolle spielt. ( ist der Rest bei der Regression von gegen .) Verwenden der univariaten Erwartungen für die StandardnormalverteilungcYX2X1

E(X4)=3, E(X2)=E(Y2)=1, EY=0

und Feststellung , dass und sind unabhängig AusbeutenXY

E(ρ2X4+(1ρ2)X2Y2+cX3Y)=3ρ2+(1ρ2)+0=1+2ρ2.

Multipliziert man dies durch gibtσ11σ22

E(X12X22)=σ11σ22+2σ122.

Die gleiche Methode gilt für das Ermitteln der Erwartung eines Polynoms in , da es in und zu einem Polynom wird dass, wenn er expandiert ist , ein Polynom in den unabhängigen normalverteilte Variablen und . Von(X1,X2)(X,ρX+(1ρ2)Y)XY

E(X2k)=E(Y2k)=(2k)!k!2k=π1/22kΓ(k+12)

für das Integral (wobei alle ungeraden Momente durch Symmetrie gleich Null sind) können wir ableitenk0

E(X12pX22q)=(2q)!2pqi=0qρ2i(1ρ2)qi(2p+2i)!(2i)!(p+i)!(qi)!

(mit allen anderen Erwartungen an Monome gleich Null). Dies ist proportional zu einer hypergeometrischen Funktion (fast per Definition: Die Manipulationen sind nicht tief oder lehrreich).

1π2p+q(1ρ2)qΓ(p+12)Γ(q+12)2F1(p+12,q;12;ρ2ρ21).

Die hypergeometrischen Funktionszeiten als multiplikative Korrektur für ungleich Null .ρ(1ρ2)qρ

whuber
quelle
1
Danke für die ausführliche Antwort! Ich denke auch über verwandte Fragen mit anderen Polynomen nach, daher ist dies ein wirklich hilfreicher Rahmen. Das ist eine sehr clevere Transformation, die ich vorher noch nicht gesehen hatte. Cool!
AGK
3
Um Ihre Untersuchung zu erleichtern, habe ich die Details für allgemeine Polynome angegeben. Als ich diese Antwort ursprünglich schrieb, war ich amüsiert zu erkennen, dass ich diese Transformation aus dem Lehrbuch für Grundstatistiken von Friedman, Pisani und Purves gelernt habe: Wir bringen dies Studienanfängern bei!
whuber