model { for (i in 1:K) { for (j in 1:n) { Y[i, j] ~ dnorm(eta[i, j], tauC) eta[i, j] <- phi[i, 1] / (1 + phi[i, 2] * exp(phi[i, 3] * x[j])) } phi[i, 1] <- exp(theta[i, 1]) phi[i, 2] <- exp(theta[i, 2]) - 1 phi[i, 3] <- -exp(theta[i, 3]) for (k in 1:3) { theta[i, k] ~ dnorm(mu[k], tau[k]) } } tauC ~ dgamma(1.0E-3, 1.0E-3) sigmaC <- 1 / sqrt(tauC) varC <- 1 / tauC for (k in 1:3) { mu[k] ~ dnorm(0, 1.0E-4) tau[k] ~ dgamma(1.0E-3, 1.0E-3) sigma[k] <- 1 / sqrt(tau[k]) } }