The integral equation approach for models of clines

We have described above the outline of a general recursive method for deriving the likelihood of a model. In order to implement this procedure automatically we need to develop two sets of algorithms, (i) for taking a pedigree in a simple numerical representation, and specifying the order in which we peel individuals off, a non-trivial problem for large pedigrees, and (ii) for carrying out the necessary algebraic manipulations so that L is ultimately obtained, as a function of the parameters which need to be estimated in the model. These have been developed: (i) for zero-loop pedigrees, (ii) for an autosomal locus with three genotypes and arbitrary penetrances, and programmed, and are in the process of generalization to more complex cases.