model{
   for(j in 1:Nt){
      for(i in 1:Ns){
         K[i,j]~dbern(theta[i,j])
      }
   }
   # Priors
   for(j in 1:Nt){
      for(i in 1:Ns){
         theta[i,j]<-pt[j]*ps[i]
      }
   }
   for(i in 1:Ns){
      ps[i]~dbeta(1,1)
   }
   for(j in 1:Nt){
      pt[j]~dbeta(1,1)
   }
}