2017-03-07 8 views
0

나는 두 개의 배열 (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] 

는 I 샘플의 N 숫자를 생성하고자하는 (실제로는 PDF를 더 미세하게 샘플링되고 합 1) (X, Y와 다음은

제가 가지고 무엇 from val_1, val_2).

답변

1

mcmc를 사용해야하는 이유가 있습니까?

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

뜻 연관된 확률과 pa에서 값을 선택 하였다. 하나의 calladd에서 N 개의 샘플을 얻으려면 size 속성을 설정하십시오. 이것으로 충분하지 않습니까? 당신은 배포판을하고 링크와 같은 플롯을 생성 할 경우

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

아, 그리고, 내가 사용하는 것이 matplotlib.pyplot.hist(array)