model {
 for(i in 1:N) {
  mu[i] <- beta * x[i]
  y[i] ~ dnorm(mu[i],psi)
 }
 beta ~ dnorm(0.5, 5)
 psi ~ dgamma(3, 1)
 errvar <- 1/psi
}
