#pymc3
計算された結果である確率分布$pp$からの値を取り出したいという問題。
サンプリング結果は$trace$に含まれるので、例えば$$pm.getvalue(trace['pp'], x)$$みたいに取り出せるものを探していたんだけど、まだ見つからない。
グラフ表示はされているので何か見落としているハズだと思いAPIリファレンスを眺めてみたのだが、見つからない。
仕方がないので、グラフ表示部分のソースコードを見てみたら、例えば離散分布のポアソン分布の場合、サンプリングされた結果からMathPlotのhistを使って表示していた。連続密度分布(例えば正規分布等)でもサンプリングされた結果から曲線を再構成しているようだった。内部的に!
何と言う事でしょう!
確率分布関数化するAPIは用意してないのかな?
そのような用途は無いの???
それとも、あまりに自明すぎるAPIがあるのか???
StackOverflowでpymc3関連のものを探っても、このような質問等はないみたいだし。。。
この件は引き続き調査を継続します。なるべくなら自前で作りたくないのでね(笑)