Ich neige dazu zu hören, dass normalerweise 3 größte Eigenwerte am wichtigsten sind, während diejenigen nahe Null Rauschen sind
Sie können das testen. Weitere Informationen finden Sie in dem in diesem Beitrag verlinkten Artikel. Wenn Sie sich mit finanziellen Zeitreihen befassen, möchten Sie möglicherweise zuerst die Leptokurtizität korrigieren (dh berücksichtigen Sie die Reihe der garchbereinigten Renditen, nicht die Roherträge).
Ich habe einige Forschungsarbeiten gesehen, in denen untersucht wurde, wie sich natürlich vorkommende Eigenwertverteilungen von denen unterscheiden, die aus zufälligen Korrelationsmatrizen berechnet wurden (wiederum Rauschen von Signal unterscheiden).
Edward:> Normalerweise würde man es anders herum machen: Sehen Sie sich die multivariate Verteilung der Eigenwerte (der Korrelationsmatrizen) an, die aus der gewünschten Anwendung stammen. Sobald Sie einen glaubwürdigen Kandidaten für die Verteilung von Eigenwerten identifiziert haben, sollte es ziemlich einfach sein, daraus zu generieren.
Das beste Verfahren zum Identifizieren der multivariaten Verteilung Ihrer Eigenwerte hängt davon ab, wie viele Assets Sie gleichzeitig berücksichtigen möchten (dh wie groß ist Ihre Korrelationsmatrix?). Es gibt einen ordentlichen Trick, wenn ( ist die Anzahl der Assets).pp ≤ 10p
Bearbeiten (Kommentare von Shabbychef)
Vier-Schritte-Verfahren:
- Angenommen, Sie haben Unterproben multivariater Daten. Sie benötigen für jede Teilstichprobe einen Schätzer der Varianz-Kovarianz-Matrix (Sie können den klassischen Schätzer oder eine robuste Alternative wie die schnelle MCD verwenden , die in matlab, SAS, S, gut implementiert ist). R, ...). Wenn Sie sich mit finanziellen Zeitreihen befassen, möchten Sie wie üblich die Reihe der garchbereinigten Renditen berücksichtigen, nicht die Roherträge.˜ C j jj = 1 , . . . , J.C.~jj
- Berechnen Sie für jede Unterprobe , ..., die Eigenwerte von .˜ Λ j = log ( ˜ λ j 1 ) log ( ˜ λ j p ) ˜ C jjΛ~j= Log( λ~j1)Log( λ~jp)C.~j
- Berechnen Sie , die konvexe Hülle der Matrix, deren j-ter Eintrag (dies ist wiederum in Matlab, R, ... gut implementiert). .J × p ˜ Λ jC.V.( Λ~)J.× pΛ~j
- Zeichnen Sie zufällig Punkte aus dem (dies erfolgt durch Gewichtung für jede Kante des wobei , wobei ein Draw aus einer Exponentialverteilungseinheit ist (weitere Details hier ).w i C V ( ~ Λ ) w i = γ iC.V.( Λ~)wichC.V.( Λ~) γiwich= γich∑pi = 1γichγich
Eine Einschränkung besteht darin, dass die schnelle Berechnung der konvexen Hülle einer Reihe von Punkten extrem langsam wird, wenn die Anzahl der Dimensionen größer als 10 ist.J.≥ 2