model{
   # Likelihood function
   for(i in 1:m){
      k[i]~dbin(theta[i],N[i])
   }
   # Prior
   for(i in 1:m){
      theta[i]~dbeta(1,1)
   }
   # Difference
   delta<-theta[1]-theta[2]
}