model {

   # genotype frequencies
   x[1] <- p.a*p.a + f*p.a*(1-p.a)
   x[2] <- 2*(1-f)*p.a*(1 - p.a)
   x[3] <- (1-p.a)*(1-p.a) + f*p.a*(1-p.a)

   # likelihood
   n[1:3] ~ dmulti(x[], N)

   # priors
   p.a ~ dunif(0, 1)
   f <- 0
   w ~ dunif(0, 1)

   # sample size
   N <- sum(n[])
}

list(n=c(53, 95, 38))