習作；混合正規分布モデル(Gaussian Mixture Model) by rstan

習作；混合正規分布モデル(Gaussian Mixture Model) by rstan

Stanのマニュアルに載っている簡単なモデルなんだけど，最後の「各要素の所属確率」を算出するのになぜ...概要を表示 Stanのマニュアルに載っている簡単なモデルなんだけど，最後の「各要素の所属確率」を算出するのになぜか手間取った。softmax関数を使うときの型の問題でした。ちぇ。 Stanコードはこの通り。generated quantitiesのところでsoftmax関数を使います。 [code] data{ int<lower =2> K; #number of clusters int<lower =1> N; #number of observations real X[N]; #observed data } parameters{ vector[K] mu; vector<lower =0,upper=10>[K] sig2; simplex[K] theta; } transformed parameters{ vector[K] ps[N]; for(n in 1:N){ for(k