the denominator as in the conjugate case, i.e., equation (7) with k = 2 and v = (1, 1, 1), and p(n∣H0)=n!2n12+2n11!n12!n22![IBe(0.5,2n11+n12+2,n12+2n22+1)+Be(2n11+n12+1,n12+2n22+2)−IBe(0.5,2n11+n12+1,n12+2n22+2)], where Be(a, b) = Γ(a)Γ(b)/Γ(a + b) and IBe(x,a,b)=∫0xza−1(1−z)b−1dz is the incomplete beta function (which is straightforward to calculate). Under the Dir(1,1,1) prior the probability that the inbreeding coefficient f is greater than 0 is 0.67, showing that the implications of innocuous looking priors should be carefully examined.
