因子分析步骤

1) 输入原变量x的观测数据阵并对其标准化，记之为X；

2) 求x的子样相关阵 $R=\frac{1}{N}{X}^{\prime }X$ ；

3) 求R的特征值 ${\lambda }_1\ge \cdots \ge {\lambda }_m$ 及相应的单位特征向量 $u_1,\cdots ,u_m$ ；

4) 确定公共因子的个数k值；

5) 计算初始因子载荷阵 $A=\left(a_1,\cdots ,a_k\right)$ ；

6) 计算共同度 $H=\left(h_1^2,\cdots ,h_m^2\right)$ ， $h_i^2={\displaystyle \sum a_{ij}^2}$ ， $i=1,\cdots ,m$ ；

7) 求正交(旋转)因子解—因子的正交旋转；

8) 计算因子得分，如可求 $R^{-1}$ 从而得因子得分矩阵。