Kovarianzmatrix für die Gaußsche Prozess- und Wishart-Verteilung

11

Ich lese dieses Papier über generalisierte Wishart-Prozesse (GWP) durch. Die Arbeit berechnet die Kovarianzen zwischen verschiedenen Zufallsvariablen (nach dem Gaußschen Prozess ) unter Verwendung der quadratischen exponentiellen Kovarianzfunktion, dh . Es heißt dann, dass diese Kovarianzmatrix GWP folgt.K(x,x)=exp(|(xx)|22l2)

Früher dachte ich, dass eine Kovarianzmatrix, die aus der linearen Kovarianzfunktion ( )K(x,x)=xTx berechnet wurde , der Wishart-Verteilung mit geeigneten Parametern folgt.

Meine Frage ist, wie können wir immer noch annehmen, dass die Kovarianz einer Wishart-Verteilung mit quadratischer exponentieller Kovarianzfunktion folgt? Was ist im Allgemeinen die notwendige Bedingung für eine Kovarianzfunktion, um eine verteilte Wishart-Kovarianzmatrix zu erzeugen?

Steadyfish
quelle

Antworten:

8

Was verwechselt wird, ist die Kovarianzspezifikation in Bezug auf den Umgebungsraum, auf dem der Gaußsche Prozess definiert ist, und die Operation, die eine endlich dimensionale Gaußsche Zufallsvariable transformiert, um eine Wishart-Verteilung zu erhalten.

Wenn eine p- dimensionale Gaußsche Zufallsvariable (ein Spaltenvektor) mit dem Mittelwert 0 und der Kovarianzmatrix Σ ist , ist die Verteilung von W = X X T eine Wishart-Verteilung W p ( Σ , 1 ) . Es ist zu beachten, dass W eine p × p- Matrix ist. Dies ist ein allgemeines Ergebnis darüber, wie die quadratische Form xx x T istXN(0,Σ)pΣW=XXTWp(Σ,1)Wp×p

xxxT.
wandelt eine Gaußsche Verteilung in eine Wishart-Verteilung um. Es gilt für jede Wahl einer positiven definitiven Kovarianzmatrix . Wenn Sie iid Beobachtungen X 1 , , X n haben, dann ist mit W i = X i X T i die Verteilung von W 1 + + W n eine Wishart W p ( Σ , n ) -Verteilung. Teilen durch nΣX.1,,X.nW.ich=X.ichX.ichT.
W.1++W.n
W.p(Σ,n)n wir dividieren, erhalten wir die empirische Kovarianzmatrix eine Schätzung von Σ- -Σ.

Für Gaußsche Prozesse gibt es einen Umgebungsraum, beispielsweise zur Veranschaulichung, dass es sich um , sodass die betrachteten Zufallsvariablen durch Elemente im Umgebungsraum indiziert werden. Das heißt, betrachten wir einen Prozess ( X ( x ) ) x R . Es ist Gaußsch (und der Einfachheit halber hier mit dem Mittelwert 0), wenn seine endlichen dimensionalen Randverteilungen Gaußsch sind, dh wenn X ( x 1 , , x p ) : = ( X ( x 1 ) , , X ( x)R.(X.(x))xR. für alle x 1 , ... , x pR . Die Wahl derKovarianzfunktion, wie vom OP erwähnt, bestimmt die Kovarianzmatrix, dh cov ( X ( x i ) , X ( x j ) ) =

X.(x1,,xp): =(X.(x1),,X.(xp))T.N.(0,Σ(x1,,xp))
x1,,xpR. Ohne Berücksichtigung der Wahl von K ist die Verteilung von X ( x 1 , , x p ) X ( x 1 , , x p ) T ein Wishart W p ( Σ ( x 1 , , x p )).
cov(X.(xich),X.(xj))=Σ(x1,,xp)ich,j=K.(xich,xj).
K.
X.(x1,,xp)X.(x1,,xp)T.
-Verteilung.W.p(Σ(x1,,xp),1)
NRH
quelle
Vielen Dank für die Beantwortung. Ich habe ein paar Fragen, reg. Ihre Antwort - Wenn Sie sagen, dass die Transformation, die Gaußsche dist in Wishart dist transformiert, für jede Wahl einer + ve bestimmten Cov-Matrix gilt, welche unterschiedlichen Möglichkeiten haben wir für diese Cov-Matrix? Um die durch die cov-Funktion definierte cov-Matrix zu verdeutlichen, geben i und j Elemente im Umgebungsraum des Gaußschen Prozesses an (wenn es sich beispielsweise um einen zeitlichen Prozess handelt, werden die Zeitpunkte t_1 und t_2).
Steadyfish
ichjxichxjΣ - -ΣΣ
xT.x
@steadyfish, oh, ich verstehe. Tatsächlich war ich schlampig mit den Transpositionen und ob die Vektoren Zeilen- oder Spaltenvektoren waren. Ich habe das jetzt präzisiert und ein wenig über die Beziehung zwischen der empirischen Kovarianzmatrix und der theoretischen Kovarianzmatrix hinzugefügt. Die Theorie ist nicht in Bezug auf die Beobachtungen definiert.
NRH