model {

   # likelihood
   n.ij ~ dbin(n, P.ij)

   # prior
   P.ij ~ dbeta(nu, omega)

}

list(k = 201, n = 372, nu = 1, omega = 1)