model {
  for (i in 1:Ny) {
    y[i] ~ dnorm(mu[x[i]], 1/sigma[x[i]]^2)
  }

  for(j in 1:2) {
    mu[j] ~ dnorm(meanY[j], 1/(100*sdY[j])^2)
    sigma[j] ~ dunif(sdY[j]/100, sdY[j]*10)
  }
}
