我正在运行 this 简单 Rjags 教程中的代码。请注意,参数的初始值由
model.inits <- list(tau=1, beta0=1, beta1=0, beta2=0)
指定。我想测试生成的 MCMC 链实际上是从这些初始值开始的。
因此,出于测试目的,在
jags.model()
内,我设置了n.adapt=0
,因此不会发生适应。打印 head(model.samples[[1]])
会导致第一行的 beta0
为 0.58,beta1
为 0.97,而 beta2
为 0.54。因此,在第一次迭代之后,没有一个参数接近于它们的初始化。
希望能解释一下我是否遗漏了什么。