model{
   # Likelihood function
   for(j in 1:ntest){
      for(i in 1:nsubj){
         K[i,j]~dnorm(mu[i],lambda)
      }
   }
   # Prior
   # mu
   for(i in 1:nsubj){
      mu[i]~dunif(0,300)
   }
   # lambda
   lambda~dgamma(0.001,0.001)
   sigma<-pow(lambda,-0.5)
}