分散が最大の方向を見つける 主成分分析は、データの雲が最も大きくばらつく1つの方向を求め、次にそれと直角な次の方向を求め、というように続けます。N N N 個の資産のリターンの雲は N N N 次元にあって見づらいものの、主方向は独立なリスク源を順位付けします。手順は短く、各列を中心化し、共分散行列を作り、固有ベクトルを固有値順に並べ、先頭のものにデータを射影します。
Σ = 1 T − 1 X c T X c = Q Λ Q T , Z = X c Q , Var ( z j ) = λ j , s j = λ j ∑ i λ i \Sigma=\frac{1}{T-1}X_c^{T}X_c=Q\Lambda Q^{T},\qquad Z=X_cQ,\qquad \operatorname{Var}(z_j)=\lambda_j,\qquad s_j=\frac{\lambda_j}{\sum_i\lambda_i} Σ = T − 1 1 X c T X c = Q Λ Q T , Z = X c Q , Var ( z j ) = λ j , s j = ∑ i λ i λ j (1) ここで X c X_c X c は中心化したリターンの T × N T\times N T × N 行列、Q Q Q の列は主方向(ローディング)、Z Z Z の列はスコア、つまり新しい座標でのデータです。スコアは無相関で、j j j 番目の分散は λ j \lambda_j λ j なので、固有値は全分散 tr Σ \operatorname{tr}\Sigma tr Σ をちょうど分け合います。方向の符号は任意で、q \mathbf q q と − q -\mathbf q − q は同じ直線です。注意が2つあります。中心化が重要で、しないと第1方向の一部が平均を向きます。単位も重要で、共分散行列のPCAは最も変動の大きい変数に支配されるため、尺度が異なる変数は先に標準化します。つまり相関行列で行います。
2資産を手計算で 2つの資産の共分散行列が
Σ = ( 4 2 2 3 ) \Sigma=\begin{pmatrix}4&2\\2&3\end{pmatrix} Σ = ( 4 2 2 3 ) だとします。トレースは7、行列式は8なので
λ 2 − 7 λ + 8 = 0 \lambda^2-7\lambda+8=0 λ 2 − 7 λ + 8 = 0 、
λ = ( 7 ± 17 ) / 2 \lambda=(7\pm\sqrt{17})/2 λ = ( 7 ± 17 ) /2 、すなわち
5.562 5.562 5.562 と
1.438 1.438 1.438 です。第1成分は全分散の
5.562 / 7 = 79.5 % 5.562/7=79.5\% 5.562/7 = 79.5% を説明します。方向は
( 4 − 5.562 ) q 1 + 2 q 2 = 0 (4-5.562)q_1+2q_2=0 ( 4 − 5.562 ) q 1 + 2 q 2 = 0 を解いて符号を除き
q 1 ≈ ( 0.788 , 0.615 ) \mathbf q_1\approx(0.788,0.615) q 1 ≈ ( 0.788 , 0.615 ) 。中心化した観測値
( 1 , 1 ) (1,1) ( 1 , 1 ) の第1スコアは絶対値で
0.788 + 0.615 = 1.403 0.788+0.615=1.403 0.788 + 0.615 = 1.403 です。両方のローディングが同符号なので、この成分は2資産の加重平均です。
手順 :中心化、共分散、固有分解、固有値順に並べ、上位の方向に射影。スコア :無相関で分散は λ j \lambda_j λ j 、寄与率は λ j / ∑ λ i \lambda_j/\sum\lambda_i λ j / ∑ λ i 。注意 :符号は任意。データを中心化。単位が違うなら標準化(相関行列)。実際の株式リターンが語ること 2010年以降の米大型株12銘柄の日次リターン(4,184日)では、相関行列の第1主成分が分散の約40%を説明し、すべての銘柄が同符号で負荷します。これが市場です。第2成分は約12%で、銘柄をテクノロジーとディフェンシブのグループに分けます。残りの10個が残り48%を分け合い、ノイズとの区別がずっと難しくなります。これが次のレッスンの主題です。
ラボ:採点される内容 PCAの手順 15点 第1固有値 15点 寄与率 15点 スコア 15点 標準化すべき場面 15点 スコアの事実 10点 なぜデータを中心化するのか 15点