연습을 위해 Darren Wilkinson이 블로그 게시물 Gibbs sampler in various languages (revisited)에 예제 프로그램을 작성했습니다. 아래 코드가 나타납니다. 이 코드는, SBCL 1.0.56을 사용 buildapp를 사용하여 코어 이미지를 생성하고이 타이밍이 계산 된 방법 이었기 때문에 time ./gibbs >
동적 디스패치를 피하기 위해 Barton 및 Nackman 트릭을 사용하여 클래스를 구현하려고합니다. (필자는 성능이 중요한 MCMC 코드를 작성하고 있습니다.) 저는 C++ 전문가는 아니지만 기본 트릭은 다른 곳에서 저에게 일하고 있습니다. 그러나 이제는 두 번째 파생 클래스를 템플릿 화해야하는 경우가 있습니다. 이것은 문제를 일으키는 것 같다. 내
MCMC method을 사용하여 이중 집적을 수행해야합니다. 나는 이미 정확한 결과를 가지고 romberg과 doublequad 통합을 사용하여 그것을 수행했다. 또한 MCMC 통합을 사용하여 결과를 비교해야합니다. PyMC을 이해하기가 어려웠습니다. 개요는 다음과 같습니다. 일부 timeseries 데이터가 있는데 어떤 배포판이 적합한 지 알아야합니다.
일부 MCMC 배포를 시뮬레이션하기 위해 JAGS를 사용하고 있습니다. 상당한 규모의 배포판에 대한 지원을 볼 수 있습니다. 그러나 JAGS에서 정의되지 않은 특정 배포본을 샘플로 추출하고 싶습니다. 누구든지 우리 고유의 커스텀 배포판을 작성하는 방법을 알고 있습니까? 감사합니다.