p.post.found function

computes the posterior probability of observations of a founder