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