R2WinBUGS - 如何再次將WinBUGS計算結果重新加載到bug對象

Question

我使用R中的R2WinBUGS包在WinBUGS（百萬次迭代）中運行了一個looong計算：

bugs.object <- bugs(...)

但是R墜毀了。 如何在不再運行winbugs的情況下將bugs.object重新加載到R中？ 我試過這個（我有3個鏈）：

out <- read.bugs(paste("coda", 1:3, ".txt", sep = ""))

但out數據結構與bug對象完全不同（因為它是不可用的）。 我試着用as.bugs.array轉換它：

bugs.object <- as.bugs.array(out, model.file = "ttest.txt", n.iter = 1000000, n.burnin = 300000, n.thin = 2, program = "WinBUGS")

但它不起作用。 請幫忙。 謝謝。

Answer 1

您可能正在閱讀錯誤消息，其中R內存bugs.array創建bugs.array對象。

您可以通過在bugs函數中設置codaPkg=T語句來解決此問題。 這會將CODA文件保存在指定的工作目錄中，而不是創建R2WinBUGS對象（在R崩潰之前）。 然后你可以使用coda包中的read.mcmc讀取coda文件，如果你真的想要，可以將mcmc對象轉換為bugs.array 。

如果您的MCMC太大或者沒有足夠的內存用於R，這可能不起作用。