model {
for (j in 1:J)
{
y[j] ~ dnorm (theta, tau.y[j])
tau.y[j] <- pow(sigma.y[j], -2)
}
theta ~ dnorm (0.0,1.0E-6)
}