data { int N; // the sample size int k; // the number of A_1 alleles observed } parameters { real p; // the allele frequency } model { // likelihood // k ~ binomial(N, p); // prior // p ~ uniform(0.0, 1.0); }