ON THE 2-LOCUS SAMPLING DISTRIBUTION

Two methods are discussed for evaluating the distribution of the configuration of unlabeled gametic types in a random sample of size n from the two-locus infinitely-many-neutral-alleles diffusion model at stationarity. Both involve finding systems of linear equations satisfied by the desired probabilities. The first approach, which is due to Golding, is to include additional probabilities in the system that allow some members of the sample to be specified at only one locus. The second approach, which is new, considers the joint distribution of the sample configuration and the number of recombination events since the time of the most recent common ancestor. The first approach is used for numerical computation, whereas the second approach is used to derive a two-locus version of Hoppe's urn model. The latter permits efficient simulation of the two-locus sampling distribution, provided the recombination parameter is not too large. © 1990, Springer-Verlag. All rights reserved.