Bedeutung der Quadratwurzel von Kovarianz- / Präzisionsmatrizen

8

Sagen XRnist eine Zufallsvariable mit der Kovarianz . Einträge der Kovarianzmatrix sind per Definition Kovarianzen: Es ist auch bekannt, dass Einträge mit der Genauigkeit erfüllen: wobei die rechte Seite die Kovarianz von wobei von allen anderen Variablen abhängig ist.ΣRn×n

Σij=Cov(Xi,Xj).
Σ1
Σij1=Cov(Xi,Xj|{Xk}k=1nXi,Xj}),
XiXj

Gibt es eine statistische Interpretation der Einträge einer Quadratwurzel von oder ? Durch die Quadratwurzel aus einer quadratischen Matrix I : jeder Matrix , so dass . Eine Eigenwertzerlegung dieser Matrizen ergibt meines Erachtens keine solche eintragsbezogene Interpretation.ΣΣ1AMMtM=A

Yair Daon
quelle
1
@kjetilbhalvorsen Erklärung hinzugefügt. Ich habe jedoch ein paar Details unter die Lupe genommen. Wovon ist der Wert abhängig? Da keine angegeben ist, muss sie über alle Werte gemittelt werden.
Yair Daon
1
Der Bericht über Regression, Korrelation und bedingte Verteilungen unter stats.stackexchange.com/questions/71260/… enthält explizite geometrische Konstruktionen von zwei verschiedenen Quadratwurzeln der inversen Kovarianzmatrix. Diese geometrischen Ideen verallgemeinern sich auf höhere Dimensionen und liefern dadurch mindestens zwei unterschiedliche, bekannte statistische Interpretationen (nämlich PCA und multiple Regression).
whuber

Antworten:

1

Ich werde Matrixquadratwurzeln von als schreiben , um mit der Cholesky-Zerlegung übereinzustimmen, die als wobei lowtri (unteres Dreieck) ist. Sei also ein Zufallsvektor mit und . Sei nun ein Zufallsvektor mit Erwartungsnull und Einheitskovarianzmatrix.ΣΣ=AATΣ=LLTLXEX=μVarX=ΣZ

Beachten Sie, dass es (mit Ausnahme des Skalarfalls) unendlich viele Matrixquadratwurzeln gibt. Wenn wir einer der sein lassen , können wir alle anderen als wobei eine beliebige orthogonale Matrix ist, . Dies ist als einheitliche Freiheit der Quadratwurzeln bekannt .AAOOOOT=OTO=I

Betrachten wir einige bestimmte Quadratwurzeln der Matrix.

  1. Zuerst eine symmetrische Quadratwurzel. Verwenden Sie die spektrale Zerlegung, um zu schreiben . Dann ist und dies kann als PCA (Hauptkomponentenanalyse) von interpretiert werden .Σ=UΛUT=UΛ1/2(UΛ1/2)TΣ1/2=UΛ1/2Σ

  2. Die Cholesky-ZersetzungΣ=LLT undList lowtri. Wir können vertretenX wie X=μ+LZ. Multipliziert man, um skalare Gleichungen zu erhalten, erhält man ein DreieckssystemZ, was im Zeitreihenfall als MA-Darstellung (gleitender Durchschnitt) interpretiert werden kann.

  3. Der allgemeine Fall A=LOMit den obigen Angaben können wir dies nach dem Drehen als MA-Darstellung interpretieren Z.

kjetil b halvorsen
quelle
+1 Re Quadratwurzelberechnung durch spektrale Zerlegung, wie ich gesehen habe Σ1/2=UΛ1/2UTauch. Aber das wäre anders als das, was Sie erwähnt haben, aber beide sind gültig, stimmt das?
NULL
0

Eine nxn-Matrix kann, wie Sie erwähnen, viele Quadratwurzeln haben. Eine Kovarianzmatrix muss jedoch positiv semidefinit sein, und eine positive semidefinitive Matrix hat nur eine Quadratwurzel, die ebenfalls positiv semidefinit ist. Schauen Sie sich den Wikipedia-Artikel mit dem Titel "Quadratwurzel einer Matrix" an.

Michael R. Chernick
quelle
3
Ich habe nicht nach einer symmetrischen Quadratwurzel gefragt. Zum Beispiel reicht eine Cholesky-Zerlegung für den Zweck dieser Frage aus.
Yair Daon
1
Ich habe nicht gesagt, dass ein Begriff in der Kovarianzmatrix nicht negativ sein kann. Ich spreche von einer anderen Eigenschaft, dass eine Kovarianzmatrix zumindest positiv semidefinit ist.
Michael R. Chernick
0

Manchmal sind Menschen daran interessiert, die Positionen von Nullen in der Präzisionsmatrix aus demselben Grund zu schätzen, den Sie oben beschrieben haben. WennM ist Ihre Quadratwurzelmatrix, dh MM=Σ1, dann für Knoten ij

Σi,j1=0MiMj=0
Ich stelle mir also vor, dass Sie beim Betrachten des inneren Produkts zwischen den Spalten Ihrer geschätzten Quadratwurzelmatrix eine Zahl erhalten, die der Nähe zu bedingt unabhängigen zwei Knoten ähnelt. Nur eine Idee.

Die Quadratwurzel der Kovarianzmatrix, das ist die Skala. Ich stelle mir vor, einen normalen normalen Zufallsvektor zu simulieren und ihn dann mit der Quadratwurzelmatrix vormultiplizieren. Wenn diese Matrix ein unteres Dreieck hat, stelle ich mir immer vor, all die kleinen Multiplikationen und Additionen vorzunehmen.

Es gibt auch Einzelfälle, in denen bestimmte Elemente der Quadratwurzelmatrix nur Quadratwurzeln einzelner Elemente der Kovarianzmatrix sind. Das ist allerdings nicht so interessant, also denke ich, dass Sie bereits darüber nachgedacht haben.

Taylor
quelle