「PAL8-CNV 〜ぱるはちこんぶ〜」より
「XYZ」のほうのTSVを読み込む | myakari3d0=read.table("clipboard",h=0) |
| 平地を除く | myakari3d <- rbind(subset(myakari3d0,V3> 0),subset(myakari3d0,V3< -3)) |
…しすてむこーる・
…いんすとおーる！ | install.packages("mclust")
install.packages("ggvoronoi")
# これは最初に1回だけ必要な操作です |
| どーん | mycmpakari3d <- prcomp(myakari3d, scale=TRUE)
mykmakari9c3 <- kmeans(mycmpakari3d$x, 9, nstart=50)
library(mclust)
clPairs(myakari3d, cl=mykmakari9c3$cluster)
# プロットされた図を右クリックしてコピーや保存をします |
| ばーん | original <- mykmakari9c3$centers %*% t(mycmpakari3d$rotation)
original <- scale(original, center = FALSE, scale = 1 / mycmpakari3d$scale)
original <- scale(original, center = -mycmpakari3d$center, scale = FALSE)
myakaricenters9co <- as.table(original)
write.table(myakaricenters9co, file="clipboard", sep="\t") |
| x,yの座標のリストを読み込む | myakaridd=read.table("clipboard",h=0)
names(myakaridd) <- c("x", "y") |
| |
| ボロノイ図を描く | library(ggvoronoi)
ggplot(myakaridd,aes(myakaridd$x, myakaridd$y)) + stat_voronoi(geom="pa |
