很想知道是否有办法禁用在MCMC()
包中运行adaptMCMC
函数期间打印的输出消息。我正在多次运行这个函数,这些打印语句正在填满我的磁盘。
以下是它的外观(直接来自manual):
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采样器比产生大量标准输出的采样更精细。