naoya_t@hatenablog

いわゆるチラシノウラであります

PRML復々習レーン#14（再）

PRML Python NumPy

http://connpass.com/event/3529/

前回、ワルプルギスの夜が来た為に延期（というか嵐のハッカソン）になった復々習レーン。
今日は９章。k-meansと混合ガウスとEMアルゴリズム。

sleepy_yoshiさんがsklearnとか使って10行ぐらいでさくっとコードを書いている間に、matplotlibと格闘しながら混合ガウス分布のEMアルゴリズムをごりごり実装しておりました。
図9.8（下巻p.153）的なものの再現です。

データセットはold faithful間欠泉データです。

実装はgithubに上げてあります：
https://github.com/naoyat/PRMLrevenge/tree/5632b99f79c7833c99faf3257cefaccb9a84e204/chap9

次回（11/9 or 11/23）は9章の残りと10章を読みます。
10章（近似推論法）はみんなでextreme readingします。頑張ろう！

おまけ