如何做主成分分析?
stats 包中的 princomp 函数。
(pc. cr <- princomp(USArrests , cor = TRUE))
plot(pc. cr , type = "lines" # or ” b a r p l o t ”
) # or s c r e e p l o t
loadings(pc. cr)
princomp() 中的参数 cor = TRUE 表示使用样本相关矩阵作主成分分析,反之使用样本协方差矩
阵。loadings() 返回因子荷载。screeplot() 绘制碎石图。