model {
for(i in 1:T) {
   for(j in 1:N) {
      y[i,j] ~ dnorm(mu[i],tau[i])
   }
   mu[i] ~ dnorm(0.0,0.001)
   tau[i] ~ dgamma(0.001,0.001)
}
   delta <- mu[2]-mu[1]
}