玻璃数据的r-谱聚类

Question

我正在尝试在UCI存储库的玻璃数据集上进行光谱聚类。如果你知道玻璃数据有9个属性-- "RI"，"Na"，"Mg"，"Al"，"Si"，"K"，"Ca"，"Ba"和"Fe"。我使用的"kknn"包具有谱聚类"specClust"，如下所示：

library(kknn)
data(glass)
specClust(glass[, 1:9], centers=2, nn=50, iter.max=30)  

但是，我得到了两个属性的集群方法，如下所示：

     [,1]       [,2]
1  0.9096567 -0.3331058
2  0.5101274  0.8273862

它应该返回所有9个属性的集群均值，而不仅仅是2个属性。当我运行K-均值集群时，R console将返回9个属性的集群均值，如下所示：

       RI      Na       Mg       Al       Si         K        Ca     Ba
1  1.519226 13.7124 0.210400 1.837800 72.87120 0.4910000 10.160400 0.56580000
2  1.518103 13.3150 3.438841 1.325122 72.58378 0.4989024  8.590061 0.05591463
      Fe
1  0.0480000
2  0.0597561

我想知道为什么"kknn"包中的谱聚类只返回两个属性的集群方法。

Answer 1

您看了两个集群的总体统计数据。

这应该会给你你想要的东西，还有更多：

pacman::p_load(kknn)
data(glass)
cl   <- specClust(glass[, 1:9], centers=2, nn=50, iter.max=30) 
pcol <- as.character(glass$Type)
pairs(glass[2:10], pch = pcol, col = c("green", "red")[cl$cluster])
table(glass$Type, cl$cluster)

aggregate(glass[,2:10],by=list(cl$cluster),mean)

组1. RI Na Mg Al Si K Ca Ba Fe 1 1.519061 13.81638 0.9115942 0.9115942 1.746232 0.4652174 9.729130 0.534202899 0.04318841 2 1.518035 13.21345 3.5282069 1.301517 72.6460 0.5122069 8.589517 0.004137931 0.06358621

​

​