model { for (i in 1:n) { y[i] ~ dnorm(muy[i],tau); muy[i] <- m * x[i] + c; } m ~ dnorm(0.0, 1.0E-4); c ~ dnorm(0.0, 1.0E-4); tau ~ dgamma(1.0E-3, 1.0E-3); sigma <- 1.0/sqrt(tau); # prediction muy1 <- m * x1 + c; y1 ~ dnorm(muy1,tau); }