# Example from WinBUGS the movie
model {
  for (i in 1:N) {
    y[i] ~ dnorm(mu,tau)
    }
  mu ~ dnorm(0,.001)
  tau ~ dgamma(.0001, .0001)
  sigma <- sqrt(1/tau)
 }

#data
list(N=10, y=c(1, 2, 2, 3, 3, 3, 3, 4, 4, 5)) 

#initial values
list(mu=0, tau=1)