2016-04-28 10 views
0

正規化定数がわからない分布からランダム偏差を抜き出したいと思います。この分布は、未知の形状およびスケール= 1のガンマ尤度の前の共役である。 PDFは、このPDFファイルの不可欠のための閉じた形はありませんScipy Stats rv_continuous正規化定数なし

enter image description here

として、ここでhttps://en.wikipedia.org/wiki/Conjugate_prior#Table_of_conjugate_distributionsウィキペディアに与えられています。ここでscipy.stats.rv_continuousを使用していますか?そうでない場合は、任意の人が任意のpdfから逸脱を描くために推奨できるものは何ですか?

答えて

1

数値正規化(scipy.integrate.quad)は常にオプションです。

+0

このチップをありがとう!少しはそのスピードを心配した。ウィキペディアによれば、いくつかのアルゴリズムを使用して正規化定数をサンプルすることを実際に知る必要はありません。私は整数を数値的に計算するのを避ける方が速いと思います。 https://en.wikipedia.org/wiki/Rejection_sampling –

+0

はい、これもオプションです:-)。あなたは良いメジャントが必要です。 –