各位,大家新年好!在做课题时,遇到一个问题,建立的的模型总是UPDATE后总是出现TRAP66,不知道什么 原因,求帮助。
数据共1660条记录,6个变量。不知道什么原因,请指教。
模型如下:
/*非时空模型*/
model
{
for (i in 1:m)
{
for (t in 1:T)
{
#Binomial likelihood for observed counts
newsp[i,t]~dbin(p[i,t],pop[i,t])
logit(p[i,t])<-a[1]+a[2]*gdp[i,t]+a[3]*sickbed[i,t]+a[4]*ncjmrjcsr[i,t]+a[5]*ncrk[i,t]
}
}
#Other priors
for (i in 1:5){
a[i]~dnorm(0.01,1.0E-2)
}
}