The IBD process along four chromosomes

In this note, we present the continuous-time Markov rate matrix that models identity by descent (ibd) patterns among four chromosomes in a population. The equilibrium distribution of this Markov process along a chromosome is the set of 4-gene state probabilities given by the Ewens sampling formula. This model will facilitate inference of identity by descent among the four chromosomes of a pair of individuals, using data at dense SNP loci among which there may be linkage disequilibrium.