我有一个大数据集“nm”:
dim(nm)
[1] 3871012 13
我正在尝试在该数据上绘制自组织地图(SOM):
library(kohonen)
som_grid = somgrid(xdim = 20, ydim=20 , topo="hexagonal")
zx = as.matrix(nm)
som_model = som(zx,
grid=som_grid,
rlen=100,
alpha=c(0.05,0.01),
keep.data = TRUE)
som_cluster = cutree(hclust(dist(as.numeric(unlist(som_model$codes)))), 4)
plot(som_model, type="codes", bgcol=rainbow(4)[som_cluster])
这个代码工作正常,但是当我添加群集边界时会出现以下错误:
add.cluster.boundaries(som_model, som_cluster)
Error in nb[i, 1] : subscript out of bounds
我不确定如何纠正这个错误或我做错了什么。
请提供解决方案。谢谢!