我使用这两个包来进行贝叶斯分析,但有一些我不明白的差异:
首先包rjags允许适应阶段,带jags.model函数,而包r2jags没有这个阶段,带函数jags(或jags.parallel)从后验分布开始采样。自适应阶段是否包含在该功能中,或者包r2jags没有考虑它?
其次,在rjags,我可以说这两段代码是相似的吗?
jmod <- jags.model(file="somefile.txt", data = data, n.chains=3)
update(jmod,100)
jsample <- coda.samples(jmod, n.iter=100, variable.names=par)
Run Code Online (Sandbox Code Playgroud)
和
jmod <- jags.model(file="somefile.txt", data = data, n.chains=3)
jsample <- coda.samples(jmod, n.iter=200,n.burnin=100, variable.names=par)
Run Code Online (Sandbox Code Playgroud)
也就是说,带update功能的老化阶段也可以在coda.samples功能上完成?谢谢你。