model{
   for(i in 1:m){
      k[i]~dbin(theta,N[i])
   }
   # Prior
   for(i in 1:m){
      N[i]~dcat(pi)
   }
   for(i in 1:n){
      pi[i]<-1/n
   }
   theta~dbeta(1,1)
}