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.87,2.33,2.81,2.83,3.11,1.82,1.37,2.91,1.78,0.52,1.4,-0.18,2.18,0.79,1.17,0.78,0.48), 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)