mod1 <- "model
{
for (i in 1:numGroups)
{
y[i] ~ dbin(p[i],n[i])
logit(p[i]) <- eta[i]
eta[i] ~ dnorm(gamma[ind[i]], tau1)
ind[i] <- step(theta[i]) + 1
}
for (i in 1:numGroups)
{
theta[i] ~ dnorm(0, tau2)
}
#mu1<-logit(phi1)
#mu2<-logit(phi2)
gamma[1] <-logit(phi1)
gamma[2] <-logit(phi2)
}"
以上是一个用BRugs包计算后验的模型,想用R2jags实现,跪求模型需要怎么改呢?