model{ for(i in 1:n){ y[i]~dnorm(mu[i],phi) mu[i] <- alpha + beta*x[i] } phi~dgamma(0.001,0.001) alpha~dnorm(0,0.001) beta~dnorm(0,0.001) }