library(rjags) jags_code = "model{ for (i in 1:n){ Y[i] ~ dnorm(mu[i], tau) mu[i] <- beta0 + beta1*x[i]} beta0 ~ dnorm(0,0.001) beta1 ~ dnorm(0,0.001) eta &