Estimating the Lambda measure in multiple-merger coalescents

This article is devoted to the construction of a non-parametric estimator of the density of Λ that is based on the observation at a single time of the so-called Site Frequency Spectrum (SFS), which describes the allelic frequencies in a present population sample. First, we produce estimates of the multiple-merger rates by solving a linear system, whose coefficients are obtained by appropriately subsampling the SFS. Then, we use a technique that aggregates the information extracted from the previous step through a kernel type of re-construction to give a non-parametric estimation of the measure Λ. We give a consistency result of this estimator under mild conditions on the behavior of Λ around 0. We also show some numerical examples of how our method performs.PMID:37742787 | DOI:10.1016/j.tpb.2023.09.002

Source: Theoretical Population Biology - September 24, 2023 Category: Biology Authors: Ver ónica Miró Pina Émilien Joly Arno Siri-J égousse Source Type: research