model{
   for(i in 1:n){
      k[i]~dnorm(mu,lambda)
   }
   # Prior
   mu~dnorm(0,0.001)
   lambda~dgamma(0.001,0.001)
   sigma<-pow(lambda,-0.5)
}