想知道是否有任何方法可以禁用在包中的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 采样器比导致大量标准输出的采样器复杂得多。