「PAL8-CNV 〜ぱるはちこんぶ〜」より
「XYZ」のほうのTSVを読み込む | myakari3d0=read.table("clipboard",h=0) |
| -80〜20mの地点のみ | myakari3dc <- subset(myakari3d0,V3< 30) |
…しすてむこーる・
…いんすとおーる！ | install.packages("mclust")
# これは最初に1回だけ必要な操作です |
| どーん | mycmpakari3dc <- prcomp(myakari3dc, scale=TRUE)
mykmakaricbd3c3 <- kmeans(mycmpakari3dc$x, 3, nstart=50)
library(mclust)
clPairs(myakari3dc, cl=mykmakaricbd3c3$cluster)
# プロットされた図を右クリックしてコピーや保存をします |
| ゆんゆん | clPairs(mycmpakari3dc$x, cl=mykmakaricbd3c3$cluster)
# プロットされた図を右クリックしてコピーや保存をします |
| ばーん | original <- mykmakaricbd3c3$centers %*% t(mycmpakari3dc$rotation)
original <- scale(original, center = FALSE, scale = 1 / mycmpakari3dc$scale)
original <- scale(original, center = -mycmpakari3dc$center, scale = FALSE)
myakaricbd3centers3co <- as.table(original)
write.table(myakaricbd3centers3co, file="clipboard", sep="\t") |
| すたぱーん | myakari2d0 <- myakari3d0[ ,c(1, 2)]
plot(myakari2d0, type="n")
text(myakaricbd3centers3co)
rect(0, 0, 256, -256, border = "red")
# プロットされた図を右クリックしてコピーや保存をします |
| …すちゃ。 | table(mykmakaricbd3c3$cluster)
table(myakari3dc$V3, mykmakaricbd3c3$cluster)
# 表示された集計表を見て内容を確かめます
write.table(table(myakari3dc$V3, mykmakaricbd3c3$cluster), file="clipboard", sep="\t")
# 表をクリップボードにコピーします |
