在adaptMCMC包中控制MCMC功能的输出消息

时间:2014-09-18 05:07:16

标签: r mcmc

很想知道是否有办法禁用在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采样器比产生大量标准输出的采样更精细。

0 个答案:

没有答案