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