model {

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

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

}

# MN allele counts from Table 2.3 (p. 44)
list(k = 201, n = 372, nu = 1, omega = 1)

# data from this moring
list(k = 1312, n = 4300, nu = 1, omega = 1)