主成分分析碎石图_科学网—主成分分析之R篇 - 林元震的博文

R中psych包可以进行主成分分析，其分析的步骤为：

(1) 判断主成分的个数；

(2) 提取主成分；

(3) 获取主成分得分；

(4) 列出主成分方程，解释主成分意义。

【例子】 测定了20株杨树树叶，每个叶子测定了4个变量(叶长x1，2/3处叶宽x2，1/3处叶宽x2，1/2处叶宽x2)，测定结果如表4-52所示。试进行本样本的主成分分析。

psych包中的fa.parallel()函数可以判断主成分的个数，其使用格式为：

fa.parallel(x, fa = , n.iter =)

其中，x为待研究的数据集或相关系数矩阵，fa为主成分分析(fa= "pc")或者因子分析(fa = "fa")，n.iter指定随机数据模拟的平行分析的次数。分析代码如下：

运行上述代码，得到结果如下：

上图中，直线与x符号链接的曲线为碎石图，1.0水平线为1准则的特征值，虚线为100次随机数据模拟的平行分析。碎石图画出了特征值与主成分分数的图形。结果表明，选择2个主成分即可保留样本中的大量分信息。

第二步，提取主成分。psych包中的principal( )函数可以根据原始数据或相关系数矩阵做主成分分析，其使用格式为：

principal(x, nfactors =, rotate