我是 STAN 和贝叶斯建模的新手,我正在尝试使用brms
包使用 COM-泊松分布来拟合模型,以获取一些分散的计数数据。但是,我正在努力应对 STAN 的发散和崩溃。警告包括
Divergent transistions after warmup
fit contains no samples
Bulk ESS low
Tail ESS low
Largest Rhat was 4 indicating chains have not mixed
数据:数据集
我的模型是
brm(EstAgeInYrs ~ poly(KnwnAgeInYrs,3) + (1|ReaderID),
cores = 4, iter = 1000
data = DataHI, family = "com_poisson")
我通过使用glmmTMB
.
有什么建议可以在这方面取得进展吗?我已经确定第一个年龄 = 1 和年龄 = 2 var / mean > 1,并且我知道brm( ... family = com_poisson)
与过度分散的斗争,但其余数据明显低于分散(<0.5)。
任何帮助将不胜感激