model{
   for(i in 1:s){
      k[i]~dnorm(mu,lambda)
   }
   # Prior
   mu~dunif(20,200)
   sigma~dunif(0,30)
   lambda<-pow(sigma,-2)
}