CまたはPythonでのバイモーダル分布

Question

CまたはPythonでのバイモーダル分布に従ってランダムな値を生成する最も簡単な方法は何ですか？

私はZigguratアルゴリズムやBox-Muller変換のようなものを実装することができますが、すぐに使えるライブラリがあるか、わかりにくい簡単なアルゴリズムがあれば、それが良いでしょう。

Answer 1

2つのモーダルディストリビューションのいずれかの値を選んでいませんか？

http://docs.python.org/library/random.html#random.triangular

あなただけの三角形への通話のためのパラメータの2つのセットの間で前後に切り替えるような音。

def bimodal(low1, high1, mode1, low2, high2, mode2): 
    toss = random.choice((1, 2)) 
    if toss == 1: 
     return random.triangular(low1, high1, mode1) 
    else: 
     return random.triangular(low2, high2, mode2) 

これは、必要なすべてを行う可能性があります。

Answer 2

古くなったストレートフォワードaccept-reject algorithmが常にあります。それがジョニーフォンノイマンのために十分だったなら、それはあなたのために十分であるはずです;-)。