model
{
for(i in 1:n){
   Y[i] ~ dnorm(mu[i],tau)
   mu[i] <- beta[1] + beta[2] * x[i]
 }

for(j in 1:2){
   beta[j] ~ dnorm(0,0.0001)
 }
tau ~ dgamma(0.001,0.001)
sigma <- 1 / sqrt(tau)
}