model { for( i in 1 : N ) { for( j in 1 : T ) { Y[i , j] ~ dnorm(mu[i , j],tau) mu[i , j] <- alpha[i] + beta[i] * (t[i,j] - 6.5) + gamma * (y0[i] - mean(y0[])) } alpha[i] ~ dnorm(alpha0,tau.alpha) beta[i] ~ dnorm(beta0,tau.beta) } tau ~ dgamma(0.001,0.001) sigma <- 1 / sqrt(tau) alpha0 ~ dnorm(0.0,1.0E-6) tau.alpha ~ dgamma(0.001,0.001) beta0 ~ dnorm(0.0,1.0E-6) tau.beta ~ dgamma(0.001,0.001) gamma ~ dnorm(0.0,1.0E-6) }