R中JAGS问题：MCMC链未从设定的初始值附近启动

核心原因：初始值是MCMC链的起点，而非第一个保存的样本 #

你看到的head(model.samples[[1]])是第一次MCMC抽样后的结果（迭代1），而你设置的初始值是链的起始状态（迭代0）。即使关闭适配（n.adapt=0），RJAGS也会从初始值出发，执行一次MCMC采样步骤生成第一个保存的样本，而非直接保留初始值。

验证初始值的正确方法 #

要确认链是否从指定初始值启动，不要查看抽样样本，而是直接查看模型的初始状态：

# 初始化模型后直接查看状态（迭代0）
model <- jags.model("model.txt", data = data, inits = model.inits, n.adapt = 0)
print(model$state())

运行这段代码会输出你设置的初始值：tau=1, beta0=1, beta1=0, beta2=0，证明链确实从指定值启动。

为什么第一次抽样参数变化大？ #

RJAGS默认使用Gibbs抽样，每个参数会根据其他参数的当前值，从自身的条件后验分布中采样。如果你的初始值和参数的真实后验分布差异较大，第一次采样就会生成远离初始值的参数值，这是MCMC算法的正常行为——算法会逐步收敛到后验分布，初始步的波动是预期内的。

内容的提问来源于stack exchange，提问作者bob