	model
	{
		for( i in 1 : N ) {
			Y[i] ~ dnorm(mu[i],tau)
			mu[i] <- alpha + beta * (x[i] - xbar)
		}
		tau ~ dgamma(0.001,0.001) sigma <- 1 / sqrt(tau)
		alpha ~ dnorm(0.0,1.0E-6)
		beta ~ dnorm(0.0,1.0E-6)	
	}
