R语言实现PCA
一次来自男神MR.来的作业
PCA
读取数据
#读取数据,数据名称为“PCA上机数据”,且对数据进行整理
setwd("/Users/auzzer_pang")
data_PCA = read.table("PCA上机数据.txt")
data = data_PCA[,1:7]
names(data) =
c('100','200','400','800','1500','3000','Marathon')
row.names(data) = data_PCA[,8]
PCA主体部分
#由于单位不同所以先对数据标准化
cor(data)
#利用princomp函数进行主成分分析
PCA=princomp(data, cor=T)
#cor = T这个选项可以帮助我们吧矩阵变为协方差阵
summary(PCA, loadings=T)
输出PCA的主要结果,列出主成分分析分析结果
利用碎石图来展示PCA结果
接下来我们利用碎石图来客观表现我们的主成分分析结果