On Computing the Coalescence Time Density in an Isolation-With-Migration Model With Few Samples

IN a recent article, [Wang and Hey (2009)][1] consider estimation of the parameters in an isolation-with-migration (IM) model for two species. For each locus, the data X consist of two samples, and therefore the probability of the data depends only on the time to the most recent common ancestor (