가역 점프 MCMC 코드 (Matlab 또는 R)


14

가역 점프 MCMC를 위해 잘 작성된 코드 (Matlab 또는 R)를 알고 있습니까? 바람직하게는 주제에 대한 논문을 보완하기위한 간단한 데모 애플리케이션으로, 프로세스를 이해하는 데 유용합니다.


BUGS Book의 288 페이지에 OpenBUGS에 리버시블 점프 MCMC를 추가하는 것에 대한 논의가 있었고 WinBUGS는 점프 모듈을 가지고 있다고 생각합니다. JAGS에 대한 유사한 개발에 대해 들어 본 적이 있습니까? 다양한 차원을 수행하는 효과를 얻기 위해 JAGS 코드를 작성하는 방법이 있습니까? 단순히 필요한 것보다 더 큰 치수를 정의하는 것입니까?
Jan Galkowski

답변:


12

RJMCMC는 인용 논문 인 1995 년 논문 에서 Peter Green 에 의해 소개되었습니다 . 그는라는 포트란 프로그램 작성 AutoRJ 자동 RJMCMC를 들어; 그의 페이지는 David Hastie의 C 프로그램 AutoMix와 연결되어 있습니다. Scott Sisson2005 년 논문 의 표 1에는 다양한 RJMCMC 알고리즘에 사용할 수있는 무료 소프트웨어 목록이 있습니다 . 구글 검색은 또한 글래스고 대학교 (University of Glasgow)의 그룹에서 의사 코드를 직접 프로그래밍하려는 경우 원리를 이해하는 데 유용한 의사 코드를 찾습니다 .


7

King et al.의 인구 생태학대한 베이지안 분석 책 . 인구 생태의 맥락에서 RJMCMC를 설명합니다. 나는 매우 명확한 설명을 발견했으며 부록에 R 코드를 제공합니다.

이 책에는 관련 웹 페이지 가 있지만이 책에서 찾은 코드 중 일부는 웹 사이트에 없습니다.


RJMCMC에 WinBUGS 코드를 제공합니까? WinBUGS는 RJMCMC가 가능합니까?
Curious

1
@Tomas 누군가가 내 대답을 편집하고 R을 WinBUGS로 변경했습니다. 내가 기억하는 것에서,이 책은 R 코드를 제공합니다.
csgillespie

감사합니다. (+1) Btw, RJMCMC가 인구 생태에 유용하다고 생각 했습니까? 나는 그것이 유용 할 수있는 예를 아직 찾지 못했다. 나는 표준 MCMC와 같은 데이터 증강을 항상 사용할 수 있고 그것은 나를 위해 일한다. 내가 책을 가지고 있지 않다.
Curious

@Tomas이 책은 모델 선택을 위해 RJMCMC를 사용합니다. 나는 설명이 아주 훌륭하다는 것을 알았다. 그래도 분노에 RJMCMC를 사용해야했습니다.
csgillespie

"분노"?? 무슨 소리 야?
Curious

1

@onestop의 답변에 한 가지 세부 사항을 추가하십시오 .Olivier Cappé (CT / RJ MCMC)에서 릴리스 한 C 소프트웨어는 가역 점프 MCMC 알고리즘을 이해하는 데 매우 도움이됩니다 (특히 출생 및 사망 확률을 설계하는 방법 병합 이동). 소스 코드에 대한 링크는 다음과 같습니다. http://perso.telecom-paristech.fr/~cappe/Code/CTRJ_mix/About/


1

Jailin Ai는 Leeds에서 석사 논문의 일부로 R 코드와 함께 RJ MCMC를 상당히 훌륭하게 보여줍니다 (그린의 원본 논문에 매우 가깝지만). 또한 Green의 1995 년 논문에 포함 된 변경점 문제에 대한 심도있는 예를 제공합니다.

논문 및 코드를 여기에서 찾으십시오.

http://www1.maths.leeds.ac.uk/~voss/projects/2011-RJMCMC/


0

Nando de Freitas는 신경망 매개 변수 추정을위한 가역 점프 MCMC 알고리즘 사용에 대한 데모를 제공합니다. 이 모델은 뉴런 수, 모델 매개 변수, 정규화 매개 변수 및 노이즈 매개 변수를 추정 할 랜덤 변수로 처리합니다.

코드와 작성은 http://www.cs.ubc.ca/~nando/software.html에서 볼 수 있습니다.

당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.