PythonのMCMC(Markov Chain Monte Carlo)ライブラリであるPyMC3は,試された方はご存知の通り,数値処理ライブラリTheanoをベースとして作られている.この組み合わせのおかげで,複雑な問題も扱える柔軟性と十分な計算処理能力持つと言われているが,いかんせんTheanoは素人には相当難しい.シンボルの扱いを先に定義し,後から数値をバインドして計算するコンセプトが理解を難しくさせている.ここでは,PyMC3 simulation作業で問題が生じたときにTheanoの機能に触らずににデバッグを行うための小技を紹介する. PyMC3 simulation作業の流れ 一般的なsimulationの作業の流れは次のようなものである. 1. データセットの用意(ファイル読み込み). 2. シミュレーションモデルの記述.(事前分布はこれとこれ,尤度関数はこれ,など.) 3
![Theano素人がPyMC3で頑張るための小技 - Qiita](https://cdn-ak-scissors.b.st-hatena.com/image/square/5215aa2435239212fe38f090bf93365d32bf0360/height=288;version=1;width=512/https%3A%2F%2Fqiita-user-contents.imgix.net%2Fhttps%253A%252F%252Fcdn.qiita.com%252Fassets%252Fpublic%252Farticle-ogp-background-9f5428127621718a910c8b63951390ad.png%3Fixlib%3Drb-4.0.0%26w%3D1200%26mark64%3DaHR0cHM6Ly9xaWl0YS11c2VyLWNvbnRlbnRzLmltZ2l4Lm5ldC9-dGV4dD9peGxpYj1yYi00LjAuMCZ3PTkxNiZoPTMzNiZ0eHQ9VGhlYW5vJUU3JUI0JUEwJUU0JUJBJUJBJUUzJTgxJThDUHlNQzMlRTMlODElQTclRTklQTAlOTElRTUlQkMlQjUlRTMlODIlOEIlRTMlODElOUYlRTMlODIlODElRTMlODElQUUlRTUlQjAlOEYlRTYlOEElODAmdHh0LWNvbG9yPSUyMzIxMjEyMSZ0eHQtZm9udD1IaXJhZ2lubyUyMFNhbnMlMjBXNiZ0eHQtc2l6ZT01NiZ0eHQtY2xpcD1lbGxpcHNpcyZ0eHQtYWxpZ249bGVmdCUyQ3RvcCZzPTE5NGYzZGQ4NzE1YzQwYzA1MjRkYWVlNWM3ZDg4ZjM2%26mark-x%3D142%26mark-y%3D112%26blend64%3DaHR0cHM6Ly9xaWl0YS11c2VyLWNvbnRlbnRzLmltZ2l4Lm5ldC9-dGV4dD9peGxpYj1yYi00LjAuMCZ3PTYxNiZ0eHQ9JTQwVG9tb2tJc2hpaSZ0eHQtY29sb3I9JTIzMjEyMTIxJnR4dC1mb250PUhpcmFnaW5vJTIwU2FucyUyMFc2JnR4dC1zaXplPTM2JnR4dC1hbGlnbj1sZWZ0JTJDdG9wJnM9OTY0M2M4ZTQzZGJjMjFkNTY2YzFhYjEwZmY0ZDUxMWU%26blend-x%3D142%26blend-y%3D491%26blend-mode%3Dnormal%26s%3D9326834ea1e7206494ae8f0a023dd3d1)