1.主成分分析
princomp()
princomp(formula, data = NULL, subset, na.action, ...)
princomp(x, cor = FALSE, scores = TRUE, covmat = NULL, subset = rep(TRUE, nrow(as.matrix(x))), ...)
x为数据集;cor默认为FALSE,表示使用样本的协方差阵作主成分分析,若为TRUE表示使用相关系数矩阵;如果不用x指定数据集,可使用covmat指定协方差阵。
prcomp()
prcomp(x, retx = TRUE, center = TRUE, scale. = FALSE, tol = NULL,...)
retx表示是否旋转数据,即中心化和标准化,默认为TRUE;center表示中心化;scale.表示用相关系数矩阵计算主成分,默认使用协方差矩阵。
绘制碎石图:由陡峭变平坦时,选择最佳主成分个数
screeplot()
screeplot(x, npcs = min(10, length(x$sdev)), type &#