一挙放送
寝かせないぞっ☆ | myakari3d0=read.table("clipboard",h=0) # XYZのほうのTSVデータ
myakari3d <- rbind(subset(myakari3d0,V3> 10),subset(myakari3d0,V3< -10))
myakari3dg <- subset(myakari3d0,V3 == 0)
myakari3dyama <- subset(myakari3d0,V3> 10)
myakari3dumi <- subset(myakari3d0,V3< -10)
myakari2dg <- myakari3dg[ ,c(1, 2)]
colnames(myakari2dg) <- c("x", "y")
mykmakari3dumi100 <- kmeans(myakari3dumi, 100, nstart=50)
myakari2dumi100 <- as.data.frame(mykmakari3dumi100$centers[, c(1,2)])
plot(myakari2dumi100)
# プロットされた図を右クリックしてコピーや保存をします
mykmakari3dyama100 <- kmeans(myakari3dyama, 100, nstart=50)
myakari2dyama100 <- as.data.frame(mykmakari3dyama100$centers[, c(1,2)])
plot(myakari2dyama100)
# プロットされた図を右クリックしてコピーや保存をします
myfunc5 <- function(x, y) cbind(x, y, myakari2dumi100[which.min(head(dist(rbind(c(x, y), myakari2dumi100), diag=FALSE, upper=FALSE), n=nrow(myakari2dumi100))), ])
plot(0, 0, xlim=c(0, 255), ylim=c(-255, 0), type="n", xlab="x", ylab="y")
for (i in 1:nrow(myakari2dyama100)) {
g <- myfunc5(myakari2dyama100[i, 1], myakari2dyama100[i, 2])
arrows(g[,1], g[,2], g[,3], g[,4], length = 0.05, col = "deepskyblue")
}
rect(0, -255, 255, 0, border = "blue")
# プロットされた図を右クリックしてコピーや保存をします |
