hei3d <- subset(myakari3d0,V3 == 0) hei2d <- hei3d[ ,c(1, 2)] plot(hei2d) |
library(mclust) |
mykmakari9 <- kmeans(hei2d, 9, nstart=50) clPairs(hei2d, cl=mykmakari9$cluster) |
mykmakari17 <- kmeans(hei2d, 17, nstart=50) clPairs(hei2d, cl=mykmakari17$cluster) |
plot(hei2d, type="n") text(mykmakari17$centers) rect(0, 0, 256, -256, border = "red") |