4

問題

2つの変数に依存する確率密度関数を計算しました。この多変量分布を使用して、PDFに比例する確率で発生する乱数を生成したいと思います。

どうやら、SciPyは現在単変量分布のみをサポートしているようです。2D配布を可能にする簡単な方法や使いやすいパッケージはありますか?

回避策として、関心のあるドメインで乱数を作成して破棄するか、PDFに関連する可能性を残しておくようにすることができますが、それでも他のオプションがある可能性があります。乱数の生成は高速である必要はありません。

ご協力いただきありがとうございます!

考えられる解決策は次のとおりです

答えに基づいて(どうもありがとう!)、私はあなたがこの要点で見つけるかもしれないいくつかのコードをハッキングしました。この例をsin^2 * Gauss PDFで実行すると、特定の条件を満たす(円の内側にある)2000個のランダムなランダム変量がPDF上にプロットされます。多分それは他の人にも役立つでしょう。

4

2 に答える 2

2

それで、あなたはPDFを持っていて、このPDFに従ってF(x,y)ペアを生成し、x配布したいですか?y

多変量バージョンの反転手法(wiki)を使用できない限り、棄却サンプリングが最適です。

于 2012-06-18T14:14:24.933 に答える
0

変数XとYの場合、Xの独立分布でxを生成し、xが与えられたYの分布でayを生成するだけで、2つの単変量分布のサンプリングに分離できませんか?

于 2012-06-18T14:14:28.160 に答える