model
		{
		   for( i in 1 : N ) {
		      y[i] ~ dnorm(mu,tau)
		   }
		mu ~ dnorm(xi,kappa)
		tau ~ dgamma(alpha,beta)
		sigma <- 1/sqrt(tau)
		}
