model{
   k1~dbin(theta1,N)
   k2~dbin(theta2,N)
   # Prior
   theta1~dbeta(1,1)
   theta2~dbeta(1,1)
   # Difference between two thetas
   delta<-theta1-theta2
}