我正在尝试在 JAGS 中运行二项式 beta 模型（请参阅下面的示例代码）。我不断收到错误消息：错误：尝试运行 JAGS 模型时遇到以下错误：
Error in node a0
Slicer stuck at value with infinite density
我正在努力理解。我想也许初始条件是将 beta 分布发送到参数空间的无限区域，但经过一些调查似乎并非如此。关于此错误意味着什么或如何调整代码以适应它的任何想法？
我把我的代码和一些编造的样本数据放在下面。这是我在数据集中可能期望的数据类型。
#Generate some sample data
counts = c(80,37,10,43,55,23,53,100,7,11)
n = c(100,57,25,78,55,79,65,100,9,11)
consp = c(1.00, 0.57, 0.25, 0.78, 0.55, 0.79, 0.65, 1.00, 0.09, 0.11)
treat = c(0.5,0.5,0.2,0.9,0.5,0.2,0.5,0.9,0.5,0.2)
#Model spec
model1.string <-"model{
for (i in 1:length(counts)){
counts[i] ~ dbin(p[i],n[i])
p[i] ~ dbeta( ( mu[i] * theta[i]) , ((1-mu[i])*theta[i]))
mu[i] <- ilogit(m0 + m1*consp[i] + m2*treat[i])
theta[i] <- exp(n0 + n1*consp[i])
} …