model{
   for(i in 1:N){
      K[i,1]~dt(mu[K[i,2]],1/sigma[K[i,2]]^2,tau)
   }
   # Priors
   for(j in 1:2){
     mu[j]~dnorm(mean[j],1/(100*sd[j])^2)
     sigma[j]~dunif(sd[j]/1000,1000*sd[j])
   }
   tau<-nutau+1
   nutau~dexp(1/29)
}