model       {for (i in 1:17) 
	{y[i] ~ dnorm(mu[i],tau.y)
                     mu[i] <- a + b*(x[i]-mean(x[]))
	}
# priors
	a ~ dnorm(0.854,.091)
	b ~ dnorm(0.278,0.036)
	tau.y ~ dgamma(4.274,3.759)
	sig2.y <- 1/tau.y
                  }

Data
list(x=c(3.76, 2.22, 2.70, 2.72, 3.00, 1.71, 1.27, 2.80, 1.67, 0.41, 1.29, -0.29, 2.07, 0.69, 1.06, 0.67, 0.37), y=c(2.51,2.5,2.49,2.35,2.29,2.08,2.01,1.9,1.9,1.72,1.68,1.68,1.67,1.62,1.34,1.15,.85))

Inits
list(tau.y=1)