あたらしい 新逢坂山田前 だぞっ☆ | library(nnclust) myakari2d0 <- myakari3d0[ ,c(1, 2)] mycmpakari0elv <- prcomp(myakari3dc, scale=TRUE) mykmakari90elvc3 <- kmeans(mycmpakari0elv$x, as.integer(nrow(myakari3dc) / 41625 * 90), nstart=50) original <- mykmakari90elvc3$centers %*% t(mycmpakari0elv$rotation) original <- scale(original, center = FALSE, scale = 1 / mycmpakari0elv$scale) original <- scale(original, center = -mycmpakari0elv$center, scale = FALSE) mykmakari90elvc3o <- as.table(original) mykmakari90to30 <- kmeans(original, as.integer(nrow(myakari3dc) / 41625 * 30), nstart=50) #write.table(mykmakari90elvc3o, file="clipboard", sep="\t") #write.table(mykmakari90to30$centers, file="clipboard", sep="\t")
# いくつ? nrow(mykmakari90elvc3o) nrow(mykmakari90to30$centers) |
154レール! (1341じゃないよ154だよ) | par(cex=0.8) # 文字を小さく Y <- matrix( mykmakari90to30$centers, ncol=3 ) mst <- mst(Y) plot(myakari2d0, type="n") segments( Y[mst$from,1], Y[mst$from,2], Y[mst$to,1], Y[mst$to,2], col="green") text(mykmakari90to30$centers) rect(0, 0, 256, -256, border = "red") # プロットされた図を右クリックしてコピーや保存をします |
472バス! (5931じゃないよ472だよ) | par(cex=0.6) # 文字を小さく X <- matrix( mykmakari90elvc3o, ncol=3 ) mst <- mst(X) plot(myakari2d0, type="n") text(mykmakari90elvc3o) segments( X[mst$from,1], X[mst$from,2], X[mst$to,1], X[mst$to,2], col="black") rect(0, 0, 256, -256, border = "red") # プロットされた図を右クリックしてコピーや保存をします |