MCMCpack 抑制 MCMCmetrop1R 功能输出

Question

函数 MCMCmetrop1R 可以使用选项verbose=FALSE 或verbose=0 抑制其输出到屏幕。但是，当函数完成其操作时，这不会停止函数报告以下内容。 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@ Metropolis 录取率为 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，这个也不会出现。