| R は、自由なソフトウェアであり、「完全に無保証」です。 |
| 一定の条件に従えば、自由にこれを再配布することができます。 |
| 配布条件の詳細に関しては、'license()' あるいは 'licence()' と入力してください。 |
# 「xyz」のほうのTSVデータをクリップボードにコピーしておきます
myakari3d0=read.table("clipboard",h=0) |
|
myakari3dn0 <- subset(myakari3d0,V3 != 0)
mykmakari17edge <- kmeans(myakari3dn0, 17, nstart=50)
library(mclust)
clPairs(myakari3dn0, cl=mykmakari17edge$cluster)
# プロットされた図を右クリックしてコピーや保存をします |
mycmpakari17edge <- prcomp(myakari3dn0, scale=TRUE)
clPairs(mycmpakari17edge$x, cl=mykmakari17edge$cluster)
# プロットされた図を右クリックしてコピーや保存をします |
mykmakari9n0c3 <- kmeans(mycmpakari17edge$x, 9, nstart=50)
clPairs(myakari3dn0, cl=mykmakari9n0c3$cluster)
# プロットされた図を右クリックしてコピーや保存をします |
original <- mykmakari9n0c3$centers %*% t(mycmpakari17edge$rotation)
original <- scale(original, center = FALSE, scale = 1 / mycmpakari17edge$scale)
original <- scale(original, center = -mycmpakari17edge$center, scale = FALSE)
clPairs(myakari3dn0[, c(1, 2)], cl=mykmakari9n0c3$cluster)
# または plot(myakari3dn0[, c(1, 2)], type="n")
text(original)
rect(0, 0, 256, -256, border = "red")
# プロットされた図を右クリックしてコピーや保存をします |
