"omne trium perfectum"또는 "threes로 오는 모든 것이 완벽합니다."라틴어가 약간 있습니다. 나는 3 개의 게시물을 연속해서 쓰려고하지 않았다. 마르코프 체인 몬테 카를로 (MCMC), 그러나 때때로 별들은 그 이야기가 계속해서 쓰여지는 방식으로 정렬됩니다.
그러나 첫째로, 농담 :
Bayesian과 Frequentist는 같은 일을 신청합니다. 면접관은 Frequentist에게 "죄송합니다. 당신은 훌륭한 후보자이지만 베이지안과 함께하기로 결정했습니다. "Frequentist는 실망하고 있지만"나는 이해합니다. 그러나 균형을 깨뜨린 것이 무엇인지 말해 줄 수 있습니까? "면접관은"종이에, 너는 둘 다 똑같은 자격 조건을 가지고있어. 그러나 베이지안은 사전 경험이 있습니다. "
이 게시물은 실제로 Google이 제공 한 의견 중 일부입니다. JMP WinBUGS를 사용하여 자주 모델에 맞는 사용자. 그는 MCMC 진단 추가 기능을 사용하여 후방 샘플을 탐색 할 때 WinBUGS 출력을 적절한 형식으로 가져 오는 것이 어려웠습니다. Wolfinger (1998)는 항공기 부품에 사용 된 복합 소재의 인장 강도 측정을 분석합니다 (표 4의 데이터는 방겔 (Vangel, 1992) ) 8 개의 매개 변수가있는 모델을 사용합니다. 단일 체인에 대해 11,000 개의 사후 샘플 (처음 1,000 개가 던져 질 번인 샘플이라고 가정)을 요청하면 WinBUGS에서 2 개의 출력 파일을 생성합니다. 첫 번째 index.txt는 시작 및 종료 관찰 번호가있는 매개 변수의 목록입니다. 이 예에서 index.txt의 내용은 다음과 같습니다.
a [1] 1 11000
a [2] 11001 22000
a [3] 22001 33000
a [4] 33001 44000
a [5] 44001 55000
뮤 55001 66000
sigma2.a 66001 77000
sigma2.e 77001 88000
chain 1.txt라는 두 번째 파일에는 두 개의 열이 있습니다. 첫 번째 열은 반복 횟수 1 ~ 11,000을 8 번 반복 한 것이고 두 번째 열은 8 개의 매개 변수에 대한 샘플 값을 포함합니다. 본질적으로 88,000 x 2 매트릭스. 이러한 출력 파일을 JMP MCMC 진단 추가 기능 (체인 내의 샘플 및 열이 매개 변수라고 가정)에서 어떻게 사용할 수 있습니까? 무료로 사용할 수 있습니다. WinBUGS to JMP Conversion 추가 기능 ! (다운로드하려면 무료 SAS 프로파일이 필요합니다.) 추가 기능은 WinBUGS 출력 파일 (index.txt 및 chain 1.txt, chain 2.txt 등)의 디렉토리 위치를 묻고 MCMC 샘플이라고하는 JMP 테이블을 생성합니다 새로운 이름으로). 그림 1을 참조하십시오.
그림 1. WinBUGS 출력 파일에서 변환 된 MCMC 샘플 JMP 테이블
몇 가지 사항을 고려해야합니다. 첫째, 그림 1의 데이터 테이블에는 번인으로 의도 된 처음 1,000 개의 샘플이 포함되어 있습니다. 테이블에서 행을 삭제하거나 처음 1000 행을 선택하고 마우스 오른쪽 버튼을 클릭하고 '제외 / 제외 해제'를 선택하여 분석에서 행을 제외 할 수 있습니다. 이로 인해 MCMC 진단 추가 기능이없는 것처럼 이러한 행을 무시합니다. 번인 샘플을 자세히 조사하고, 즉 플롯에서의 수렴 속도를 평가하고 싶다면 언 클리어 할 수 있습니다. 둘째, 매개 변수 a가 a_i로 이름이 바뀌 었습니다. 변수 이름이 변환 프로세스에 의해 배열로 해석되지 않도록 의도적이었습니다. WinBugs 내에서 매개 변수 이름을 지정할 때이 점을 기록하십시오. 변환으로 인해 이름이 중복 된 열이 생길 수 있습니다.
참고
Wolfinger RD. (1998). 베이지안 시뮬레이션을 이용한 분산 구성 요소 모델의 허용 오차 간격. 저널 품질 기술 30 : 18-32.
