library(psych)
mydata <- iris[,1:4]
# 计算最优主成分个数
mydata.screePlotsModel <- fa.parallel(mydata, fa="pc", n.iter=100, show.legend=F, main="崖底碎石图");
# 主成分个数
n=mydata.screePlotsModel$ncomp;
# 主成分分析函数 模型
mydata.pr <- princomp(mydata, scores=T, cor = TRUE)
#勾画碎石头
screeplot(mydata.pr,type="line",mian="碎石图",lwd=2)
summary(mydata.pr,loadings = TRUE)
mydata <- iris[,1:4]
# 计算最优主成分个数
mydata.screePlotsModel <- fa.parallel(mydata, fa="pc", n.iter=100, show.legend=F, main="崖底碎石图");
# 主成分个数
n=mydata.screePlotsModel$ncomp;
# 主成分分析函数 模型
mydata.pr <- princomp(mydata, scores=T, cor = TRUE)
#勾画碎石头
screeplot(mydata.pr,type="line",mian="碎石图",lwd=2)
summary(mydata.pr,loadings = TRUE)