0

想知道是否有任何方法可以禁用在包中的MCMC()函数运行期间打印的输出消息。adaptMCMC我多次运行这个函数,这些打印语句填满了我的磁盘。

这是它的外观(直接来自手册):

library(adaptMCMC)
#log-pdf to sample from
p.log <- function(x) {
B <- 0.03 # controls 'bananacity'
-x[1]^2/200 - 1/2*(x[2]+B*x[1]^2-100*B)^2
}

samp <- MCMC(p.log, n=200, init=c(0, 1), scale=c(1, 0.1),
adapt=TRUE, acc.rate=0.234)

  generate 200 samples 

输出消息是:

  |========================================================================================================================================================| 100%

我的 gibbs 采样器比导致大量标准输出的采样器复杂得多。

4

0 回答 0