とある技術者の徒然草

生産技術者の適当な日記(統計言語Rに関するメモがメイン)

【R言語】x-means法でクラスタリング

x-means法でクラスタリング

k-means法はクラスター数を指定する必要がある。
クラスター数を指定しないx-means法を下記を参考(マルパクリ)に試してみた。

aaaazzzz036.hatenablog.com

#http://aaaazzzz036.hatenablog.com/entry/2013/11/27/210355
#s-kmeans.R xmeans_sub.r xmeans.r mkmeans.R をDL
#スクリプトとして作業ディレクトリに保存

source("xmeans_sub.r")
source("xmeans.r")
source("s-kmeans.R")

#x-meansでクラスタリングする
xkm  <- xmeans$new ()
xkm$ishioka_xmeans (tsne_plot) #-tsne_plotはデータフレーム
tsne_plotX <- cbind(tsne_plot,xkm$Cluster)
tsne_plotX <- rename(tsne_plotX,"class"="xkm$Cluster")
tsne_plotX$class <- as.factor(tsne_plotX$class) #ファクター化する

mkm_tsneplot_X <- ggplot(tsne_plotX) + geom_point(aes(x=x, y=y,col=class))


f:id:M_taka072:20190210141804p:plain
xmeans

62までクラスター数が分かれた。
うまく分類できていない。。。