for (i in 1:3) {
a.female[i] ~ dexp(1)
x.female[i] <- a.female[i]/a.female.sum
}
a.female.sum <- sum(a.female[])
and use the analogous one for the genotype frequency in males.
p.bar <- mean(p[]) p.var <- sd(p[])*sd(p[])
sd(p[]) is the standard deviation (the square root of the variance) of p[].