Old Faithful + variational_gaussian_mixture

http://d.hatena.ne.jp/n_shuyo/20100412/prml で話題になってるのを実装してみた。
ただRの奴をそのまま移植しても勉強にならないので、わからなくなるまでpythonで書き続けた。

結果縮退はしてるけど、潰れすぎてる感もある。というか必ずある点に来てる気がする。入力値を正規化してないのが問題かしら？

あとディガンマ関数がわからなくて結局polygamma(0,x)を使ったけど、あとでしましまさんにscipy.special.psiがあることを教えてもらった。psiなのか・・

(mk)
[1.9517697309034685e-69, 3.9674122478717695e-68]
[1.9517697309034685e-69, 3.9674122478717695e-68]
[3.3374566183149512, 67.84133422396188]
[1.9517697309034685e-69, 3.9674122478717695e-68]
[3.336836379708978, 67.828726475194429]
[3.194198653198653, 64.929292929292927]