主成分分析应用之聚类

数据的读取，通过嵌套的for循环把R矩阵填满，即得到一个对称矩阵。

 

princomp函数covmat参数的含义，通过R语言的帮助文档可知。

 

 

 

当输入矩阵为协方差矩阵，可以通过此参数的输入，从而使princomp函数能求出其主成分。

 

Loadings(pr），为主成分分析输出结果中的负荷因子矩阵

 

 

从summmary(pr),结果可知，当主成分个数m=3时，累计贡献率为70.52%，大于70%。即取三个主成分。

 

从而求得第k个主成分在上的因子负荷量矩阵x，按照最近距离进行聚类。

首先通过：d <- dist(scale(X)) 中心标准化数据之后，求得距离矩阵。

 

hc1 <- hclust(d,method = "single") ，以最近距离，进行聚类

 

 

 

plot(hc1,hang= -1) 输出聚类图

 

 

 

结合聚类分析的一般还有以下函数

 

 

 

 

本人大三，不是985，不是211，简单普通的本科生，这些也是最基础的知识。勿喷，谢谢。还有对于这些函数，最好就先把主要的知识点原理搞明白。我大一大二总觉得学原理没用，不愿花时间。现在错了。会命令不懂原理，就好像你得到一个爱的人，得不到她的心，除非你只是玩玩。