昨天跟今天在研究主成分分析,想要做一个主成分回归的模型,在这过程中,发现一些细节的东西,稍微不注意就出错的,请大家指导!
- rm( list = ls ( all = TRUE))
- shuju<-read.table("C:/Users/Administrator/Desktop/sj1.txt",header=T);
- data<-as.matrix(shuju);
- ##先求矩阵标准化,再求相关系数,最后主成分分析。
- data1<-scale(data[,1:3],scale=T)
- data3<-cor(data1)
- pca.pr<-princomp(data3)
- summary(pca.pr,loadings=TRUE)###结果与princomp()不同
- ##应该用covmat参数直接指定协方差矩阵
- pca.pr<-princomp(covmat=data3)
- summary(pca.pr,loadings=TRUE)
- #直接用数据来分析
- data.princomp<-princomp(data[,1:3],cor=T);
- summary(data.princomp,loadings=T);
复制代码数据
- X1 X2 X3 Y
- 1 149.3 4.2 108.1 15.9
- 2 161.2 4.1 114.8 16.4
- 3 171.5 3.1 123.2 19.0
- 4 175.5 3.1 126.9 19.1
- 5 180.5 1.1 132.1 18.8
- 6 190.7 2.2 137.7 20.4
- 7 202.1 2.1 146.0 22.7
- 8 212.4 5.6 154.1 26.5
- 9 226.1 5.0 162.3 28.1
- 10 231.9 5.1 164.3 27.6
- 11 239.0 0.7 167.6 26.3
复制代码