有一些拉丁文表示“omne trium perfectum”或“三分之一是完美的”。我還沒打算連續寫三個帖子馬爾可夫鏈蒙特卡洛 （MCMC），但有時候星星的對齊方式使故事繼續寫下來。
但首先，一個笑話：
貝葉斯和Frequentist申請同樣的工作。採訪者告訴Frequentist，“很抱歉告訴你這件事。你是一個很棒的候選人，但是我們決定和貝葉斯一起去。“常旅客很失望，但是說，”我明白了，但是你能告訴我這對她有利的平衡嗎？“採訪者說：”在紙面上，你們都有相同的資格。然而，貝葉斯先前有過經驗。“
這篇文章實際上是我們收到的一些反饋的結果 JMP 經常使用WinBUGS適合模型的用戶。雖然他找到了 MCMC診斷有助於探索後驗樣本的加載項，他發現很難將WinBUGS輸出轉換為適當的形式。他向我發送了以下示例：Wolfinger（1998）分析了飛機部件中使用的複合材料的拉伸強度測量結果（表4中的數據） 範格爾（1992） ）使用具有8個參數的模型。對於單個鏈，請求11,000個後驗樣本（假設前1000個是將被拋出的老化樣本）在WinBUGS中生成兩個輸出文件。第一個是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的文件包含2列。第一列是重複8次的迭代次數1到11,000，而第二列包含8個參數的採樣值;基本上是88,000 x 2矩陣。這些輸出文件如何被JMP MCMC Diagnostic插件使用（假設行是鏈中的樣本，列是參數）？只需申請免費提供 WinBUGS到JMP轉換加載項 ！（下載需要免費的SAS配置文件。）外接程序要求WinBUGS輸出文件的目錄位置（index.txt和鏈1.txt，鏈2.txt等）並創建一個名為MCMC Samples的JMP表（可以保存）用一個新的名字）。見圖1。
圖1。從WinBUGS輸出文件轉換的MCMC樣本的JMP表
要考慮幾點。首先，圖1中的數據表包括預計為老化的前1,000個樣本。我可以從表中刪除行，也可以通過選擇前1000行來排除分析中的行，右鍵單擊並選擇“排除/取消排除”。這會導致MCMC診斷加載項忽略這些行，就好像它們不存在一樣。如果我想更詳細地檢查老化樣本，比如說，以評估跟踪圖中的收斂速度，我總是可以在以後取消。二，注意參數如何a 被重命名為a_i。 這是故意的，因此轉換過程不會將變量名稱解釋為數組。在WinBugs中命名參數時請記下這一點，因為轉換可能會導致列具有重複的名稱。
參考
Wolfinger RD。（1998）。使用貝葉斯模擬的方差分量模型的容差區間。質量技術雜誌 30：18-32。
