想问楼主可以解答模型转化的问题嘛?我在用R2jags,想把以下的模型改成R2jags能识别的语法不知道怎么改,跪求解答!
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)
}"
|