静粛に、只今統計勉強中

仕事でデータ分析をすることになったバリバリ文系アラフィフのおっさんが、独学で統計の勉強を始めました。

Excel VBAで『マンガでわかる統計学 因子分析編』の主成分分析の手順をマクロにしてみた3

前回に続き、『マンガでわかる統計学 因子分析編』で説明される分析の流れを、作ったマクロのソースコードで追っていきます。

主成分分析の流れ(承前)

Step6 各観測の主成分得点を求める。

f:id:cyclo-commuter:20180711125025p:plain
f:id:cyclo-commuter:20180711125440p:plain
こうなります。
本では後のステップになりますが、寄与率と累積寄与率もここで計算しています。
f:id:cyclo-commuter:20180711125532p:plain
f:id:cyclo-commuter:20180711125551p:plain

今回、『マンガでわかる統計学 因子分析編』を読んで勉強になった!
と実感できる瞬間がありました。

本の中では特に説明がなかったのに、主成分得点を標準化得点と固有ベクトルの行列の積で求めればよいのだ、と直感できたときです。

{\displaystyle \begin{pmatrix}-0.7 & 0.4 & 1.5 \\-1.4 & 1.3 & -1.4 \\1.4 & -0.5 & 0.8 \\-0.7 & -1.4 & 0.1 \\0.0 & 1.3 & 1.5 \\0.7 & -0.5 & -0.7 \\0.7 & 0.4 & 0.1 \\-1.4 & -1.4 & -1.4 \\0.0 & -0.5 & -0.7 \\1.4 & 1.3 & 0.1\end{pmatrix}\times\begin{pmatrix}0.57157 & -0.60447 \\0.52212 & 0.789607\\0.633055 & -0.10553\end{pmatrix}}*1

第3章のたいく…(ゲフンゲフン)行列の計算に真面目に取り組んだのは無駄じゃなかった~。 

Step7 主成分得点を元に散布図を描く。

f:id:cyclo-commuter:20180711140753p:plain
こうなります。
f:id:cyclo-commuter:20180711140820p:plain

*1:数値は、前々回からの記事に貼り付けた画像の見た目で表現しています。実際の値は見た目と異なります。