关于rjags中的以下尾波样本,
samples<-coda.samples(jags, c('B', 'A'), 5000,thin=5)
迭代设置如下 5000 以及 thin 参数设置为 5 在贝叶斯分析中,是否有任何机制或规则来设置这两个参数。
5000
thin
5
你的意思是“有没有办法为这些参数定义‘正确’或‘最佳’值? 如果是这样,那么答案是“否”。不是R或任何语言。您需要查看模型拟合的诊断(如轨迹图、Gelman-Rubin、Raftery-Lewis和Heidelberger-Welch统计),并调整MCMC过程的参数,直到获得“可接受”的结果。这是一门艺术,而不是一种确定性算法。
您还应该查看您使用的链的数量,以及每条链的老化、变薄和其他参数。