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结果

接下来我们利用碎石图来客观表现我们的主成分分析结果

screeplot (PCA, type="lines") 
# 画碎石图，用直线图类型
abline (v = 2 ,  
        col = "#900021",
        lwd = 2,
        lty = 2)
#添加划分