ICAからモデルベースのクラスタリング
昨日気になったので、とりあえずどんな感じか確かめようと。
> library("fastICA")
> library("mclust")
> data("iris")
> b <- fastICA(iris[1:4],2)
> mc <- Mclust(b$S,G=3,modelNames="VVV")
> plot(mc,b$S)
混合分布はクラスタリングの本にも載っているが、PRMLの下巻のほうがK-meansからEM,混合モデルの導入の流れがスムーズで分かりやすかったかな。
Rで学ぶクラスタ解析
パターン認識と機械学習 下 - ベイズ理論による統計的予測
