2017-03-07 6 views
0

私は2つの配列(pdf_#)を持っており、それぞれがパラメータ(val_#)の確率を含んでいます。どのようにしてこれらの分布からサンプルを採取して、共同事後分布を作成することができますか(つまりcorner plot hereのように)?私は、配列ではなく、サンプリングされる関数を渡すだけのように見えます。ここ確率分布関数をサンプリングするにはemceeを使用しますか?

は、私が持っているものである。

pdf_1 = [.1, .1, .25, .3, .15] 
pdf_2 = [.25, .3, .2, .1, .4] 

val_1 = [2, 3, 4, 5, 6] 
val_2 = [1, 2, 3 ,4 ,5] 

(実際にはPDFファイルをより細かくサンプリングされ、1に合計)

Iは、Xと(サンプルNの数を生成するために、Yを希望しますval_1、val_2から)をpdfsの後に置きます。

答えて

1

mcmcを使用する理由はありますか?

numpy.random.choice(a=val_1, size=N, p=pdf_1) 

意志が関連する確率paから値を選択しました。 1つのcalladdのN個のサンプルについては、size属性を設定します。これで足りないの?

https://docs.scipy.org/doc/numpy/reference/generated/numpy.random.choice.html#numpy.random.choice

ああ、あなたがディストリビューションを持っているし、あなたのリンクのようにプロットを生成する場合、私はmatplotlib.pyplot.hist(array)

を使用します
関連する問題