MCMCpack抑制MCMCmetrop1R功能輸出

Question

函數MCMCmetrop1R可以選擇使用選項verbose=FALSE或verbose=0禁止輸出到屏幕。 但是，當函數完成其操作時，這並不會阻止函數報告以下內容。 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@大都會的接受率為0.15862 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@

我嘗試在MEMCmetrop1R函數周圍使用invisible()函數，但這沒有任何區別。 還有其他建議停止報告接受率的功能嗎？

Answer 1

當我今天遇到相同的問題時，這是我的解決方案：您可以使用R.utils包中的capture.output或captureOutput() 。 據說后者更快（請參閱http://www.jottr.org/2014/05/captureOutput.html ）。

R.utils::captureOutput(expr={
  post.sample <- MCMCmetrop1R(loglik, theta.init = c(0, 1, 1), ...)
})

注意，這也會影響所有其他輸出。 因此，如果您設置verbose > 0 ，則也不會顯示。