Sagen ist 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.
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.
Antworten:
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 Σ=LLT L X EX=μ 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 .A AO O OOT=OTO=I
Betrachten wir einige bestimmte Quadratwurzeln der Matrix.
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 Σ
Die Cholesky-ZersetzungΣ=LLT undL ist 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.
Der allgemeine FallA=LO Mit den obigen Angaben können wir dies nach dem Drehen als MA-Darstellung interpretieren Z .
quelle
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.
quelle
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 M′M=Σ−1 , dann für Knoten i≠j
Σ−1i,j=0⟺M′iMj=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.
quelle