我正在尝试在JAGS中运行我的第一个混合模型,但是遇到了问题。我已经修改了代码,以获得连续的结果来代替泊松使用。根据其他帖子,我显然将JAGS与模型中某处的反陈述相混淆。
运行时错误:第4行的编译错误。Y的子集表达式中的尺寸不匹配
下面是我的模型-如果需要添加其他信息,请告诉我。
for(i in 1:N) {
for (j in 1:k[i]) {
Y[i,j] ~ dpois(lambda[i,j])
log(lambda[i,j]) <- true.y[i,j,z[i]]
true.y[i,1] <- m2 * (X[i,j]) + a1[i]
true.y[i,2]<- m2 *(X[i,j] ) + c +a2[i]
}
z[i]~dcat(pi[1:2])
a1[i]~dnorm(0,tau.a[1])
a2[i]~dnorm(0,tau.a[2])
}
pi[1:2] ~ ddirch(alpha[])
m2 ~ dnorm(0,0.0001)
c ~ dnorm(0,0.0001);
tau.a[1] ~ dgamma(0.001,0.001)
tau.a[2] ~ dgamma(0.001,0.001)
}
'''